Abstract
To control the roof during gob-side entry retaining by roof cutting in inclined coal seams, the retained gob-side roadway is zoned based on the mechanical principle and technological process of no-pillar mining with gob-entry retention. A simplified mechanical model for surrounding rocks in different subzones was established by using theoretical analysis and numerical simulation to attain the demand for the support resistance and deformation of the roof in different subzones. According to load and deformation characteristics of the roof and mechanical characteristics of NPR cables, single props, and a sliding-type gangue-retaining structure formed by U-shaped steel inserts, the supporting systems for roadways in different subzones and the constitutive model thereof were established. On this basis, the action of the supporting system was analysed and a field test was performed. The results show that the supporting system undergoes three stages of behaviour, i.e., pressure growth, yielding under constant pressure, and stabilisation during whole entry retention. It can guarantee the collaborative deformation of the supporting systems with the roof on the premise of constant support resistance, thus satisfying the requirement for roadway protection. The roadway 150 m back from the working face is stable, and the final convergence between the roof and floor of the retained entry is 257 mm, showing a favourable entry-retention effect.
1. Introduction
Coal is among the most important basic sources of energy, accounting for approximately 30% of global energy consumption [1, 2]. China, as a large energy consumer, consumes coal in an amount sufficient to account for over 50% of all primary energy demand by 2030 [3, 4]. Statistics show that about 55% of the recoverable coal reserves occur in inclined coal seams in China [5]. Therefore, safe and efficient exploitation of inclined coal seams is important when trying to relieve the pressure on the nation in terms of the demand for energy. At present, traditional longwall mining techniques are mostly used for mining inclined coal seams, which leads to the average mine recovery rate being only about 50%, as many coal pillars acting as shields are left underground and cannot be reclaimed [6, 7]. Therefore, it is necessary to develop no-pillar mining techniques for use in inclined coal seam to increase resource recovery.
The gob-side entry-retention technique based on pressure relief by roof cutting, one of the several no-pillar mining techniques, can actively change the structural shapes of the roadway roof and the roof of goaf by pressure relief by roof cutting. Furthermore, it retains the gob-side entry and no-pillar mining by roof caving, and bulking of caved materials in goaf occurs [8]. The technique overcomes various problems (including waste of resources, difficulty in subsequent mining operations, and accumulation of gas in corners) induced by remaining coal pillars, improves the stress environment in roadways, and increases the safety and economic benefit of mines. The technique has been regarded as one of major directions guaranteeing the sustainable development of China’s coal industry. To develop such nonpillar mining techniques, many researchers have investigated gob-side entry retention based on pressure relief by roof cutting. Compared with the traditional gob-side entry retention based on roadside packing, the gob-side entry-retaining technique based on roof cutting should be cut off the floor; the structure of surrounding rock and the stress distribution in the roof change after cutting; thus, it is necessary to investigate the roof structure and control techniques used when cutting a roadway. Based on the mechanism of pressure relief by roof cutting, He et al. [9] established several “surrounding rock structure-roadside support body” mechanical models, using surrounding rock control techniques such as constant-resistance large-deformation bolts and roadside props; Ma et al. [10] established mechanical models of the coal support area, dynamic pressure zones, and roadway table areas in gob-side entry retention formed by roof cutting, and the support system was constructed using NPR cables, pier-type unit supports, and I-beams. Gao et al. [11] theoretically studied changes in roof structures during the entire entry-retaining process, and the mechanical models of the roof during the entire entry-retaining process were established; the constant-resistance and large-deformation anchor cables, gangue prevention structures, and pier-beam unit supports were designed as entry-in support bodies. Fan et al. [12, 13] constructed a mechanical model of “surrounding rock‐supporting body” for bolt-grouting controlled roof cutting for gob-side entry retaining; a grouting cable, a single prop, and a steel beam were used for roadway retention. Han et al. [14] established a lateral cantilever fractured structural mechanical model of gob-side entry retaining and found that the use of an unduly short cantilever would result in insufficient pressure relief and an unduly long cantilever would generate greater stresses in the roof, resulting in intense deformation. Guo et al. [15] constructed a mechanical model of the surrounding rock in a gob-side entry under conditions of roof cutting and without roof cutting: the surrounding rock control mechanism using blasting fracturing in the roof was analysed. Based on the principle of gob-side entry retention by roof cutting, Sun and Zhang et al. [16, 17] established a simplified mechanical model of the roof structure formed in the gob-side entry: the roadway was reinforced by large-deformation anchor cables, a lattice beam and a single prop, I-beam, and steel mesh.
In addition, scholars have conducted research into surrounding rock structures and roof control during gob-side entry retention based on roadside packing. Ning et al. [18] established a mechanical model of the support structure for gob-side entry retention using cement-based filling beside a roadway in inclined coal seams; a roadside backfill-truss support structure comprising three anchor bolts and a steel strip was adopted. Based on the gob-side entry-retention method involving the use of backfill material for wall construction, Du et al. [19] established a gob-side entry-retaining mechanical model and calculated the side support resistance; the gob-side entry-retaining approach by backfilling was proposed. Yan et al. [20] proposed a roof-cutting mechanical model of gob-side entry retention with a roadside backfill body; the roof-cutting resistance and the width of the roadside backfill body were calculated. Luan et al. [21] analysed the surrounding rock structure formed in gob-side entry retention by building a concrete wall beside the roadway; bolt, cable, steel pipe, and angle irons were used to reinforce the surrounding rock. Wu et al. [22] presented a dynamic-equilibrium mechanical model of the main roof structure; the mechanisms of formation of different types of short cantilever rock beam structures were clarified. This review of the literature elucidates the significant contributions made to the understanding of surrounding rock structure and control techniques of gob-side entry retention: however, previous research into the roof structure was mainly aimed at near-horizontal coal seams; there are few studies on the roof structure of the retained roadway formed by roof cutting in inclined coal seams. The research into control techniques applicable to retained roadways was mainly based on engineering test data; there remains a lack of theoretical research on supporting mechanisms and techniques suited to retained roadways formed by roof cutting in inclined coal seams.
In view of this, the aims of this study were fourfold: (1) a mechanical model for surrounding rock structures in different subzones of a retained roadway in inclined coal seams was established, and a method of calculation of the support resistance in the roadway was determined. (2) The deformation characteristics of roadway roof in the process of retaining the roadway under the condition of roof cutting were studied. (3) The yielding supporting systems in different subzones and their constitutive models were constructed to ascertain the working principles and efficacy of the supporting systems. (4) The effect of the support system was verified by in situ testing. The results enrich the understanding on the mechanism of retained roadway formed by roof cutting and provide a reference and guidance for the technology’s further promotion and application in inclined coal seams.
2. Engineering Background
The test was conducted based on the 110,301 working face of Taiping Coal Mine located in Xuanwei, Yunnan Province, China. The working face is at a burial depth of 50 to 150 m, with a strike length of 600 m and a dip length of 125 m (it lies adjacent to the 110,303 working face). The working face is mined by using fully mechanised coal mining along the strike, with the retreating longwall moving by 4 to 5 m/d. The coal seam in the working face has a thickness of 1.2 to 1.8 m, with an average thickness of 1.5 m; the dip angle of the coal seam is between 25° and 32°, with an average dip of 28°. The roof is siltstone with a thickness of 13.0 m while the floor of the coal seam in the working face is a silty mudstone of 18.7 m thickness. The basic conditions of the working face are shown in Figure 1.

(a)

(b)

(c)
3. Materials and Methods
3.1. Research Framework
Figure 2 depicts the framework of this study. Firstly, combined with the geological features and the principle of roadway retention as a result of roof cutting, the mechanical characteristics of the rock surrounding the retained roadway in different subzones were analysed; the deformation characteristics of the retained roadway roof were calculated by numerical simulation; and the mechanical characteristics of the support materials were also studied.

Secondly, combined with the mechanical characteristics of the retained roadway, the deformation characteristics of the retained roadway roof, and the mechanical characteristics of support materials, the supporting systems for roadways in different subzones and the constitutive model thereof were established, and the working principles and efficacy of the supporting system were analysed.
Thirdly, a full-scale engineering test was performed, and the support effect was assessed.
3.2. The Mechanical Analysis of Surrounding Rocks of the Retained Roadway
3.2.1. The Principle of Gob-Side Entry Retention through Pressure Relief by Roof Cutting
At the core of the gob-side entry-retention technique through pressure relief by roof cutting is the directional presplitting of the roof of the retained roadway along the direction of advance direction of the working face through bilaterally cumulative explosion. This can cut off the mechanical connection between the roof of the retained roadway and the goaf; additionally, the roadway roof is strengthened with negative Poisson’s ratio (NPR) cables to enable the controllable deformation of the mined working face during entry retention. As the working face is advanced, the roof of the goaf is caved timeously along the cutting face and then bulked to support the overlying strata; finally, the roadway roof gradually stabilises, as shown in Figure 3 [23].

3.2.2. The Mechanical Analysis of Surrounding Rocks in Different Subzones
According to the principle of gob-side entry retention through pressure relief by roof cutting, gangues in goaf experience dynamic evolution such as caving, compaction, and stabilisation with the mining of the working face after presplitting cutting of the roof ahead of the working face. In this case, the corresponding surrounding rock structures in the roadway are also changed, leading to the difference of stress conditions. According to different stress conditions on surrounding rocks, the gob-side roadway retained by roof cutting can be partitioned into an advanced zone, a dynamic pressure zone, and a stable zone (Figure 4).

The advanced zone of the retained roadway is located within a certain range ahead of the working face. In this case, although a cutting face has been formed in the roadway roof under the effect of cumulative explosion, the cutting face delivers a certain included angle with the vertical direction and coal can still provide the support force for the roadway roof. Therefore, the two ends of the roadway roof in the zone are supported by coal and are insignificantly deformed. The surrounding rock structure is demonstrated in Figure 4(a); the dynamic pressure zone of the retained roadway is located within a certain range at back of the working face. In this context, one end of the roadway roof is supported by coal while the other end is suspended. The gangues in the goaf within the zone are in constant motion, and thus they are greatly influenced by dynamic pressure. The surrounding rock structure is shown in Figure 4(b); the stable zone of the retained roadway is situated at back of the dynamic pressure zone, where the caved gangues in the goaf are compacted and the roadway roof does not subside any further. Furthermore, a stable structure with one end supported by coal and the other supported by caved gangues is formed. The surrounding rock structure is shown in Figure 4(c).
The simplified mechanical models for surrounding rocks in different subzones were established based on the characteristics of the surrounding rock structures. The deformation of the roof and the demand for support resistance in the corresponding subzones were calculated.
3.3. Numerical Simulation of Deformation Characteristics of the Roof during Entry Retention
Numerical simulation and laboratory test are popular methods used to study rock mechanics and engineering problems [24–26]. Continuum modelling approaches based on numerical simulation are appropriate with regard to large-scale engineering problems [27]; therefore, the continuous-medium software FLAC3D was selected to simulate the deformation characteristics of a roadway roof in the process of roadway retention.
According to geological conditions in the Taiping Coal Mine, the conveyor gateways of the 110,301 working face therein were simulated by using FLAC3D. The model measures 200 × 100 × 160 m (length × width × height). The simulated working face has a dip length of 110 m and a strike length of 100 m. Two conveyor gateways have a height of 2.55 m, with the widths of the roof and floor being 3.6 and 4.4 m, respectively. Two conditions (with and without roof cutting) were simulated, as shown in Figure 5. The constitutive model used was a Mohr–Coulomb model, the roof cutting was simulated by applying the null unit, and the goaf near the position of roof cutting was backfilled based on the strain-hardening model. The working face was excavated for 60 m, and the model was loaded in 3000 steps. The physicomechanical parameters of various strata are listed in Table 1.

To compare the deformation characteristics of the retained roadway roof under the two conditions (with and without roof cutting), displacement-monitoring points were distributed in the middle part of the roadway roof at intervals of 10 m within the first 60 m behind the working face. On this basis, data pertaining to the change in vertical displacement of the roadway roof throughout the loading process under the two conditions were provided. The vertical displacement of the roof of the retained roadway was extracted from the model after completion of the calculations at positions shown in Figure 6.

3.4. Supporting Material and Its Mechanical Characteristics
Cable, bolt, single props, steel, and bolt-grouting reinforcement are commonly used methods to support the fractured surrounding rock in rock engineering [28, 29]. The support material of the roadway formed by roof cutting involves three types: active support material such as negative Poisson’s ratio (NPR) anchor cable, temporary support material in roadways such as single props, and support material with the gangue-retaining structure.
3.4.1. Active Support of the NPR Cables
The NPR cables are anchor cables which can be axially extended with increasing rock deformation to slowly absorb the produced energy during surrounding rock deformation [30]. The NPR cables have the properties of high constant resistance force, large elongation, and strong energy absorption [31]. Since the development of the NPR cable, it has been successfully applied to the stability control of many coal mines. The results are remarkable.
In terms of the core technology, the constant-resistance device is added to an ordinary anchor cable; consequently, the NPR cables maintain a constant resistance (preset artificially) and are deformed through relative slip in the interaction with surrounding rocks. The NPR cables are mainly subjected to elastic deformation when the load on the NPR cables is lower than the constant resistance; after the load reaches a constant resistance, the NPR cables start to slide and yield on the premise of keeping the resistance unchanged until the cables are stable again. The constitutive model of the NPR cables can be represented by using a Saint-Venant body (S body), as shown in Figure 7 [10].

The NPR cables are deformed under the effect of the tension σ, whose total strain is equal to the sum of the strains in various elements. The stress-strain relationship of the NPR cables satisfies the following equation:
According to Figure 7 and (1), it can be found that the tension on the anchor cables increases with increasing ε2 when the friction disk in the S body does not slide; in the case that the friction disk in the S body slides, σ remains unchanged as ε1 rises and ε2 remains constant; therefore, the NPR cables show the following mechanical characteristics: the steel strands are subjected to the elastic deformation, and the support resistance provided by anchor cables grows with the increase of the deformation on condition that the tension is lower than the constant resistance; when the tension exceeds the constant resistance of the anchor cables, the anchor cables are able to adapt to surrounding rock deformation by slip of the constant-resistance device which takes over the role of the stretching of the steel strand; moreover, the resistance remains the same during the deformation of the anchor cables.
3.4.2. Temporary Support Using Single Props
The temporary support using single props mainly relies on the single hydraulic prop, which is composed of various parts including the cylinder block, movable column, and valves. In the case that the load applied to the prop exceeds the rated working resistance thereof, the high-pressure liquid in the internal cavity of the prop drives opens the safety valve, thus relieving the pressure. The constitutive model of the single hydraulic prop can be formed by connecting the S body with a damping element (N body) in parallel, as shown in Figure 8.

The single hydraulic prop is deformed under the effect of the pressure σ, whose total stress equals the sum of stresses of various elements and the total strain is equivalent to the sum of strains of various elements connected in parallel. The stress-strain relationship of the single hydraulic prop satisfies the following equation:
It can be found from Figure 8 and (2) that when the pressure applied on the single props exceeds the rated working resistance, the support force provided by the single props is closely related to the rate of deformation. When the rate of deformation is large and the liquid fails to escape timeously through the safety valve, the total support force increases. Therefore, the single props exhibit the following mechanical characteristics: before reaching the rated working resistance, the support resistance provided by single props rises with the growth of the external load; after reaching the rated working resistance, single props realise yielding by leakage of liquid through the safety valve. The support resistance remains unchanged if the liquid can be timeously leaked through the safety valve; otherwise, the support resistance increases. This indicates that single props are also able to control the rate of roof subsidence to prevent the occurrence of impact failure.
3.4.3. The Sliding-Type Gangue-Retaining Support Formed by U-Shaped Steel Inserts
The sliding-type gangue-retaining support structure is composed of two ordinary U-shaped steel members. Fastened by the clamping device, the two U-shaped steel members can undergo relative slip, thus realising collaborative deformation with the roadway roof. The frictional resistance between the two U-shaped steel members can be adjusted by controlling the torque applied to the clamping nuts. When the roof pressure is lower than the frictional resistance between the U-shaped steel members, they are mainly subjected to the elastic deformation and do not slide; the U-shaped members slide and yield when the roof pressure exceeds the frictional resistance between them. The constitutive model of the U-shaped steels can be expressed by use of an S body, as shown in Figure 9.

The stress-strain relationship of the sliding-type U-shaped steel members satisfies (3)); that is, they allow the support resistance to be the same during deformation, thus realising yielding under constant resistance. This guarantees that the roadway roof rotates and subsides in a controllable manner.
As shown in Figure 9 and (3), the stress σ remains unchanged as ε5 rises and ε6 remains constant when the friction disk in the S body slides. Thus, the sliding-type gangue-retaining structure formed by U-shaped steel inserts delivers the following mechanical characteristics: when the pressure is less than the static friction between U-shaped members, the members are subjected to elastic deformation, and the support resistance provided by the gangue-retaining structure increases with the elastic deformation of the U-shaped steel members; when the pressure exceeds the static friction between the U-shaped steel members, the gangue-retaining structure is able to realise yielding through the sliding occurring between the two members, thus keeping the resistance unchanged during yielding.
3.5. Engineering Test
According to the geological conditions of the working face and the deformation characteristics of the roadway roof, the asymmetric support form in the roadway is determined to strengthen the roof in the side near the cutting. The specific parameters are shown as follows (Figure 10): (1) presplitting cutting to a depth of 7 m occurs 200 mm away from the upper corner of the roadway in the mining side; it shows the included angle of 15° with the vertical direction and is inclined to the side of the mining working face (Figure 10(a)). (2) The spacing between blast holes is set to 500 mm, with an explosive charging size of φ32 × 330 mm and a charge structure of 2.5 + 2 + 2 + 1(Figure 10(b)). (3) Two columns of NPR cables are distributed at the roof, with the length and diameter of anchor cables of 9.3 m and 21.6 mm, respectively. The first column of NPR cables is 300 mm away from the roof-cutting line, with a row spacing of 1 m; these are connected with φ16 mm reinforced bar ladders along the strike of the roadway. The second column of NPR cables is 1300 mm from the first column, with a row spacing of 1800 mm. (4) The temporary support is formed by using single props in combination with π-shaped beams with a length of 2.8 m, with each beam having four columns. (5) The gangue-retaining structure is made of two 29 U-shaped steel inserts, and the spacing between adjacent supports is set to 500 mm.

(a)

(b)
4. Results
4.1. The Structural Model and Mechanical Analysis of Surrounding Rocks of the Retained Roadway Formed by Roof Cutting in Inclined Coal Seams
4.1.1. Advanced Zone of the Retained Roadway
The roadway roof in the zone is synchronously supported by coal on two sides and the supporting structure in the roadway, whose simplified mechanical model is shown in Figure 11.

After presplitting cutting of the roof, the cutting face makes a certain included angle to the vertical. The coal still presents a certain supporting effect on the roadway roof. In this condition, the roadway roof can be regarded as a beam with one end clamped and the other end simply supported. Accordingly, the deflection of the roadway roof is given bywhere E and I represent the elastic modulus (MPa) and the sectional moment of inertia (m4) of the beam, respectively.
It can be seen from (4) that the deflection of the roof is maximised at x = a/2, as shown in the following equation:
In this case, the anchor bolts and anchor cables are taken as the main supporting structures of the roadway roof. Accordingly, it is feasible to calculate the support resistance Pi in the roadway:
In the mining of the working face, the roadway roof not only bears the deadweight of the overlying strata but also suffers from the dual influence of the abutment pressure and the mining-induced dynamic pressure. Therefore, it is necessary to support the roadway roof within a certain range ahead of the working face to prevent damage under the influence of the dynamic pressure during the mining of the working face.
4.1.2. The Dynamic Pressure Area in the Retained Roadway
The gangues in the goaf within the area are in a state of near-caving and they have not, as yet, effectively supported the roadway roof (Figure 12).

In this case, the roadway roof can be considered as a cantilever, whose deflection is given by
When x = a, the defection of the roof reaches a maximum:
The lateral fracture span L of the strata in the main roof and the horizontal side thrusts on rock blocks A (TA) and B (TB) are separately expressed as follows [32]:where L, S, and l refer to the lateral fracture span (m) of the main roof, the length (m) of the working face, and the step distance (m) of periodic caving of the working face, respectively.where q′, h, and ΔSA separately denote the heavy load (kN) of the main roof within the unit length, the thickness (m) of the main roof, and the subsidence (m) of rock block A at position B′.
The width x0 of the limit equilibrium zone within the coal sides and the support force σ of the plastic zone on the coal are separately expressed as follows [33]:where c0 and φ separately refer to the cohesion (MPa) and internal friction angle (°) of the interfaces between the coal seams with the roof and floor strata; hd, A, k, γ, H, and px denote the thickness (m) of the coal seams, the lateral pressure coefficient, the maximum stress concentration factor, the average bulk density (kN/m3) of strata, the mining depth (m), and the support strength (MPa) of the coal sides, respectively.
In this context, it is feasible to perform stress analysis on rock blocks A and B by using the static equilibrium method to calculate the support resistance required to reach roadway stability.
As for rock block B, the static equilibrium equation can be attained according to , , and .
According to (13) to (15), the shear forces on rock blocks A (NA) and B (NB) can be calculated:
As for rock block A, the following equation can be obtained according to :
By simultaneously solving (10) to (18), the support resistance Pi required to reach roadway stability can be calculated:
It can be seen from (8) that the roof of the retained roadway presents the maximum vertical displacement in the side of the roof cutting. Therefore, more attention must be paid to strengthening of the support on the roof in the side of the roof cutting of the roadway. As shown in Figure 11, the roadway roof is mainly supported by coal sides and the supporting structure in the roadway within the dynamic pressure area of the retained roadway. Thus, it is necessary to strengthen the support in the roadway based on the primary support. The resultant force on the support required for reach roadway stability can be calculated according to (19).
4.1.3. The Stable Zone within the Retained Roadway
The gangues in the goaf in the stable zone of the retained roadway have been completely caved and compacted; moreover, the rotation of the main roof and its overlying strata ceased. In this case, the broken gangues can support the rock beam; furthermore, the overburden strata movement stops, and the roadway roof tends to stabilise. In this condition, the simplified mechanical model of the roadway is shown in Figure 13.

As for rock block B, the following equation can be obtained according to (8):
As for rock block A, the following equation can be attained based on :
According to (20) to (22), the support resistance in the roadway within the stable zone of the retained roadway can be calculated:
In (23), the support force imparted by gangues in the goaf on the roof can weaken the support resistance in the roadway. Therefore, it is essential to maximise the bearing capacity of the broken stone sides in the stable zone of the retained roadway. By comparing (23) with (19), it can be found that it is feasible to remove the temporary supports in the roadway progressively and in accordance with stress-monitoring results from the roadway: if , all temporary supports in the roadway are removed.
It can be seen from (6), (19), and (23) that the demands for support resistance in the roadway within different subzones all gradually reduce with increasing dip angle. The maximum deflection of the roof in the dynamic pressure area of the retained roadway is found near the mining side. Therefore, strengthening of the roof near the mining side warrants greater attention, thus reducing the support cost on the premise of guaranteeing roof stability.
4.2. Analysis of Deformation Characteristics of the Roof during Entry Retention by Roof Cutting
Figure 14 shows changes at different positions of the roof of the retained roadway under load (under both conditions). It can be seen from the figure that the roadway roof is always deformed in conditions without roof cutting. The average subsidence of the roadway roof reaches 1.52 m at the end of the loading stage; after pressure relief by roof cutting, the roadway roof gradually stabilises after 2000 loading steps. The average subsidence of the stabilised roadway roof is 0.55 m, which decreases by 64% compared with the result before roof cutting.

(a)

(b)
The nephogram of the vertical displacement of the roof is drawn, as shown in Figure 15. According to the spatial characteristics of the roof subsidence, the positions with the maximum subsidence of the roadway roof under the two conditions are both on the side near the goaf (roof cutting), which conforms to the result deduced from (25). After pressure relief by roof cutting, the average roof subsidence in the side near roof cutting within the range of 60 m behind the working face is 0.55 m; at the same position, the average subsidence of the roadway roof without roof cutting reaches 2.78 m, which is five times as large as that in former case. In this condition, this means that failure and collapse have occurred in the roadway.

(a)

(b)
4.3. Establishment of a Roof-Support System in Different Subzones
4.3.1. The Advanced Zone of the Retained Roadway
The rock mechanical properties will gradually become weaker under the blasting and construction disturbances [34]. The roadway roof weakens slightly under the influence of the top cutting blasting and mining-induced dynamic pressure. The supporting system in the zone can be constructed by applying NPR cables and single props, as shown in Figure 16. According to analysis in Section 3.4, the constitutive model of the supporting system is formed by connecting the S body and a Coulomb body in parallel (known as the S-N-S model), as shown in Figure 16.

(a)

(b)
Under these parallel conditions, the total stress on the supporting system is the sum of the stresses on various parallel elements, and the total strain is equivalent to the total of the strains in various parallel elements, shown as follows:
The constitutive equation of the model is expressed as follows:
On condition that the roof bears only a low pressure, the supporting system is not subjected to yielding deformation (ε1 and ε4 equal zero). The support resistance of the supporting system is mainly offered by the NPR cables and single props in the form of the elastic deformation, which is calculated as follows:
When the roof pressure is high, the supporting system is subjected to yielding deformation (ε1 and ε4 are nonzero). The support resistance of the supporting system is mainly provided by the constant resistance σ1 of the NPR cables and the rated working resistance σ4 of single props, which is calculated as follows:
According to (27) and (28), it can be found that the support system for the advanced zone of the retained roadway shows the following mechanical characteristics: the support resistance increases with increasing external load. After reaching the maximum support resistance of the support units, the system enables yielding deformation and maintains constant stiffness during deformation.
4.3.2. The Dynamic Pressure Area of the Retained Roadway
Under the disturbance of the caving of gangues in the goaf and overburden strata movement, the surrounding rock of roadway is prone to dynamic disaster, which influences the long-time stability of rock structures and rock engineering [35]. This requires the supporting system in the roadway to provide sufficient resistance to control the rate of roof subsidence, so that the rock beams can make contact with the gangues in a slow and controlled manner; therefore, it is necessary to establish the supporting system by utilising the NPR cables, the single props, and the sliding-type gangue-retaining structure formed by U-shaped steel members. In this case, the combined constitutive model of the three support units is formed by connecting three S bodies and a Coulomb body in parallel (known as the S-N-S-S model), as shown in Figure 17.

(a)

(b)
The stress on various units is equal to the sum of stresses on various parallel elements, and the total strain is equivalent to the strains in various parallel elements, as given by the following equations:
The constitutive equation for the model is expressed as follows:
The support resistance is given by
When the roof pressure is low, the supporting system is not subjected to yielding deformation (ε1, ε4, and ε6 are zero). The support resistance of the supporting system is mainly offered by the elastic deformation of the NPR cables, the single props, and the gangue-retaining structure formed by U-shaped steel members, as given by
On condition that the roof bears high pressure, the supporting system undergoes yielding deformation (ε1, ε4, and ε6 are nonzero). The support resistance of the supporting system is mainly provided by the constant resistance σ1 of the NPR cables, the rated working resistance σ4 of single props, and the maximum static friction σ5 between U-shaped steel members in the gangue-retaining structure, as given by
According to (32) and (33), the mechanical characteristics of the supporting system for the dynamic pressure area of the retained roadway are shown as follows: different support units can be collaboratively deformed. When various support units do not reach their maximum support resistance, the support resistance provided by the supporting system gradually increases with increasing roof deformation. After the support resistances of various support units all reach the maximum, the maximum resistance provided by the supporting system is related to the rate and amount of compression of the single props. Therefore, in (33) can be regarded as a fixed value; afterwards, σ also becomes a fixed value with the change of ε. This implies that the supporting system can realise yielding under constant resistance after reaching the maximum support resistance, thus realising collaborative deformation between the roof and the supporting structure. Furthermore, the roof pressure is transferred to the gangues in the goaf. If the support strength is not sufficient, becomes extremely large, bursting or bending the single props, thus risking a collapse.
4.3.3. The Stable Zone of the Retained Roadway
The mining-induced dynamic pressure disappears, and the gangues in the goaf can support the roadway roof and its overlying strata. Thus, the roadway roof tends to stabilise. On this condition, the rock beam of the main roof is supported by coal and gangues to form a stable structure. The immediate roof is suspended below the main roof with the NPR cables. Therefore, it is feasible to withdraw the single props and retain the NPR cables and the sliding-type gangue-retaining structure to support the roof and retain gangues. In this case, the combined constitutive model for the two support units is formed by using two S bodies connected in parallel (known as the S-S model), as shown in Figure 18.

(a)

(b)
On this condition, the relationships between the total stress and strain of the supporting system with stresses and strains of various parallel elements satisfy the following equations:
The constitutive equation of the model is shown as follows:
The support resistance provided by the supporting system is calculated as follows:
In this stage, the roadway roof tends to stabilise, so it can be thought that the total deformation ε of the supporting system stabilises; thus, both ε1 and ε6 are zero in (36). The support force offered by the supporting system is calculated as follows:
Therefore, the supporting system for the stable zone of the retained roadway can still provide the constant support resistance for the roadway roof. On this condition, the roadway roof remains stable under the combined support of the gangues in the goaf, the NPR cables, and the sliding-type gangue-retaining structure formed by the U-shaped steel members. Without being disturbed, the pressure on the supporting system does not vary.
Therefore, the established supporting systems in different subzones can adapt to the deformation characteristics of the roof in different subzones. On the premise of playing the strongest supporting role in the system, they are able to realise yielding under constant resistance and guarantee controllable deformation of the roof, thus finally ensuring the stability of the retained roadway.
4.4. Engineering Test Effect
The convergence between the roadway roof and floor and the stress on single props during entry retaining are monitored, as shown in Figures 19 and 20: the deformation of the roadway roof and floor is mainly between 0 and 150 m behind the working face. Within the first 55 m, the roof and floor move faster, where the roof subsidence reaches 47 mm, accounting for 61% of the total deformation; the convergence between the roof and floor is 147 mm, which accounts for 57% of the total deformation; in the range of 55 to 150 m behind the working face, movement of the roadway roof decelerates and gradually stabilises beyond 150 m. Eventually, the roof subsidence, floor heave, and convergence between the roadway roof and floor reach 77, 180, and 257 mm, respectively.

(a)

(b)

According to the change in the stress on each single prop, it can be seen that the whole entry-retaining process can be divided into three stages: (I) Pressure-growth stage: in the range of 0 to 20 m behind the working face, the pressure on each single prop constantly increases. In this stage, single props have not reached their rated working resistance, and the support resistance linearly rises with increasing roof subsidence. (II) Yielding stage under constant pressure: in the range of 20 to 150 m behind the working face, the deformation and rate of deformation of the roof are large, and therefore the support resistance of single props increases rapidly during compression due to a failure to leak liquid timeously via the pressure-relief valves. Thus, the resistance is found to fluctuate. (III) Stable stage: beyond 150 m behind the working face, the roof no longer subsides, and the pressure on each single prop gradually stabilises. This indicates that the broken rock has been able to provide stable support to the overlying strata of the main roof. In this case, single props can be progressively withdrawn, leaving only the anchor bolts (cables) and the sliding-type gangue-retaining structure formed by U-shaped steel members in the roadway to support the roof.
By comprehensively analysing the deformation characteristics of the roadway and the stress characteristics of single props, the supporting system is seen to undergo three stages, i.e., pressure growth, yielding under constant resistance, and then stabilisation, under the influence of the roof movement during entry retention. It can adapt to the deformation characteristics and satisfy the demand for support to the roof of the retained roadway. After withdrawing the single props (Figure 21), the final convergence between the roof and floor of the roadway reaches 257 mm, showing a favourable effect of entry retention.

5. Discussion
5.1. The Structural Model and Mechanical Characteristics of Surrounding Rocks of the Retained Roadway
A simplified mechanical model of the surrounding rock structure in different subzones of the retained roadway formed by roof cutting in an inclined coal seam was established; the formulae predicting the support resistance in the roadway were attained; this was used to calculate the characteristics of the roof deformation and the support requirements in the roadway under the influence of the dip angle of the coal seam. This is different from previous studies of horizontal coal seams [10–16]. However, the influence of mining dynamic pressure was not reflected in the formula; the influence of mining dynamic pressure should be fully considered in the construction of support system.
Due to the influence of the advanced dynamic pressure generated during mining, the roof of the advanced zone of the retained roadway is necessary to increase the support strength in the roadway to resist the advanced dynamic pressure. The roof in the dynamic pressure area of the retained roadway can be regarded as a cantilever beam with one end supported by coal and the other end suspended. Influenced by overburden strata movement and caving of gangues, the roof subsides. In this case, the support resistance required in the zone is maximised. The roof in the stable zone of the retained roadway can be taken as a simply supported beam, with one end supported by coal and the other end supported by gangues. The influence of the mining-induced dynamic pressure vanishes: in this condition, it is feasible to weaken the support resistance in the zone and maintain roadway stability by relying on coal and broken gangues.
5.2. Deformation Characteristics of the Roof during Entry Retention by Roof Cutting
We have obtained the deformation characteristics of the roadway roof under the conditions of roof cutting and without roof cutting. Through comparative analysis, it is found that when relieving the pressure by roof cutting, the roadway roof gradually stabilises after 2000 loading steps, suggesting that the gangues in the goaf have been compacted; in this context, the roadway roof is synchronously supported by the coal sides and gangues in the goaf and eventually tends to stabilise; without pressure relief by roof cutting, gangues fail to provide the effective support force for the roadway roof and its overlying strata owing to their caving over a certain height. As a result, the roadway roof suffers significant subsidence and then collapses. It can be found that, given a sufficient height when cutting the roof, the key to realising self-stabilisation of the roadway is that the support in the roadway can guarantee a controllable rotation and subsidence of the roof to acquire support from gangues.
5.3. Deformation Characteristics of the Roof-Support System in Different Subzones
According to mechanical characteristics of supporting materials, the rheological mechanical models and constitutive equations of supporting systems for roadways in different subzones were established for the first time.
We found that the roof in the advanced zone of the retained roadway is slowly deformed. In this condition, the support resistance provided by the supporting system rises with the increased deformation of the roof, thus inhibiting roof subsidence. This is because the NPR cable has not reached the constant resistance, and the single props have not reached the rated working resistance; the support force delivered by the supporting system is mainly provided through the elastic deformation of the steel strands of the NPR cables and single props.
The roof in the dynamic pressure area of the retained roadway is significantly deformed. The supporting system can keep the resistance constant after resistances of various support units all reach a maximum, thus realising collaborative deformation with the roof. This is because the support units have reached their maximum support resistance, the NPR cables realise yielding by slip of the constant-resistance device, and the single props realise yielding by leakage of liquid through the safety valve. Thus, the roadway roof and its overlying strata are controlled in their descent to contact with gangues.
As for the stable zone of the retained roadway, the supporting system is no longer deformed. This is because the main pressure on the roof has been transferred to gangues in the goaf, on the premise of a constant support resistance, and the support system works in tandem with the gangues to ensure the stability of the roadway roof.
5.4. Advantages and Limitations
The technique overcomes the waste of resources induced by the retention and abandonment of the supporting coal pillars, and it has been widely used in near-horizontal, thin, medium-thick, and thick coal seams [10, 11, 36, 37]. In this study, the technique of roadway retention using roof cutting was successfully applied in an inclined coal seam, thus addressing the problem of the waste of resource during mining and providing guidance for future projects relying on roadway retention formed by roof cutting under similar geological conditions.
It should be pointed out that this study did not fully consider the influences of coal seam thickness and lithological changes in the roof on the retained roadway formed by roof cutting in inclined coal seams. In future studies, research into roadway retention by roof cutting in coal seams of different thickness and under different roof conditions should be conducted.
6. Conclusions
By exploring the structural evolution of the deformation of the rock surrounding a retained roadway formed by roof cutting in inclined coal seams, support systems for the roof in different subzones of the retained roadway were established and applied in practice. The following conclusions are drawn:(1)The mechanical models for surrounding rock structures in the advanced zone, dynamic pressure area, and stable zone of the retained roadway formed by roof cutting in inclined coal seams were separately established. On this basis, the deformation of the roof in the corresponding subzones and formulae predicting the support resistance in the roadway were attained. These provide a theoretical basis for the parameter design of the retained roadway by roof cutting in inclined coal seams.(2)The results obtained through numerical simulation show that the roof of the retained roadway by roof cutting tends to stabilise after subsiding by some 0.55 m in the entry-retaining process; the roadway roof without roof cutting always subsides and is deformed. This implies that the roof of the retained roadway formed by roof cutting can be effectively supported by gangues in the goaf after subsiding by a certain distance, thus realising self-stabilisation of the roadway.(3)Aiming at the stress characteristics of different subzones and mechanical characteristics of the self-stabilisation of the roadway, the supporting systems in corresponding subzones were established by applying NPR cables, single props, and a sliding-type gangue-retaining structure formed using U-shaped steel members. Furthermore, three rheological mechanical models (S-N-S, S-N-S-S, and S-S) corresponding to the support systems and their constitutive equations were proposed.(4)The working principles of the support systems in different subzones were analysed from the perspective of rheological mechanics. The results show that the established supporting systems can realise yielding under constant resistance on the premise of providing the maximum support resistance among various units. Furthermore, the roof pressure is slowly transferred to gangues in the goaf, thus finally realising the controlled stabilisation of the roadway.(5)The results obtained through in situ testing and measurement indicate that the supporting system is best placed in the pressure-growth stage of the gob-side entry-retention process within the first 20 m behind the working face; the system mainly realises yielding under constant pressure in the range of 20 to 150 m behind the working face; the roof tends to stabilise beyond 150 m behind the working face. The final convergence between the roof and floor of the roadway reaches 257 mm, showing a favourable entry-retention effect.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research was funded by the National Natural Science Foundation of China (no. 52074300); Science and Technology Department Joint Fund Project of Guizhou Province, China (no. LH[2017]7036); and Youth Science and Technology Talent Growth Project of Guizhou Provincial Department of Education (no. ky[2018]414).