Abstract

This paper is aimed at the inverted arch support instability of track roadway with mining level +1100 m in Liuyuanzi Coal Mine. By means of field investigation, theoretical analysis, numerical calculation, and engineering practice, the instability reasons of inverted arch structure are expounded, the mechanical mechanism of instability of inverted arch structure is revealed, and the “sliding-rotating beam” for the instability of inverted arch structure is put forward. Based on Fenner’s formula and mechanics principle, the equilibrium equation of “sliding-rotating beam” is given. The results show the following: Firstly, the insufficient stiffness at the joint of the inverted arch structure and the U-shaped steel support on the floor is the key reason for the floor instability. Secondly, when the action stress of the “sliding-rotating beam” is less than the critical value, three kinds of instability modes of the inverted arch structure may occur, that is, sliding upward, rotating upward, or sliding-rotating upward. Each instability criterion and critical value are also different. Considering the axisymmetric condition, the critical value calculation formula of the three modes can be simplified into one formula. Thirdly, the equivalent friction factor restricts the stability of the “sliding-rotating beam,” and there is a “breaking point.” The relationship between the equivalent friction factor and the action stress of the “sliding-rotating beam” is “class hyperbola.” When the equivalent friction factor is greater than the “breaking point value,” the “sliding-rotating beam” may remain stable. Moreover, with the increase of equivalent friction factor, the action stress required for the stability of the “sliding-rotating beam” is smaller, and it tends to be more stable. The breaking point value of equivalent friction factor is 18.6. Finally, 36U-shaped steel round frame with bolt-mesh-shotcrete-combined support is applied to improve the equivalent friction factor and the foot stiffness of U-shaped support in roadway. After two months of on-site implementation, the floor heave was reduced by 69.1%. In conclusion, the theoretical analysis is correct and the control method is effective.

1. Introduction

Floor heave is one of the focuses concerned by scholars at home and abroad [17]. A series of dynamic phenomena [811] such as roadway instability and floor water inrush are caused by floor heave, which not only affect mine ventilation, transportation, and pedestrians but also increase a lot of maintenance workload and mining costs [12, 13].

Due to the complexity of roadway occurrence environment, scholars have different opinions on the mechanism of roadway floor heave [1417]. According to different floor heave mechanism, the control technology is also different [11, 1821]. Mo [22] introduced some of the main floor heave events in the development of the Glencore Bulga Underground plant. The previously mentioned study indicates that the high horizontal stress of the roadway surrounding rocks and certain types of floor lithology configuration are the reasons for the failures of floor strata. Malkowski and Ostrowski [23] carried out the research of numerical modeling based on convergence monitoring records. The convergence measurements were carried out on three stations in a selected main gate in a coal mine from Upper Silesia Coal Basin near Jastrzebie-Zdroj. The width, the height, and the heave of the floor of the gate were measured. The authors highlight that, during the measurements period, the horizontal stress was 1.45 to 1.61 times greater than the concurrent vertical stress. The so high horizontal stress causes heave of unsupported gallery floor, which is commonly observed in the mines in Silesia. Mo et al. [24] studied especially the buckling failure mechanism in the floor and relevant analysis techniques on the coal mine floor. This paper introduces a new floor classification system, the Coal Mine Floor Rating, based on the underlying failure mechanisms. The new floor mass rating system quantifies the stability of floor strata by considering two main factors: uniaxial compressive strength and discontinuity spacing of the floor units. The floor heave and nonfloor heave cases are determined depending on the Coal Mine Floor Rating and the Horizontal Stress Rating. Lai et al. [25] studied the mechanism of roadway floor heave by field geological investigation, theoretical analysis, and numerical simulation. They thought that the surrounding rock of the roadway presents asymmetrical shrinkage under the original support condition, and it is the extrusion flow type floor heave. The bottom without support and influence of mining are the important causes of floor heave. Jia et al. [26] analyzed the stress environment, deformation, and failure characteristics of mining roadway in the fully mechanized caving face. The results revealed that the intrinsic cause of the large asymmetrical floor deformation in the mining roadway is the asymmetrical phenomenon of the surrounding rock’s stress environment caused by mining. This also results in the nonuniform distribution of the mining roadway floor’s plastic zone. Zhang et al. [27] analyzed the deformation and failure characteristics of the mining roadway’s floor and elucidated the floor heave mechanisms of shallow bulging and deep sinking. Moreover, the displacement law of the floor rock at different depths was further analyzed, and the characteristics of “two dots and three areas” of the floor heave were determined. Kang et al. [28] analyzed the influence of mining and side pressure coefficient on the roadway floor heave. The relationships between the thickness of the soft rock stratum and the mechanical properties of the surrounding rock and floor heave were analyzed. Qin et al. [29] proposed and applied a reinforcement scheme for deep dynamic soft rock roadway. The results indicate the following: by increasing the supporting strength of the internal bearing structure, cohesion, and internal friction angle of the surrounding rock and by reducing the influence of mining, making the external bearing structure close to the roadway, and reducing the thickness of the bearing structure, the bearing capacity of the shallow surrounding rock in the roadway can be improved. Peng et al. [30] designed a U-type steel inverted arch and bolts support scheme to control the floor heave of the roadway. After calculation through the mechanical model of the inverted arch, its bearing capacity is 0.56 MPa. Through comparison and analysis of various supporting schemes, the “U-shaped steel + inverted arch + pouring concrete + backwall grouting” technology was selected, and the engineering practice showed that the supporting scheme can effectively improve surrounding rock stability. Shreedharan and Kulatilake [31] presented stability studies on two tunnels, a horseshoe-shaped tunnel and an inverted arch-shaped tunnel, in a deep coal mine in China, performed using the 3DEC distinct element code. They illustrated the importance of incorporating stress relaxation before support installation through the increased support factor of safety and reduced grout failures. Finally, a comparison between the two differently shaped tunnels established that the inverted arch-shaped tunnel may be more efficient in reducing roof sag and floor heave for the existing geo-mining conditions. Chun et al. [32] designed a yielding inverted arch of composite structures to control serious floor heave of tunnel. They studied constitution of the inverted arch to present mechanism of the composite structure to restrain floor heave. By field trial of the yielding inverted arch with composite structures, it is an effective method to control floor heave, and its important function is to absorb main elastic deformation coming from bottom strata.

The above-mentioned scholars have done a lot of successful research and application on the mechanism and control of floor heave. Because of the complex and changeable environment of coal mine roadway or the fact that the key points of floor heave control are different, the researches on the mechanism of floor instability under inverted arch support are relatively few or only a part of the influencing factors are studied. Therefore, further study was needed to solve the instability and control of the antiarch support of the roadway with mining level +1100 m in Liuyuanzi Coal Mine, especially the failure mechanism of the roadway floor under the antiarch support. According to the failure characteristics of the roadway floor, the U-shaped steel inverted arch support was regarded as a beam that can slide (or rotate) at both ends. The stability conditions and influencing factors of the beam were analyzed by using Fenner’s formula and mechanical principle. The causes of floor heave under the support condition were given, and the feasible and effective control scheme was put forward. The method and technology can provide theoretical basis and control measures for the normal use and safe operation of roadway.

2. Project Overview

The coal-bearing strata are Yan’an Formation of Jurassic Mesozoic in Liuyuanzi Coal Mine. The coal seams No.4-1 and No.5-1 are divided into upper coal group, and the coal seams No.7-1 and No.8-3 are divided into lower coal group. According to the borehole of No. L503, the track roadway, conveyor roadway, and return air roadway are arranged horizontally at +1100 m (i.e., at the lowest point of coal seam No.5-1). The track roadway is buried at a depth of 455 m and lied in mudstone, fine sandstone, medium sandstone, and sandy mudstone.

The average RQD index of the rock mass is from 39% to 87% around the track roadway. The experimental results indicate that the compressive strength of the mudstone in 5-1 coal seam floor and 7-1 coal seam roof is from 3.94 MPa to 22.2 MPa in natural state, while it is from 0.16 MPa to 10.5 MPa in saturated state, and the softening coefficient is from 0.04 to 0.47. The characteristics of the rocks are easy to weather, soften, and swell. The original support design of the track roadway in Liuyuanzi Coal Mine is as follows: the bolt-mesh support is used for the first support, the shotcreting and 36U-shaped steel are used for the second support, shotcreting thickness is 200 mm, the net section is 15.7 square meter, and the driving width and height, respectively, are 4900 mm and 4500 mm. All these are shown in Figure 1.

3. Characteristics and Influence Factors of Floor Heave in Roadway

3.1. Characteristics of Floor Heave in Roadway

According to the on-site investigation, the main characteristics of the floor heave in roadway are as follows:(1)Large deformation: According to the field investigation, the floor heave of the roadway with mining level +1100 m is generally 0.5 m to 2.2 m. The floor heave can still reach 0.5 m after the inverted arch support is applied to the local position.(2)Floor crack: The floor of the roadway is strongly deformed, and floor crack is particularly serious when encountering water.(3)High deformation speed: According to the on-site measurement, the floor heave speed of the track roadway can reach 10 mm to 20 mm/d, which causes floor dinting every month; otherwise it will affect transportation.

3.2. Factors Affecting Deformation
3.2.1. Improper Support

Although the antiarch support is applied to track roadway, because of the insufficient rigidity of the contact between the antiarch structure and U-shaped steel support, which leads to the failure of the supporting structure of the floor, floor heave occurs under the action of the floor pressure.

3.2.2. Strata Structure

The track roadway is located near the syncline axis, and the axis of the track roadway is not perpendicular to the syncline axis, so that the track roadway is subjected to a certain level of tectonic stress, forcing the floor rock to extrude into the tunnel, thus causing the extrusion of the antiarch structure and intensifying the floor heave.

3.2.3. Properties and Hydraulic Characteristics of Surrounding Rock

The coal-bearing strata are Yan’an Formation of Jurassic Mesozoic in Liuyuanzi Coal Mine, with well-developed bedding and joints, extremely low strength, easy weathering and fragmentation, easy softening in water, and certain expansibility. Under the influence of Cretaceous aquifer water, the roof water of No. 4 coal, the roof water of No. 5 coal, floor water, and structural water, the water inflow can reach 396.54 m3/d. The floor mudstone of the track roadway will become muddy and swell when it meets the water, which intensifies the plastic deformation of the floor. This is an important reason for the drastic and long duration of floor heave.

4. Mechanical Environment and Instability Analysis of the Floor

4.1. The Stress Field in Roadway

For convenience of the analysis and consideration of engineering practice, the roadway is simplified as the circular hole problem of axisymmetric plane strain in elastic-plastic mechanics in which both the load and the structure are axisymmetric. The mechanical model is shown in Figure 2. is original rock stress; is roadway radius; is plastic zone radius; is support strength; this paper refers to the stress acting on surrounding rock of roadway or other supporting structures, which is called the action stress; and are polar coordinates of any point in plastic zone.

In polar coordinates, according to Fenner’s formula and Mohr-Coulomb criterion [3335], the stress calculation formula of plastic zone of surrounding rock of circular tunnel can be obtained:where is radial stress at any point, MPa; is tangential stress at any point, MPa; is shear stress at any point, MPa; is the same as the previous one, MPa; C is rock cohesion, MPa; φ is angle of internal friction for rock, °; is radius of roadway, m; the rest are as before.

Formula (1) in polar coordinates is transformed into a stress equation in a rectangular coordinate system by means of the following equations:

That is,

Substituting (1) into (3), the stress pattern in the Cartesian coordinates near the roadway adjacent to the plastic zone is

In general, the shear stress can be neglected when considering the rock mass fracture at the bottom corner of the roadway after excavation. Formula (4) is simplified to

4.2. The Instability of Inverted Arch Structure

The key point to control the floor heave is the joint of U-shaped steel support and inverted arch structure. The friction force is used to ensure the stiffness of the structure. Therefore, the mechanical properties of the joint determine the performance of the whole supporting structure. According to the mechanical behavior of the joint, the floor heave is divided into three types: “sliding upward”, “rotating upward,” and “sliding-rotating upward,” as shown in Figure 3.

(1)“Sliding upward” type: As shown in Figure 3(a), when the floor stress is symmetrically distributed and the friction force between the inverted arch and the U-shaped steel support and its load is insufficient to balance the floor pressure, the inverted arch structure will rise, resulting in the overall uplift of the floor. This type was called “sliding upward” floor heave. Generally, this kind of floor is risen upward on the whole, which is not affected by asymmetric external forces and is not easy to cause floor crack.(2)“Rotating upward” type: As shown in Figure 3(b), when the stress on the floor is asymmetrically distributed, the floor relatively keeps balance at the narrower end of the concentrated stress zone (called the hinge end), while it will rotate upward around the hinge end at the other end of the concentrated stress zone (called the active end). This type was called “rotating upward” floor heave. Because the floor is subjected to asymmetric external forces after rotation, the bending resistance of concrete is smaller than U-shaped steel, which is easy to cause floor crack.(3)“Sliding-rotating upward” type: On the one hand, as shown in Figure 3(a), when the inverted arch structure is sliding upward, the stress on the floor may be distributed asymmetrically for some reason, resulting in the change of the floor from sliding upward to rotating upward; on the other hand, as shown in Figure 3(c), when the inverted arch structure is rotating upward, the contact area between the hinged end and the U-shaped steel support decreases, resulting in a great reduction in the frictional force at the joint. The hinge ends cause sliding upward due to mechanical unbalance. All of the above were collectively called “sliding-rotating upward” floor heave. As mentioned earlier, this situation is liable to cause floor crack.

5. “Sliding-Rotating Beam” Model and Its Mechanical Behavior

5.1. Model of Sliding-Rotating Beam

As shown in Figure 1, the U-shaped steel support strictly limits the transverse displacement of the inverted arch structure, but its longitudinal constraint is relatively weak. At this point, the inverted arch structure can be assumed as a “board.” To simplify the analysis, the following assumptions are made for the “board” structure:(1)The inverted arch and the overlying concrete are a continuous integral structure along the axial direction of the roadway. In order to simplify the calculation, the curved beam is regarded as a straight beam.(2)The two ends of the “board” are movable supports, which can move up and down but are controlled by friction factors at the contact.(3)Because of the high bending strength of U-shaped steel antiarch support, the bending of “board” is not considered, and the stress distribution of the floor is consistent with that of Figure 3.

Therefore, the inverted arch structure of the floor can be regarded as a “beam” supported transversely by U-shaped steel at both ends, which is called a “sliding-rotating beam,” which can either slide vertically or rotate around an end, as shown in Figure 4(a). The force of “sliding-rotating beam” is shown in Figure 4(b), where L is the length of “sliding-rotating beam” and L = a + b + c, m; a and c are the width of concentrated stress action zone, m; b is the width of low stress action zone, m; d is the thickness of the beam, m; and are the concentrated stress of the floor rock on the “sliding-rotating beam” in the width zone equal to “a,” MPa; and are the concentrated stress of the floor rock on the “sliding-rotating beam” in the width zone equal to “c,” MPa; is the stress of the floor rock on the “sliding-rotating beam” in the width zone equal to “b,” MPa; is the friction stress of U-shaped steel support and antiarch support, MPa; q is the self-weight stress of concrete overlying inverted arch, MPa.

5.2. Balance and Motion of “Sliding-Rotating Beam”
5.2.1. The Equilibrium Condition of “Sliding-Rotating Beam”

As shown in Figure 4(b), the “sliding-rotating beam” keeps static balance in the “x” direction, taking the length in units at the direction perpendicular to the “xoy” plane. The static balance of the “y” direction and the moment balance (using “o” as the center of moment) of the “sliding-rotating beam” can be expressed aswhere ; ; ; f is the equivalent friction factor; γ′ is the body force of concrete, N/m3; h is the thickness of the concrete, m; the others are the same as before.

In this paper, the equivalent friction factor is a mechanical parameter to measure the difficulty of reverse arch supporting structure moving to roadway space. The equivalent friction factor here mainly considers the interaction of contact mode (connection mode) between inverted arch structure and roadway wall (or supporting structure). It can be calculated by the ratio of the resistance between the surfaces of a structure or component to the compressive stress applied to the surface. Under the general surface contact and friction conditions, the equivalent friction factor is equal to the static friction coefficient, and its value is less than 1. The equivalent friction factor will be greater than 1 or even greater if the resistance increasing accessories are set between the reverse arch supporting structure and the U-shaped steel support foot of the roadway. Therefore, the value range of equivalent friction factor in this paper is 0 to ∞.

The stress on the “sliding-rotating beam” around the roadway can be calculated according to formula (4); the radius of the roadway in the a, b, and c zones of the “sliding-rotating beam” is , , and , respectively.

As shown in Figures 2 and 4(b), replace in (5) with ; then substituting (5) into (6) and (7), respectively, they can be written aswhere and are the action stress of the floor when the “sliding-rotating beam” is in critical state of sliding and rotating, which is called critical value, MPa; ; .

In fact, is positive. If the result calculated by formula (8) or formula (9) is negative, it means that the roadway floor is in instability state under the current conditions. In conclusion, the condition of “sliding-rotating beam” motion is as follows:(1)If , the floor will maintain balance; otherwise, the “sliding upward” type floor heave will be produced.(2)If , the floor will maintain balance; otherwise, the “rotating upward” type floor heave will be produced.(3)If and , the “sliding-rotating upward” type floor heave will be produced.

5.2.2. The Influence of f on Ps

According to the analysis in Section 4.2, the main reason for the failure of reverse arch support in this paper is the contact relationship and connection mode at the foot of support. Therefore, this section mainly discusses the effect of equivalent friction factor on .

As shown in Figure 4(b), considering the symmetry distribution of the concentrated stress on both sides of the roadway floor, the thickness of the “sliding-rotating beam” is very small compared with the roadway radius. Let a = c, and ; formulas (8) and (9) can be written aswhere is the action stress of the “sliding-rotating beam” under the condition of symmetrical stress distribution on the floor, also known as the critical value, MPa.

We consider the value of equivalent friction factor when the roadway floor is in the limit of self-stability, that is, when the roadway floor does not need to be supported. Let in formula (10); we can get

The calculation parameters of track roadway with mining level +1100 m in Liuyuanzi Coal Mine are shown in Table 1. By substituting the data in Table 1 into (11), we can draw the relationship between the equivalent friction factor and , as shown in Figure 5. If the equivalent friction factor can be arbitrarily increased, the equivalent friction coefficient will have a “breaking point value,” which divides the diagram into “meaningless area” and “actual engineering area.” When the equivalent friction factor is less than the “breaking point value,” that is, is negative, it is indicated that the reverse arch support is unstable. Only when the equivalent friction factor is greater than the “breaking point value,” is positive, it means that the reverse arch support of roadway is likely to be stable. As can be seen from the “actual engineering area” in Figure 5, the relationship between the equivalent friction factor and the action stress of the “sliding-rotating beam” is “class hyperbola.” We mainly analyze the “actual engineering area.”

It can be seen from Figure 5 that the support of roadway floor is meaningful only when the equivalent friction factor of the roadway support structure exceeds the “breaking point value” of 18.6 under other unchanged conditions. With the increase of the equivalent friction factor, the action stress decreases gradually. This shows that increasing the equivalent friction factor is helpful to the stability of the reverse arch supporting structure. However, the friction factor of reverse arch support is only 0.15. Obviously, the floor heave of roadway must occur under this condition.

According to formula (11), the equivalent friction factor of roadway without support is 56.3. This value is far beyond the range of static friction coefficient. This confirmed the inevitability of floor heave once again. In order to increase the equivalent friction factor of the footings, other resistance-increasing accessories or connecting components must be used in the reverse arch support of the mine roadway. Of course, if considering the influence of lateral force on the reverse arch support, the resistance-increasing accessories and connecting members are still in an unstable state. Therefore, in this case, the U-shaped steel support with full closure is recommended.

6. Design and Implementation of Floor Heave Control in Roadway

6.1. Design of Floor Heave Control in Roadway

Aiming at the original support scheme, bolt-grouting support and inverted arch support were used to repair the roadway for many times, while the floor heave had reached 0.5 m in three months. At the same time, the maintenance cost increases gradually, which has seriously affected the normal use of track roadway. From the previous theoretical analysis, it can be seen that the floor heave is mainly caused by the insufficient stiffness between the inverted arch and the foot of U-shaped steel support. In order to enhance the support stiffness at the bottom corner of the roadway, 36U-shaped steel round frame with bolt-mesh-shotcrete-combined support is proposed to control the floor heave. The support design of roadway is shown in Figure 6.

The construction sequence of the support scheme is as follows: firstly, bolt-mesh support is adopted; then 36U-shaped steel round frame is used to support roadway, and finally concrete is poured behind the 36U-shaped steel round frame. The bolt is left-hand thread bolt with a diameter of 20 mm and a length of 2300 mm. The row spacing between bolts is 900 mm × 900 mm. The Z2335 resin capsule is used as anchoring agent. The metal mesh is welded with concrete irons with diameter of 6 mm, and the mesh size is 80 mm × 80 mm. The diameter of the round frame is 4.6 m, which is composed of six 36U-shaped steel round frames and the overlapping length is 500 mm. The spacing of the round frame is 600 mm. 200 mm C20 concrete is poured behind the round frame, and 800 mm thick plain concrete is laid on the bottom of the round frame.

6.2. Feasibility Analysis of Support Design
6.2.1. Numerical Model

In order to accurately analyze the feasibility and effectiveness of the new support design, the three-dimensional models of semicircular arch section (original support design) and circular section (new support design) are established and analyzed by FLAC3D, as shown in Figure 7.

6.2.2. Computing Mode and Parameters

To accurately analyze the problem, combining with the surrounding rock environment of the roadway to restore the actual roadway engineering, Mohr-Coulomb strength criterion, Isotropic-Fluid-Flow constitutive model, and Fluid-Mechanical Interaction calculation model are used in the models. The physical and mechanical parameters used in the calculation are shown in Table 2. The bolts and the U-shaped steel supports are generated by cable units and shell units, respectively. The specific mechanical parameters are shown in Tables 3 and 4.

6.3. Simulation Results

As shown in Figure 8, when excavating 6 m, the floor stress of the semicircular arch roadway shifts to the bottom corner, the compressive stress is 8.7 MPa, and the center zone of the floor bears tensile stress for about 0.6 MPa, and the former is 14.5 times of the latter. The floor stress of the circular roadway shifts to the two sides of the rock mass, the compressive stress is 7.0 MPa, and the center zone of the floor bears compressive stress for about 2.0 MPa, and the former is 3.5 times of the latter. So the circular section roadway support scheme may change the floor from tension to compression. Moreover, to a certain extent, it can homogenize the stress distribution of the floor, which is conducive to the stability of the floor.

As shown in Figure 9, the plastic depths of floor of semicircular arch tunnel and circular tunnel are 2.5 m and 1.0 m, respectively, when excavating 6 m, while the latter is only 40% of the former, which shows that the circular section roadway support scheme can effectively restrain the plastic expansion of the floor.

As shown in Figure 10, the floor displacements of semicircular arch roadway and circular roadway are 12.3 mm and 7.8 mm, respectively, and the latter is only 63% of the former when excavating 6 m. Thus, circular roadway support scheme can effectively control floor heave.

By the above analysis, 36U-shaped steel round frame with bolt-mesh-shotcrete-combined support is feasible and reliable. After the implementation of the scheme, the roadway waistline is demarcated again, and the roadway floor heave is calculated based on it. After two months of monitoring, the floor heave of roadway is reduced from 246 mm to 76 mm, and the floor heave is reduced by 69.1%.

7. Conclusion

(1)Insufficient strength of the connections between U-shaped steel and inverted arch structure is the key to instability, and the roadway is located in the lower part of synclinal structure coal seam, the floor rock mass is broken, and the floor heave is aggravated by water softening.(2)The mechanical model of “sliding-rotating beam” is proposed for the instability of the inverted arch floor structure, and the equilibrium equation of the “sliding-rotating beam” is deduced. If , “sliding-rotating beam” will keep stability. If , three instability modes will occur in the “sliding-rotating beam”: “sliding upward,” “rotating upward,” and “sliding-rotating upward.” Under normal conditions, the critical value calculation formulas of the three modes can be simplified into one formula.(3)The equivalent friction factor affects the stability of reverse arch supporting structure. By increasing the equivalent friction factor between the inverted arch structure and the U-shaped steel support, the “sliding-rotating beam” can keep balance, but the equivalent friction factor has a “breaking point value”; when the equivalent friction factor is less than the “breaking point value,” the “sliding-rotating beam” cannot keep balance, and the floor heave must occur. If the equivalent friction factor is greater than the “breaking point value,” then a sufficiently large equivalent friction factor can keep the “sliding-rotating beam” stable. The “breaking point value” of the foot of reverse arch support is 18.6, which is much larger than the one of the reverse arch supporting structure. So the floor heave must happen.(4)The 36U-shaped steel round frame with bolt-mesh-shotcrete-combined support is applied to improve the coupling relationship of bottom corner support. After the implementation of the scheme, the results show that the floor heave is reduced by 69.1%. Therefore, the fully closed yieldable metal support is more suitable for roadway floor heave control.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant no. 51964037) and the Key Projects of the National Natural Science Foundation of China (Grant no. 51634007).