Quantitative Assessment of Bed-Separation Dynamic Development Caused by Inclined Coal Seam Longwall Mining

: Coal mining under the Quaternary thick loose layer a ﬀ ects key strata breakage, Bed-sep-arations development, ground subsidence, and other studies. This paper presents a method for solving the de ﬂ ection of a large-de ﬂ ection inclined thin plate under a thick loose-layer cover with additional lateral loads and midplane forces. The methods presented are based on the principle of large-de ﬂ ection of thin-plate, energy method, and fracture mechanics theory. The 7225 work face in Anhui Province, China, was studied. Combined with the large-de ﬂ ection inclined thin plate model, the initial breakage distance within the main roof plate was calculated to be 33 m with the initial break-age angle of 61.2°, and the period breakage distance was calculated to be 21 m with the period breakage angle of 55.4°. The distribution range of “Vertical Three Zones” from 7225 working face to the ground, including the height of the caved zone is 38.07 m, the height of the fractured zone is 41.13 m, and the height of the curved zone with the thick loose layer removed is 187.56 m. During the dynamic development of the principal key strata (PKS), the de ﬂ ection value develops from 0 mm to 2714 mm with 7225 working face mining, and the maximum value of the spatial volume is 56,485 m 3 , which is veri ﬁ ed by Three-dimensional Discrete Element Code (3DEC) numerical simulation. The dynamic development of Bed-separation within the overlying strata, with a maximum development height of 545.2 mm and a maximum volume of 11,228.1 m 3 of the Bed-separation cavity. The dynamic development of the Bed-separation height and the cavity under di ﬀ erent mining length and width conditions of the working face are also discussed. The large-de ﬂ ection inclined thin plate model proposed in this paper e ﬀ ectively explores the dynamic de ﬂ ection and fragmentation law of the overlying strata induced by the inclined working face of Longwall mining and provides a theoretical basis and computational model for quantitatively evaluating the dynamic development of the Bed-separation cavity.


Introduction
China's Anhui, Shanxi, and Shandong Provinces extensively distribute quaternary thick loose mining layers [1].Within the mining areas of Huaibei in Anhui Province, there is a specific location where mining operations occur at a depth of 600 m.The 200 m thin rock strata encumber the working face, followed by a 400 m thick loose layer.This results in general movement damage to the rock strata.This necessitates developing and implementing a highly accurate model that adheres to the inherent principles of the overburden rock movement for practical use in production.The thick, loose layers mainly comprise sparsely scattered sediments that have yet to undergo the process of solidification and hardening into rock [2].These sediments exhibit similar physical features.The underlying rock layers, located beneath the broad, loose strata, consist of progressively deposited mudstone, fine sandstone, siltstone, and other stones with various physical characteristics [3].Qian Minggao proposed the Key Strata hypothesis to explain how thick rock layers affect the mobility of different rock layers, especially those with different levels of softness and hardness beneath thick, loose layers.Subsidence deformation in the key strata leads to simultaneous deformation of the overloaded rock layers, entirely or partially, along with the non-key strata.The underlying non-key layer experiences asynchronous deformation, causing it to separate from the bed beneath the key stratum [4,5].Research on bed separation plays a crucial role in gas extraction above coal seams [6], examination of water risks related to bed separation [7], and the implementation of separation grout Injection mining.
The deterioration of rock layers due to excessive pressure and the development of gaps between beds caused by underground mining is a complex process driven by various spatial and temporal factors [8,9].Prior research [10,11] conducted experiments with comparable materials and augmented their results with numerical simulations [12][13][14].Nevertheless, it is crucial to acknowledge that these studies have a comparatively less robust theoretical basis for their approaches and mostly concentrate on visualizing the process of structural deterioration in overloaded rock layers.Theoretical analyses of elastic foundation support beams [15,16] and plates [17][18][19] provide the basis for comprehending the structural deterioration process of heavily loaded rock layers.The purpose of the elastic foundation support beam and plate concept is to specifically analyze the primary stratum and investigate the structural failure of the principal strata and the occurrence of the Bedseparation phenomenon beneath it.Within the framework of the Winkler foundation, it is hypothesized that thick, unconsolidated layers function as the weight on top.In contrast, other rock layers are considered the weight on top of the main layer.The settlement deformation of the primary layers shows a slight variation in thickness, and it is estimated by utilizing a thin plate with negligible bending.It is hypothesized that the working face is positioned horizontally to the rock layer above it, and the load is delivered at a right angle to the center of the plate.The Winkler foundation is based on the idea that it comprises several independent and non-interacting springs.The force exerted by the foundation depends only on the amount it is deflected.
When a coal seam is sloped, it is necessary to use Longwall mining techniques that have prolonged lengths along the direction of the seam's slope and strike.As a result, there is a significant difference in the vertical distance between the boundaries of the coal seam opening and closing and the mining depths along the uphill and downhill boundaries of the coal seam [20,21].The assumption that the overburden weight is equal or perpendicular to the midface of the key strata is not applicable in the context of inclined working face mining.This mining technique causes the overlying material to apply forces on the central plane of the targeted rock layers [22].The mining process entails the progression of slanted excavations, resulting in changes in the mining depth.As a result, this causes changes in the load distribution of each essential stratum on the top of the bedrock.If the seam is mined to a greater thickness, the deformation of the critical strata towards the edge of the slab can be approximated to the deformation of a thin plate [23,24].However, it must be recognized that the magnitude of deformation is roughly equal to the plate thickness.Therefore, if we use a slightly deformed thin plate model and ignore the strain in the middle plane, the prediction accuracy of the critical strata deformation will be reduced.Although it is possible to obtain an analytical solution for the deflection of a thin plate by using the Navier solution for a simply supported thin plate with small deflection, it is difficult to obtain the Navier solution for a large-deflection inclined thin plate with an additional mid-plane force, and the accuracy of the results is lower than that of the energy method [25][26][27].The main concept of the energy method is to describe the process of deflection of a thin plate in which the strain energy of the plate is dynamically balanced with the work performed by the external loads [28].The strain energy method has the advantages of a simple method and easy programming for calculating the deflection coefficients of large-deflection inclined thin plates, while the calculation results can effectively meet the engineering requirements [29,30].
This study thoroughly examines the angle of inclination of the coal seam.It is assumed that the rock layers above the working face are parallel to the working face and have the same angle as the coal seam.A thin plate model with large deflection is generated using the Key Strata Theory.The strata located above the load maintain a vertical orientation.The item provided can be divided into two loads: one exerting force at a right angle to the inclined plate and another exerting force parallel to the mid-plane of the inclined plate.Theoretical improvement of the Winkler foundation discontinuity problem entails adding constraints to individual springs in the Winkler model, creating a two-parameter foundation model.This work suggests utilizing the inclined large deflection thin plate model to improve the accuracy of predicting the deflection curve of significant strata.Furthermore, its objective is to ascertain the location, form, and beginning and ending times of bed separation formation under dynamic stress resulting from inclined coal seam mining.

Calculation of Overlying Load Based on Stress Transfer
As the rock strata subside, consolidation, extrusion, and shear produce opposite stresses.The direct use of the calculation method for the product of capacity and thickness will make the load calculation value large.The stress transfer method can calculate the load above the bedrock for soft and fragile rock strata or thick, loose strata soils.Since the rock strata have a certain cohesion, their strength obeys the Mohr-Coulomb theorem [31][32][33], i.e., =c+ tan where  is the soil shear strength, c is the soil cohesion,  is the soil tensile strength, and  is the angle of friction within the soil.
The bedrock internal cohesion is considerable, the surrounding rock on both sides is stable, and the subsidence area is limited to the upper part of the top plate.First, a schematic diagram of the calculation of the loads applied to the roof by the overlying bedrock is established.
As shown in Figure 1, the shear forces on both sides of a differentially thin bedrock of depth z and thickness dz are as follows: where z is the depth of a differential thin strata in bedrock, h  is the differential thin strata horizontal stress at depth z, v  is the differential thin strata vertical stress at depth z,  is the soil shear strength, c is the soil cohesion,  is the soil internal friction an- gle, and  is the lateral pressure coefficient.
Diagram of stress transfer from the overlying bedrock to the roof plate with both sides stable surrounding rocks.
Assuming a uniform load p at the surface, the equilibrium equation for the differ- ential thin strata d z in the vertical direction is Calculations obtained tan ( ) Solving the differential Equation obtains where  is the rock weight, a is the span of the subsidence area, v  is the differen- tial thin vertical stress at depth z, and  is the lateral pressure coefficient.
According to the surface boundary conditions, when Now, the vertical stress v  at the bedrock depth z can be found as follows: When z H  , the vertical load transferred to the top plate during bedrock subsidence is obtained as

Calculation of Overlying Load Based on Stress Transfer
The Winkler model is an assumption made by the Czech engineer Winkler in 1867 when calculating railroad tracks that the pressure at any point on the surface of the foundation is proportional to the displacement at this point [34,35].
where k is called the foundation reaction coefficient.
According to this assumption, the displacement of a point on the foundation is independent of the other points.The foundation comprises many independent springs that do not affect each other, and the Winkler foundation model is illustrated in Figure 2. The model assumes that the foundation deformation occurs only within the footing area and that there is no deformation outside.The Winkler foundation model cannot diffuse strain and deformation and thus has severe defects in the theory.Although the elastic half-space model is perfect in theory, it has more difficulties in mathematical calculation.For this purpose, some foundation calculation models between the two are cited, among which the two-parameter foundation model uses two independent parameters to characterize the foundation properties, which theoretically improves the disadvantages of foundation discontinuity in the Winkler foundation model and is more straightforward than the elastic half-space model in mathematical calculations.
Since the plate is tilted in space, the Winkler foundation model considers no interaction between adjacent units irrelevant.A two-parameter model is chosen, as shown in Figure 3.The feature is to connect a series of springs in a Winkler foundation with a film subjected to a constant value of tensile force T. The balance between the film force and the spring system is considered.The relationship between the lead hammer load P and the displacement of the foundation in the three-dimensional case is For the inclined plate on the foundation, the film tension N can be approximated as a vector sum of the gravity of the inclined plate and the force of the overlying load along the inclined plate.

Establishment of Thin Plate Model
The coordinate axes x and y are taken in the plate's mid-plane, and the z-axis is perpendicular to the mid-plane.Two pairs parallel to the xz and yz planes are used to intercept the cell from the plate to establish the bending moment MX, MY, the torque MXY and the lateral vertical shear QX, QY.The results are shown in Figure 4. Considering the load distribution on the upper surface of the plate, the load intensity is denoted by q, and the load acting on the unit body is qdxdy.The equilibrium equation is obtained by projecting the force acting on the unit body in the z-direction.
The force acting on the unitary body is de-mobilized in the x-direction to obtain the equilibrium equation [36].
When ignoring the effect of shear QX,QY on the bending, as well as the compressive stress z  generated by the load q.At this point, the deflection is expressed by Equation [24]: , the strain energy increment per unit volume is For the Kirchhoff assumption for thin plates, the above Equation can be simplified as follows: The strain and stress components are replaced by the generalized deformation component , , x y xy    and the inner torque, respectively, and integrated along the thick- ness direction to obtain the deformation energy increment per unit area of the middle surface expressed by the generalized deformation increment: Under the condition of linear elasticity, the deformation energy density is a positive definite quadratic function of the generalized deformation: Dividing the above Equation along the median plane and expressing the generalized deformation component by the deflection w, the deformation energy for bending of thin plates is obtained as follows: When all the boundaries of the plate are rigidly fixed, or the deflection of the plate is zero on all sides, the strain energy of thin plate bending can be simplified as The strain energy resulting from the bending of the plate is mainly caused by the strain energy caused by the elongation of the mid-plane and the additional strain energy resulting from the bending of the plate by the lateral load.
Assuming that at each point on the mid-plane of the plate, the unit length force Thus, the strain energy caused by the elongation of the mid-plane is as follows: After the lateral load is applied, it bends the plate and causes additional strain on the mid-plane.When the plate is bent, the small displacements generated at a point in the mid-plane in the x, y, z direction components are represented by ux, uy, w.Now consider a mid-plane linear cell in the x-axis; from Figure 5, it can be seen that the elongation produced in the cell because of displacement ux is , compare the length of cell A1B1 and its projection length in the x-axis, the elongation produced in the cell due to displacement Additional elongation due to bending of mid-plane by lateral load Therefore, the total unit elongation generated in the x-axis of the unit taken in the mid-plane is Similarly, the total unit elongation produced in the y-axis is Due to the mid-plane shear strain generated by the plate bending, the shear strain generated by displacement ux,uy is  Summing the shear strains resulting from the various displacements obtains.
The additional strain energy caused by the mid-plane bending strain is By adding U1, U2, the strain energy generated by the bending of the plate due to the action of the mid-plane force is obtained.
Due to the thin plate theory, the displacements parallel to the mid-plane within the mid-plane of the thin plate are neglected in Kirchhoff's assumption, i.e., the displacements ux,uy along the x,y direction are neglected.
To the theory of large deflection of thin plates, the strain energy expression generated by the face force in bending of inclined plates satisfying Kirchhoff's assumption condition is

Solving the Key Strata Breakage Distance and Angle Based on the First Strength Theory
Due to coal seam mining, the main roof is prone to tensile damage when the main roof first shows some inclined cracks at the edge of the tensile area on both sides of the embedded end, causing the rock under the main roof to collapse.The angle of inclination of fracture inversions and fracture traces in the bedrock on the main roof plays a decisive role in the extent of damage to key strata in the overburden.The larger the fracture inversion step of each strata in bedrock, or the smaller the angle between the fracture trace and horizontal plane, the smaller the breakage damage of key strata in overlying rock, and the smaller the influence range of surface subsidence [37].Therefore, the analysis of the formation of bedrock fracture inversion steps and the inclination angle of each fracture trace is necessary to accurately predict the internal movement boundary of the overlying rock strata and the extent of surface subsidence [31], and the results are shown in Figure 7.According to the theory of maximum principal stress, it is known that rocks are usually resistant to compression but not to tension.The main influencing factor for the fracture of the rock formation is the maximum principal stress, and the rock formation will break in tension along the direction perpendicular to the principal stress.According to the distribution regulation of the primary tensile stress, the fracture trace of the threshold rock formation can be determined [38].
According to the relationship between the principal stress and the stress components, the principal tensile stress 1  and the principal compressive stress 3  at any point can be obtained: For the thin plate model, the theoretical geometric equations of the thin plate are derived.
Derivation of the theoretical physical equations for thin plates leads to Combining the physical equations of the thin plate with the expressions of principal tensile stress 1  and principal compressive stresses 3  , the distribution of principal tensile 1  and principal compressive stresses 3  in the thin plate during the initial and periodic fracture of the top plate is obtained.Since the roof plate was not destroyed before the initial fracture, the thin plate was the surrounding solid support model.The slat profile along the orientation of the thin top plate is plotted to obtain the first fracture principal tensile stress (blue line).The principal compressive stress (red dashed line) is symmetrically distributed about the neutral axis of the slat.The principal tensile and principal compressive stresses are perpendicular to each other at any point.The main compressive stress is near the upper edge of the slat, the main tensile stress is near the lower edge of the slat, and the main tensile stress shows a "funnel-shaped" inclined distribution on both sides of the slat.The results are shown in Figure 8.After the initial fracture of the roof plate occurs, with the mining of the working face, a thin plate model with three sides solidly supported and one side supported is established.Meanwhile, the slat profile along the strike of the thin top plate is drawn to obtain a symmetric distribution of the principal tensile stress 1  (blue line) and the principal compressive stress 3  (red dashed line) in the top plate at cyclic fracture about the neu- tral axis of the slat.The main tensile stress is near the upper edge of the middle of the slat, the main compressive stress is near the lower edge of the middle of the slat, and the writing cracks are formed by the concentrated distribution of the main tensile stress in the slat, showing a "funnel-shaped" inclined distribution.The results are shown in Figure 9.The angle between the principal tensile stress and the horizontal plane in the initial and periodic fractures is The angle between the fracture trace and the horizontal plane can be expressed as

Engineering Geology Background
This  As a result of the coal seam's slant, the bedrock above the opening and closing line experiences a loose strata load with a variation in thickness of 102.7 m following the extraction of the coal seam over a considerable mining distance, with the results shown in Figure 11.If we disregard the impact of the coal seam inclination angle and consider the coal seam and overlaying rock strata as a horizontal plate, applying an approximately equal load throughout the entire plate will increase inaccuracy.The overburden load on the lower seam cannot be accurately estimated as perpendicular to the thin plate due to the degree of inclination of the coal seam.The angle of inclination of the coal seam will correspond to a similar angle in the system.The overburden load can be divided into two components: stresses perpendicular to the thin plate and stresses parallel to the thin plate's central surface.The subsequent sections provide an analysis of the impacts associated with inclined coal seam mining and horizontal coal seam mining.Analyze the rock distribution and physical properties of rock strata in a drill hole data above the working face and determine the distribution of key strata of bedrock above the working face.The information of each rock strata number, name, thickness, depth from the surface, rock strata classification schematic, and key strata distribution were obtained from drilling, as shown in Figure 12.
Concerning the physical properties of each rock strata in the actual borehole measured by the mine and related studies, the parameters of each rock stratum in the above borehole, such as density, tensile modulus, shear modulus, elastic modulus, cohesion, internal friction angle, and Poisson's ratio, were obtained as shown in Table 1.

Deflection Function Based on Minimum Potential Energy Inclined Key Strata Thin Plate
The Equation of static equilibrium by elastic mechanics can be derived with the help of the principle of minimum potential energy.The expression of minimum potential energy is: where U denotes the strain energy of the thin plate, and W denotes the work by the external force on the thin plate.
In this paper, the strain energy of the inclined critical laminate is equal to the total strain energy under the combined action of lateral load and mid-plane force, which can be expressed as When the inclined key strata of the thin plate under the elastic foundation, consider the foundation reaction force work, then the external force on the thin plate work is equal to the lateral load and flexible foundation combined force in the vertical direction with the plate work and the work of the troops in the mid-plane along the horizontal path of the plate.
The work by the lateral load and the elastic foundation combined force can be expressed as The work by the mid-plane force along the mid-plane can be expressed as The principle of minimum potential energy states that the actual displacement always makes the total potential energy take the minimum value among all displacements satisfying the displacement boundary.Therefore, the actual displacement should make the total potential energy zero.
3.2.1.Four-Sided Solid-Supported Inclined Key Strata Plate Deflection before Initial Fracture Assume that the deflection function for the bending of a solidly supported thin plate around the perimeter is Strain energy by combined lateral load and mid-plane force in an inclined key strata thin plate model with four-sided solid-supported before initial fracture.
The work by the overlying load transfer model to calculate the work by the overlying load during the bending deformation of the inclined key strata.
The inclined key strata lamellae are subjected to reaction forces in the elastic foundation to do work as The total potential energy in the four-sided solid-supported inclined key strata model is By the principle of minimum potential energy, when , the deflection coefficient of the solidly supported inclined plate deflection is 3.2.2.Three-Side Solid-Supported and One-Side Simple-Supported Inclined Key Strata Plate Deflection before Periodic Fracture Assume that the deflection function of a three-sided solid-supported and one-sided simple-supported thin plate deflection is Strain energy by combined lateral load and mid-plane force in a three-sided solidsupported and one-sided simple-supported thin plate model before periodic fracture.
The overlying load transfer model does the work to calculate the work by the overlying load during the bending deformation of the inclined key strata.
a b mn x Strike y Dip q q q wdxdy a Strike b Dip A ab a Strike q q a Strike q q b Dip q q a S t r ik e b D ip The inclined key strata lamellae are subjected to reaction forces in the elastic foundation to do work as The total potential energy in the three-sided solid-supported one-sided simple-supported inclined key strata thin plate model is )) 3 By the principle of minimum potential energy, when , the deflection coefficient of three-sided solid-supported and one-sided simple-supported is

Calculation of Inclined Key Strata Breaking Distance and Fracture Trace Angle
The direct top collapse is performed without foundation support as the coal seam below the main roof is mined.Thus, the elastic foundation reaction force integral term is removed, and the angle between the main roof's initial and periodic breakage distance and the fracture trace is calculated.Regarding the development of the key strata under the Bed-separation and no foundation support, the overlying key strata's limit break distance and the fracture trace's angle are calculated.By the principle of minimum potential energy, assume that the mining length of the working face is _ L broken , and the overlying load on the inclined plate is The deflection coefficient of the four-sided solid-supported inclined key strata thin plate model is The deflection coefficient of the three-sided solid-supported one-sided simple-supported inclined key strata thin plate model is When the length of the coal seam mining area reaches the breaking distance limit, the main roof's principal stress is equal to the measured ultimate tensile stress of the main roof borehole.From the calculation formula of the principal tensile stress and the deflection coefficient of the thin plate, it is known that After the breakage distance of initial breakage and periodic breakage is known, the breakage angle can be calculated by substituting the formula.
As the working face is mining, the angle between the fracture trace and the horizontal plane, when the solid-supported inclined thin plate is initially broken around, is When the working face continues mining, the three-sided solid-supported one-sided simple-supported periodically fracture; the angle between the periodic fracture trace and the horizontal plane is Because of the rock layers' complex composition, each key layer contains mudstone, coal seam, fine sandstone, siltstone, etc.The angle between the fracture traces of each rock strata and the horizontal plane is different.Each key strata is examined, and the angle between the fracture traces within the rock strata between the key strata and the averages simplifies the working surface.Assuming that the two key strata contain several n strata, the breakage distance and the angle between the fracture line and the working surface of each stratum have been calculated by Equations ( 54)-(56).The average angle between the fracture line and the working surface of the rock strata with the number of n layers between the two adjacent key strata is Referring to the construction process of the large-deflection inclined thin plate model, the flow chart of the model is shown in Figure 13.

Large-Deflection Inclined Thin-Plate Model to Calculate the Dynamic Development of Key Strata Deflections and Bed-Separation
The primary strata of interest within the extracted rock formations are the hard rock strata, which significantly influence the alteration of stress fields in each rock stratum due to coal seam mining and rock movement.When the primary layers experience displacement, the adjacent softer rock layers undergo deformation in tandem.Similarly, the over-lying, softer rocks may incur complete or partial damage when the primary layers fracture.The strata of primary importance are organized hierarchically, extending from the work face to the surface.The uppermost strata in this hierarchy are called the principal key strata, while the remaining are classified as key strata.Coal seam mining operations often result in varying degrees of deflection, which can be attributed to the disparities in the elastic modulus, load-bearing capacity, and foundation support conditions of individual important strata.
Consequently, this leads to a bed separation underneath each stratum.A study on mining inclined coal seams established a large-deflection inclined thin-plate model for each key strata.The objective was to investigate these key strata's dynamic growth pattern of deflection and bed separation.
Set the corner point of 7225 work face in the direction of the open cutting line and inclined uphill as the origin of the regional coordinate system, the direction of work face mining advance as the x-axis, the tip pointing downhill from the work face as the y-axis, and the direction perpendicular to the xOy plane and pointing to the surface as the z-axis.From the 7225 working face inclination parameters, it is known that compared with the 1954 Beijing coordinate system, the angle between the x-axis direction of the regional coordinate system and the horizontal plane is the working face strike inclination angle of 7.2°.The angle between the y-axis direction of the regional coordinate system and the horizontal plane is the working face inclination downhill slope of 8°.

Breakage of the Main Roof and the "Vertical Three Zones" Height Development
Liu Tianquan put forward the theory of "vertical three belts" of rock movement from horizontal and vertical directions and initially elaborated on the law of overall breakage and movement of rock strata.The vertical direction divides the overlying rock strata of the mining area into, from bottom to top: caved zone, fractured zone, and continuous bending zone [39].Similarly, McNally et al. divided the overburden of the coal seam into "three zones" after mining, replacing the continuous bending zone with the elastic zone [40].Forster divided the overlying rock strata into "four zones," which are caved zone, fractured zone, mobile continuous bending zone, and topsoil layer.Holla et al. divided the overlying rocks into "four zones," called the complete caved zone, fractured zone, expansion zone, and mobile continuous bending zone [41].Forster estimates the height of the entire caved zone as 5 times the mining thickness and the development of the fractured zone as 21 to 33 times the mining thickness.Chinese scholars give the empirical formula for calculating the height of the caved zone, a function related to the mining thickness and the rock expansion coefficient of the caved zone.The procedure for calculating the height of the caved zone in practice is as follows: The caved zone and the fissured zone are called two zones, also called the waterconducting fissured zone, and the water in the two zones can flow into the mining area with the fissures.It is essential to study the development height of the hydraulic fractured zone when carrying out Bed-separation grouting and filling mining, which must be designed to avoid slurry flow from the fractured zone into the mining area.The formula for calculating the fractured zone height in practice is During the initial stages of mining operations, the primary roof undergoes tensile stress and surpasses its maximum tensile strength, resulting in tensile failure in the main roof.This event marks the advancement of the working face up to the point of initial pressure exertion on the main roof.The main roof periodic pressure step refers to advancing the working surface distance following the first pressure on the main roof period, accompanied by stretch breaks.The main roof's initial and recurrent pressure steps can be calculated by considering the initial and periodic braking distances of the large-deflection inclined thin-plate model of the main roof.Due to the predominance of soft rocks or thin layers of hard rocks in the uppermost strata of the primary roof, the tensile strength is diminished, resulting in a reduced braking distance compared to that of the main roof.
Consequently, the two roofs are prone to collapse simultaneously.The siltstone strata located in the 15th layer of the drill hole are situated at a distance of 49 m from the extraction region.These strata are found within the fractured zone, as determined by using the empirical formula for computation.Determining the height of fractured zone development involves the computation of the breakage distance of the primary layers, the fracture trace of the 15th siltstone stratum, and the average fracture trace of the overlaying rock strata.The utilization of the large deflection achieves this inclined thin-plate model.Table 2 presents the data regarding the braking distance, fracture trace, and coal seam angle of the primary strata and additional soft rock seams located in the bubble fall zone and fracture zone above the coal seam resulting from mining activities at the 7225 work face.The findings indicate that the initial separation distance between the key strata and the soft rocks exceeds the recurring distance.Additionally, the angle formed by the initial fracture trace and the coal seam is larger than that formed by the recurring fracture trace and the coal seam.The angular deviation between the initial and periodic fracture traces in soft rock formations is larger than the corresponding deviation in hard key strata.According to the data in Table 2 and Figure 14, the results of the analysis are shown in Table 3, and the revealed height of the collapse zone in the overburden movement is 38.07 m after mining 7225 working face.Additionally, the fractured zone exhibits a development height of 41.13 m.The thick, loose layer was assumed to be laterally loaded, and the height of the overall bending zone was assumed to be the remaining height of the bedrock thickness minus the caved and fractured zones.The "Vertical Three Zones" heights calculated for the large-deflection inclined thin plate are included in the interval estimated by the empirical formulas, indicating that the model's calculated values are highly accurate and meet the engineering requirements.The initial and periodic pressure steps of the main roof demonstrate a strong alignment with the measured pressure steps seen in the coal mine.This demonstrates that the model exhibits a significant level of precision and may be effectively utilized in examining primary roof pressure and the vertical extent of the "upper three zones".As the inclined working face descended, the rock strata within the collapsed zone experienced breakage and yielded.In contrast, the rock strata within the fractured zone underwent fragmentation, forming fractured inverted steps.The angle between the fracture trace resulting from the fracture inversion step and the inclined coal seam is determined by referencing Table 2. Figure 14 illustrates the schematic representation of the regional partitioning of the collapsed, fractured, and continuous bending zones above the mining area.Additionally, it depicts the progression of the deviated strata beneath the primary key strata within the continuous bending zone.The second layer of mudstone is situated directly above the coal seam, while the third layer consists of fine sandstone and serves as the primary roof over the coal seam.
Additionally, a section of the soft rocks between the fourth and fourteenth layers is encompassed within the collapsed zone.The fractured zone encompasses a portion of the soft rocks spanning from the 4th to the 14th strata, the siltstone found in the 15th stratum, and a section of the soft rocks extending from the 16th to the 23rd strata.The continuous bending zone encompasses several strata contributing to the overall surface area.Notably, the 24th strata consists of siltstone and is the third key.Additionally, the 39th strata comprises fine sandstone and is considered the principal key layer.Lastly, the 51st strata consists of thick, loose strata, further contributing to the overall composition of the region.
Within this context, it is observed that the fine sandstone present in the collapsed zone's primary layer has bubbled.Similarly, the siltstone found in the primary layer of the fractured zone has experienced fracturing.Furthermore, the area beneath the bed separation zone has undergone compaction.The siltstone found in the third key strata, and the fine sandstone in the 39th principal key strata inside the continuous bending zone exhibit the presence of bed separation cavities beneath the principal key strata.Furthermore, these bed separation cavities below the principal key strata demonstrate a greater height and volume of development.Beneath the primary roof structure lies the secondary roof and the mining zone, which has seen a loss of support from the underlying elastic foundation.The Bed-separation voids inside the cracked zone undergo prompt compaction after their formation.When employing the energy technique to analyze the key strata inside the fractured zone, it is necessary to consider the contribution of the elastic foundation term to calculate the work performed and predict the deflection of the key strata.The elastic foundation work term can be eliminated using two essential techniques: creating Bed-separation cavities beneath each critical stratum inside the continuous bending zone and utilizing large-deflection inclined thin plate modeling.
. With the known height of the caved zone and fractured zone and the angle between the fracture trace and the coal seam, it is concluded that after 96 m of mining, the fracture trace of the inverted step extends to the continuous bending zone, and the overall deflection of the continuous bending zone begins to appear.The Bed-separation cavity below the principal key strata is caused by the non-synchronous bending of the principal and third key strata.
Continuous bending zone key strata starting deflection position is equal to the quotient of the height of each key strata and soft rocks in the caved zone and the fractured zone with the tangent of the angle between the corresponding initial fracture trace and the coal seam: Continuous bending zone key strata ending deflection position is equal to the quotient of the height of each key strata and soft rocks in the caved zone and fractured zone with the tangent of the angle between the corresponding periodic fracture trace and the coal seam: The rock layers within the collapsed region are undergoing subsidence.In contrast, in the fractured region, the offsets are present for a shorter duration due to the fracturing and compression of the critical rock layers.In the context of work face mining, it is observed that the continuous bending zone above the coal seam consists of two significant strata: the 24th strata siltstone and the 39th strata fine sandstone.Among them, the 39th strata fine sandstone is the key stratum.The 24th strata siltstone, the third significant strata, is situated within the continuous bending zone and can be chosen as the primary strata for grout filling mining.
Nevertheless, the proximity of its lower section to the fractured zone raises concerns regarding the potential impact of high-intensity grouting activities on the soft rocks beneath the 24th stratum.This could lead to the inadvertent flow of slurry into the fractured zone and subsequently into the mining region, disrupting regular mining operations.
Hence, this section primarily addresses the height and volume of the bed separation cavity beneath the key strata that undergo dynamic development during the mining process.
In order to investigate the correlation between the deflection of principal key strata, the dynamic development of bed separation, and the surface subsidence, we referred to the comprehensive observation data of surface observation points at work face 7225 and the corresponding mining line location.The research utilized the dates of 15 September 2015, 27 October 2015, 2 December 2015, 26 December 2015, and 15 April 2016, along with the matching working face mining distances.Table 4 presents the research date, which includes the mining distance of the working face, the location of the start and finish of deflection, the maximum deflection value, and the results of the deflection space volume of the primary essential strata.The maximum deflection value and the significance of the deflection space of the primary key strata exhibited a continual increase.After mining, the most excellent deflection value reached a magnitude of 2714 mm, creating a voluminous space measuring 56,485 m 3 .Based on the deflection function of the key strata, the mining work face at depths of 135 m, 300 m, 414 m, 535 m, and 651 m exhibits deflection in the significant vital strata.This deflection is observed in the three-dimensional expansion of the regional coordinate system.The dynamic representation of these results is depicted in Figure 15, which showcases the time series of the deflection.When the distance between the working face and the mining area is 135 m, the primary key strata exhibit a maximum deflection value of 9.6 mm.The primary key strata do not exhibit discernible deflection.As the distance of work face mining rises, there is a corresponding increase in the deflection range of the primary key stratum, which expands rapidly in the downward direction.
Additionally, the maximum deflection value experiences a rapid increase.After excavating, the working face reaches a depth of 535 m; the primary key strata exhibit a maximum deflection value of 2683.2 mm.Furthermore, it is seen that the maximum deflection value experiences a gradual increase over time.At a mining depth of 651 m, when the working face intersects the closure line, the primary key strata exhibit a maximum deflection value of 2714 mm, thereby attaining its peak deflection.The deflection curve of the Bed-separation below the principal key strata in the continuous bending zone is equal to the difference between the deflection curve of the third key strata and the deflection curve of the principal key strata: The volume of the Bed-separation cavity below the principal key strata in the continuous bending zone can be obtained: The work face has a width of 7225 units and a length of 149 m, with a mining distance ranging from 0 m to 96 m.Within this range, the key strata exhibit continuous bending without any deflection.However, beyond a mining distance of 96 m, the key strata in the continuous bending zone experience deflection deformation.Figure 16 displays the maximum development height of the bed separation cavity underneath the primary key strata within the overlying rock of the 7225 workings to the mining distance.The absence of bed separation below the primary key strata is observed when the working face is advanced from 0 m to 96 m.The increase rate in the bed separation's maximum height is significant, as it reaches a maximum value of 545.2 mm when the working face progresses from 96 m to 579 m.The vertical extent of the bed separation diminishes gradually as the working face progresses from a distance of 579 m to 651 m.

Verification of the Principal Key Strata Deflection and Lower Delamination Accuracy Based on 3DEC
The rock consists of a fundamental continuum of rock blocks and non-continuous geological structure surfaces.In a non-equilibrium stress field, rocks can exhibit mechanical behavior like a continuous medium.The structural surfaces between individual rock pieces facilitate this mechanical interaction.When the structural surface experiences stress that exceeds its bearing capacity, the rock blocks exhibit vertical misalignment with each other.Three-dimensional Discrete Element Code (3DEC, v7.0) is a computational software for the numerical analysis of discontinuous models using the discrete element method, corresponding to the Mathematics Subject Classification (MSC) code number "65-04".
The Mohr-Coulomb model has been selected as the intrinsic model for the rock blocks.The physical and mechanical properties used as input are listed in Table 1.These properties primarily consist of bulk weight, elastic modulus, tensile modulus, shear modulus, cohesion, internal friction angle, Poisson's ratio, and others.The Mohr-Coulomb slippage model was chosen for the nodal ontogenetic model of the surface structure between the rock blocks.A discrete element model with three-dimensional characteristics is constructed, spanning from the uppermost plate of the 7225 working face to the surface.The model boundary is defined as the spatial extent of the basin undergoing surface movement.It is postulated that the vertical movement of the lower plate is negligible on the zaxis, and the force of gravity influences the surface.The velocity in the x and y-directions is zero, indicating the absence of any regular movement.The inside zone corresponds to the mining area encompassing the 7225 work face, where grid encryption measures are used.On the other hand, the exterior region refers to the non-mining area characterized by a somewhat lower density of grid placement.The transparent bottom plate was implemented to enhance the visibility of the coal seam excavation region, as depicted in Figure 18.Similar to the depiction in Figure 16, the mining distance for the 7225 work face was intentionally set at distances of 300 m, 414 m, 535 m, and 651 m.Subsequently, a numerical simulation using the 3DEC software was conducted to assess the deflection of the primary critical strata.The findings are depicted in Figure 19, wherein the red hue represents the downward deflection of the key strata, the blue hue represents the upward uplift of the key strata, and the black rectangle box denotes the designated mining range of work face 7225.After excavating 300 m, the primary strata of interest located above the mining area are predominantly observed along the x-axis within the range of 0 m to 280 m.Notably, higher deflection values are observed near the open cut line, while two minor deflections are observed between 280 and 651 m.The primary strata above the active mining area are predominantly observed along the x-axis within the interval of 0 m to 349 m, following a mining depth of 414 m.
Nevertheless, these surfaces are linked to a pair of deflection structures close to the anticipated termination point of the primary geological layers above the active mining area.Following the extraction of 535 m, the primary strata above the active mining area exhibited a collective deflection, with the deflection surface extending farther downhill towards the 7225 work face.Upon extracting a total of 651 m, it has been observed that the primary strata above the 7225 working face exhibit a fully developed deflection surface.Furthermore, the boundary of this deflection surface extends towards the working face, displaying an inclination in the downhill direction.The extent of the deflection surface in the x-axis direction is limited to the boundary between the open cut line and the finish line.The 3DEC program simulates the deflection range of the major key strata above the 7225 working face in mining operations.The simulated deflection range closely aligns with the deflection value determined by the big deflection inclined thin plate model.The occurrence of inexplicable uplift outside the designated mining area for the 7225 work face can be attributed to the correctness of the deformation zone setup.

Influence of Work Face Mining Length on the Development Height and Cavity Volume of the Bed-Separation
The dynamically evolving nature of bed separation beneath the PKS horizons within the continuous bending zone leads to its non-development during the initial phase.Subsequently, the vertical distance between the bed separation experiences a notable escalation till it attains its furthest magnitude, coinciding with the extraction of the mining front.Following a specific time frame, the progressive expansion of mining distance continually increases the load exerted on the primary key strata.Additionally, the height of the lower bed separation gradually declines until it reaches a point of closure.
This study investigates the correlation between the maximum height of the bed separation development and the mining distance.It is anticipated that the mining operations will continue on the 7225 working face, and the mining distance of this working face will be extended to 3000 m. Figure 20   Given the hypothetical scenario where the working face persists in its mining activities and decides to expand the mining distance to 3000 m, it is imperative to investigate the comprehensive correlation between the volume of the bed separation cavity and the mining distance.The findings are depicted in Figure 21.The volume of the bed separation cavity exhibits a progressive increase until it reaches its maximum value of 21,379.4m 3 when the mining distance of the working face is less than 1917 m.Subsequently, when the mining distance extends, there is a gradual reduction in the volume of the bed separation cavity.Specifically, at a mining distance of 3000 m, the volume of the bed separation cavity is observed to be 18,339.8m 3 , which is equivalent to the volume observed at a mining distance of 1170 m.

Influence of Working Face Width on the Development Height and Cavity Volume of the Bed-Separation
The utilization of various mining widths on the 7225 working face can lead to adequate and inadequate mining tendencies.Consequently, this will result in variations in the starting and ending positions of the main roof's initial and periodic braking distances and the deflection caused by the thin plate of the principal key strata.These factors, in turn, will impact the height and volume of the cavity in the overlying bed separation that develops beneath the principal key strata.The design mining width of the 7225 working face is 150 m.It is hypothesized that the mining width might be 120 m, 180 m, 210 m, or 240 m.Table 5 displays the parameters encompassing the initial breaking distance, periodic breaking distance of the main roof, starting and ending position of the deflection of the principal key strata thin plate, as well as the height and cavity volume of the overlying Bed-separation of the principal key strata.Based on the assumption that the mining length of the 7225 working face remains constant at 651 m, the correlation between the development height of the overlying Bedseparation under the principal key strata and various design mining widths is established using the large-deflection inclined thin plate model.The outcomes of this analysis are depicted in Figure 22.As the mining breadth expands to either 120 m or 150 m, the overlaying bed separation beneath the primary key strata undergoes continuous development as the mining length increases.When the mining width is reduced to 120 m, the growth rate of the overlaying stratum is seen to be slower compared to a mining width of 150 m.The maximum height of bed separation development is smaller in the former case.When the design mining width is set at 180 m, 210 m, or 240 m, the bed separation quickly develops and achieves its maximum height within a short distance after the major key strata above the 7225 working face deflect.Subsequently, when the displacement of the underlying key strata is sustained, the underlying rocks progressively support the primary key strata.The progressive obstruction and subsequent contraction of the bed separation ultimately culminate in its total closure underneath the PKS.Based on the assumption that the mining length of the 7225 working face remains constant at 651 m, the correlation between the volume of the cavity in the overlying Bedseparation under the principal key strata and various design mining widths is established using the large-deflection inclined thin plate model.The findings are visually presented in Figure 23.As the width of the design reaches 120 m, an increase in the volume of the bed separation cavity located above the primary key strata is observed as the mining length progresses.However, the maximum volume is one-third of the volume achieved when the mining breadth is 150 m.When the design mining width is set at 180 m, 210 m, or 240 m, the development of the bed separation cavity occurs swiftly once the major key strata above the 7225 working face deflect.Furthermore, the volume of the cavity exceeds the mining width of 150 m until the mining length reaches 280 m.Nevertheless, after the delaminated cavity reached its maximum volume, the volume of the bed separation cavity rapidly decreased to zero due to the swift closure of the bed separation.

Conclusions
This study proposes a large-deflection inclined thin plate mechanics model for calculating PKS deflection and Bed-separation induced by an inclined working face mining activity covered by a thick loose layer.The model was applied to quantify and evaluate the dynamic development of Bed-separation during Longwall mining in the 7225 inclined workings of a mine in Anhui Province, China.As previously mentioned, this paper provides the following findings: 1. Based on the drilling data obtained above the 7225 working face, each rock stratum's physical and mechanical properties from the coal seam to the ground were determined.The thick, loose layer overlying the 7225 working face was determined to contain four key strata, numbered 3, 15, 24, and 39.It was also determined that the main roof plate was No.3, 1.47 m above the 7225 work face, and the PKS was No.39, 171.99 m above the 7225 work face; 2. A model of a large-deflection inclined thin plate with one unsupported side and three supported sides was created by mining the 7225 inclined face.Then, the initial breakage distance of the main roof plate was 33 m, and the initial breakage angle was 61.2°; the periodic breakage distance was 21 m, and the periodic breakage angle was 55.4°.
The distribution range of the "Vertical Three Zones" from the 7225 face to the ground surface included the height of the caved zone of 38.07 m, the height of the fissured zone of 41.13 m, the height of the bending zone of 187.56 m excluding the thick loose layer.It also analyses the relationship between the spatial location of each key strata's starting and ending deformation points within the "Vertical Three Zones" and the mining activities' distance; In the future, the model will be applied more to the dynamic assessment of bed-separation development caused by horizontal or inclined coal seam mining, and it will provide a theoretical basis for the rational use of Bedseparation cavities, such as grouting backfill.

U2
The additional strain energy caused by the mid-plane bending strain k the foundation reaction coefficient

Figure 3 .
Figure 3. Diagram of the two-parameter foundation model.

Figure 4 .
Figure 4. Diagram of distribution of bending moment and shear force in a thin plate.

. 3 . 3 .
Strain Energy Generated by the Bending of the Plate by the Action of the Mid-Plane Stresses N is known, the strain component on the mid-plane of the plate can be known by Hooke's law as ( The results are shown in Figure6.

Figure 6 .
Figure 6.Additional shear strain generated by lateral load bending the mid-plane.

Figure 8 .
Figure 8. Diagram of the distribution of the main stress traces and fracture shape of the initial bedrock fracture.

2 Figure 9 .
Figure 9. Diagram of the distribution of main stress traces and fracture shape of bedrock periodic fracture.
study selects the 7225 working face in Huaibei, Anhui Province, with a strike length of 651 m, inclination length of 150 m, coal seam strike angle of 7.2°, an inclination angle of 8.1°, and coal seam mining thickness of 3.7 m.The elevation of the corner point of the working face cutting line and the uphill boundary is −410.8m, the corner point of the finishing line and the downhill boundary line is −513.5 m, and the ground elevation is 28 m.The thickness of bedrock is 266.76 m, and the thickness range of thick loose strata is 172.04 m to 274.74 m.The spatial distribution and geological structure of the working face are shown in Figure 10.

Figure 10 .
Figure 10.Diagram of the geographical location and geological structure of the study area.

Figure 11 .
Figure 11.Analysis of different stratigraphic models and rock seam forces by coal seam inclination.

Figure 14 .
Figure 14.Bending and breaking of key strata caused by mining, the regional division of " Vertical Three Zones", dynamic development of Bed-separation.

4. 2 .
Deflection of the Principal Key Strata and the Dynamics of the Height and Volume of Bed-Separation Development With mining the 7225 work face, the initial fracture trace and the periodic fracture trace extend upward from the open cut line position from the main roof position to above the fractured zone.The starting position of each key strata in the continuous bending zone is _ L Start Coordinatsystem from the origin, and the ending position of each key strata in the continuous bending zone is _ min tan _ L ing dis ce L End Coordinatsystem  .After the mining distance _ min tan L ing dis ce of the work face, as the 7225 work face takes Longwall mining, the bending width of the principal key strata sheet of the continuous bending zone and the inclined width of the work face is approximated as b.The actual bending length of the principal key

Figure 15 .
Figure 15.Diagram of the three-dimensional expansion of the principal key strata deflection face in the regional coordinate system (working face mining distance 135 m, 300 m, 414 m, 535 m, 651 m)

Figure 16 .
Figure 16.Variation of the maximum height of Bed-separation development below the principal key strata with the advancing distance of the 7225 work face.The dynamic variation results of the volume of the off-bed cavity below the principal key strata in the overlying rock strata caused by the 7225 working face mining distance are shown in Figure 17.When the working face advances from 0 m to 96 m, the volume of the Bed-separation cavity is 0 m 3 because the Bed-separation is not developed below the principal key strata.As the mining distance increases from 96 m to 651 m, the volume of the Bed-separation cavity increases and reaches the maximum value of 11,228.1 m 3 when the mining distance is 651 m.

Figure 17 .
Figure 17.Volume of the Bed-separation cavity under the principal key strata with the advancing distance of the 7225 working face.

Figure 18 .
Figure 18.3DEC program 7225 working face excavation deformation area design model.

Figure 19 .
Figure 19.3DEC numerical simulation program calculates the deflection value of the principal key strata (work face mining distance 300 m, 414 m, 535 m, 651 m).
displays the outcomes of the diagram illustrating the fluctuations in the maximum height of the bed separation development to the primary key strata, as influenced by the mining distance.Once the mining distance of the working face surpassed 579 m, there was a gradual decrease in the development height of the bed separation under the primary important strata.The height of the bed separation development experienced a reduction to 169.4 mm upon reaching a mining distance of 3000 m.

Figure 20 .
Figure 20.Diagram of the maximum height variation of Bed-separation developed under the principal key strata assuming a mining distance of 3000 m.

Figure 21 .
Figure 21.Diagram of the assumed mining distance of 3000 m for the volume change of the Bedseparation cavity under the principal key strata.

Figure 22 .
Figure 22.Influence of different mining widths on the development height of the overlying Bedseparation under the principal key strata.

Figure 23 .
Figure 23.Impact of different mining widths on the volume of the overlying Bed-separation cavity below the principal critical strata.

3 .
The large-deflection inclined thin plate model was applied to the 7225 inclined working face mining process.The PKS deflection value develops from 0 mm to 2714 mm with the mining of the working face, and the maximum value of the space volume is 56,485 m 3 , verified by 3DEC numerical simulation.The dynamic development of bedseparation below the PKS within the overburden strata, the maximum development height is 545.2 mm, and the maximum value of Bed-separation cavity volume is 11,228.1 m 3 .Finally, the dynamic development of bed-separation height and cavity under the working face's different mining length and width conditions is discussed.As the length of the working face increases, the Bed-separation height increases and then decreases, and the maximum value of the development height is 545.2 mm, corresponding to the mining distance of 579 m.The volume of the Bed-separation cavity increases until it starts to decrease after 1917 m.The Bed-separation height and cavity volume show an increasing trend when the working face width increases within 150 m.After the working face width exceeds 180 m, the Bed-separation cavity develops and closes rapidly; 4. The large-deflection inclined thin plate model proposed in this paper can effectively explore the dynamic development of PKS deflection and Bed-separation induced by Longwall mining inclined workings overlying thick loose layers.Of course, the largedeflection inclined thin plate model has limitations in its application to the analysis of overlying stratigraphic deflection caused by open-pit mining and partially sharp inclined working face mining.

1  3 WWThe minimum potential energy 1 qKey strata deflection coefficient 2 qNN
Key strata principal tensile stress N A vector sum of the gravity along the inclined plate Key strata principal compressive stresses MX Thin plate bending moment in the x-axis  The angle between the principal tensile stress and the horizontal plane MY Thin plate bending moment in the y-axis  The angle between the fracture trace and the horizontal plane MXY Thin plate torque W The work by the external force on the thin plate QX Lateral vertical shear in the x-axis v The work by the lateral load and the elastic foundation combined force QY Lateral vertical shear of the y-axis h The work by the mid-plane force along the mid-plane w Thin plate deflection p E Lateral load at the junction point between the uphill direction and the cut line mn A Lateral load at the junction point between the downhill direction and the finishing line x Unit length mid-plane tension in the x-axis tensile  Key strata ultimate tensile stress y Unit length mid-plane tension in the y-axis _ initial  The angle between the key strata initial fracture trace and the horizontal plane xy N Unit length mid-plane shear _ periodic  The angle between the key strata periodic fracture trace and the horizontal plane ux Small displacements of mid-plane in the xaxis _ initial average  The average angle of the initial fracture trace between the two adjacent key strata uy Small displacements of mid-plane in the yaxis _ periodic average  The average angle of periodic fracture trace between the two adjacent key strata ' x Total unit elongation in the x-axis _ h caved zone The height of the caved zone ' Key strata tensile stress in the x-axis _ h fractured zone The height of the fractured zone y Key strata tensile stress in the y-axis _ min tan L ing dis ce The mining distance of the work face xy  Key strata shear stress _ L Start Coordinatsystem Continuous bending zone key strata starting deflection position D Bending stiffness of the thin plate _ L End Coordinatsystem Continuous bending zone key strata ending deflection position v Poisson's ratio of thin plate _ Bed separation w The deflection curve of the Bed-separation a Length of plate _ Bed separation V The volume of the Bed-separation cavity b Width of plate

Table 1 .
Reference table of physical properties of each rock strata in drill holes.Figure 12. Work face drilling data and key strata judgment.

Table 2 .
Break distance and fracture trace angle between key strata and other rock strata in the caved zone and fractured zone.

Table 3 .
Comparison of empirical formulas and modeling results for the development height of the "Vertical Three Zones".

Table 4 .
Mining distance of the working face corresponding to the research date and deflection parameters of the principal key strata.

Table 5 .
Study date corresponding working face mining distance and surface subsidence basin inscription parameters.