Journal of Applied Mathematics

Journal of Applied Mathematics / 2013 / Article
Special Issue

Mathematical and Numerical Modeling in Geotechnical Engineering

View this Special Issue

Research Article | Open Access

Volume 2013 |Article ID 859803 |

María-Belén Prendes-Gero, José Alcalde-Gonzalo, Pedro Ramírez-Oyanguren, Francisco-José Suárez-Domínguez, Martina-Inmaculada Álvarez-Fernández, "Longwall Mining Stability in Take-Off Phase", Journal of Applied Mathematics, vol. 2013, Article ID 859803, 12 pages, 2013.

Longwall Mining Stability in Take-Off Phase

Academic Editor: Ga Zhang
Received06 Jun 2013
Accepted23 Jul 2013
Published28 Aug 2013


Mechanised longwall mining is one of the more commonly employed exploitation methods in underground mines in the north of Spain as well as in the rest of the world. It is continuously changing and evolving, with new techniques, technology, equipment, and face management practices and systems appearing for the purposes of improving aspects such as operational and financial performances and, above all, the safety of the miners. Despite its importance, there are no regulations for the mining of longwall coal seams. This work aims to contribute to an advance in the design and optimisation of the roof support in longwall mining, analysing the stability of the roof using a method based on the resistance of materials, which considers the characteristics of the properties of the roof materials. The influence of not only the individual elements of support but also the coalface, which is considered one more supporting element, is investigated. The longitudinal and transverse spacings of the support and the number of walkways constituting the exploitation panel are analysed. The proposed formulation is validated by information gathered in a mine located in the region of Castilla-Leon.

1. Introduction

Important reserves of coal exist in Spain, which can contribute to reducing the energetic dependence on other countries [1]. Though their importance has diminished gradually, the National Plan of Strategic Reservation of Coal 2006–2012 and the New Model of Integral and Sustainable Development of the Mining Regions [2] established a minimal level of production to keep open the possibility of relying on the coal in case of crisis or a considerable increase in the price of crude oil. For Spain, this level of production was established as 9.12 million tons in 2012, and though at present the new National Plan of Coal (2013–2018) is being developed, the principal aim is to allow the continuity of the Spanish developments.

Obtaining the existing reserves requires using levels of mechanisation of labour that allow economic competitiveness as well as safety within the context of the regulation RGNBSM (General Regulation of Basic Procedure of Mining safety, 1985) [3]. In this regard, methods of excavation in longwall mining using individual hydraulic support elements (mechanical props) have special importance (Figure 1). This mining method is commonly used in the coal-bearing deposits of Castilla-Leon in the north of Spain [4].

But longwall mining is not a new approach to coal mining. In fact, the basic principles of longwall mining have been traced back to the latter part of the 17th century, to Shropshire and other counties in England [5]. In the United States, longwall mining is actually 50 percent of the underground coal production [68], with 49 operating longwalls producing over 175 million tonnes per year [7]. Besides USA, among the most productive longwalls are those in Australia and the Shenhua coalfield in China. Shenhua was developed from a green-field site in 1992 and now has more than seven major underground mines producing over 200 million tonnes per year. Australia has 29 operating faces with a total longwall production of 47.5 million tonnes.

2. Engineering Background

Longwall mining [9, 10], for panels of coal, is an exploitation method where an extended wall of coal is mined in a sole slice. The wall or longwall is around 150–300 m long and 1-2 m thick, and the slab of coal that is being mined is around 1000–3500 m long and 150–300 m broad. The coal is removed from the wall by shearing machinery, which travels back and forth across the coalface, and load the coal onto a conveyor belt that carries it out of the mine. The area immediately in front of the coalface is supported by a series of hydraulic roof supports, which temporarily hold up the roof strata and provide a working space for the shearing machinery, the face conveyor and of course the mines. After each slice of coal is removed, the hydraulic roof supports, the face conveyor, and the shearing machinery are moved forward.

From that moment, the roof immediately above the seam is allowed to collapse into the void that is left as the face retreats. Miners working along the coalface, operating the machinery, are shielded from the collapsing strata by the canopy of the hydraulic roof supports. As the roof collapses into the goaf behind the roof supports, the fracturing and settlement of the rocks progresses through the overlying strata and results in sagging and bending of the near surface rocks and in some cases subsidence of the ground above [11].

Mechanised longwall mining is ever changing and evolving with new techniques, technology, equipment, face management practices, and systems appearing as a direct means to continually improve all aspects of operational and financial performances [12].

Nevertheless, and in spite of the importance of this method of exploitation, no legal regulation exists regarding the types of support and thickness and characteristics of the rock mass surrounding the excavation, which would guarantee the safe functioning of these developments.

Calculation of the pressure that the working roof wall is to exert on hydraulic props is essential for support design, both to ensure working global stability and to avoid prop punching on gables [13]. These two concepts are analysed in multiple papers, as are the behaviour of the overlying rock strata and the performance of the support. R. Singh and T. N. Singh [14] verified the influence that the additional load on the chock shield was exercising by the broken rock mass in thick seams, with values of limiting span given by the clamped and cantilever beam equations. Ju and Xu [15] analysed the strata behaviour during the operation of great mining height, studying the structural characteristics of key strata (cantilever beam and/or voussoir beam) as well as the movement law. Bilim and Özkan [16] analysed the effect of excavation schedules on the overlying rock strata and the supports. González-Nicieza et al. [17] and Juárez-Ferreras et al. [18, 19] analysed the maximum pressure that the coal hanging wall and footwall are capable of supporting, as well as the density of the props and the conditions of the supports they are resting on so that penetration of the props does not occur.

This work tries to contribute to an advance in the design and optimisation of the roof support in longwall mining, analysing the stability of the roof with a method based on the resistance of materials [20] from the characteristics of the roof materials.

In addition, two very important aspects of the safety and the productivity of the exploitation have been considered: the dimensions and the number of walkways of the panel. The width of the rear walkway at waist level needs to be considered from an ergonomic aspect to allow for lamps and self-rescuers worn by underground personnel. The horizontal dimension across an average person’s waist is between 620 and 700 mm, so a rear walkway width of ≥750 mm at waist height is required to permit ergonomically effective passage of people along the face so their productivity is not unnecessarily impeded [12]. The number of walkways impacts the safety of the panels and their productivity due to the fact that it produces variations in the bending moment and deflection laws.

3. Exposition of the Problem

It is considered as a longwall exploitation (Figure 2) if the two dimensions according to the and axes are very large in comparison with the height or dimension along the axis. For simplicity, in this paper we consider a length of the workshop to be equal to 1 m in the direction. The supporting hydraulic elements divide the panel into sections or spans; in each span, the following characteristics are considered constants:(i)length: (m);(ii)Young’s modulus: (Pa);(iii)thickness: (m);(iv)load: (N/m).

In addition, it is considered that the stiffness of the supporting elements can be different and it is possible that there exist other supporting elements with stiffness , for example, keys of wood, the coalface, or the protection rock mass between two panels. So, each span, besides the properties previously indicated, presents the following: (i)stiffness of the uniform support: (Pa);(ii)stiffness of the props: (Pa).

Along with these assumptions is imposed the restriction that within each span the roof rock is perfectly elastic, homogeneous, and isotropic and the neutral axis coincides with the centre line of the thickness.

Considering each roof of a span as a beam subjected to a uniformly distributed load, , acting in the principal plane of the symmetric cross-section (Figure 3), the deflection of this roof is described through a differential equation of fourth degree [20]: where is the deflection in each point and is a constant given by the expression . with the moment of inertia.

The solution of the differential equation (1) is given by where and are given by the expressions , and , respectively, and , , , and , are constants that must be determined using the boundary conditions, which constitute the unknown quantities of the problem.

Once the deflection is obtained, the angle of rotation at any point of the panel is given by the first derivative of the deflection (): The bending moment () and the shear force () at every point of the panel are given by (4) and (5), respectively: As mentioned previously, once the equations that define the problem have been established, it is necessary to determine the boundary conditions (null deflections and null rotations at the ends of the panel) and the conditions of compatibility (equal deflections, rotations, and bending moments in the points between spans) in order to know the constants , , , and .

For this purpose, the panel has been analysed in the phase of take-off by two configurations commonly used in the mines of Castilla-León, with two walkways (panel type 1) and with one walkway (panel type 2).(i)Panel type 1 has two walkways, three elements of support, and six spans (Figure 4), the first and last spans being supported elastically by the coalface;(ii)Panel type 2 has one walkway, two elements of support, and five spans (Figure 5), and as in the previous case, the first and last spans are supported elastically by the coalface.

The application of the boundary and compatibility conditions gives rise to the system of (6); where will be 6 for panel type 1 and 5 the panel type 2. where Replacing the values of the deflection and its derivatives in (6) and expressing the system in matrix form, the following is obtained: where the transposed matrix of is given by where Therefore, the not null elements of the matrices and are obtained (the Appendix) for the panel type 1. In panel type 2, the matrixes and are equal to the first type in its first 14 rows, with the last ones also being equal with the exception of the following changes of indexes: indexes 5 and 6 of the problem type 1 transform into the indexes 4 and 5, respectively, of the problem type 2.

4. Calculating the Factor of Safety of the Workshop

Once the efforts have been calculated, it is possible to calculate the factor of safety (FS) of the panel dividing the tensile strength of each section for the normal stress (14). If this coefficient is bigger than one, the roof of the panel is capable of supporting the stresses to which it is subjected, and, therefore, the work is realised in safe conditions.

To calculate this FS in each section of the roof, the shear stress (11) and the bending stress (12) are calculated as follows: From (11) and (12), the normal stress is calculated as Once the stresses are calculated, the FS is given by

5. Practical Case

A panel in the take-off phase is analysed with two configurations commonly used in the mines of Castilla-León and more specifically in the Feixolin mine [17], with two walkways (panel type 1) and with one walkway (panel type 2). Figures 4 and 5 show that in the phase of take-off, the first and last spans rest elastically on the coal, whereas the halfway sections rest on siltstone with an apparent density of 2.75 t/m3. The hydraulic props are of the type EA 25 manufactured by Salzgitter and reach a maximum extend length of 2.5 m.

To calculate the stiffness of the siltstone, the pressures obtained in the penetration test are considered, and to calculate the stiffness of the props [21], the load-deformation curves from the load plate test are used. All of these calculations were for the Feixolin mine [17]. Table 1 shows the properties used in the analysis.

Section 1
Section 2
Section 3
Section 4
Section 5
Section 6

Length of span (m)
Thickness (m)
Young’s modulus (MPa)355035503550355035503550
Load (kN/m)172001720017200172001720017200
Supporting stiffness (MPa/m)47.847.847.847.847.847.8
Tensile strength (MPa)
Prop number123
Prop stiffness (MPa/m)150015001500
Spacing of props (m)0.550.550.55

One of the most important parameters from the safety point of view as well as in terms of accessibility is the deflection of the roof. In this case, it is considered that negative values of the deflection are equivalent to an increase of the deflection and vice versa.

As it is observed in Figure 6 for both configurations analysed, there is a symmetrical deflection from the centre with a shape very similar to that of a bifixed beam, though in this case the curve does not increase progressively to reach a maximum deflection in the centre. On the contrary, several segments are distinguished in the curve. These segments coincide with the changes of spans of the panel. The first and last spans (1 and 6 in panel type 1 and 1 and 5 in panel type 2) present an increase of the deflection reaching the maximum at the edge of the span, that is, at the edge of the elastic support. These maximum values of deflection are −3.98 mm for panel type 1 and −4.39 mm for panel type 2. At this point, the presence of hydraulic props produces a decrease in the value of the deflection, even from negative values that indicate that the roof turns down to positive values at spans 3 and 4 of panel type 1. These positive values are due to the pressure in the opposite direction to the deflection, which is realised by the supports. The deflection descends to 0.31 mm being null in the middle of the workshop. On the contrary, in panel type 2, the decrease of the deformation is constant until the minimum is reached in the middle point of span 3, with 0.96 mm.

In any case, the deflection of the roof can be considered low and barely affects working conditions.

The fall of the roof, that is, the increase of its deflection, depends directly on the stiffness of the hydraulic props. As it is observed in Figure 7, with a stiffness of the props 100 times less, similar curves are obtained. Nevertheless, the values of deflection in the ends of the first spans are major (−5.35 mm). On the other hand, the decrease of this deflection is minor not reaching positive values at any time. This is due to the fact that in these conditions the capacity of support is for the roof and not for the hydraulic props.

On the contrary, considering an increase in the stiffness of the hydraulic props, it does not suppose a change over the results obtained in the initial conditions (Figure 6). This is because the stiffness of the props initially is so big compared to the stiffness of the rest of the materials that an increase of its value does not produce any effect.

The bending moment (Figure 8) shows symmetrical curves to the deflection ones. In this case, the maximum values for panels type 1 and type 2 are 5.05 kNm and 5.38 kNm, respectively. The maximum values are at the edges of the elastic supports, whereas the minimum values occur at the points of placement of the hydraulic props.

In this case, a decrease of the stiffness of the props does not produce a change in the shape of the curves and only produces small variations in the maximum values obtained (panel type 1: 5.24 kNm, and panel type 2: 5.50 kNm).

A change in the materials of the roof does not alter the shape of the curve, but it modifies the maximum and minimum values, which does not happen with the deflection. This behaviour is due to how the bending moment is obtained (4), multiplying the second derivative of the deflection by a constant for each section, . For spans of equal thickness, is directly proportional to Young's modulus of the material of the span.

The representation of the rotation angle of the spans produces symmetrical curves from the centre of the panel, as much in the case of type 1 as type 2. In Figure 9, it is observed that the use of props produces fluctuations in the rotation, but always reaching lower values to those presented for zones close to the extremes. The maximum rotation is around 4.79 radians in panel type 1 and 5 radians in panel type 2.

The shear force (Figure 10) presents jumps of values in those spans where the hydraulic props are placed. These maximum values (in absolute value) are in the second and last and the values reached as much in panel type 1 as in panel type 2 are very similar and around 20 kN.

Although most of the parameters shown in Table 1 depend on the type of rock, there are some of them that can vary. The type of support, hence its stiffness (Figure 6), as well as the length of the spans, is a compromise between the safety and the productivity of the panel. While the depth of the panel, that is, the load per unit of length, will increase over time, the exploitation advances.

A decrease in the length of the spans (Figure 11) decreases the value of all the parameters analysed: deflection, rotation angle, bending moment, and shear force. However, decreasing this value means, on the one hand, placing a bigger number of props, thus increasing the cost of production, and on the other reduce the space step for miners. In any case, this length cannot be less than 0.75 m considering all the equipment that the miners wear around their body.

An increase in the length of the spans increases the value of all parameters and specifically the value of the deflection. So it is necessary to find a compromise between length and safety/productivity.

With the advance in the exploitation, the depth of the panels grows and therefore the load on spans increases. This increase results in a greater deflection of the roof (Figure 12) and also an increase of the bending moment and of the shear force.

The analysed parameters, and specifically the shear force and the bending moment, let us know one of the most important points in the design of a panel of longwall: the FS. In the two examples with the properties shown in Table 1, the values of the FS are bigger than 5. That is to say, the studied mine is safe for this design. Nevertheless, it is possible to determine in what spans and in what type of panel a minor FS is produced.

From (14), it is possible to deduce that the FS is directly proportional to the tensile strength and inversely proportional to the sum of the bending moment and the shear force. The variation of the values of the bending moment and shear force (Figures 8 and 10) indicates that the most critical points, in the analysis of the safety, are at the edges of the elastic supports (the borders of the first and last span).

The influence of the distribution of the props, their stiffness, or the number of walkways in the panel over the FS has been analysed. Nevertheless, unless extreme values were used, the FS scarcely varied in its value. On the contrary, changes in the properties of the materials, and especially changes in the Young's modulus, produce great variations in the FS.

6. Conclusions

(i)The calculation of the stability of roofs in longwall mining can be resolved by employing the classic resistance of materials. (ii)In addition, because the calculation process is very fast, it is possible to design a more appropriate roof support for a specific longwall mining workshop, to know the behaviour of the roof as a support element, and to analyse the influence and capacity of the individual support elements and the effect of their longitudinal and transverse spacing.(iii)While the most influential factor in the deflection of the roof is the stiffness of the props, the bending moment depends directly on the properties of the materials and specifically on their Young's modulus.(iv)In all cases, the disposition of two walkways in the panel against one walkway reduces the maximum values of the parameters analysed with the exception of the shear force. In any case, these maximum values are placed, for both configurations analysed, in the second and last but one sections. So, these sections are critical in the analysis of the FS.(v)Finally, the variables that have most influence in the FS, are undoubtedly the depth of the panel and above all the changes in the properties of the rock mass of the roof.

7. Highlights

(i)Analyse the influence that the dimensions and the number of walkways of the panel have on the stability of the roof.(ii)Study the influence of not only the individual elements of support but also the coalface, which is considered one more supporting element.(iii)Analyse the configuration of the support (longitudinal and transverse spacing) in longwall mining.(iv)Know in a simple way the safety factor for a workshop.


The not null elements of the matrices and of the problem type 1 are as follows.

Row 1.Row 2. Row 3.Row 4.Row 5.Row 6.Row 7.Row 8.Row 9.Row 10.Row 11.Row 12.Row 13.Row 14.Row 15.Row 16.Row 17.Row 18.Row 19.Row 20.Row 21.Row 22.Row 23.Row 24.


  1. I. Rosal-Fernandez, “How costly is the maintenance of the coal-mining jobs in Europe? The Spanish case 1989–1995,” Energy Policy, vol. 28, no. 8, pp. 537–547, 2000. View at: Publisher Site | Google Scholar
  2. Ministry of Industry, Energy and Tourism. Institute for Restructuring of Coal Mining and Alternative Development of Mining Regions. Report of the national strategic coal reserve 2006–2012 and new model for comprehensive and sustainable development of mining regions, 2013,
  3. RGNBSM, General Regulation of Basic Procedure of Mining Safety. Spain, Ministry of Industry and Energy, 2006, (Spanish).
  4. J. Gómez de las Heras, Provisions and design support in unions panel-gallery on pit-horizontal layers machined. Madrid, Geo-Miner Technologic Institute of Spain, Ministry of Industry, Trade and Tourism, 1988, (Spanish).
  5. R. Lindsay Galloway, A History of Coal Mining in Great Britain, Augustus M. Kelley, New York, NY, USA, 1969.
  6. Energy Information Administration, Office of Coal, Nuclear, Electric and Alternate Fuels, Longwall Mining, U. S. Department of Energy, Washington, DC, USA, 1995.
  7. Medium-term coal market report 2012, OECD Publishing. IEA, 2012.
  8. Coal mining. World coal. World Coal Institute, 2009.
  9. P. Darling, SME Mining Engineering Handbook, Society for Mining, Metallurgy, and Exploration (SME), 3rd edition, 2011.
  10. A. A. Waddington and D. R. Kay, “Recent development in the incremental profile method of predicting subsidence, tilt and strain over a series of longwall Panels,” in Proccedings of the International Conference on Geomechanics/Ground Control in Mining and Underground Construction, Wollongong, Australia, July 1998. View at: Google Scholar
  11. A. García-Bello, J. Ordieres-Meré, and C. González-Nicieza, “Generalization of the influence function method in mining subsidence,” International Journal of Surface Mining, Reclamation and Environment, vol. 10, no. 4, pp. 195–202, 1996. View at: Publisher Site | Google Scholar
  12. G. W. Mitchell, “Longwall mining,” in Australasian Coal Mining Practice, R. J. Kininmonth and E. Y. Baafi, Eds., pp. 340–373, The Australasian Institute of Mining and Metallurgy, Melbourne, Australia, 2010. View at: Google Scholar
  13. M. S. Venkata Ramayya and L. Sudhakar, “Selection of powered roof support for weak coal roof,” Journal of Mines, Metals and Fuels, vol. 50, no. 4, pp. 114–117, 2002. View at: Google Scholar
  14. R. Singh and T. N. Singh, “Investigation into the behaviour of a support system and roof strata during sublevel caving of a thick coal seam,” Geotechnical and Geological Engineering, vol. 17, no. 1, pp. 21–35, 1999. View at: Publisher Site | Google Scholar
  15. J. Ju and J. Xu, “Structural characteristics of key strata and strata behaviour of a fully mechanized longwall face with 7.0 m height chocks,” International Journal of Rock Mechanics and Mining Sciences, vol. 58, pp. 46–54, 2013. View at: Publisher Site | Google Scholar
  16. N. Bilim and I. Özkan, “Determination of the effect of roof pressure on coal hardness and excavation productivity: an example from a Çayırhan lignite mine, Ankara, Central Turkey,” International Journal of Coal Geology, vol. 75, no. 2, pp. 113–118, 2008. View at: Publisher Site | Google Scholar
  17. C. González-Nicieza, A. Menéndez-Díaz, A. E. Álvarez-Vigil, and M. I. Álvarez-Fernández, “Analysis of support by hydraulic props in a longwall working,” International Journal of Coal Geology, vol. 74, no. 1, pp. 67–92, 2008. View at: Publisher Site | Google Scholar
  18. R. Juárez-Ferreras, C. González-Nicieza, A. Menéndez-Díaz, A. E. Álvarez-Vigil, and M. I. Álvarez-Fernández, “Forensic analysis of hydraulic props in longwall workings,” Engineering Failure Analysis, vol. 16, no. 7, pp. 2357–2370, 2009. View at: Publisher Site | Google Scholar
  19. R. Juárez-Ferreras, C. González-Nicieza, A. Menéndez-Díaz, A. E. Álvarez-Vigil, and M. I. Álvarez-Fernández, “Measurement and analysis of the roof pressure on hydraulic props in longwall,” International Journal of Coal Geology, vol. 75, no. 1, pp. 49–62, 2008. View at: Publisher Site | Google Scholar
  20. S. P. Timoshenko and J. N. Goodier, Theory of Elasticity, McGraw- Hill, New York, NY, USA, 1970.
  21. V. I. Klishin and T. M. Tarasik, “Stand tests of hydraulic props in dynamic loads,” Journal of Mining Science, vol. 37, no. 1, pp. 77–84, 2001. View at: Publisher Site | Google Scholar

Copyright © 2013 María-Belén Prendes-Gero et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.