Soft Roof Failure Mechanism and Supporting Method for Gob-Side Entry Retaining

To study the soft roof failure mechanism and the supporting method for a gateway in a gently inclined coal seam with a dip angle of 16 ̋ kept for gob-side entry retaining, and through the methodology of field investigation and numerical and analytical modeling, this paper analyzed the stress evolution law of roof strata at the working face end and determined that the sharp horizontal stress unloading phenomenon along the coal wall side did not appear after the working face advanced. Conversely, the horizontal stress along the gob side instantly decreased and the tensile stress produced, and the vertical stress in the central part of the roof had a higher reduction magnitude as well. An in-depth study indicates that the soft roof of the working face end subsided and seriously separated due to the effect of the front abutment pressure and the roof hanging length above the gob line, as well as certain other factors, including the rapid unloading of the lateral stress, tension and shear on the lower roof rock layer and dynamic disturbance. Those influencing factors also led to rapid crack propagation on a large scale and serious fracturing in the soft roof of the working face end. However, in the gob stress stabilized zone, the soft roof in the gob-side entry retaining has a shearing failure along the filling wall inside affected by the overburden pressure, rock bulking pressure, and roof gravity. To maintain the roof integrity, decrease the roof deformation, and enable the control of the working face end soft roof and the stabilization of the gob-side entry retaining roof, this study suggests that the preferred bolt installation angle for the soft roof situation is 70 ̋ based on the rock bolt extrusion strengthening theory.


Introduction
The retaining of the gob-side entry is to maintain the head entry of current mining panel behind the working face to be reused for the next panel as the tail entry.This technology can effectively increase the coal recovery rate, reduce the roadway development rate, and mitigate the outburst risk, as no pillar is needed for the retained entry, but rather an artificial filling wall is constructed on the gob side with a special support for entry stability [1].
Since the 1950s, pillarless roadways have been constructed worldwide, and extensive studies have been carried out on support resistance [2,3].Previous studies showed that the gateway roof rock mass failure mechanism and supporting methods have a significant impact on gob-side entry retaining, determining the procedure's success.Many researchers categorized five failure types of a gateway roof rock mass, including the compressive stress failure type, tensile stress failure type, shear stress failure type, squeezing and fluidity failure type, and geological structure failure type [4][5][6].Based on the above failure mechanism of gateway roof strata, many active support models were suggested, and a combined supporting technology using a bolt, mesh, and cable was found to be effective for ground control.It was found that a combined support is unable to control the main roof in roadway retaining according to the theory of "given deformation", but it can stabilize the immediate roof with the main roof to maintain the retained gob-side entry very well [7,8].
Rock bolts have become a popular technique for reinforcing rock masses all over the world.Rock bolts are installed to reinforce a fractured rock mass by resisting dilation or shear movement along the fractures.Nemcik et al. [9] determined that the non-linear bond-slip relationship used in the FLAC model for bolting accurately matches the experimental data reported by Ma et al. [10].Yang [11] divided the anchored force evolution process into three stages and found that a surface structure-like bearing plate produces a maximum anchorage force acting directly on the surface of rock that can provide the greatest degree of support.This not only improved the rock mass stress state but also increased the thickness of the reinforcement structure formed by the anchoring force.Zheng and Zhang [12] determined the shear stress and pressure stress distribution equation for the anchored segment and upon studying the stress distribution rule of anchored segment, found that the effect of the pre-tensioned bolt in the soft rock was superior to that in hard rock.Cao [13] showed that the integrity of the supporting system that prevents local failure of surrounding rock from progressing into overall failure is important in rock bolting, so that reinforcing measures should be taken when necessary.Fan [14] designed a roadway heterogeneity controlling technology on a siltstone roof.Hua [15] made use of bolt support and anchor cable reinforcement support technologies inside and beside a retained roadway, respectively, and maintained the retained gob-side entry very well.Yan [16] used pre-tensioned bolts on a roadway roof with medium and fine-grained sandstones by an extrusion lockset to limit its vertical deformation for entry integrity and also analyzed the mechanical mechanism, technical principle and technical characteristics of bolt and cable coupling supports at a mine site.Chen and Bai [17] used bolts with high pre-tensioning, high strength, and a large elongation rate and cables as the basic supporting element, a single hydraulic prop with a metal hinge top beam as a reinforcement support inside the roadway, and a high-water rapid-solidifying material to support the side of the roadway to effectively control the roof rock deformation with medium and fine-grained sandstones in gob-side entry retaining.
In the past experimental studies and supporting theory [4,[14][15][16][17], the bolt support is mainly aimed at the control of the medium-hard roof rock mass, which are generally medium and fine-grained sandstones in the horizontal and flat seams.The supporting theory also provides supporting references for the similar conditions of the roadway in the initial excavation stage and retaining stage.However, the working face end roof and retained roadway roof seriously deformed in the soft roof rock mass and gentle inclined coal seam with a soft roof rock after referencing the existed supporting theory.Therefore, to further develop our understanding of the soft roof failure mechanism and provide a supporting method, the following works were conducted in this article.First, the roof failure phenomenon of a gently inclined coal seam working face end and retained roadway were investigated at a coal mine in the southwest of China.Second, the roof failure mechanism was analyzed through the field investigation and stress evolution law obtained by three-dimensional distinct element code (3DEC).Finally, a roof support was designed based on results of this investigation to maintain roadway stability.

Survey of Study Site
The target coal mine is located in the southwest of China.The mining coal seam is #C19, with an average dip angle of 16 ˝and a mining thickness of 1.1 m.The mining method is the fully mechanized longwall, and gob side entry retaining is used, as shown in Figure 1a.The investigated head entry is buried at 315 m, with a cross-section of 3.8 m in width and 2.1 m in height of the short rib, where resin-anchored rock bolts (Φ20 ˆ2200 mm), cables (Φ12.54 ˆ6300 mm), and reinforcement mesh were used, and three rows of single hydraulic supports with top hinged beams were used as advanced strengthening support, as shown in Figure 1b.The roof rock bolts are non-fully anchored (anchor length is 1.4 m) with a pre-tension of 90 kN, inter-row spacing of 800 ˆ800 mm and bolt angles of 90 ˝, 80 ˝, 70 ˝, and 60 ˝for different bolts, and the extra cables have a pre-tension of 200 kN and an inter-row spacing of 3200 ˆ1200 mm.According to a geological survey, the gateway roof strata typically consist of weak rocks.The immediate roof is approximately 1.0 m in thickness and is made up of sandy mudstone and mudstone.The main roof consists of more competent siltstones and silty mudstones containing siderite nodules.In other words, the roof strata have low strength and poor stability and exhibit clear stratification, thus resulting in an immediate roof caving with the working face advancing as well as part of the main roof.Also, the floor is made up of sandy mudstone.It must be stressed that, at the working face end, the gateway soft roof was often affected by dynamic loading for the following reasons: (1) caving of the overburden rock mass in the gob roof, and the gangue movement, collision; (2) periodic weighting of the main roof; and (3) blasting operations of adjacent mining and excavating faces.

The Survey Results of Soft Roof Failure Characteristics
The soft roof failure characteristics of the front abutment pressure zone, the stress decreasing zone, and the stress stability zone are shown in Figure 2. According to a filed investigation of the head entry, there was a serious roof deformation due to the impact of front abutment pressure and a few cracks with less opening over a small area of the roof surface, as shown in Figure 2a.With the working face advancing, the roof overhanging length increased, and cracks formed and joined together in the roof rock mass, as shown in Figure 2b, resulting in a decrease of the roof strength, loss of self-bearing capacity, separation from the main roof, and roof deformation and becoming unload rock body, especially for the rock mass indicated by the red line.In addition, two bolts and one cable (B1, B2 and C1) on the gob side lost their support ability.After building the filling wall and removing the supports in the stress stability zone, the roof fractured along the inside of the filling wall, indicated as the "fracture line" in Figure 2c.
Field observations found that there were a large number of cracks in the vertical direction of roof near the end of the working face, so the roof caving occurred generally within a height range of 1.5 to 2 m, as shown in Figure 3.According to a geological survey, the gateway roof strata typically consist of weak rocks.The immediate roof is approximately 1.0 m in thickness and is made up of sandy mudstone and mudstone.The main roof consists of more competent siltstones and silty mudstones containing siderite nodules.In other words, the roof strata have low strength and poor stability and exhibit clear stratification, thus resulting in an immediate roof caving with the working face advancing as well as part of the main roof.Also, the floor is made up of sandy mudstone.It must be stressed that, at the working face end, the gateway soft roof was often affected by dynamic loading for the following reasons: (1) caving of the overburden rock mass in the gob roof, and the gangue movement, collision; (2) periodic weighting of the main roof; and (3) blasting operations of adjacent mining and excavating faces.

The Survey Results of Soft Roof Failure Characteristics
The soft roof failure characteristics of the front abutment pressure zone, the stress decreasing zone, and the stress stability zone are shown in Figure 2. According to a filed investigation of the head entry, there was a serious roof deformation due to the impact of front abutment pressure and a few cracks with less opening over a small area of the roof surface, as shown in Figure 2a.With the working face advancing, the roof overhanging length increased, and cracks formed and joined together in the roof rock mass, as shown in Figure 2b, resulting in a decrease of the roof strength, loss of self-bearing capacity, separation from the main roof, and roof deformation and becoming unload rock body, especially for the rock mass indicated by the red line.In addition, two bolts and one cable (B1, B2 and C1) on the gob side lost their support ability.After building the filling wall and removing the supports in the stress stability zone, the roof fractured along the inside of the filling wall, indicated as the "fracture line" in Figure 2c.
Field observations found that there were a large number of cracks in the vertical direction of roof near the end of the working face, so the roof caving occurred generally within a height range of 1.5 to 2 m, as shown in Figure 3.

Stress Evolution Law in Roof
To analyze the stress evolution law in the roof, the numerical modeling was adopted.

Numerical Simulation Model
Because of the coal seam condition and as the rock mass was discontinuous [18,19], 3DEC of ITASCA (Minneapolis, MN, USA) [20,21] was used to study the stress evolution law in the roof.The hexahedral model has a length, width, and height of 230 m, 100 m, and 200 m, respectively, includes coal seams and rock strata, with a total of 14 layers, as shown in Figure 4, in accordance with the geological conditions of the investigated mine.The Mohr-Coulomb yield criterion for the materials and the Coulomb slip model for contact were used.The mechanical and physical properties of all the layers and the contacts between every two layers are described in [13,22], respectively.
According to the reference [23], the minimum coefficient of lateral pressure is very close to 0 and the maximum can be up to 6 due to tectonic movement.For this case, we consider the coefficient of lateral pressure is 0.5, as the gateway is placed at shallow depth and close to the anticline axis and seam outcrop, which resulting in a much low horizontal stress.So, the state of the in situ stresses is σx = σy = 4.25 MPa and σz = 8.5 MPa, with σy parallel to the longwall advance direction and σx perpendicular, as shown in Figure 4b.A vertical pressure of 8.5 MPa is applied on the top surface, and the velocity of the bottom surface was restricted in all three directions, vx = vy = vz = 0 m/s.The velocity of the other four surfaces were restricted in the normal direction vn = 0 m/s.
The shape and size of the head entry are introduced in Section 2.1.After an initial equilibrium calculation, rock bolts were installed, as shown in Figure 1b, once the head was entry excavated.The rock bolts were represented as built-in "cable" elements.For resin-grouted rock bolts, the stiffness (Kbond) and the cohesive strength (Sbond) of the grout are the two key properties that govern the anchor characteristics [22]; Kbond = 3.06 × 10 9 N/m/m and Sbond = 2.3 × 10 5 kN/m were adopted in this study.In addition, a cross-sectional area of 3.142 × 10 −4 m 2 , an elastic modulus of 200 GPa, and a tensile yield

Stress Evolution Law in Roof
To analyze the stress evolution law in the roof, the numerical modeling was adopted.

Numerical Simulation Model
Because of the coal seam condition and as the rock mass was discontinuous [18,19], 3DEC of ITASCA (Minneapolis, MN, USA) [20,21] was used to study the stress evolution law in the roof.The hexahedral model has a length, width, and height of 230 m, 100 m, and 200 m, respectively, includes coal seams and rock strata, with a total of 14 layers, as shown in Figure 4, in accordance with the geological conditions of the investigated mine.The Mohr-Coulomb yield criterion for the materials and the Coulomb slip model for contact were used.The mechanical and physical properties of all the layers and the contacts between every two layers are described in [13,22], respectively.
According to the reference [23], the minimum coefficient of lateral pressure is very close to 0 and the maximum can be up to 6 due to tectonic movement.For this case, we consider the coefficient of lateral pressure is 0.5, as the gateway is placed at shallow depth and close to the anticline axis and seam outcrop, which resulting in a much low horizontal stress.So, the state of the in situ stresses is σx = σy = 4.25 MPa and σz = 8.5 MPa, with σy parallel to the longwall advance direction and σx perpendicular, as shown in Figure 4b.A vertical pressure of 8.5 MPa is applied on the top surface, and the velocity of the bottom surface was restricted in all three directions, vx = vy = vz = 0 m/s.The velocity of the other four surfaces were restricted in the normal direction vn = 0 m/s.
The shape and size of the head entry are introduced in Section 2.1.After an initial equilibrium calculation, rock bolts were installed, as shown in Figure 1b, once the head was entry excavated.The rock bolts were represented as built-in "cable" elements.For resin-grouted rock bolts, the stiffness (Kbond) and the cohesive strength (Sbond) of the grout are the two key properties that govern the anchor characteristics [22]; Kbond = 3.06 × 10 9 N/m/m and Sbond = 2.3 × 10 5 kN/m were adopted in this study.In addition, a cross-sectional area of 3.142 × 10 −4 m 2 , an elastic modulus of 200 GPa, and a tensile yield

Stress Evolution Law in Roof
To analyze the stress evolution law in the roof, the numerical modeling was adopted.

Numerical Simulation Model
Because of the coal seam condition and as the rock mass was discontinuous [18,19], 3DEC of ITASCA (Minneapolis, MN, USA) [20,21] was used to study the stress evolution law in the roof.The hexahedral model has a length, width, and height of 230 m, 100 m, and 200 m, respectively, includes coal seams and rock strata, with a total of 14 layers, as shown in Figure 4, in accordance with the geological conditions of the investigated mine.The Mohr-Coulomb yield criterion for the materials and the Coulomb slip model for contact were used.The mechanical and physical properties of all the layers and the contacts between every two layers are described in [13,22], respectively.
According to the reference [23], the minimum coefficient of lateral pressure is very close to 0 and the maximum can be up to 6 due to tectonic movement.For this case, we consider the coefficient of lateral pressure is 0.5, as the gateway is placed at shallow depth and close to the anticline axis and seam outcrop, which resulting in a much low horizontal stress.So, the state of the in situ stresses is σ x = σ y = 4.25 MPa and σ z = 8.5 MPa, with σ y parallel to the longwall advance direction and σ x perpendicular, as shown in Figure 4b.A vertical pressure of 8.5 MPa is applied on the top surface, and the velocity of the bottom surface was restricted in all three directions, v x = v y = v z = 0 m/s.The velocity of the other four surfaces were restricted in the normal direction v n = 0 m/s.
The shape and size of the head entry are introduced in Section 2.1.After an initial equilibrium calculation, rock bolts were installed, as shown in Figure 1b, once the head was entry excavated.The rock bolts were represented as built-in "cable" elements.For resin-grouted rock bolts, the stiffness (K bond ) and the cohesive strength (S bond ) of the grout are the two key properties that govern the anchor characteristics [22]; K bond = 3.06 ˆ10 9 N/m/m and S bond = 2.3 ˆ10 5 kN/m were adopted in this study.In addition, a cross-sectional area of 3.142 ˆ10 ´4 m 2 , an elastic modulus of 200 GPa, and a tensile yield strength of 165 kN were assigned to the "cable" element.A stepwise excavation in the y-direction was adopted to simulate the working face advancing by deleting blocks in five steps of 10 m each, as shown in Figure 4c.
strength of 165 kN were assigned to the "cable" element.A stepwise excavation in the y-direction was adopted to simulate the working face advancing by deleting blocks in five steps of 10 m each, as shown in Figure 4c.

Numerical Results
To understand the stress evolution law, the horizontal stresses in the x-direction of points A1 and A3 and the vertical stress of point A2 in the roof at y = 45 m, as shown in Figure 4d, were monitored during the working face advancing, and the results are shown in Figure 5. Especially, it can be found that from A to B Zone, the horizontal stress on the coal side (A3) increases slightly in magnitude and then reduces to approximately 5.9 MPa, while the horizontal stress on the gob side (A1) drops instantaneously and appears tensile, and the vertical stress at the gateway central (A2) reduces significantly to 4.5 MPa as soon as the working face advances beyond the monitor points.The stresses in the gateway roof changed over the whole process from the beginning of the caving to the roof stabilization, until a tensile state was reached, which will influence the stability of the gob-side entry retaining.

Force States of Roof Rock Mass
Here, combining with the numerical results and field investigation results, the force states of roof rock mass was analyzed to facilitate the failure mechanism analysis of roof rock mass and later theoretical analysis of supporting.
The stress evolution law showed in Figure 5 experienced the front abutment pressure zone (section A-A in Figure 1) and stress decreasing zone (section B-B in Figure 1).In section A-A, there are lateral forces F1a and F1b parallel to the rock strata and an overburden pressure P1 resembling the stress

Numerical Results
To understand the stress evolution law, the horizontal stresses in the x-direction of points A 1 and A 3 and the vertical stress of point A 2 in the roof at y = 45 m, as shown in Figure 4d, were monitored during the working face advancing, and the results are shown in Figure 5. Especially, it can be found that from A to B Zone, the horizontal stress on the coal side (A 3 ) increases slightly in magnitude and then reduces to approximately 5.9 MPa, while the horizontal stress on the gob side (A 1 ) drops instantaneously and appears tensile, and the vertical stress at the gateway central (A 2 ) reduces significantly to 4.5 MPa as soon as the working face advances beyond the monitor points.The stresses in the gateway roof changed over the whole process from the beginning of the caving to the roof stabilization, until a tensile state was reached, which will influence the stability of the gob-side entry retaining.
Minerals 2015, 5, page-page 5 strength of 165 kN were assigned to the "cable" element.A stepwise excavation in the y-direction was adopted to simulate the working face advancing by deleting blocks in five steps of 10 m each, as shown in Figure 4c.

Numerical Results
To understand the stress evolution law, the horizontal stresses in the x-direction of points A1 and A3 and the vertical stress of point A2 in the roof at y = 45 m, as shown in Figure 4d, were monitored during the working face advancing, and the results are shown in Figure 5. Especially, it can be found that from A to B Zone, the horizontal stress on the coal side (A3) increases slightly in magnitude and then reduces to approximately 5.9 MPa, while the horizontal stress on the gob side (A1) drops instantaneously and appears tensile, and the vertical stress at the gateway central (A2) reduces significantly to 4.5 MPa as soon as the working face advances beyond the monitor points.The stresses in the gateway roof changed over the whole process from the beginning of the caving to the roof stabilization, until a tensile state was reached, which will influence the stability of the gob-side entry retaining.

Force States of Roof Rock Mass
Here, combining with the numerical results and field investigation results, the force states of roof rock mass was analyzed to facilitate the failure mechanism analysis of roof rock mass and later theoretical analysis of supporting.
The stress evolution law showed in Figure 5 experienced the front abutment pressure zone (section A-A in Figure 1) and stress decreasing zone (section B-B in Figure 1).In section A-A, there are lateral forces F1a and F1b parallel to the rock strata and an overburden pressure P1 resembling the stress

Force States of Roof Rock Mass
Here, combining with the numerical results and field investigation results, the force states of roof rock mass was analyzed to facilitate the failure mechanism analysis of roof rock mass and later theoretical analysis of supporting.
The stress evolution law showed in Figure 5 experienced the front abutment pressure zone (section A-A in Figure 1) and stress decreasing zone (section B-B in Figure 1).In section A-A, there are lateral forces F 1a and F 1b parallel to the rock strata and an overburden pressure P 1 resembling the stress (A 1 , A 2 and A 3 ) in A Zone in Figure 5.In addition, field investigation showed that, there are bolt, mesh, and cable support force T 1 , reinforced supporting force R 1 in advance, shear forces Q 1a , Q 1b on the two sides and gravity G.The force state was shown in Figure 6a.
In section B-B, the gob side lateral force disappears along the strata strike direction resembling the stress (A 1 ) in B Zone in Figure 5, because of the roof rock mass caving near the gob side and part of the roof rock mass fracturing, but the coal side lateral force is F 2b resembling the stress (A 3 ) in B Zone in Figure 5. On the other hand, the overburden of pressure is assumed to become 0 resembling the stress (A 2 ) in B Zone in Figure 5, due to the roof separating from main roof.There are also a bolt, mesh, and cable support force T 2 , a strengthen support force R 2 at the working face end and a shear force Q 2b in field, also the gravity G.The force state was shown in Figure 6b.
Furthermore, the roof rock mass stress in cross section C-C (in Figure 1) has experienced the stress states of sections A-A and B-B, but, field investigation showed that it will change after the construction of the artificial filling wall and the main roof rotary sinking, and the lateral force F a on the coal side and the overburden pressure P will restore, the lateral force on the gob side will change to F b .In addition, the bolt, mesh and cable support force become to T, the shear forces at the two ends become to Q 1a and Q 1b and gravity remain the same G.The force state was shown in Figure 6c.
Minerals 2015, 5, page-page 6 (A1, A2 and A3) in A Zone in Figure 5.In addition, field investigation showed that, there are bolt, mesh, and cable support force T1, reinforced supporting force R1 in advance, shear forces Q1a, Q1b on the two sides and gravity G.The force state was shown in Figure 6a.
In section B-B, the gob side lateral force disappears along the strata strike direction resembling the stress (A1) in B Zone in Figure 5, because of the roof rock mass caving near the gob side and part of the roof rock mass fracturing, but the coal side lateral force is F2b resembling the stress (A3) in B Zone in Figure 5. On the other hand, the overburden of pressure is assumed to become 0 resembling the stress (A2) in B Zone in Figure 5, due to the roof separating from main roof.There are also a bolt, mesh, and cable support force T2, a strengthen support force R2 at the working face end and a shear force Q2b in field, also the gravity G.The force state was shown in Figure 6b.
Furthermore, the roof rock mass stress in cross section C-C (in Figure 1) has experienced the stress states of sections A-A and B-B, but, field investigation showed that it will change after the construction of the artificial filling wall and the main roof rotary sinking, and the lateral force Fa on the coal side and the overburden pressure P will restore, the lateral force on the gob side will change to Fb.In addition, the bolt, mesh and cable support force become to T, the shear forces at the two ends become to Q1a and Q1b and gravity remain the same G.The force state was shown in Figure 6c.

Mechanism for Failure of Roof Rock Mass
Roof failure was affected by many factors, so we will combine them with the failure characteristics obtained by filed investigation, stress evolution law obtained by numerical modeling and the force states to study the roof failure mechanism.

Unloading Effect of the Lateral Stress
Affected by the abutment pressure in the mining and excavation process, the roof rock mass is in the yield state, as shown in cross section A-A in Figure 1.The pre-tensioned bolt-cable-mesh supporting system improved the strength, significantly increased the yield strength, and altered the deformation characteristics of the roof rock mass.At the same time, the support system exerted a pressure stress on the rock mass, so the compressive zone stress state had to be altered, which can offset some of the tensile stress and friction and enhance the shear capacity.In addition, the axial and lateral anchored forces increased the shear strength of the weak structural plane, preventing the roof rock mass from moving and sliding along the block structure plane.The pre-tensioned support system controlled the expansion deformation and destruction, preventing roof separation, sliding, fracture opening, and new crack generation in the anchorage zones, not only maintaining the integrity but also forming a pre-tensioned bearing structure with a large stiffness [24,25].

Mechanism for Failure of Roof Rock Mass
Roof failure was affected by many factors, so we will combine them with the failure characteristics obtained by filed investigation, stress evolution law obtained by numerical modeling and the force states to study the roof failure mechanism.

Unloading Effect of the Lateral Stress
Affected by the abutment pressure in the mining and excavation process, the roof rock mass is in the yield state, as shown in cross section A-A in Figure 1.The pre-tensioned bolt-cable-mesh supporting system improved the strength, significantly increased the yield strength, and altered the deformation characteristics of the roof rock mass.At the same time, the support system exerted a pressure stress on the rock mass, so the compressive zone stress state had to be altered, which can offset some of the tensile stress and friction and enhance the shear capacity.In addition, the axial and lateral anchored forces increased the shear strength of the weak structural plane, preventing the roof rock mass from moving and sliding along the block structure plane.The pre-tensioned support system controlled the expansion deformation and destruction, preventing roof separation, sliding, fracture opening, and new crack generation in the anchorage zones, not only maintaining the integrity but also forming a pre-tensioned bearing structure with a large stiffness [24,25].
Influenced by the rotary sinking of the overlying roof and the superposition of the overburden pressure, the vertical displacement of the anchor bolt and cable increased significantly, leading to an increase of T 1 .The stress increment of the roadway roof rock mass mainly comes from the bulking deformation.When the bolt, mesh, and cable apply pressure on the broken rock mass and adds an anchoring force, the rock mass volume or volume rise rate decrease, and a bulking force is produced by the broken rock mass that would work the bolts and surrounding rock at the same time and put the roof rock mass in an extrusion state.As shown in cross section A-A in Figure 1, it is precisely because of the interaction of the "support-surrounding rock", that the roof rock mass presented a failure state in Figure 2a.
Although the single hydraulic props support the roof strata with pressure, the extrusion state rock mass above the hydraulic support started loosening on the working face end.In cross section B-B in Figure 1, the loosening roof rock mass caved behind the support body on the gob side, causing a sharp reduction of the lateral pressure, especially for the accumulated bulking force.Meanwhile, the overburden pressure decreased almost to 0 due to the separation of the soft roof from the hard roof above.
It can be observed that, accompanied by the gob roof caving, the lateral pressure unloaded and the bulking force decreased in the broken rock mass of the roadway roof, leading to its stress state changing from a three-dimensional stress state to a two-dimensional or uniaxial stress state, leading more easily to failure.At this point, together with the tensile stress effect, the rock mass volume "elastic" expansion in the layer strike direction caused the generation and propagation of a large number of cracks and the further failure of the rock mass.

Tensile and Shear Failure of the Working Face End Lower Layer Rock Mass
Through the field investigation of the roof support body layout and its work process and based on the features of the roof rock mass deformation and failure, it is found that the working face end soft roof presents the following phenomenon: (1) An uneven supporting at the working face end roof.A strong mine ground pressure appeared, causing the significant sinking of the working face end roof in the working face advancing process, so a single hydraulic prop with an articulated roof beam was usually set as a reinforcement support.There is uneven pressure on the roof surface due to the low strength and stiffness of the sandy mudstone and mudstone and the higher strength and stiffness strengthening the support body, resulting in part of the support body inserting the roof and leading to a roof rock shear failure with shear stress q, as shown in Figure 7.In addition, the passive supporting force is too large and does not couple support with the bolts (cables), causing part of the bolts (cables) to be loosened, affecting the roof rock mass local stress state and resulting in ultimate rock mass failure.(2) Effect of local moment.Treating the reinforcement body (individual hydraulic prop and articulated roof beam, etc.) as the fulcrum, the lower roof strata produced a local bending on a small scale, but when the fulcrum bending moment (M) is too large, the layered roof strata cause tensile and compressive damage to the upper surface and lower surface, respectively, additionally increased the possibility of the support body inserting the upper roof and leading to rock mass shear failure, as shown in Figure 7.
The two phenomena above are shown in the Figure 1 A-A and B-B cross sections and may occur simultaneously.

Dynamic Failure
The working face end roof rock mass is often affected by vibration as described above.Here, 3DEC was also applied to study the dynamic failure characteristics of the roof rock mass, the coal and rock mass physical and mechanical parameters and the boundary conditions used in the numerical model shown in Figure 4, additionally adding the model viscous boundary conditions for all boundaries to make the stress wave propagate or be absorbed to simulate the infinite foundation environment.However, we alter the length to 1.6 m in the y direction and treat it as a plane strain model.As the location in cross section A-A in Figure 1, under the action of a triangle stress wave to simulate the influence of a dynamic disturbance to the roof rock mass.In the simulation process, the stress peak value of the triangle stress wave is 8 MPa rise time is 1 ms and the decrease time is 7 ms [26][27][28][29].Roof bolt (like a 1 , a 2 , a 3 , a 4 , a 5 , a 6 in Figure 4c) end displacement and the surrounding rock plastic zone were monitored, and the results are shown in Figures 8 and 9 respectively.Minerals 2015, 5, page-page 8 environment.However, we alter the length to 1.6 m in the y direction and treat it as a plane strain model.As the location in cross section A-A in Figure 1, under the action of a triangle stress wave to simulate the influence of a dynamic disturbance to the roof rock mass.In the simulation process, the stress peak value of the triangle stress wave is 8 MPa rise time is 1 ms and the decrease time is 7 ms [26][27][28][29].Roof bolt (like a1, a2, a3, a4, a5, a6 in Figure 4c) end displacement and the surrounding rock plastic zone were monitored, and the results are shown in Figures 8 and 9, respectively.It can be observed from Figure 8 that the displacement of all the monitoring points had a small increase in magnitude after the disturbance.
Figure 9a,b show the plastic distribution zone before and after the dynamic disturbance of the gateway surrounding rock mass, respectively.The stress wave had countless reflections and refractions on the roof crack surface according to the stress wave propagation theory, and the tensile stress would environment.However, we alter the length to 1.6 m in the y direction and treat it as a plane strain model.As the location in cross section A-A in Figure 1, under the action of a triangle stress wave to simulate the influence of a dynamic disturbance to the roof rock mass.In the simulation process, the stress peak value of the triangle stress wave is 8 MPa rise time is 1 ms and the decrease time is 7 ms [26][27][28][29].Roof bolt (like a1, a2, a3, a4, a5, a6 in Figure 4c) end displacement and the surrounding rock plastic zone were monitored, and the results are shown in Figures 8 and 9, respectively.It can be observed from Figure 8 that the displacement of all the monitoring points had a small increase in magnitude after the disturbance.
Figure 9a,b show the plastic distribution zone before and after the dynamic disturbance of the gateway surrounding rock mass, respectively.The stress wave had countless reflections and refractions on the roof crack surface according to the stress wave propagation theory, and the tensile stress would environment.However, we alter the length to 1.6 m in the y direction and treat it as a plane strain model.As the location in cross section A-A in Figure 1, under the action of a triangle stress wave to simulate the influence of a dynamic disturbance to the roof rock mass.In the simulation process, the stress peak value of the triangle stress wave is 8 MPa rise time is 1 ms and the decrease time is 7 ms [26][27][28][29].Roof bolt (like a1, a2, a3, a4, a5, a6 in Figure 4c) end displacement and the surrounding rock plastic zone were monitored, and the results are shown in Figures 8 and 9, respectively.It can be observed from Figure 8 that the displacement of all the monitoring points had a small increase in magnitude after the disturbance.
Figure 9a,b show the plastic distribution zone before and after the dynamic disturbance of the gateway surrounding rock mass, respectively.The stress wave had countless reflections and refractions on the roof crack surface according to the stress wave propagation theory, and the tensile stress would It can be observed from Figure 8 that the displacement of all the monitoring points had a small increase in magnitude after the disturbance.
Figure 9a,b show the plastic distribution zone before and after the dynamic disturbance of the gateway surrounding rock mass, respectively.The stress wave had countless reflections and refractions on the roof crack surface according to the stress wave propagation theory, and the tensile stress would cause changes to the roof rock mass stress state and result in a combination of shear failure on existing joints/weakness horizons, an extension of critically oriented joints and propagation of new fractures through previously intact rock so that the integrity of the roof became poor, the plastic zone expanding into the higher strata, and the ultimate widespread destruction of the roof rock mass [7].At the same time, the two gateway sides and floor rock mass also produced plastic damage under the stress wave action, but had small extended failure zones.
The above results show that the increment of the displacement and plastic zone was small, but due to the constant increase of the dynamic load, it can be presumed that the working face end soft roof rock deformation and failure magnitude would significantly increase under the repeated dynamic load action, and the more cracks are formed, the more damage caused to the rock mass.To minimize the power damage, there is a need to improve the soft roof rock mass stress state and reduce the development degree of cracks.

Mechanism for Failure of the Roof Rock Mass in Retaining Roadway
As shown in Figure 2c, cross section C-C in the stabilized zone in retaining roadway experience a fracture deformation process at the working face end that can be influenced by the main roof rotary deformation and stress recovery process beside the working face, resulting in its larger damage and deformation degree, causing a weakening of its integrity.After strengthening the support of the roof rock mass in the retaining roadway and constructing the artificial filling wall body, their integrity and stability were improved and the bulking force was restored to a certain degree.For a certain retaining distance, after the strengthened support body was recycled, due to the artificial filling wall's higher strength and stiffness, the bulking force would release along its inside.Therefore, the soft roof of retained entry caved at the man-made constructed wall (see shear failure line in Figure 10) due to forces such as the overburden pressure (P), the bulking stress (P s ), rock mass gravity (G) and the bolt, mesh, and cable support force (T), which formed a "similar cantilever beam" structure, as shown in Figure 10.
Minerals 2015, 5, page-page 9 cause changes to the roof rock mass stress state and result in a combination of shear failure on existing joints/weakness horizons, an extension of critically oriented joints and propagation of new fractures through previously intact rock so that the integrity of the roof became poor, the plastic zone expanding into the higher strata, and the ultimate widespread destruction of the roof rock mass [7].At the same time, the two gateway sides and floor rock mass also produced plastic damage under the stress wave action, but had small extended failure zones.
The above results show that the increment of the displacement and plastic zone was small, but due to the constant increase of the dynamic load, it can be presumed that the working face end soft roof rock deformation and failure magnitude would significantly increase under the repeated dynamic load action, and the more cracks are formed, the more damage caused to the rock mass.To minimize the power damage, there is a need to improve the soft roof rock mass stress state and reduce the development degree of cracks.

Mechanism for Failure of the Roof Rock Mass in Retaining Roadway
As shown in Figure 2c, cross section C-C in the stabilized zone in retaining roadway experience a fracture deformation process at the working face end that can be influenced by the main roof rotary deformation and stress recovery process beside the working face, resulting in its larger damage and deformation degree, causing a weakening of its integrity.After strengthening the support of the roof rock mass in the retaining roadway and constructing the artificial filling wall body, their integrity and stability were improved and the bulking force was restored to a certain degree.For a certain retaining distance, after the strengthened support body was recycled, due to the artificial filling wall's higher strength and stiffness, the bulking force would release along its inside.Therefore, the soft roof of retained entry caved at the man-made constructed wall (see shear failure line in Figure 10) due to forces such as the overburden pressure (P), the bulking stress (Ps), rock mass gravity (G) and the bolt, mesh, and cable support force (T), which formed a "similar cantilever beam" structure, as shown in Figure 10.

Roof Support Countermeasures
To develop a supporting method for the roof rock mass at the working face end and retained gateway, the deformation and force in the roof rock mass must be analyzed.The rock bolt is the basic supporting tool.First, we calculated the "cable" element to represent the rock bolts on the working face end with several different bolt installation angles in 3DEC model to find the bolt end displacement change characteristics.Then, we analyzed the bolt limit equilibrium tension force change features in the retained gateway roof rock mass with the change of the bolt installation angle though the force balance equation.

Roof Support Countermeasures
To develop a supporting method for the roof rock mass at the working face end and retained gateway, the deformation and force in the roof rock mass must be analyzed.The rock bolt is the basic supporting tool.First, we calculated the "cable" element to represent the rock bolts on the working face end with several different bolt installation angles in 3DEC model to find the bolt end displacement change characteristics.Then, we analyzed the bolt limit equilibrium tension force change features in the retained gateway roof rock mass with the change of the bolt installation angle though the force balance equation.

Deformation Analysis of Working Face End Roof
From the soft roof failure process and failure characteristics, the following mechanism for the soft roof support can be obtained: (1) Providing the roof rock mass extrusion stress and changing the rock mass stress state to improve their strength [11]; (2) Preventing cracks from generating and propagating to increase both the strength of the affected roof strata and the stiffness of the whole bolted strata to reduce roof deformation and dynamic damage; and (3) Reducing the broken rock mass bulking force and maintaining the stability of the roadway roof by supporting in time.According to the above points and the "bolt extruding reinforcement theory" by pre-tensioned bolts providing a compressive zone in the axial direction [11], the working face end soft roof mechanical structure should be similar to that in Figure 11 with bolts supporting, so the installation angle is suggested to be 0 ˝< α < 90 The distributed axial stress of the pre-tensioned bolts in the extension direction and the lateral stress vertical axial direction can alter the rock stress state, especially for the thin and weak immediate roof conditions [13].Here, we use the same numerical model and the physical and mechanical parameters of the rock mass in Figure 4 to analyze the effect of the bolt installation angle α on the stability of the roof rock mass.Bolt installation angles of 50 ˝, 60 ˝, 70 ˝, and 80 ˝(in Figure 12) were applied in the numerical model and the monitoring bolt end displacement was at y = 45 m with the same excavation process as shown in Figure 4 conducted, and the results are shown in Figure 13.As actual field bolt installation at an angle α in the mine roof close to 90 ˝, similar to in Figure 1b, it was taken as 90 ˝to facilitate the analysis, and the results are shown in Section 3.1.

Deformation Analysis of Working Face End Roof
From the soft roof failure process and failure characteristics, the following mechanism for the soft roof support can be obtained: (1) Providing the roof rock mass extrusion stress and changing the rock mass stress state to improve their strength [11]; (2) Preventing cracks from generating and propagating to increase both the strength of the affected roof strata and the stiffness of the whole bolted strata to reduce roof deformation and dynamic damage; and (3) Reducing the broken rock mass bulking force and maintaining the stability of the roadway roof by supporting in time.According to the above points and the "bolt extruding reinforcement theory" by pre-tensioned bolts providing a compressive zone in the axial direction [11], the working face end soft roof mechanical structure should be similar to that in Figure 11 with bolts supporting, so the installation angle is suggested to be 0° < α < 90°.
The distributed axial stress of the pre-tensioned bolts in the extension direction and the lateral stress vertical axial direction can alter the rock stress state, especially for the thin and weak immediate roof conditions [13].Here, we use the same numerical model and the physical and mechanical parameters of the rock mass in Figure 4 to analyze the effect of the bolt installation angle α on the stability of the roof rock mass.Bolt installation angles of 50°, 60°, 70°, and 80° (in Figure 12) were applied in the numerical model and the monitoring bolt end displacement was at y = 45 m with the same excavation process as shown in Figure 4 conducted, and the results are shown in Figure 13.As actual field bolt installation at an angle α in the mine roof close to 90°, similar to in Figure 1b, it was taken as 90° to facilitate the analysis, and the results are shown in Section 3.1.Figure 13 shows that: (1) The monitored displacement had a good regularity and the displacement curve showed a trend of a "concave" type with α increasing as a whole; the vertical displacement of roof was minimal when α = 60°, and the displacement increased when α = 70° and 80°, but it had a smaller increment; (2) For a certain α, the point displacement increased as the distance to the gob decreased, that is to say, the roof rotated and sank; and (3) The displacement is larger under the current anchor installation angle (α = 90°) condition; its

Deformation Analysis of Working Face End Roof
From the soft roof failure process and failure characteristics, the following mechanism for the soft roof support can be obtained: (1) Providing the roof rock mass extrusion stress and changing the rock mass stress state to improve their strength [11]; (2) Preventing cracks from generating and propagating to increase both the strength of the affected roof strata and the stiffness of the whole bolted strata to reduce roof deformation and dynamic damage; and (3) Reducing the broken rock mass bulking force and maintaining the stability of the roadway roof by supporting in time.According to the above points and the "bolt extruding reinforcement theory" by pre-tensioned bolts providing a compressive zone in the axial direction [11], the working face end soft roof mechanical structure should be similar to that in Figure 11 with bolts supporting, so the installation angle is suggested to be 0° < α < 90°.
The distributed axial stress of the pre-tensioned bolts in the extension direction and the lateral stress vertical axial direction can alter the rock stress state, especially for the thin and weak immediate roof conditions [13].Here, we use the same numerical model and the physical and mechanical parameters of the rock mass in Figure 4 to analyze the effect of the bolt installation angle α on the stability of the roof rock mass.Bolt installation angles of 50°, 60°, 70°, and 80° (in Figure 12) were applied in the numerical model and the monitoring bolt end displacement was at y = 45 m with the same excavation process as shown in Figure 4 conducted, and the results are shown in Figure 13.As actual field bolt installation at an angle α in the mine roof close to 90°, similar to in Figure 1b, it was taken as 90° to facilitate the analysis, and the results are shown in Section 3.1.Figure 13 shows that: (1) The monitored displacement had a good regularity and the displacement curve showed a trend of a "concave" type with α increasing as a whole; the vertical displacement of roof was minimal when α = 60°, and the displacement increased when α = 70° and 80°, but it had a smaller increment; (2) For a certain α, the point displacement increased as the distance to the gob decreased, that is to say, the roof rotated and sank; and (3) The displacement is larger under the current anchor installation angle (α = 90°) condition; its Figure 13 shows that: (1) The monitored displacement had a good regularity and the displacement curve showed a trend of a "concave" type with α increasing as a whole; the vertical displacement of roof was minimal when α = 60 ˝, and the displacement increased when α = 70 ˝and 80 ˝, but it had a smaller increment; (2) For a certain α, the point displacement increased as the distance to the gob decreased, that is to say, the roof rotated and sank; and (3) The displacement is larger under the current anchor installation angle (α = 90 ˝) condition; its differences were 155 and 321 mm from the minimum displacement.

Bolt Limit Equilibrium Tension Force
It is known that the working principle of bolts is to maintain the roof rock mass stability in the early stage and to control roof deformation in the later stage.Moreover, carrying out gob-side entry retaining is a systematic project and the roof deformation in the later retaining stage, when maintaining the working face end roof rock mass stability, should be considered.After the construction of the artificial filling wall and recycled reinforcement support body, the roof shear failure (see Figure 2c) is mainly influenced by the bolt-cable-mesh force T, gravity G, shear forces Qa and Qb, and the lateral forces Fa and Fb (reference in Figure 6c), regardless of the overburden pressure because of the roof separation, as shown in Figure 14.
At the moment, the roadway roof rock mass undergoes serious failure and its integrity is poor, so it loses its rigid body properties.Therefore, there is distance h1 between rotating point A and the upper endpoint, so it cannot have a certainty value under the action of torque and 0 ≤ h1 ≤ h.On the other hand, the roof had been cut along the edge of the gob near the artificial filling wall, so the lateral pressure and shear force can be ignored.In the simplified calculation, it can be set that Fa = 0, Qa = 0.The stress balance relationship in the layer strike direction and vertical layer direction are established, as well as the moment balance relationship established by point A, as shown in Equations ( 1)-( 3), respectively.

Bolt Limit Equilibrium Tension Force
It is known that the working principle of bolts is to maintain the roof rock mass stability in the early stage and to control roof deformation in the later stage.Moreover, carrying out gob-side entry retaining is a systematic project and the roof deformation in the later retaining stage, when maintaining the working face end roof rock mass stability, should be considered.After the construction of the artificial filling wall and recycled reinforcement support body, the roof shear failure (see Figure 2c) is mainly influenced by the bolt-cable-mesh force T, gravity G, shear forces Q a and Q b , and the lateral forces F a and F b (reference in Figure 6c), regardless of the overburden pressure because of the roof separation, as shown in Figure 14.
At the moment, the roadway roof rock mass undergoes serious failure and its integrity is poor, so it loses its rigid body properties.Therefore, there is distance h 1 between rotating point A and the upper endpoint, so it cannot have a certainty value under the action of torque and 0 ď h 1 ď h.

Bolt Limit Equilibrium Tension Force
It is known that the working principle of bolts is to maintain the roof rock mass stability in the early stage and to control roof deformation in the later stage.Moreover, carrying out gob-side entry retaining is a systematic project and the roof deformation in the later retaining stage, when maintaining the working face end roof rock mass stability, should be considered.After the construction of the artificial filling wall and recycled reinforcement support body, the roof shear failure (see Figure 2c) is mainly influenced by the bolt-cable-mesh force T, gravity G, shear forces Qa and Qb, and the lateral forces Fa and Fb (reference in Figure 6c), regardless of the overburden pressure because of the roof separation, as shown in Figure 14.
At the moment, the roadway roof rock mass undergoes serious failure and its integrity is poor, so it loses its rigid body properties.Therefore, there is distance h1 between rotating point A and the upper endpoint, so it cannot have a certainty value under the action of torque and 0 ≤ h1 ≤ h.On the other hand, the roof had been cut along the edge of the gob near the artificial filling wall, so the lateral pressure and shear force can be ignored.In the simplified calculation, it can be set that Fa = 0, Qa = 0.The stress balance relationship in the layer strike direction and vertical layer direction are established, as well as the moment balance relationship established by point A, as shown in Equations ( 1)-( 3), respectively.On the other hand, the roof had been cut along the edge of the gob near the artificial filling wall, so the lateral pressure and shear force can be ignored.In the simplified calculation, it can be (2) When α takes the values of 50 ˝, 60 ˝, 70 ˝, 80 ˝, and 90 ˝, the relation curve of the bolt tension T under roof rock mass ultimate stability conditions with an installation angle α is shown in Figure 15a, and the change trend is similar to the curve in Figure 12.
Here, it can be explained that when α is kept constant, the smaller the T value, the smaller the force required for the roof rock mass limit equilibrium under the action of the torque; when the actual bolt tension values remain unchanged, the roof rock mass have a minimum displacement, which requires a minimal bolt tension force (T min ) for limit equilibrium.For example, in Figure 15a, the largest displacement or the worst stability is when α = 50 and the minimum displacement or the best stability is when α = 70 ˝of the roof rock mass.Hence, the roof displacement is proportional to T. Combining Figures 13 and 15a, it can be observed that when α = 70°, it can provide a failure rock mass with extrusion pressure and alter the stress state, and it can also compress joints and fissures to reduce its opening and the dynamic damage at the same time.Hence, it is advantageous to control the working face end roof and maintain the roof stability in the gob-side entry retaining.

Discussion
The thickness of the soft roof stratum in a gently inclined coal seam is small and its strength is low, so mining activities can cause the roof rock mass to serious break and extremely easily cave in, producing a great threat to normal production activities and personnel safety.The selection of a pre-tensioned anchoring technique and bolt arrangement style can improve the bulking roof rock mass stress state by preventing cracks from expanding and reducing the roof separation.To meet the requirements of the retained head entry, a bolt arrangement form is proposed.Furthermore, the application of Equation ( 6) can be used for the discussion of the supporting measures under similar coal seam geological conditions of, for instance, different coal seam dip angles, gateway widths, and soft roof thicknesses.Because researchers provided more insights into the soft roof failure mechanism of horizontal coal seam by retaining gateway to set α = 90°, thereby assuming the caving height h = 2 m, coal seam dip angle θ = 0°, other conditions remain unchanged as above, thereby determining the proper relationship concerning α and T, as shown in Figure 15b.It shows that T is much large when α = 90° under the high thickness soft roof and horizontal coal seam condition.
On the one hand, the lateral shear forces of Qa and Qb in the theoretical calculation process of Equation ( 7) are not considered here, but the calculation results and the numerical simulation results are very consistent.Hence, it is feasible to treat the theoretical calculation results as a bolt support reference in the field.On the other hand, only the bolts support is mimicked in the numerical simulation, without incorporating the action of the individual hydraulic props, anchor mesh, and anchor cable, implying that the findings present here are conservative.
In addition, being restricted by current existing technology, it may be difficult to apply this theory practically, as there will be a high degree of drilling difficulty to set the anchor installation angle to α = 70° near the gob edge, as recommended by this paper.However, bolts with a theoretical value α can be installed on the coal side, and for the gob side bolts, α can be set as close as possible Combining Figures 13 and 15a, it can be observed that when α = 70 ˝, it can provide a failure rock mass with extrusion pressure and alter the stress state, and it can also compress joints and fissures to reduce its opening and the dynamic damage at the same time.Hence, it is advantageous to control the working face end roof and maintain the roof stability in the gob-side entry retaining.

Discussion
The thickness of the soft roof stratum in a gently inclined coal seam is small and its strength is low, so mining activities can cause the roof rock mass to serious break and extremely easily cave in, producing a great threat to normal production activities and personnel safety.The selection of a pre-tensioned anchoring technique and bolt arrangement style can improve the bulking roof rock mass stress state by preventing cracks from expanding and reducing the roof separation.To meet the requirements of the retained head entry, a bolt arrangement form is proposed.Furthermore, the application of Equation ( 6) can be used for the discussion of the supporting measures under similar coal seam geological conditions of, for instance, different coal seam dip angles, gateway widths, and soft roof thicknesses.Because researchers provided more insights into the soft roof failure mechanism of horizontal coal seam by retaining gateway to set α = 90 ˝, thereby assuming the caving height h = 2 m, coal seam dip angle θ = 0 ˝, other conditions remain unchanged as above, thereby determining the proper relationship concerning α and T, as shown in Figure 15b.It shows that T is much large when α = 90 ˝under the high thickness soft roof and horizontal coal seam condition.
On the one hand, the lateral shear forces of Q a and Q b in the theoretical calculation process of Equation ( 7) are not considered here, but the calculation results and the numerical simulation results are very consistent.Hence, it is feasible to treat the theoretical calculation results as a bolt support reference in the field.On the other hand, only the bolts support is mimicked in the numerical simulation, without incorporating the action of the individual hydraulic props, anchor mesh, and anchor cable, implying that the findings present here are conservative.
In addition, being restricted by current existing technology, it may be difficult to apply this theory practically, as there will be a high degree of drilling difficulty to set the anchor installation angle to α = 70 ˝near the gob edge, as recommended by this paper.However, bolts with a theoretical value α can be installed on the coal side, and for the gob side bolts, α can be set as close as possible to the theoretical value so that they will not be influenced by the roof caving in near the gob, as shown in Figure 2, and to make the most use of the bolts resource by improving their force.The point of the roof failure angle decreased as the horizontal stress level increased, indicating that failure tends to occur around the entry corners when the horizontal stress was low [30].However, the corner bolts (like e and f in Figure 12) that are installed tilted to the coal side can increase the horizontal stress, so they can prevent shear failure around the entry corners very well.

Conclusions
This paper analyzed the roof rock mass failure characteristics and their failure mechanism in the working face end and the gob-side entry retaining, and discussed its support technology according to the special gently inclined seam occurrence conditions.The following were reached from the analysis process: (1) After the working face advances, it is found that the horizontal stress of the soft rock mass at the working face end does not exhibit a large magnitude unloading on the coal wall side, but the horizontal stress momentarily fell, and a tensile stress appeared on the gob side.The vertical stress in the gateway central dropped significantly, almost down to zero.(2) The sinking and separation of the soft roof rock mass in the gently inclined coal seam working face end is affected by the front abutment pressure and the hanging roof on gob side.The initiation and propagation of cracks and the fractures of the rock mass are produced by the actions of the lateral stress unloading loose, tensile, and shear stresses in the low layer caused by uneven support and no coupling support and dynamic disturbances.(3) The roof rock mass failed in shear mode along the inside of man-made constructed wall in the stability zone of the retained gateway, due to the overburden pressure, bulking force, roof gravity, and combined supporting force.The failed roof forms a "similar cantilever beam" structure.(4) The equation of the bolt ultimate equilibrium tension force, a function of the seam inclination, gateway width, soft roof thickness, and bolt installation angle, was established according to the stress balance analysis of the roof rock.(5) To prevent the working face end soft roof rock mass from increasing its deformation and becoming significantly fractured, and also to maintain the gob-side entry retaining roof's stability, it is suggested that the gateway roof bolt installation angle be 70 ˝to provide an extrusion stress, change the rock mass stress state, and improve their strength for better entry maintenance.

Minerals 2015, 5 ,
page-page used, and three rows of single hydraulic supports with top hinged beams were used as advanced strengthening support, as shown in Figure1b.The roof rock bolts are non-fully anchored (anchor length is 1.4 m) with a pre-tension of 90 kN, inter-row spacing of 800 × 800 mm and bolt angles of 90°, 80°, 70°, and 60° for different bolts, and the extra cables have a pre-tension of 200 kN and an interrow spacing of 3200 × 1200 mm.

Figure 2 .
Figure 2. Schematic of soft roof rock mass failure; (a) In front abutment pressure zone; (b) In working face end; (c) In stress stability zone of retaining gateway.

Figure 3 .
Figure 3. Caving roof near the working face end.

Figure 2 . 4 Figure 2 .
Figure 2. Schematic of soft roof rock mass failure; (a) In front abutment pressure zone; (b) In working face end; (c) In stress stability zone of retaining gateway.

Figure 3 .
Figure 3. Caving roof near the working face end.

Figure 3 .
Figure 3. Caving roof near the working face end.

Figure 4 .
Figure 4. (a) Numerical model; (b) Coal seam mining model; (c) Working face advancing; (d) Immediate roof grid, bolt and monitoring points; (e) Gateway cross section of model.

Figure 5 .
Figure 5. Roof rock mass stress evolution law in working face end.

Figure 4 .
Figure 4. (a) Numerical model; (b) Coal seam mining model; (c) Working face advancing; (d) Immediate roof grid, bolt and monitoring points; (e) Gateway cross section of model.

Figure 4 .
Figure 4. (a) Numerical model; (b) Coal seam mining model; (c) Working face advancing; (d) Immediate roof grid, bolt and monitoring points; (e) Gateway cross section of model.

Figure 5 .
Figure 5. Roof rock mass stress evolution law in working face end.

Figure 5 .
Figure 5. Roof rock mass stress evolution law in working face end.

Figure 6 .
Figure 6.Roof rock mass force; (a) Force of cross section A-A; (b) Force of cross section B-B; (c) Force of cross section C-C.

Figure 6 .
Figure 6.Roof rock mass force; (a) Force of cross section A-A; (b) Force of cross section B-B; (c) Force of cross section C-C.

Figure 7 .
Figure 7. Schematic of support body damage to the roof rock mass.

Figure 8 .
Figure 8. Displacement variation before and after the dynamic disturbance.

Figure 9 .
Figure 9. Plastic zone of surrounding rock: (a) Plastic zone before disturbance; (b) Plastic zone after disturbance.

Figure 7 .
Figure 7. Schematic of support body damage to the roof rock mass.

Figure 7 .
Figure 7. Schematic of support body damage to the roof rock mass.

Figure 8 .
Figure 8. Displacement variation before and after the dynamic disturbance.

Figure 9 .
Figure 9. Plastic zone of surrounding rock: (a) Plastic zone before disturbance; (b) Plastic zone after disturbance.

Figure 8 .
Figure 8. Displacement variation before and after the dynamic disturbance.

Figure 7 .
Figure 7. Schematic of support body damage to the roof rock mass.

Figure 8 .
Figure 8. Displacement variation before and after the dynamic disturbance.

Figure 9 .
Figure 9. Plastic zone of surrounding rock: (a) Plastic zone before disturbance; (b) Plastic zone after disturbance.

Figure 9 .
Figure 9. Plastic zone of surrounding rock: (a) Plastic zone before disturbance; (b) Plastic zone after disturbance.

Figure 10 .
Figure 10.Roof rock mass stress model of the retaining roadway.

Figure 10 .
Figure 10.Roof rock mass stress model of the retaining roadway.

Figure 11 .
Figure 11.Stress state diagram of soft top supporting effect.

Figure
Figure Stress state diagram of soft top supporting effect.

Figure 11 .
Figure 11.Stress state diagram of soft top supporting effect.