Analysis on Catastrophe Theory during First Weighting Sliding Instability and Support Crushing of Main Roof with Large Mining Height in Shallow Coal Seam

Roof sliding and instability along the coal wall usually occur in the working face at large mining heights during the process of the first weighting, which causes roof cutting and support crushing. A mechanical model consists of the main roof, immediate roof, and support based on the nonlinear characteristics of the failure and instability of the immediate roof under the abutment pressure, which we constructed to study the step sinking of the main roof, as well as to assign the reasonable value of the support resistance during the first weighting. The instability mechanism of the system was studied by the catastrophe theory and the principle of energy conservation. A conclusion was drawn that the combined cantilever beam structure for the immediate roof will form with the increase of the mining height, and the instability of the immediate roof causes the catastrophic instability of the system. The system instability was found to be related to the stiffness ratio K, material parameters, the load Q, and the first weighting interval of the main roof by analyzing the necessary and sufficient conditions for system instability. The influence degree of each parameter on the stiffness ratio K was as follows: elastic modulus E > support stiffness k1 > cross-section area a > immediate roof thickness H. The calculation equations of support resistance and subsidence of roof step were obtained. The method of judging the roof instability using catastrophe theory was proved as reasonable on the basis of the monitoring example of no. 12401 working face in Shendong mining area, China. On this basis, a reasonable value of support resistance was further calculated, and the working face was maintained safely when the support resistance exceeded 19,232 kN.


Introduction
Shallow coal seams are widely distributed in the northwestern territory of China, which is characterized by shallow burial depth, thin bedrock, and an overlying thick and loose sand layer [1]. The mining height of coal seam is increasing, reaching a maximum of 8.8m at present with the continuous improvement of fully mechanized mining degrees, which makes the goaf space formed by coal seam mining increase constantly. The original main roof is broken into a cantilever beam structure with the increasing of goaf space, which forms a new immediate roof, and increases the thickness of the immediate roof. As a result, a new main roof is formed in the upper rock stratum and a cantilever Appl. Sci. 2020, 10, 5408; doi:10.3390/app10165408 www.mdpi.com/journal/applsci Appl. Sci. 2020, 10, 5408 2 of 23 beam-masonry beam structure is formed, which makes the strata behavior violent and the dynamic load phenomenon obvious [2,3]. The strata behavior is the most violent and difficult to control during the first weighting. Roof sliding and instability along the coal wall usually occurs, resulting in support crushing or water/sand burst, which brings many security risks to coal mining. Therefore, it is of great significance to study the relationship between the support and the surrounding rock during the first weighting of working face, as well as to analyze the mechanical mechanism and key disaster-causing factors of the main roof step sinking along the coal wall then to provide the determination method of the support resistance, which has an important role in ensuring the safety of the working face. Many scholars have already done much work on the characteristics of mine pressure and the relationship between support and surrounding rock in working face with large mining height by various methods. Ju et al. [4][5][6] found that the breaking of the upper sub-key strata of the working face will cause the break of the lower sub-key strata, with the interval of step and the intensity of weighting in working face being changed alternately. Zhao et al. [7] found that the structure instability of the cracked roof rock beam depends on the interaction between the vertical load and the horizontal thrust. Wang et al. [8] established three kinds of mechanical models in the mining field-these being symmetrical pressure-arch, stepped pressure-arch, and rotating-squeezed pressure-arch-and derived the instability criterion of each pressure-arch structure. Szurgacz and Brodny [9,10] developed an effective method for determining the position (geometry) of the support section during the operation process and analyzed the influence of dynamic load on the work parameters of a powered roof support's hydraulic leg via tests. The research results were conducive to the safe and efficient mining of coal. Yang and Liu [11] pointed out that mining intensity, single roof structure, easy sliding, and instability were the main reasons for the violent appearance of strata behavior. Huang and Zhou [12,13] proposed the definition of the equivalent immediate roof of large mining height and analyzed the pressure characteristics of the oblique step voussior beam structure of the working face, giving the support resistance of the working face without slipping and instability. Pang et al. [14] proposed the cantilever beam + masonry beam structure and stability control technology for ultra-height mining height working face on the basis of 8.0m ultra-height mining height mining practice. Xing et al. [15] established a fracture mechanical model of key roof and derived the theoretical calculation formula for the fracture step of key roof and support resistance. Sun et al. [16] proposed a new model to describe the movement and subsidence of the inner burden during shallow coal seam mining. Feng et al. [17] found that the mining height had an exponential function relationship with the average first weighting interval and had a binomial functional relationship with average periodic weighting interval. Chinese scholars focused on the characteristics of roof breaking, pressure mechanism, and the adaptability of the support during the first weighting. Huang et al. [18] analyzed the roof control effect and working resistance overrun percentage and put forward an inversion method for working resistance optimization. Zhao et al. [19] established the calculation method for determining the dynamic load factor when the slipping instability of main roof structure occurs. Wang and Wang [20] pointed out that the increase of height to length ratio of key block causes the increase of the rotation angle required for the hinged structure to reach equilibrium. This usually results in sliding and instability of the main roof. The cut down of the entire thickness of the overlying strata along the working face is a dynamic evolutionary process during the first weighting of main roof [21]. Yang [22] derived the formulae for breaking convergence of the main roof and built the criterion of step convergence after the main roof breaking during the first weighting by initial post-buckling theory and the cusp catastrophic model. Zuo et al. [23] analyzed the blasting mechanism of shallow depth mining face with a large cutting height by fracture mechanics and calculated the initial fracturing distance of main roof after blasting.
The research on the immediate roof started earlier. Scholars mainly studied the bearing, deformation, and failure characteristics during the weighting process of the immediate roof. Qian et al. [24][25][26] pointed out that the relationship between the support resistance and the amount of roof subsidence was caused by the action of the support and the immediate roof under a given deformation of the main roof. Gong and Jin [27] carried out the in situ measurement and the simulation Appl. Sci. 2020, 10, 5408 3 of 23 experiment on the structural characteristics of roofs with large mining height, dividing the immediate roof into three types and finding that the stability control of the type III immediate roof mainly considers the thickness, location, and engineering mechanical characteristics of the key layer of the immediate roof. Du and Bai [28] divided the immediate roof into two types and found that the instability of the type II immediate roof usually causes periodic weighting on the working face. Yuan et al. [29] found that the key strata in immediate roof (KSIF) controls the caving of its upper immediate roof strata and has a direct dynamic impact on its lower strata when the KSIF breaks. Yang et al. [30][31][32] gave a method to determine the range of the immediate roof and deduced the calculation equation of roof pressure on the basis of the principle of conservation of energy. Liu et al. [33] found that a relative thick immediate roof tends to form a stable voussoir beam, while a relative thin immediate roof tends to form a cantilever beam. Their research mainly focused on the failure characteristics of the immediate roof and the determination of the support resistance, which is less involved in the role of the immediate roof in the relationship between the support and surrounding rocks and the fracture and instability mechanism of the immediate roof.
At present, the main research methods of roof instability include numerical simulation, elastic plastic theory, elastic damage theory, and catastrophe theory. As a new subject, the theory of catastrophe theory has been widely applied, not only in the field of rock system instability, but also in the field of coal and rock system instability. Wang et al. [34] analyzed the mutation instability rules of the room mining residual pillars in the shallow depth seam. Xu et al. [35] deduced the expression of the limit thickness of the goaf roof under the filling body by building the cusp catastrophe model. Mu et al. [36] revealed that the stability of layered rock mass was controlled by the external force, the modulus of the elasticity, the severity of rock mass, the inter layer shear stress, and the geometry based on the cusp catastrophe theory. Xu et al. [37,38] established the cusp catastrophe model on the basis of the energy theory, pointing out that the instability of the support system was caused by the incompatibility of energy release. The results showed that the catastrophe theory can explain the nonlinear dynamic characteristics of the instability of a coal-rock system. The deformation of the immediate roof accumulates elastic energy under the load of the overlying rock layer, and when the energy accumulates to a certain extent, it will be released instantly, resulting in the step sinking of the main roof. The crushing of the support and the ribbing of the coal wall have obvious instantaneity and abruptness. Therefore, the catastrophe theory is introduced into the stability analysis and parameter calculation during the first weighting of the working face to solve the problem of sudden disaster in coal mining.
In the system of main roof-immediate roof-support, the deformation of the main roof is generally given, and the influence degree of the support on the main roof and the effect degree of the main roof on the support depend on the mechanical properties (stiffness) of the immediate roof medium and the yield ability of the support. The key point of roof control is the control of the immediate roof, and the role of immediate roof in the relationship between the support and surrounding rock is of vital importance [21,22,39]. According to the unstable collapse process and pressure characteristics of the main roof structure of large mining height, the main roof will break and rotate, compressing the immediate roof and support when the immediate roof and support system are not enough to support the main roof. When the deformation of the immediate roof rock exceeds its deformation critical value, the immediate roof will be crushed, and the main roof will slide along the coal wall, resulting in a wide range of roof cutting, support crushing, and step sinking ( Figure 1) [40]. Therefore, the immediate roof rock mass controls the instability of the main roof, and the key point to control the large-scale pressure of the main roof is to control the instability of the immediate roof under the abutment pressure.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 4 of 23 large-scale pressure of the main roof is to control the instability of the immediate roof under the abutment pressure. Figure 1. The main roof is cut off and the immediate roof is crushed [40].
Therefore, the main objective of this study was to explore the general process and weighting characteristics of the main roof along the coal wall during the first weighting of large mining height in shallow coal seam, to analyze the mechanical mechanism of the main roof instability, to reveal the disaster-causing factors, and to explore the determination method of support resistance to ensure the stability of working face. In this study, according to the structural characteristics of the main roof cutting and instability along the coal wall during the first weighting, the system mechanical model consists of the main roof, immediate roof and support is established. This model is based on catastrophe theory and the energy conservation principle analysis method, as well as through considering the nonlinear characteristics of the failure and instability of the immediate roof. A new analysis method was proposed to study the mechanical mechanism of roof cutting along a coal wall during the first weighting. The instability process was analyzed theoretically, and the sufficient and necessary conditions for the sudden instability of the system were successfully obtained. The main influencing factors of the system catastrophe, such as stiffness ratio K, rock material parameters, external load Q, and weighting interval l are discussed. We also compared and analyzed the effects of elastic modulus E, support stiffness k1, cross-sectional area a, and immediate roof thickness H on the stiffness ratio K. After further derivation, the calculation equation of support resistance was obtained, and the main influencing factors of step sinking were analyzed. Finally, the rationality of the theoretical analysis results were verified by an engineering example, and the reasonable value of the support is given. This paper is organized as follows: the structural characteristics and cutting process of the main roof during the first weighting in the working face of large mining height in the shallow coal seam are described in Section 2. A catastrophe theoretical analysis model consisting of main roof, immediate roof, and support is given in Section 3. The catastrophe instability condition of the immediate roof, the calculation method of the support resistance, and the step sinking are deduced by catastrophe theory in Section 4. Furthermore, the theoretical results are verified in Section 5. Conclusions are drawn in Section 6.

Working Face Overview and Strata Behavior Monitoring
To analyze the roof structural characteristics of thin bedrock under thick loose layer during shallow coal seam mining, we took the no. 12401 panel of 5 −2 coal seam in Shendong mining area as the engineering background. The length of working face was 260 m with a dip angle of 1°-3°, and the average burial depth was 120 m. The thickness of the bedrock was 70 m on average. The thickness of the loose layer was 50 m on average. The thickness of the bedrock fluctuated in a small range. The roof was mostly mudstone, fine-grained sandstone, and siltstone, with an average compressive Therefore, the main objective of this study was to explore the general process and weighting characteristics of the main roof along the coal wall during the first weighting of large mining height in shallow coal seam, to analyze the mechanical mechanism of the main roof instability, to reveal the disaster-causing factors, and to explore the determination method of support resistance to ensure the stability of working face. In this study, according to the structural characteristics of the main roof cutting and instability along the coal wall during the first weighting, the system mechanical model consists of the main roof, immediate roof and support is established. This model is based on catastrophe theory and the energy conservation principle analysis method, as well as through considering the nonlinear characteristics of the failure and instability of the immediate roof. A new analysis method was proposed to study the mechanical mechanism of roof cutting along a coal wall during the first weighting. The instability process was analyzed theoretically, and the sufficient and necessary conditions for the sudden instability of the system were successfully obtained. The main influencing factors of the system catastrophe, such as stiffness ratio K, rock material parameters, external load Q, and weighting interval l are discussed. We also compared and analyzed the effects of elastic modulus E, support stiffness k 1 , cross-sectional area a, and immediate roof thickness H on the stiffness ratio K. After further derivation, the calculation equation of support resistance was obtained, and the main influencing factors of step sinking were analyzed. Finally, the rationality of the theoretical analysis results were verified by an engineering example, and the reasonable value of the support is given. This paper is organized as follows: the structural characteristics and cutting process of the main roof during the first weighting in the working face of large mining height in the shallow coal seam are described in Section 2. A catastrophe theoretical analysis model consisting of main roof, immediate roof, and support is given in Section 3. The catastrophe instability condition of the immediate roof, the calculation method of the support resistance, and the step sinking are deduced by catastrophe theory in Section 4. Furthermore, the theoretical results are verified in Section 5. Conclusions are drawn in Section 6.

Working Face Overview and Strata Behavior Monitoring
To analyze the roof structural characteristics of thin bedrock under thick loose layer during shallow coal seam mining, we took the no. 12401 panel of 5 −2 coal seam in Shendong mining area as the engineering background. The length of working face was 260 m with a dip angle of 1-3 • , and the average burial depth was 120 m. The thickness of the bedrock was 70 m on average. The thickness of the loose layer was 50 m on average. The thickness of the bedrock fluctuated in a small range. The roof was mostly mudstone, fine-grained sandstone, and siltstone, with an average compressive strength of 23.15 MPa. The floor was mainly siltstone, and the rock mass was complete with an average  Table 1. The working face adopted 150 sets of ZY18000/32/70D electro-hydraulic control shield-type hydraulic supports. The measured first weighting interval of the main roof was about 54 m, and the average of periodic weighting interval was 15.8 m. 5-2 coal seam 6.2 120 Figure 2 shows the monitoring curve of support resistance during the first weighting of working face, and Figure 3 shows the variation curve of working resistance of supports during the first weighting. The working resistance of the support was small and the change range was not large, most of which was 25.2-30.0 MPa, found by analyzing the curves in Figures 2 and 3. The support resistance of no. 55-120 in the middle of the working face increased rapidly when the working face advanced about 53 m, and the maximum value reached 47.3 MPa (red area in Figure 4). There was frequent opening of the safety valve, rapid roof sinking, and great danger of support crushing. The rib fall of  Table 1. The working face adopted 150 sets of ZY18000/32/70D electro-hydraulic control shieldtype hydraulic supports. The measured first weighting interval of the main roof was about 54 m, and the average of periodic weighting interval was 15.8 m. 5-2 coal seam 6.2 120 Figure 2 shows the monitoring curve of support resistance during the first weighting of working face, and Figure 3 shows the variation curve of working resistance of supports during the first weighting. The working resistance of the support was small and the change range was not large, most of which was 25.2-30.0 MPa, found by analyzing the curves in Figures 2 and 3. The support resistance of no. 55-120 in the middle of the working face increased rapidly when the working face advanced about 53 m, and the maximum value reached 47.3 MPa (red area in Figure 4). There was frequent opening of the safety valve, rapid roof sinking, and great danger of support crushing. The rib fall of  Table 1. The working face adopted 150 sets of ZY18000/32/70D electro-hydraulic control shieldtype hydraulic supports. The measured first weighting interval of the main roof was about 54 m, and the average of periodic weighting interval was 15.8 m. 5-2 coal seam 6.2 120 Figure 2 shows the monitoring curve of support resistance during the first weighting of working face, and Figure 3 shows the variation curve of working resistance of supports during the first weighting. The working resistance of the support was small and the change range was not large, most of which was 25.2-30.0 MPa, found by analyzing the curves in Figures 2 and 3. The support resistance of no. 55-120 in the middle of the working face increased rapidly when the working face advanced about 53 m, and the maximum value reached 47.3 MPa (red area in Figure 4). There was frequent opening of the safety valve, rapid roof sinking, and great danger of support crushing. The rib fall of  Table 1. The working face adopted 150 sets of ZY18000/32/70D electro-hydraulic control shieldtype hydraulic supports. The measured first weighting interval of the main roof was about 54 m, and the average of periodic weighting interval was 15.8 m. 5-2 coal seam 6.2 120 Figure 2 shows the monitoring curve of support resistance during the first weighting of working face, and Figure 3 shows the variation curve of working resistance of supports during the first weighting. The working resistance of the support was small and the change range was not large, most of which was 25.2-30.0 MPa, found by analyzing the curves in Figures 2 and 3. The support resistance of no. 55-120 in the middle of the working face increased rapidly when the working face advanced about 53 m, and the maximum value reached 47.3 MPa (red area in Figure 4). There was frequent opening of the safety valve, rapid roof sinking, and great danger of support crushing. The rib fall of  Table 1. The working face adopted 150 sets of ZY18000/32/70D electro-hydraulic control shieldtype hydraulic supports. The measured first weighting interval of the main roof was about 54 m, and the average of periodic weighting interval was 15.8 m. 5-2 coal seam 6.2 120 Figure 2 shows the monitoring curve of support resistance during the first weighting of working face, and Figure 3 shows the variation curve of working resistance of supports during the first weighting. The working resistance of the support was small and the change range was not large, most of which was 25.2-30.0 MPa, found by analyzing the curves in Figures 2 and 3. The support resistance of no. 55-120 in the middle of the working face increased rapidly when the working face advanced about 53 m, and the maximum value reached 47.3 MPa (red area in Figure 4). There was frequent opening of the safety valve, rapid roof sinking, and great danger of support crushing. The rib fall of  Table 1. The working face adopted 150 sets of ZY18000/32/70D electro-hydraulic control shieldtype hydraulic supports. The measured first weighting interval of the main roof was about 54 m, and the average of periodic weighting interval was 15.8 m. 5-2 coal seam 6.2 120 Figure 2 shows the monitoring curve of support resistance during the first weighting of working face, and Figure 3 shows the variation curve of working resistance of supports during the first weighting. The working resistance of the support was small and the change range was not large, most of which was 25.2-30.0 MPa, found by analyzing the curves in Figures 2 and 3. The support resistance of no. 55-120 in the middle of the working face increased rapidly when the working face advanced about 53 m, and the maximum value reached 47.3 MPa (red area in Figure 4). There was frequent opening of the safety valve, rapid roof sinking, and great danger of support crushing. The rib fall of  Table 1. The working face adopted 150 sets of ZY18000/32/70D electro-hydraulic control shieldtype hydraulic supports. The measured first weighting interval of the main roof was about 54 m, and the average of periodic weighting interval was 15.8 m. 5-2 coal seam 6.2 120 Figure 2 shows the monitoring curve of support resistance during the first weighting of working face, and Figure 3 shows the variation curve of working resistance of supports during the first weighting. The working resistance of the support was small and the change range was not large, most of which was 25.2-30.0 MPa, found by analyzing the curves in Figures 2 and 3. The support resistance of no. 55-120 in the middle of the working face increased rapidly when the working face advanced about 53 m, and the maximum value reached 47.3 MPa (red area in Figure 4). There was frequent opening of the safety valve, rapid roof sinking, and great danger of support crushing. The rib fall of  Table 1. The working face adopted 150 sets of ZY18000/32/70D electro-hydraulic control shieldtype hydraulic supports. The measured first weighting interval of the main roof was about 54 m, and the average of periodic weighting interval was 15.8 m. 5-2 coal seam 6.2 120 Figure 2 shows the monitoring curve of support resistance during the first weighting of working face, and Figure 3 shows the variation curve of working resistance of supports during the first weighting. The working resistance of the support was small and the change range was not large, most of which was 25.2-30.0 MPa, found by analyzing the curves in Figures 2 and 3. The support resistance of no. 55-120 in the middle of the working face increased rapidly when the working face advanced about 53 m, and the maximum value reached 47.3 MPa (red area in Figure 4). There was frequent opening of the safety valve, rapid roof sinking, and great danger of support crushing. The rib fall of 4.69 96.63 Mudstone average of 7.2 m, and the designed mining height was 7.0 m. The comprehensive histogram is shown in Table 1. The working face adopted 150 sets of ZY18000/32/70D electro-hydraulic control shieldtype hydraulic supports. The measured first weighting interval of the main roof was about 54 m, and the average of periodic weighting interval was 15.8 m. 5-2 coal seam 6.2 120 Figure 2 shows the monitoring curve of support resistance during the first weighting of working face, and Figure 3 shows the variation curve of working resistance of supports during the first weighting. The working resistance of the support was small and the change range was not large, most of which was 25.2-30.0 MPa, found by analyzing the curves in Figures 2 and 3. The support resistance of no. 55-120 in the middle of the working face increased rapidly when the working face advanced about 53 m, and the maximum value reached 47.3 MPa (red area in Figure 4). There was frequent opening of the safety valve, rapid roof sinking, and great danger of support crushing. The rib fall of 3.68 103.31 Fine sandstone in Table 1. The working face adopted 150 sets of ZY18000/32/70D electro-hydraulic control shieldtype hydraulic supports. The measured first weighting interval of the main roof was about 54 m, and the average of periodic weighting interval was 15.8 m. 5-2 coal seam 6.2 120 Figure 2 shows the monitoring curve of support resistance during the first weighting of working face, and Figure 3 shows the variation curve of working resistance of supports during the first weighting. The working resistance of the support was small and the change range was not large, most of which was 25.2-30.0 MPa, found by analyzing the curves in Figures 2 and 3. The support resistance of no. 55-120 in the middle of the working face increased rapidly when the working face advanced about 53 m, and the maximum value reached 47.3 MPa (red area in Figure 4). There was frequent opening of the safety valve, rapid roof sinking, and great danger of support crushing. The rib fall of 4.72 108.03 Mudstone in Table 1. The working face adopted 150 sets of ZY18000/32/70D electro-hydraulic control shieldtype hydraulic supports. The measured first weighting interval of the main roof was about 54 m, and the average of periodic weighting interval was 15.8 m. 5-2 coal seam 6.2 120 Figure 2 shows the monitoring curve of support resistance during the first weighting of working face, and Figure 3 shows the variation curve of working resistance of supports during the first weighting. The working resistance of the support was small and the change range was not large, most of which was 25.2-30.0 MPa, found by analyzing the curves in Figures 2 and 3. The support resistance of no. 55-120 in the middle of the working face increased rapidly when the working face advanced about 53 m, and the maximum value reached 47.3 MPa (red area in Figure 4). There was frequent opening of the safety valve, rapid roof sinking, and great danger of support crushing. The rib fall of 5.86 113.9 5-2 coal seam type hydraulic supports. The measured first weighting interval of the main roof was about 54 m, and the average of periodic weighting interval was 15.8 m. 5-2 coal seam 6.2 120 Figure 2 shows the monitoring curve of support resistance during the first weighting of working face, and Figure 3 shows the variation curve of working resistance of supports during the first weighting. The working resistance of the support was small and the change range was not large, most of which was 25.2-30.0 MPa, found by analyzing the curves in Figures 2 and 3. The support resistance of no. 55-120 in the middle of the working face increased rapidly when the working face advanced about 53 m, and the maximum value reached 47.3 MPa (red area in Figure 4). There was frequent opening of the safety valve, rapid roof sinking, and great danger of support crushing. The rib fall of 6.2 120 Figure 2 shows the monitoring curve of support resistance during the first weighting of working face, and Figure 3 shows the variation curve of working resistance of supports during the first weighting. The working resistance of the support was small and the change range was not large, most of which was 25.2-30.0 MPa, found by analyzing the curves in Figures 2 and 3. The support resistance of no. 55-120 in the middle of the working face increased rapidly when the working face advanced about 53 m, and the maximum value reached 47.3 MPa (red area in Figure 4). There was frequent opening of the safety valve, rapid roof sinking, and great danger of support crushing. The rib fall of the coal wall was seriously cracked, and the average rib fall depth was 200-400 mm. After that, the weighting entered a continuous stage.
of which was 25.2-30.0 MPa, found by analyzing the curves in Figures 2 and 3. The support resistance of no. 55-120 in the middle of the working face increased rapidly when the working face advanced about 53 m, and the maximum value reached 47.3 MPa (red area in Figure 4). There was frequent opening of the safety valve, rapid roof sinking, and great danger of support crushing. The rib fall of the coal wall was seriously cracked, and the average rib fall depth was 200-400mm. After that, the weighting entered a continuous stage.   The measured first weighting interval of the main roof was about 54 m, which was larger than that of the general mining height. The support resistance increased sharply during first weighting, and the weighting in the middle of the working face was strong. The dynamic loading phenomenon was obvious, with the average dynamic load coefficient being 1.81. The roof moved violently, causing the safety valve of the support to open frequently, and there was great danger of support crushing. Therefore, it was necessary to analyze the applicability of the hydraulic support with a rated working resistance of 18,000 kN under the conditions of large mining height.   The measured first weighting interval of the main roof was about 54 m, which was larger than that of the general mining height. The support resistance increased sharply during first weighting, and the weighting in the middle of the working face was strong. The dynamic loading phenomenon was obvious, with the average dynamic load coefficient being 1.81. The roof moved violently, causing the safety valve of the support to open frequently, and there was great danger of support crushing. Therefore, it was necessary to analyze the applicability of the hydraulic support with a rated working resistance of 18,000 kN under the conditions of large mining height. Figure 5 shows the first weighting process of working face under the condition of 7 m mining height of the shallow coal seam obtained from the similar material model test. The analysis of the simulation results showed that as the space of the goaf increased continuously with the advancement of the working face, the immediate roof had less filling to goaf (Figure 5a), causing the mininginduced fracture to extend continuously to the upper part of the goaf. The upper part of the original main roof formed the separation layer and working face pressure, and the previously broken rock mass could not form an articulated structure with the broken lower main roof, and thus the lower main roof formed a cantilever beam structure (Figure 5b). The upper rock layer broke continuously under the action of mining unloading with the advancement of the working face, and the thickness of the cantilever beam structure increased (Figure 5c) until the upper main roof formed a stable  Figure 4. The subsidence was between 210 mm and 235 mm, and the average was 224 mm. The sinking amount was between 210 mm and 235 mm, with an average of 224 mm. The monitoring results showed that the roof damage caused by the first weighting directly affected the ground surface, and there were obvious step subsidence and dynamic load phenomenon in the process of weighting.

Characteristics of Roof First Weighting
The measured first weighting interval of the main roof was about 54 m, which was larger than that of the general mining height. The support resistance increased sharply during first weighting, and the weighting in the middle of the working face was strong. The dynamic loading phenomenon was obvious, with the average dynamic load coefficient being 1.81. The roof moved violently, causing the safety valve of the support to open frequently, and there was great danger of support crushing. Therefore, it was necessary to analyze the applicability of the hydraulic support with a rated working resistance of 18,000 kN under the conditions of large mining height. Figure 5 shows the first weighting process of working face under the condition of 7 m mining height of the shallow coal seam obtained from the similar material model test. The analysis of the simulation results showed that as the space of the goaf increased continuously with the advancement of the working face, the immediate roof had less filling to goaf (Figure 5a), causing the mining-induced fracture to extend continuously to the upper part of the goaf. The upper part of the original main roof formed the separation layer and working face pressure, and the previously broken rock mass could not form an articulated structure with the broken lower main roof, and thus the lower main roof formed a cantilever beam structure (Figure 5b). The upper rock layer broke continuously under the action of mining unloading with the advancement of the working face, and the thickness of the cantilever beam structure increased (Figure 5c) until the upper main roof formed a stable articulated structure due to rotation or sliding instability. In this process, a new immediate roof was formed, which was different from the original immediate roof, being composed of multiple layers of rock, as is shown in Figure 5d.

Characteristics of Roof First Weighting
Appl. Sci. 2020, 10, x FOR PEER REVIEW 7 of 23 formed, which was different from the original immediate roof, being composed of multiple layers of rock, as is shown in Figure 5d. In order to analyze the characteristics of the first weighting in the shallow coal seam with large mining height, we constructed a numerical analysis model with a mining height of 4-7m [41][42][43]. The simulation results are shown in Figure 6, with there being different degrees of sliding instability and violent strata behavior in the working face of different mining heights. Similar to the physical simulation, the roof successively appeared as immediate roof (the mining height shown in Figure 6a is 4 m), cantilever beam structure (the mining height shown in Figure 6b is 5 m), and composite cantilever beam structure (the mining height in Figure 6c,d is 6 m and 6 m, respectively) with the increase of mining height. The space of the goaf was insufficient in terms of the increase of mining height, and thus the collapsed main roof could not form an articulated structure, continuing to collapse into the goaf until the goaf was filled to support the main roof. The collapsed main roof and the original immediate roof formed a new immediate roof of composite cantilever beam. In the process of the pressure of the upper main roof, the load formed by the sliding instability of the main roof along the coal wall acted on the immediate roof. In this process, the support and the immediate roof constituted a system that jointly bore the main roof loads. Considering the nonlinear deformation characteristics of the immediate roof, if the immediate roof broke under the pressure and lost its load-Immediate roof cut off Up immediate roof

Cantilever beam
Composite cantilever beam In order to analyze the characteristics of the first weighting in the shallow coal seam with large mining height, we constructed a numerical analysis model with a mining height of 4-7m [41][42][43]. The simulation results are shown in Figure 6, with there being different degrees of sliding instability and violent strata behavior in the working face of different mining heights. Similar to the physical simulation, the roof successively appeared as immediate roof (the mining height shown in Figure 6a is 4 m), cantilever beam structure (the mining height shown in Figure 6b is 5 m), and composite cantilever beam structure (the mining height in Figure 6c,d is 6 m and 6 m, respectively) with the increase of mining height. The space of the goaf was insufficient in terms of the increase of mining height, and thus the collapsed main roof could not form an articulated structure, continuing to collapse into the goaf until the goaf was filled to support the main roof. The collapsed main roof and the original immediate Appl. Sci. 2020, 10, 5408 8 of 23 roof formed a new immediate roof of composite cantilever beam. In the process of the pressure of the upper main roof, the load formed by the sliding instability of the main roof along the coal wall acted on the immediate roof. In this process, the support and the immediate roof constituted a system that jointly bore the main roof loads. Considering the nonlinear deformation characteristics of the immediate roof, if the immediate roof broke under the pressure and lost its load-bearing effect, the force borne by the support would suddenly increase, wherein it would be easy to induce the roof cutting (Figure 6c,d). Therefore, the stability of the immediate roof plays an important role in the safety of the working face.

Mechanical Model Analysis
In the mining process, the sliding instability of the main roof was the biggest threat to the working face during the first weighting. The first weighting process is shown in Figure 7. The dead weight of the main roof and the load of overlying thick loose layer Q were transmitted to the immediate roof and support through the main roof, resulting in the compression deformation of the immediate roof and support (as shown in Figure 7b). If the system composed of the immediate roof and support could not form an effective support for the main roof, the main roof would rotate and deform to compress the immediate roof rock and support. The immediate roof would be crushed if the deformation increased to a certain critical value, causing the main roof to slide along the coal wall, and step sinking would occur at the main roof and the overlying thick loose layer (as shown in Figure  7c), resulting in coal mine accidents. Therefore, it was important to analyze the deformation and instability of the main roof immediate roof support system and provide measures to control the sliding instability of the coal wall when the main roof is broken and unstable, so that the main roof can rotate to the goaf for normal pressure.
Step sinking

Composite cantilever beam
Coal wall cutting

Cantilever beam
Cantilever beam

Sliding instability
Composite cantilever

Mechanical Model Analysis
In the mining process, the sliding instability of the main roof was the biggest threat to the working face during the first weighting. The first weighting process is shown in Figure 7. The dead weight of the main roof and the load of overlying thick loose layer Q were transmitted to the immediate roof and support through the main roof, resulting in the compression deformation of the immediate roof and support (as shown in Figure 7b). If the system composed of the immediate roof and support could not form an effective support for the main roof, the main roof would rotate and deform to compress the immediate roof rock and support. The immediate roof would be crushed if the deformation increased to a certain critical value, causing the main roof to slide along the coal wall, and step sinking would occur at the main roof and the overlying thick loose layer (as shown in Figure 7c), resulting in coal mine accidents. Therefore, it was important to analyze the deformation and instability of the main roof immediate roof support system and provide measures to control the sliding instability of the coal wall when the main roof is broken and unstable, so that the main roof can rotate to the goaf for normal pressure. From the perspective of mechanical theory research, this paper studied the mechanical mechanism of the main roof sliding instability, establishing a mechanical system composed of the main roof, immediate roof, and support. Generally, the elastomer can be used to represent the support, and the immediate roof can be used as a variable medium. Because of the fracture at both ends of the main roof, it can be simplified as the hinged connection with the front and rear rock mass. During the process of the main roof sliding instability, the immediate roof and support bore the weight of the roof and the load of the loose layer Q, with the subsidence of the instability rock block A being set as Δh. As the subsidence of the main roof compressed the immediate roof and the support, the step sinking of the main roof was thus equal to the total displacement of the support and immediate roof. The coal gangue in the goaf was scattered, and thus the absorbed energy could be ignored. In the process of the main roof and follow-up rock block sinking, the gravitational potential energy decreased, the potential energy of the support and the immediate roof increased, and the friction of the hinge surface worked with energy conservation.
Combined with the above analysis, the mechanical analysis model was established, as shown in Figure 8, where qs is the weight of overlying loose layer and main roof, lA is the length of rock block A, u is the compression displacement of immediate roof rock mass, w is the compression displacement of support, a is the total displacement of the support-immediate roof, and k1 is the stiffness of support.  From the perspective of mechanical theory research, this paper studied the mechanical mechanism of the main roof sliding instability, establishing a mechanical system composed of the main roof, immediate roof, and support. Generally, the elastomer can be used to represent the support, and the immediate roof can be used as a variable medium. Because of the fracture at both ends of the main roof, it can be simplified as the hinged connection with the front and rear rock mass. During the process of the main roof sliding instability, the immediate roof and support bore the weight of the roof and the load of the loose layer Q, with the subsidence of the instability rock block A being set as ∆h.
As the subsidence of the main roof compressed the immediate roof and the support, the step sinking of the main roof was thus equal to the total displacement of the support and immediate roof. The coal gangue in the goaf was scattered, and thus the absorbed energy could be ignored. In the process of the main roof and follow-up rock block sinking, the gravitational potential energy decreased, the potential energy of the support and the immediate roof increased, and the friction of the hinge surface worked with energy conservation.
Combined with the above analysis, the mechanical analysis model was established, as shown in Figure 8, where q s is the weight of overlying loose layer and main roof, l A is the length of rock block A, u is the compression displacement of immediate roof rock mass, w is the compression displacement of support, a is the total displacement of the support-immediate roof, and k 1 is the stiffness of support. From the perspective of mechanical theory research, this paper studied the mechanical mechanism of the main roof sliding instability, establishing a mechanical system composed of the main roof, immediate roof, and support. Generally, the elastomer can be used to represent the support, and the immediate roof can be used as a variable medium. Because of the fracture at both ends of the main roof, it can be simplified as the hinged connection with the front and rear rock mass. During the process of the main roof sliding instability, the immediate roof and support bore the weight of the roof and the load of the loose layer Q, with the subsidence of the instability rock block A being set as Δh. As the subsidence of the main roof compressed the immediate roof and the support, the step sinking of the main roof was thus equal to the total displacement of the support and immediate roof. The coal gangue in the goaf was scattered, and thus the absorbed energy could be ignored. In the process of the main roof and follow-up rock block sinking, the gravitational potential energy decreased, the potential energy of the support and the immediate roof increased, and the friction of the hinge surface worked with energy conservation.
Combined with the above analysis, the mechanical analysis model was established, as shown in Figure 8, where qs is the weight of overlying loose layer and main roof, lA is the length of rock block A, u is the compression displacement of immediate roof rock mass, w is the compression displacement of support, a is the total displacement of the support-immediate roof, and k1 is the stiffness of support.

Constitutive Relationship of Immediate Roof Rock Mass
The constitutive relation of the immediate roof rock mass is a nonlinear relationship with a softening property [30,44]. The stress-strain relationship of rock materials is discussed in relevant literature [45,46], giving the relation between the stress σ and strain ε.
where E is the initial value of the elastic modulus, and the integral ε 0 ϕ(t)dt is the damage parameter that is related to the distribution density of defects in rock materials. When the defects in rock conform to Poisson distribution, the nonlinear constitutive relation can be expressed as For the immediate roof rock mass with section A and height H, Equation (2) can be used as an expression of the relationship between the load R and deformation u (as shown in Figure 9) where λ = EA/H is the initial stiffness of the rock mass and u 0 is the strain value corresponding to the peak load. The nonlinear curve of Equation (3) has an inflection point at the strain u 1 = 2u 0 , and the absolute value of the corresponding slope is λ 1 = λe −2 .

Constitutive Relationship of Immediate Roof Rock Mass
The constitutive relation of the immediate roof rock mass is a nonlinear relationship with a softening property [30,44]. The stress-strain relationship of rock materials is discussed in relevant literature [45,46], giving the relation between the stress σ and strain ε.
where E is the initial value of the elastic modulus, and the integral dt  ) t ( 0   is the damage parameter that is related to the distribution density of defects in rock materials. When the defects in rock conform to Poisson distribution, the nonlinear constitutive relation can be expressed as For the immediate roof rock mass with section A and height H, Equation (2) can be used as an expression of the relationship between the load R and deformation u (as shown in Figure 9) where λ = EA/H is the initial stiffness of the rock mass and u0 is the strain value corresponding to the peak load. The nonlinear curve of Equation (3) has an inflection point at the strain u1 = 2u0, and the absolute value of the corresponding slope is λ1 = λe −2 .

System Potential Function
The total potential energy of the mechanical system composed of the main roof, immediate roof, and support can be expressed as where WJ is the work done by external force Q, since gravity and friction work have the same path; the resultant force Q can be expressed as

System Potential Function
The total potential energy of the mechanical system composed of the main roof, immediate roof, and support can be expressed as where W J is the work done by external force Q, since gravity and friction work have the same path; the resultant force Q can be expressed as Q = Q A + Q s − F C , where Q A is the gravity of the main roof rock block, Q A = h j l A ρg; l A is the length of the broken block A; h j is the main roof height, m; ρg is the main roof volume weight, kN/m 3 ; and Q s is the load of overlying thick loose layer. According to the calculation principle of Terzaghi earth pressure, the weight of load layer is Q s = ρ 1 gl 2 2λ z tan ϕ , where ρ 1 g is the average volume weight of the load layer, kN/m 3 ; ϕ is the internal friction angle of the load layer, • ; and λ z is the lateral stress ratio of the load layer, λ z = 1-sinϕ. F C is the friction caused by sliding instability, and the horizontal force on the rock contact hinge is T = (1+2i tan θ 1max ) tan ϕ 2(1−tan ϕ tan θ 1max ) Q A , where θ 1max is the maximum value of rock block rotating angle, tan ϕ is the friction coefficient between rock block Appl. Sci. 2020, 10, 5408 11 of 23 and gangue, and i is the main roof rock block size. The main roof slipped and lost stability during the first weighting, which made it become a balance structure again after the roof cutting; therefore, it can be considered that the main roof was a horizontal structure. At this point θ 1max = 0, the friction shear stress on the contact hinge can be expressed as F C = T tan ϕ = tan ϕ tan ϕ 2 Q A , kN [47,48]. According to the conservation of energy, the work done by the external force during the catastrophe instability of immediate roof can be expressed as where ∆h is the step sinking of roof, m. The strain energy of the immediate roof is U z = u 0 λue −u/u 0 du, the compression capacity of support is w = a − u, and the elastic potential energy of the support is U s = 1 2 k 1 (a − u) 2 . Because the step sinking of the roof is equal to the total displacement of the support-immediate roof system, ∆h = w + u. Therefore, combining W J , U Z , and U S into (4), the total potential energy of the system can be expressed as

Catastrophe Analysis
Taking the compression displacement u of immediate roof rock mass as the state variable, according to the cusp catastrophe theory, obtaining of the equilibrium surface M can be expressed as Then, the singular point set equation can be expressed as The equilibrium surface M satisfies V = 0 at the cusp, and the cusp equation can be obtained as At the cusp, It was found that the cusp was the inflection point of the constitutive curve of rock mass. In order to transform the cusp catastrophe model to a standard form, we carried out a Taylor series expansion at the cusp, and the first three terms were truncated as Introducing dimensionless parameters x = u−u 1 u 1 , and substitute u 1 = 2u 0 , the above equation can be simplified to obtain the equation of equilibrium surface in the standard form of cusp catastrophe, which can be expressed as The parameter K is the ratio of the support stiffness to the slope of the constitutive relation curve of the immediate roof at the inflection point, which is called the stiffness ratio. ζ is the total displacement parameter, which is related to the total displacement of the support and the immediate roof. According to Equations (13) and (14), the control variables m and n of the system are related to the stiffness ratio K, the total displacement parameter ζ, the external load Q, and the first weighting interval l of the main roof.
The bifurcation set equation of cusp catastrophe with m and n as control parameters can be expressed as The bifurcation set Equation (17) is the critical point of the system instability. Only when m ≤ 0 will the system cross the bifurcate set and cause a sudden jump of deformation, and thus m ≤ 0 is the necessary condition of the system instability. According to Equation (13), the stiffness ratio K must be less than or equal to 1, that is According to the catastrophe theory, only when m and n satisfy the bifurcation set Equation (17) will the system have a sudden jump in instability. Therefore, the bifurcation set equation is a sufficient condition for the sudden jump in instability of the system. According to the cusp catastrophe model in Figure 10, when the bifurcation set crossed from the right branch (n > 0) to the left branch (n < 0), the corresponding point was in an unstable state, and the state variable x was a sudden jump. The fracture instability occurred when the deformation value of the immediate roof increased instantly. The sufficient and necessary conditions for the catastrophic instability of the mechanical system were obtained from Equation (17), which can be expressed as The lower lobe of the equilibrium surface in Figure 11 represents the accumulation process of system energy; at this time, the elastic potential energy increased gradually and the system was at the critical point of stable state. With the continuous mining work, the deformation of support and immediate roof rock increased gradually, and the compression displacement of immediate roof increased gradually, with the stress-strain curve developing to the critical state of the middle lobe. When the upper lobe was reached, the compression displacement of immediate roof u = 2u 0 = u 1 , and the deformation increased suddenly and instantaneously, resulting in system instability, which led to the main roof cutting down along the coal wall, wherein the system reached a new equilibrium state. The whole process of the sudden jump deformation of the immediate roof was the catastrophic process of the system state. The lower lobe of the equilibrium surface in Figure 11 represents the accumulation process of system energy; at this time, the elastic potential energy increased gradually and the system was at the critical point of stable state. With the continuous mining work, the deformation of support and immediate roof rock increased gradually, and the compression displacement of immediate roof increased gradually, with the stress-strain curve developing to the critical state of the middle lobe. When the upper lobe was reached, the compression displacement of immediate roof u = 2u0 = u1, and the deformation increased suddenly and instantaneously, resulting in system instability, which led to the main roof cutting down along the coal wall, wherein the system reached a new equilibrium state. The whole process of the sudden jump deformation of the immediate roof was the catastrophic process of the system state.   The lower lobe of the equilibrium surface in Figure 11 represents the accumulation process of system energy; at this time, the elastic potential energy increased gradually and the system was at the critical point of stable state. With the continuous mining work, the deformation of support and immediate roof rock increased gradually, and the compression displacement of immediate roof increased gradually, with the stress-strain curve developing to the critical state of the middle lobe. When the upper lobe was reached, the compression displacement of immediate roof u = 2u0 = u1, and the deformation increased suddenly and instantaneously, resulting in system instability, which led to the main roof cutting down along the coal wall, wherein the system reached a new equilibrium state. The whole process of the sudden jump deformation of the immediate roof was the catastrophic process of the system state. Figure 11. The state jump when leaping over the bifurcation point sets. Figure 11. The state jump when leaping over the bifurcation point sets.
Since the parameter K was the ratio of the support stiffness k 1 to the slope λ 1 of the stress-strain curve of the immediate roof at the inflection point, and the stiffness ratio K was only related to the internal properties of the material, the internal characteristics of the material were therefore the necessary conditions for the system to have a catastrophic change. The main influence factors of the stiffness ratio K were analyzed by Equations (18) and (19). The influence of the support stiffness k 1 , the elastic modulus E, the immediate roof cross-sectional area A, and the immediate roof thickness on the stiffness ratio K are shown in Figure 12. It can be found from the curves that the stiffness ratio K was positively related to the support stiffness k 1 and the immediate roof thickness, and was negatively related to the elastic modulus E and the immediate roof cross-sectional area A. The support resistance and the supporting capacity of the roof increased with the increase of the support stiffness k 1 , which was conducive to the stability of the roof and the system (Figure 12a). Therefore, increasing the support stiffness within a reasonable range played an important role in maintaining the stability of the roof. Because the immediate roof rock mass was composed of structural plane and joint fissures, the more complete the immediate roof, the greater the elastic modulus E and the stiffness λ 1 , and the more unstable the system (Figure 12b). easy to induce the catastrophic instability of the working face, which was also the inducement of the roof cutting and sinking during the first weighting ( Figure 12c). As the combined immediate roof was generally of inverted step structure, and the more the number of immediate roof layers, the greater the length, therefore, under the condition of large mining height, the increase of A played a positive role in ensuring the stability of working face (Figure 12d). Comparing the four curves, we found that the influence degree of each parameter on the stiffness ratio was elastic modulus E > support stiffness k1 > cross-sectional area A > immediate roof thickness H. Figure 12. Analysis of influence factors of stiffness ratio K. Figure 12. Analysis of influence factors of stiffness ratio K.
In addition, the cross-sectional area A and thickness H of the immediate roof also had a great influence on the stiffness ratio. The immediate roof of the large mining height working face was composed of multiple layers of rock, which was thicker than the immediate roof of normal mining height. Because the stiffness ratio was directly proportional to the immediate roof thickness H, it was easy to induce the catastrophic instability of the working face, which was also the inducement of the roof cutting and sinking during the first weighting ( Figure 12c). As the combined immediate roof was generally of inverted step structure, and the more the number of immediate roof layers, the greater the length, therefore, under the condition of large mining height, the increase of A played a positive role in ensuring the stability of working face (Figure 12d). Comparing the four curves, we found that the influence degree of each parameter on the stiffness ratio was elastic modulus E > support stiffness k 1 > cross-sectional area A > immediate roof thickness H.
From the sufficient condition expression of the system catastrophe instability in Equation (19), we found that the system instability was also related to the external load Q and the length l A of the fracture block, i.e., the first weighting distances l. With the increase of load Q, the load on the system increased and the system stability decreased gradually. When the load reached a certain value, the catastrophe instability of the system occurred. The larger the first weighting distances l, the lower the system stability. When the first weighting interval exceeded a certain value, the sufficient conditions for the catastrophe instability were met, and the catastrophe instability of the system occurred. Therefore, the external load Q and the first weighting interval l were the external determinants of the system instability, which were large enough to change the stability of the system.
Through the analysis of the necessary and sufficient condition Equation (19), we found that the instability of the system was not only related to its internal characteristics (immediate roof elastic modulus E, thickness h, cross-sectional area A, stress-strain relationship and support stiffness k 1 ), but also related to the deadweight of the main roof, the load of the overlying thick loose layer, and the weighting distance L. For the immediate roof, the stronger the strain softening property, the larger the elastic modulus E and corresponding λ 1 ; the smaller the stiffness ratio K, the easier it was for the system lose stability. The larger the overburden load, the longer the first weighting interval of the main roof, and the easier it would be to cross the bifurcation point set, and therefore, the easier it would be for the system to lose stability.
The overlying thick loose sand layer in shallow coal seam made the load all act on the roof bedrock. During the advancement of the working face of large mining height, the main roof overlying load increased gradually, which would inevitably lead to the probability of catastrophe instability. In addition, when the support stiffness was certain, the influence of the support on the system stability and the load on the support were all determined by the mechanical properties and material parameters of the immediate roof. The lithology compositions and integrity of the immediate roof had a decisive influence on the occurrence of a catastrophe. Therefore, the catastrophe instability of the immediate roof determined the cut-off of the main roof.

Calculation of Step Sinking
When the system satisfied the necessary and sufficient conditions for catastrophe instability, we obtained three real roots of Equation (12). When crossing the set of bifurcation points, the state variable x would suddenly jump ( Figure 10). The jump amount can be expressed as The corresponding sudden jump compression of the immediate roof before and after system instability can be expressed as It can be found from Equation (21) that the sudden jump was only related to the characteristics of the immediate roof and the stiffness ratio of support. The amount of sudden jump was determined by the properties of the materials in the system, and had nothing to do with other factors.
Taking m, n, Equations (13) and (14), and λ = EA/H into Equation (14), we can formulate the total displacement a during the sudden jump of deformation as From ξ = a−u 1 u 1 and u 1 = 2u 0 (ξ takes the larger value), we can obtain The total displacement of the system when the sudden jump occurs is the sinking value of the main roof, that is, ∆h = a, and thus the sinking value of the key block can be expressed as In which ∆h is the step sinking value of the main roof when the catastrophe instability of immediate roof occurs, which is formed under the action of external conditions. From the constitutive relation curve of the rock in Figure 3, it can be found that the step sinking value was the allowable subsidence of the roof when the sliding instability occurs. If the step sinking value exceeded ∆h, the main roof would suddenly jump and lose stability, resulting in a large range of roof cutting and support crushing accidents. Therefore, the prediction of roof subsidence can be realized by calculating ∆h.

Calculation of Support Resistance
Before the sudden jump of the immediate roof, the compression displacement of the support is the difference between the displacement of the system and the immediate roof, and then the compression displacement of the support during the catastrophe instability can be expressed as When the catastrophe occurs, the compression displacement of the support is the support-surrounding rock condition needed for the sudden jump. Combined with the bifurcation set Equation (14), which is the critical point of the system instability, when m ≤ 0, the system will cross the bifurcation to generate the sudden jump of deformation. Combined with Equation (18), the support stiffness value of the catastrophic instability of the system can be obtained as Therefore, in order to keep the system stable without the sudden jump deformation of the immediate roof, the support stiffness k 1 should be greater than EA Le 2 , and the minimum support resistance required for system stability can be expressed as Considering the support efficiency, the support resistance when the catastrophe instability occurs can be expressed as where µ is the support efficiency of the support, equal to 0.9. Through the interaction between the support and the main roof, the value of the minimum support resistance can be obtained to ensure the stability of the working face, which can be used as the basis for support selection.

Analysis of Influencing Factors
In order to deeply analyze the main influence factors of roof sinking ∆h (total displacement a) during the catastrophe instability of the system, we analyzed the main parameters, combined with the engineering geological conditions of the no. 12401 fully mechanized coal mining face in Shendong mining area. Figures 13-15 show the corresponding roof subsidence curve when the overburden load q, the thickness H of the immediate roof, and the support stiffness k 1 changed.
When the thickness of the immediate roof increased to a certain value, the step sinking value remained basically stable, indicating that the larger the thickness of the rock layer, the more stable the immediate roof, as the overlying rock layer could easily form a self-supporting structure, and the influence of mining action on the roof sinking would be reduced. By comparing the three parameters, we that the influence degree on step sinking was overburden load q > immediate roof thickness H > support stiffness k1.

Engineering Example
On the basis of the above research, we analyzed the catastrophe instability condition of no. 12401 working face in Shendong mining area, calculating its reasonable support resistance to verify the theoretical analysis results. According to Equation (19), the necessary and sufficient conditions for system instability can be determined when hj = 17.
where ∑h is the thickness of an immediate roof rock stratum, m; M is the thickness of coal seam, 6.2

Engineering Example
On the basis of the above research, we analyzed the catastrophe instability condition of no. 12401 working face in Shendong mining area, calculating its reasonable support resistance to verify the theoretical analysis results. According to Equation (19), the necessary and sufficient conditions for system instability can be determined when hj = 17.  With the increase of external load q, the total displacement a of the roof increased linearly, indicating that the damage degree of system instability to the working face increased with the increase of the external load, as shown in Figure 13. It can be seen from Figures 14 and 15 that the thickness H of the immediate roof and the support stiffness k 1 were inversely proportional to the roof step sinking. With the increase of the thickness H of the immediate roof and the support stiffness k 1 , the step sinking of the roof decreased, which was conducive to the stability of the working face. When the support stiffness and the immediate roof thickness changed within a small range, the increase of the support stiffness had a greater effect on ∆h. When the support stiffness was close to or exceeded the immediate roof stiffness, the effect on the roof step sinking became smaller until it basically remained unchanged. Therefore, increasing the support stiffness (i.e., support resistance) in a certain range played a positive role in reducing the roof step sinking and ensuring the stability of the working face.
When the thickness of the immediate roof increased to a certain value, the step sinking value remained basically stable, indicating that the larger the thickness of the rock layer, the more stable the immediate roof, as the overlying rock layer could easily form a self-supporting structure, and the influence of mining action on the roof sinking would be reduced. By comparing the three parameters, we that the influence degree on step sinking was overburden load q > immediate roof thickness H > support stiffness k 1 .

Engineering Example
On the basis of the above research, we analyzed the catastrophe instability condition of no. 12401 working face in Shendong mining area, calculating its reasonable support resistance to verify the theoretical analysis results. According to Equation (19), the necessary and sufficient conditions for system instability can be determined when h j = 17.
where h is the thickness of an immediate roof rock stratum, m; M is the thickness of coal seam, 6.2 m; P is the recovery rate of the working face, 80%; and k p is the bulking coefficient of the immediate roof rock stratum, 1.15-1.30, 1.25 in this study. According to Equation (29), the thickness of the immediate roof was 31.3 m. Combined with the mining conditions of no. 12401 working face, the thickness of the immediate roof was 31.3 m when the collapsed gangue filled the goaf. Therefore, in the process of collapse, the mudstone with a thickness of 5.86 m and the fine sandstone with a thickness of 4.72 m above the coal seam will be transformed into an immediate roof and enter the goaf, forming a combined cantilever beam structure. At the time of the working face first weighting, the step sinking of the main roof was between 210 mm and 235 mm, with an average of 224 mm. The theoretical calculation results were close to the monitoring data of the step sinking, which proved the rationality of the theoretical analysis. However, there was also a difference of 0.021 m, indicating that the immediate roof had undergone a catastrophic instability, resulting in the step sinking of the main roof, which was consistent with the measured results.
(2) Judgment of system stability Necessary conditions for catastrophic instability: From λ= EA/L, λ 1 = λe −2 = E Z A H Z e −2 , and substituted into Equation (22)  It can be found from the calculation results of the catastrophe instability criterion that the mechanical system composed of the no. 12401 working face can meet the necessary and sufficient conditions of instability. Under the current support conditions, the catastrophic instability of the system is prone to occur and induce roof accidents. It can be found from the in-site monitoring that the roof had different degrees of step sinking with the average subsidence of 224 mm during the first weighting. The rib fall of coal wall is serious, which proves that the theoretical calculation is consistent with the actual situation. Therefore, it is necessary to adjust the external conditions, such as increasing the support resistance or improving the advance speed, in order to prevent the roof cutting along the coal wall.

(3) Reasonable calculation of support resistance
The support resistance of system instability can be determined by Equation (28) R Z > EA The calculated support resistance required to ensure the stability of the system was greater than the rated support resistance of 18,000 kN, and thus it is necessary to increase the support resistance to prevent roof cutting accidents and ensure the safety of the working face.

Conclusions and Summary
In this paper, the structural characteristics and instability process of the first weighting in large mining height working face were analyzed through a physical experiment and a numerical simulation, performed in order to study roof sliding and instability along the coal wall during the first weighting in the working face with large mining height in shallow coal seam. An analysis model based on catastrophe theory was developed to study roof sliding and instability along the coal wall. The necessary and sufficient conditions of system instability were obtained, and the main factors affecting the system stability were analyzed. The influence of each parameter on the stiffness ratio K was discussed and compared. The calculation equations of support resistance and roof step sinking were obtained, and the mechanical mechanism and main influence factors of the main roof step sinking were analyzed. Combined with the monitoring example, the rationality of judging roof instability by catastrophe theory was verified, and the rationality of the support resistance selected in the working face was checked. The following conclusions can be drawn through analysis of the theoretical results.
(1) The formation and weighting characteristics of composite cantilever beam structure of the first weighting in large mining height working face were analyzed through a physical test and a numerical simulation. The premise of working face stability was to ensure that the direct roof was not destroyed. (2) The sufficient and necessary conditions for system instability were obtained by analyzing the system model consisting of main roof, immediate roof, and support through catastrophe theory.
The system instability was related not only to the stiffness ratio K and material parameters of the support and immediate roof, but also to the load Q and the first weighting interval of the main roof. (3) The influence degree of each parameter on the stiffness ratio K was as follows: elastic modulus E > support stiffness k 1 > sectional area A > immediate roof thickness H. Increasing the support stiffness within a certain range had a positive effect on reducing the roof step sinking. (4) The main roof step sinking ∆h increased linearly with the external load q, and was inversely proportional to the thickness of immediate roof h and support stiffness k 1 . By comparing the three parameters, we found that the influence degree on step sinking was overburden load q > immediate roof thickness H > support stiffness k 1 . (5) The stability of no. 12401 working face was calculated and analyzed by the catastrophe theory, with the calculation results showing that the working face met the instability condition of the system, which was consistent with the monitoring results, and the support resistance should be greater than 19,232 kN to ensure the stability of the main roof.