SciELO - Scientific Electronic Library Online

vol.111 issue12'Educating the Future Mining Engineering Practitioner': Tuks Mining Engineering: 1961-2011 'Celebrating 50 years of Excellence in Mining Engineering Education'Numerical computation of average pillar stress and implications for pillar design author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Journal of the Southern African Institute of Mining and Metallurgy

On-line version ISSN 2411-9717
Print version ISSN 0038-223X

J. S. Afr. Inst. Min. Metall. vol.111 n.12 Johannesburg Dec. 2011




The design of stable pillars in the Bushveld Complex mines: a problem solved?



D.F. MalanI; J.A.L. NapierII, I

IDepartment of Mining Engineering, University of Pretoria, Pretoria, South Africa
IICSIR, Pretoria and Department of Mining Engineering, University of Pretoria, South Africa




This paper gives an overview of the difficulties associated with determining the strength of hard-rock pillars. Although a number of pillar design tools are available, pillar collapses still occur. Recent examples of large-scale pillar collapses in South Africa suggest that these were caused by weak partings that traversed the pillars. Currently two different methods are used to determine the strength of pillars, namely, empirical equations derived from back analyses of failed and stable cases, and numerical modelling tools using appropriate failure criteria. The paper illustrates that both techniques have their limitations and additional work is required to obtain a better understanding of pillar strength.
Empirical methods based on observations of pillar behaviour in a given geotechnical setting are popular and easy to use, but care should be exercised that the results are not inappropriately extrapolated beyond the environment in which they are established. An example is the Hedley and Grant formula (derived for the Canadian uranium mines), which has been used for many years in the South African platinum and chrome mines (albeit with some adaptation of the K-value). Very few collapses have been reported in South Africa for layouts designed using this formula, suggesting that in some cases it might yield estimates of pillar strength that are too conservative.
As an alternative, some engineers strongly advocate the use of numerical techniques to determine pillar strength. A close examination unfortunately reveals that these techniques also rely on many assumptions. An area where numerical modelling is invaluable, however, is in determining pillar stresses accurately and for studying specific pillar failure mechanisms, such as the influence of weak partings on pillar strength.
In conclusion, it appears that neither empirical techniques nor numerical modelling can be used solely to provide a solid basis for conducting pillar design. It is therefore recommended that both these techniques should be utilized to obtain the best possible insight into a given design problem. Owing to the uncertainties regarding pillar strength and loading stiffness, monitoring in trial mining sections and in established mining areas is also an essential tool to test the stability of pillar layouts in particular geotechnical areas.

Keywords: pillar design, bord and pillar mining, stable pillars.




Appropriate pillar design is a fundamental building block of mine design to ensure the safe and economic extraction of valuable national resources. It is therefore worthwhile to take a critical look at the tools currently available to rock engineers to conduct these designs. This paper focuses only on stable pillars and does not address the issues associated with crush/yield pillars. The design of stable pillars is currently very topical in the shallow hard-rock mining sector in South Africa, and this paper will highlight a number of examples from this country. Hard-rock pillar design nevertheless appears to be of universal interest, as shown by the examples described below.

Zipf1 describes collapses of room and pillar mines in the USA. The term catastrophic pillar failure or CPF is used to describe the mechanism whereby a few pillars fail initially; their load is then transferred to adjacent pillars, which also fail. This may result in a 'pillar run' and hundreds of pillars may fail in the process. A number of examples of CPF collapses in 'metal' mines are given in Zipf's paper, and apparently at least four such examples have occurred in the USA since 1972. One of the more recent examples is a large pillar collapse in a room and pillar base metal mine2. Figure 1 illustrates the area of the collapse. The failure began in four centrally located pillars and spread rapidly to include almost 100 pillars. The pillar width was 8.5 m and the room width was 9.7 m. The pillar height was about 12 m, resulting in a widthheight ratio of 0.70. The extraction ratio was approximately 78%. Based on the examples available to him, Zipf1 made the comment that mines experiencing CPF generally exhibit the following characteristics (quoted directly from his paper):



'(1) Extraction ratios are usually more than 60%. A high extraction ratio will put pillar stress close to peak strength and provide ample expansion room for the failed pillar material.

(2) Width-height ratios of pillars are always less than 3 for coal mine failures, usually much less than 1 in metal mine failures, and less than about 2 for nonmetal mine failures. A low width-height ratio ensures that the failed pillar material can easily expand into the surrounding openings and that the failed pillar will have little residual load-bearing capacity.

(3) The number of pillars across the panel width is always at least five and usually more than 10, which typically ensures that pillars have reached their full tributary area load. Minimum panel widths for CPF are at least 80 m.

(4) Substantial barrier pillars with width-height ratios more than 10 are absent from the mine layout.

(5) Although CPF seems more prevalent in shallow mines less than 100 m deep, this may be only a reflection of the prevalence of shallow room-and-pillar coal mines.'

An extensive database of hard rock pillar failures was compiled by Lunder and Pakalnis3, who analysed 178 case histories from hard-rock mines of which 98 were located in the Canadian Shield. Of the pillars investigated, 76 were classed as stable; 62 were classed as failed; and 40 were classed as unstable. Many of these pillars were rib or sill pillars from steeply dipping ore bodies. Lunder and Pakalnis proposed that the pillar strength can be adequately expressed by two factor-of-safety (FOS) lines. Pillars with a FOS < 1 fail, while those with a FOS > 1.4 are stable. The region 1 < FOS < 1.4 is referred to as unstable and pillars in this region are prone to spalling and slabbing, but have not completely failed. The data collected by these authors is given later in this paper.

Esterhuizen4 conducted an evaluation of the strength of slender pillars. Figure 2 illustrates published case histories of failed pillars from hard-rock metal mines. This figure illustrates that the pillar strength becomes highly variable as the width to height ratio decreases. For a ratio of 1, the pillar strength is not expected to exceed 65% of the laboratory strength (but it can also be significantly lower than this value). Esterhuizen4 suggests that a larger factor of safety is required for these pillars to account for the increased variability in strength. Numerical modelling indicated that for slender pillars, the difference between the pillar load at the onset of spalling and the ultimate pillar strength can be small. This implies that slender pillars are at or near the point of failure at the onset of brittle spalling.



Although a large number of additional papers on pillar design can be found in the literature, the objective of this paper is not to give an overview of all of these, but rather to highlight the difficulties associated with pillar design in South Africa and to emphasize the need for additional research. Some examples of pillar failures in South Africa are therefore described below, and the commonly used pillar design methodologies are critically examined.

The Coalbrook disaster and the evolution of pillar design in South Africa

Although this paper focuses on pillar design in hard rock mines, it is worthwhile to investigate the evolution of pillar design in South Africa in general. This will shed some light on the apparent lack of appropriate pillar design tools for shallow hard rock mines.

Initially, pillar dimensions and mining spans in South African coal mines were based on experience obtained through trial and error. This resulted in a number of collapses. The first report of a coal pillar collapse in South Africa was in 1904 at the Witbank colliery5. Since that date, 81 pillar collapses have been recorded in 31 collieries in the former Transvaal and the Free State. Between 1904 and 1965, there were 50 pillar collapses6. The research into coal pillar strength in South Africa gained momentum after the Coalbrook disaster. As described by Van der Merwe7, the major need for coal mine research in South Africa was identified when the Coalbrook disaster occurred on 21 January 1960. In total, 437 workers lost their lives when pillars over an area of 324 hectares collapsed (Figure 3). Following the disaster, the South African government sponsored research into coal mine safety by forming the Coal Mines Research Controlling Council (CMRCC). Research into coal pillar strength received top priority from the Council for Scientific and Industrial Research (CSIR) and the Chamber of Mines Research Organization (COMRO). A large in situ testing programme to determine pillar strength was conducted by the CSIR8 and tests were also conducted by the Chamber of Mines9.

A key study on coal pillar strength was conducted by Salamon and Munro10 and the findings as described in their paper are still being used today. It was postulated that the strength of the pillars can be expressed, in the given range of dimensions, as a power function of the height and width. Of interest is that this 'power law' was originally proposed by Greenwald, Howard, and Hartmann11. Salamon and Munro10 derived their equation by estimating the value of the constant K and the powers of width and height by the method of maximum likelihood. The derivation was empirical and was based on data from stable and collapsed cases. The analysed data included 125 cases of which 27 were collapsed pillars. It is noteworthy that they warned in the paper against the following mistake commonly made in rock engineering:

'The work described in this paper is essentially empirical, and the results, therefore, should not be extrapolated beyond the range of the data which were used to derive them.'

Although the formula is well known, it is repeated below for completeness:

where K reflects the fitted 'strength' of a metre cube of coal (7.2 MPa), w is the width of the (square) pillar, and h is the height in metres. The parameters a and b are equal to 0.46 and 0.66 respectively. The pillar volume is given by V=w2h. Defining the width:height ratio, R = w/h, Equation [1] can be expressed in the alternative form12:

It can be seen from Equation [2] that if a = b, the pillar strength is independent of the pillar volume whereas if b > a, as in Equation [1], the pillar strength is predicted to decrease as the pillar size is increased even if the pillar shape is unchanged.

Did the use of Equation [1] improve pillar designs? Wagner and Madden13 reported that since 1967, approximately 1100 million tons of coal had been mined in South Africa and they estimated that 1.2 million pillars were left underground during that time. During the same period, 13 cases of pillar collapse were reported, involving a total of about 4000 pillars. This corresponded to a probability of failure of only 0.003. In spite of the satisfactory performance of the pillar design procedure, three areas requiring further research were identified. These were:

Regional differences in coal seam strength

Effect of mining method on pillar strength

Strength of squat coal pillars.

Regarding the strength of the squat pillars, the original database used by Salamon contained no pillar with a width to height ratio greater than 3.8. Evidence collected in the field suggested that beyond a critical width to height ratio, the pillar strength exceeds that suggested by Equation [1]. Salamon14 proposed therefore that when the width to height ratio exceeds a critical ratio, the pillar strength formula should be replaced by the following:


K = the strength of a unit cube of coal

V = the pillar volume (m3)

R = the pillar width to height ratio

Ro = the critical width to height ratio

ε = rate of pillar strength increase.

From field data, no evidence was available of a collapse of a pillar with a width to height ratio greater than 4. Therefore the critical width to height ratio was selected as 5. A value of 2.5 was chosen for ε as it was considerably lower than that obtained from laboratory tests on sandstone.

Regardless of the apparent success of the coal pillar power strength formula, some criticism can be raised regarding its applicability. Virtually all laboratory and field evidence indicates that the width-height strengthening curve has a zero or positively upwards curvature15. The power formula forces downward curvature and this leads to the inelegant form of the squat pillar formula. An objection raised by Bieniawski16 is that according to the power law formulation, the cube strength (w =h) would continue to decrease indefinitely with side length. This is considered unreasonable17.

Based on these criticisms, an alternative 'linear' equation, with no volumetric size effect, was proposed which directly expresses the strengthening effect of the width-height (w/h) ratio16.

As noted above, Equation [1] implies that pillars with self-similar dimensions (the same ratio of w to h) will have different strengths and that the predicted strength decreases as the pillar volume is increased. By contrast, Equation [4] has no size effect in that the strength depends only on the ratio of w to h. York and Canbulat18 compared the relative goodness-of-fit power formula versus linear forms for both coal and hard-rock materials, and concluded that the latter behaved at least as well.

Further work regarding coal pillar design was conducted by other workers (e.g. Van der Merwe7), but it is beyond the scope of this paper and the reader is referred to this reference for additional information. In conclusion, Ryder and Jager15 state:

'The power law and its derivatives are perhaps too entrenched in coal engineering to warrant withdrawing from them at this time, but in hard rock engineering, the simpler and probably more realistic linear forms are advocated for general use.'

Regarding the gold mining industry in South Africa, the pillar problem was a completely different issue owing to the great depth of these mines and the associated seismicity. Leaving small pillars or remnants in stopes is in fact detrimental to stability as these may become a source of seismicity. This was a key motivation for the early adoption of the longwall mining method as it minimized the formation of remnants19. These longwalls remained vulnerable to seismicity and one method to reduce this risk was to leave strips of ground behind as strike stabilizing pillars. A key function of these pillars is to reduce the levels of energy release rate15. Owing to the size of the gold mining industry and the role it played in the South African economy, significant research into the behaviour of these pillars was conducted over the years. Their widths typically vary from 30 m to 40 m and at a typical stoping width of 1.5 m, this gives a width to height ratio of at least 20. These pillars are therefore considered 'indestructible'. Seismic problems are nevertheless experienced when these pillars 'punch' into the footwall, and a simple empirical formula is used during design to ensure that that the average stress on these pillars does not exceed a specified multiple of the strength of the rock in the footwall or hangingwall. So-called bracket pillars are also used to clamp seismically active geological structures in the gold mines. As the nature and function of the bracket and stabilizing pillars are different to those found in shallow hard-rock mines, these topics are not discussed in this paper and the reader is referred to Ryder and Jager15 for additional information.

In comparison to the gold and coal mining industries in South Africa, the platinum and chrome mines were the proverbial 'stepchildren' and have never received the same attention in terms of research efforts and funding. Fortunately, the mines in the Bushveld Complex have never had catastrophic pillar failures analogous to the Coalbrook disaster or rockbursts to kick-start research, and it was only with the rapidly increasing price of platinum from 1999 onwards that research activity started to receive more attention. In terms of rock engineering knowledge, much was therefore 'borrowed' from the other mining sectors and overseas research findings. A contributing factor may also have been that pillar cutting in the hard-rock industry is not as precisely controlled as in the coal industry, and it is therefore rather difficult to repeat the statistical approach followed by Salamon and Munro10. An example is shown in Figure 4, which illustrates a mapped survey of typical pillar cutting in one of the Bushveld mines, where the pillars can be seen to have a wide variety of shapes. Gay et al.20 stated that at that time:

'The design of pillar layouts in shallow to medium depth chromium and platinum mines has not reached the same advanced stage as has the design of pillars in coal mines.'

Three reasons for this were given in the paper namely:

Very little is known about the strength of small pillars composed of chromitite or pyroxenitic platinum ore

Because of the very competent and stiff strata in the hangingwall, it is difficult to determine the pillar loads accurately

The presence of and shear along near vertical faults can change the loading condition from a stiff displacement controlled system to a soft load controlled system.

A large collapse involving pillars at the No 6 Shaft, Bafokeng North mine is described by Kotze21. The collapse occurred in February 1975. Luckily, nobody was injured during this collapse as it occurred overnight. This collapse is referred to as the 'Hospital collapse' by mine personnel as the mine's hospital was located on surface above the collapsed area and some damage was sustained by the hospital building owing to the resulting subsidence. The building was repaired and it is still being used today. Extraction of the reef in this area was by room and pillar methods, with rooms approximately 25 m wide and pillars nominally 5 m square, providing an extraction of 96%. It appears as if the pillars were cut smaller than this, as Kotze21 estimated from the mine plans that the equivalent 'average' pillar width in this area was 3.6 m with an associated extraction of about 98%. Unfortunately no mention is made in the paper of the pillar height. Local hangingwall support was provided by timber sticks and mat packs. The strike span of the stope at the time of the collapse was approximately 400 m, and the average depth of the overburden was 58 m. The extent of the collapse was of the order of 350 m on strike and 350 m on dip, and the outline of the collapsed region is shown in Figure 5. The upper and lower limits of the fall were bounded by two faults. Up to the time of the collapse there were no signs of pillar scaling and no timber poles were broken, suggesting that conditions were indeed stable. This condition, however, changed overnight with the entire rock mass to surface moving down between the two faults shown in Figure 5. No doubt this observation resulted in the hypothesis given by Gay et al.20 that the presence of and shear along near vertical faults can change the loading condition from a stiff displacement controlled system to a soft load controlled system. Kotze21 used numerical modelling and stress measurements conducted by the CSIR in 1975 in pillars adjacent to the collapse area to estimate a Merensky pillar strength of 75 MPa (for a pillar with an effective width of 4.3 m). This value is of course applicable only to the particular width to height ratio of the pillars (the pillar height was given by Kotze21 as 1 m).



Empirical methods to estimate hard rock pillar strength

Martin and Maybee22 give a very good overview of the different empirical strength formulae that were developed to predict pillar strength. A comparison of these formulae is given in Figure 6. These curves were calculated for a pillar height of 5 m.



In South Africa, based on the success of the Salamon and Munro strength formula in the coal mines, it was natural for the shallow hard-rock mines to adopt a similar power-law strength formula. The necessary research to develop and calibrate a formula for local conditions was, unfortunately, never conducted. Instead, the Hedley and Grant formula23 developed for the Canadian uranium mines was adopted. Only the K-value was modified to reflect local rock strengths. This approach seemed to work well and over the years it has become firmly entrenched. It would not be incorrect to state that it is currently the 'industry accepted' method for designing pillars in shallow hard-rock mines in South Africa. Ryder and Jager15 comment that this original formulation: 'for want of local data, have subsequently been applied in South African hard-rock mines'. This is a good example of what frequently happens in rock engineering. Over time, initial assumptions and interim solutions become entrenched as common practice and the original assumptions are rarely revisited or questioned. Based on this wide acceptance of the Hedley and Grant formulation and lack of local research into more appropriate formulations, it will be worthwhile to examine the assumptions made in the original paper.

The uranium mines in the Elliot Lake area used a stope and pillar layout to mine the orebody. Narrow pillars about 250 ft (76 m) long were left along dip (see Figure 7). The pillars which were chosen for analysis were typically conducted 10-20 ft wide (3-6 m) and 8-20 ft high (2.5- 6 m). The width to height ratio of most of the pillars was close to 1 and only a very few (3 in the database) had a width to height ratio of 2.5. It is clear that this original formulation was derived for slender rib pillars and it can be questioned whether it is applicable to square pillars with a width to height ratio greater than 2.5. The original data used by Hedley and Grant is reproduced in Table I. It is immediately obvious that the dataset used was very small (28 pillars). This should be compared to the coal database of Salamon and Munro10, which included 125 pillars of which 27 were collapsed. The width to height ratio of the failed pillars varied from 1.1 to 1.5. Only 3 of these pillars were 'crushed' and 2 were 'partially crushed'. Of further concern is that it is stated in the paper that: 'The information on complete pillar crushing was obtained second-hand because it happened in mines which are closed.' This work was conducted in the days before computer-based numerical modelling could be used to determine pillar stress. The approach followed was therefore to use tributary area theory, which relates the pillar stress to the pre-mining stress and the extraction ratio by:


σP = pillar stress

So = pre-mining stress normal to the orebody

e= extraction ratio.



For workings inclined at an angle α to the horizontal, the normal stress So is a combination of the vertical stress component Sv and the horizontal stress component Sh:

The vertical stress was assumed to be simply a function of the weight of the overlying strata and the horizontal stress perpendicular to strike was assumed to be 3000 pound per square inch based on measurements in two mines. The stress given in the Table I for each pillar is therefore only a rough estimate. The methodology followed in the paper to derive the pillar strength formula can be summarized as follows:

The first step was to adopt the power law strength formulation used by Salamon and Munro10. In the notation of Hedley and Grant23 it is given as:


Qu = pillar strength (psi)

w = pillar width (ft)

h = pillar height (ft)

K = strength of 1 ft cube (psi)

a and b are constants

Hedley and Grant23 acknowledge that this equation refers to square pillars, whereas those in the uranium mines are long and narrow. Their assumption was therefore that the strength of the slender pillars will not be much greater than a square pillar of width equalling the minimum width of the long pillar.

Secondly, from extrapolation of laboratory tests, it was estimated that the value of K is 26 000 psi for a 1-ft cube.

Thirdly, appropriate values for parameters a and b had to be derived. Three different sets of values were available to them at that stage in the literature. The value for a was relatively constant at 0.5 and Hedley and Grant23 therefore decided to also adopt this value. As b varied more, a new value was computed and their approach was to focus on the three failed pillars in the database. For each of these pillars, the tributary area stress in the table was assumed to be the pillar strength. This value, as well as the K-value and a = 0.5, were substituted into Equation [7] and the value of b was solved for each pillar. The calculated values of b ranged from 0.736 to 0.768 with a mean of 0.75. This value was adopted and it resulted in the now familiar Hedley and Grant formulation.

Clearly the formulation above is based on a large number of assumptions, and the applicability of this formulation to the design of hard-rock pillars in the Bushveld Complex in South Africa becomes highly questionable. The first use of this formula in a South African mine is not clear, but Ozbay et al.25 stated that it was 'popularized by Wagner and Salamon26 as quoted by Kersten27'. Kersten27 used it to design pillars for Agnes Gold Mine. Subsequently, it has been used to design a large number of bord and pillar layouts in the country with an appropriate modification of the value of K. The rule frequently used in South Africa is to estimate K at between one-third and two-thirds of the Uniaxial Compressive Strength (UCS) of the pillar material (e.g. see Ozbay et al.25). Almost no collapses in the Bushveld Complex have been reported to date using this formulation (except where weak clay layers are present in the pillar, see sections below), so the uncomfortable question therefore remains: Why does it work and are the current designs perhaps too conservative? One hypothesis is that the 'squat' behaviour of hard-rock pillars may occur significantly earlier than the w:h ratio of 5 assumed for coal pillars15.

As an alternative to the estimate of K being one-third of the UCS of the pillar material, Stacey and Page28 and Stacey and Swart29 suggested that the value of K be chosen as the design rock mass strength (DRMS). The DRMS takes into account the rock quality, unfavourable joint orientations, and the excavation method30. Wesseloo and Swart31 conducted a study on LG6 Chromitite pillars to determine the constants K, a and b (see Equation [7]) more accurately for these pillar types. The study showed that the value of the exponents a and b suggested by Hedley and Grant23 do not apply to these chromitite pillars. They indicated in the study that values of a = 1 and b =1 may be more applicable to these pillars.

As a further attempt to develop a new formula for the South African mines, Watson et al.32 derived new values for the power law formulation given in Equation [7]. They used a maximum likelihood evaluation similar to that used by Salamon and Munro23. Their database consisted of 179 Merensky Reef pillars of which 109 were stable. The width to height ratio of the pillars in the database ranged from 1 to 8, with the majority between 3 and 6. Only one pillar had a width to height ratio of less than 1. The values derived are K = 86 MPa, a = 0.76, and b = 0.36. It is interesting to note that the b parameter is much lower than for the Hedley and Grant formula. The formula predicts pillar strengths that are much greater than the traditional Hedley and Grant formulation, with K values assumed to be a third of the UCS (see Figure 8). Unfortunately, it is not known if this formula has been tested in any trial mining sections.



Corrections for rectangular pillars

Equation [7] is applicable to square pillars. Another awkward assumption commonly made when designing pillars is to calculate the 'effective width' for rectangular shaped pillars. Holland and Gaddy33 state that only the minimum lateral dimension, w, affects the strength of a pillar, while the other dimension, L, has no effect. Wagner9 suggested that the effective width of a pillar should take the form:

where A is the cross-sectional area of the pillar and C is the perimeter. This formula is based on observations made by Wagner during large scale underground tests, namely that: 'the strength of the circumferential portions of a pillar is virtually independent of the width-to-height ratio, whereas the strength of the centre increases with increasing ratio.' The effective width as defined by Equation [8] approaches a finite value of twice the minimum pillar width for very long and narrow pillars. It is not clear if this assumption is correct and Ryder and Ozbay34 suggested an alternative shape strengthening factor of the form f= 1.0/1.1/1.2/1.3 for pillars having w 1/ w 2 ratios of 1/2/4/∞. This implies that Wagner's perimeter rule may be overestimating the strengthening effect of very long pillars.

Unfortunately no good experimental evidence is available regarding the effect of pillar shape on strength for hard-rock pillars in South Africa and this area requires further research. Somewhat concerning is that the whole empirical design philosophy rests on a huge number of unproven assumptions, and pillar strength is clearly an area that requires systematic research in future.


Recent examples of pillar failures

To illustrate the inherent dangers of using empirical design formulae for rock masses in environments for which they were not originally designed, consider three case studies of recent mine collapses in the South African Bushveld Complex.

Spencer35 reported on the failure of the pillar system and the subsequent closure of the Wonderkop Chrome Mine in May 1998 (stoping commenced only in September 1996). The mine is located close to Rustenburg and exploited the Lower Group chromitite seams, namely the LG6 and LG6A. An internal pyroxenite waste band is found between these two reefs and this results in the pillars having a multi-layered appearance (Figure 9) wherever these two reefs are mined. The Wonderkop Mine was the most easterly situated LG6 mining operation in the Rustenburg area and was situated adjacent to the Spruitfontein dome (an upfold structure which separates the Rustenburg section from the Marikana section). Its close proximity to the Spruitfontein dome has influenced the structure of the LG6 and LG6A seams and this has resulted in thick clay layers (up to 300 mm in some places) traversing the pillars in some areas (Figure 10). The position and thickness of this weak layer is highly variable (e.g. see Figure 11).







The original pillar design at the mine was conducted using the Hedley and Grant pillar strength formula. The pillar sizes were 12 m x 6 m giving an 'effective width' of 8 m according to Wagner's perimeter rule. K was assumed to be 27.3 MPa, which was a third of the laboratory strength of the rock (using samples obtained from another mining section). Owing to the complex multilayer structure of the pillars, it is not clear which of these layers were tested. The stoping width was 2 m, so the width to height ratio was at least 3 if the smallest dimension of the pillars is considered. When using the 'effective width' of the pillars as 8 m, the strength of the pillars was estimated to be 45.9 MPa.

The first underground inspection of the pillars by the consulting rock engineer was conducted in July 1997. During this visit it was noted that some joints were beginning to open at the corners of some pillars. Some joints were also opening up at the sides of the pillars and in a few cases sliding along the clay layer was noted. Figure 12 illustrates the layout and the positions where pillar failure was observed. Following these observations, steps were taken to introduce a system of barrier pillars (Figure 13) with a width to height ratio of at least 10 (it is generally believed that at this ratio the pillar becomes indestructible). During the following nine months, the condition of the pillars continued to deteriorate. To strengthen the pillars along the main dip belt and road declines, two strategies were adopted; namely waste stowing between the pillars and mesh and lacing of the pillars. The success of the mesh and lacing of the pillars appeared to be doubtful, as the drilling process introduced additional water into the clay which probably weakened the pillars further. During April 1998, the failure process accelerated and the rate of closure in some areas increased to 1.8 mm/day. Numerous falls of ground occurred and management decided to cease operations at this stage. Recently, a back analysis of this pillar failure was conducted by the authors using the TEXAN boundary element program36, which can explicitly simulate the individual pillars and calculate the stresses on these pillars relatively accurately. The results are shown in Figure 14. Two face positions were simulated; namely July 1997 (Figure 12) and July 1998 (Figure 13). From this study, the Hedley and Grant formulation was used to back-calculate the K-value for the pillars. A value as low as 6 MPa was obtained, which is in agreement with earlier back analysis studies by Spencer35, who estimated a value as low as 4.6 MPa. It should nevertheless be noted that a power law strength formula might not be applicable to these pillars owing to the presence of the clay layer. Although the mechanism of pillar failure needs further investigation, the low strength of the pillars can probably be attributed to the low friction angle of the clay layer, which allowed lateral deformation of the pillar and a reduction of confinement.

Two other large scale pillar collapses recently occurred in the Bushveld Complex. Detailed information regarding these collapses is not available in the public domain, and therefore these mines will only be referred to as Mine A and Mine B.

Mine A is a platinum mine located in the eastern portion of the Bushveld Complex. At this mine, a clay layer is also present at the hangingwall/pillar contact (Figure 15). The reef exploited in this area is the UG2. The original mine design was conducted using the Hedley and Grant pillar strength formulation with a K value of 35 MPa. The mining height was 2 m. In mid 2008, some concern was expressed regarding the stability of the pillars and a minor collapse occurred during this time. In an attempt to reinforce some of the pillars, many were supported using fibre-reinforced shotcrete. This did not stop the deterioration, however, as shown in Figure 16 with the cracked shotcrete clearly visible. During December 2008, operations were temporarily suspended at the mine when the decline was affected by the instability.





Similar to Wonderkop Mine, Mine B also exploits the LG6/LG6A chromitite reefs. The problem is essentially similar to that experienced at Wonderkop Mine with the presence of a clay layer in some of the pillars. This resulted in collapses in parts of the mine. Typical pillar failure at the mine is shown in Figure 17. Experience has indicated that increasing the pillar sizes does not necessarily work in these cases.



From these studies, the drawbacks of using empirical pillar strength formulae are obvious. The failure in all three cases is caused by the presence of clay layers which substantially weaken the pillars. The original empirical formulae were developed for different rock types and the application of these formulae outside the limits for which they were developed resulted in the large mine collapses described here.


Numerical modelling to estimate pillar strength

From the sections above, it is clear that the applicability of the current empirical formulas in hard-rock mines is highly uncertain and additional verification and calibration work is required. The orebodies in the Bushveld Complex are very different to that for which the Hedley and Grant formulation23 was derived. Furthermore, the restricted range of slender width to height ratios used when deriving the equation is unfortunate and the choice of the appropriate K value is undefined and highly uncertain.

The alternative to the empirical approach is to use numerical modelling with appropriate failure criteria to determine pillar strength. A vast amount of literature is available on attempts to simulate pillar failure, and it is not the objective of this paper to summarize all these findings. Focus will rather be placed on recent work that is applicable to the pillars in the Bushveld Complex.

Joughin et.al38 presented a risk based approach to the design of chromitite pillars. Numerical modelling was conducted using the finite element program PHASE2 and a Hoek and Brown failure criterion. The variability of the strength properties, pillar dimensions, and spans were taken into account by applying the Point Estimate Method. Many numerical analyses were conducted with varying input parameters to obtain statistical distributions of pillar safety factors. Probabilities of failure were then determined from the distributions. Back analysis of an area with a large collapse and another in which only a few pillars had failed indicated good agreement with the actual number of pillar failures.

Day and Godden39 presented a paper describing the design of panel pillars on Lonmin's platinum mines. They state that the validity of their approach is supported by extensive underground surveys and by computer backanalysis studies. According to the authors, over 300 pillars per month were cut at Lonmin at that stage and these pillars behaved as expected. The method apparently works well up to width to height ratios of 5.5, but not at greater values owing to the onset of squat pillar effects. This seems rather disappointing as the expectation is that a numerical method with an appropriate constitutive model will 'automatically' take care of the onset of squat pillar behaviour. Pillar strength was estimated by two-dimensional FLAC modelling using the original Hoek and Brown40 failure criterion:

where σc = uniaxial compressive strength (MPa) and m and s are constants that depend on the properties of the rock. The constant m was determined from laboratory testing, and s equals to 1 for intact specimens. In situ values for the constants m and s were derived by application of rock mass quality ratings using the equations of Priest and Brown41 for undisturbed rock masses:

The authors39 derived in situ values for m and s for the UG2 and Merensky reefs at Lonmin's Mines. Typical values used in the modelling are as follows: UG2; m = 25.83, s = 0.51 for a RMR of 94, Merensky (Type B); m = 8.7, s = 0.57. The resulting simulated pillar strengths seem plausible when the pillar width to height ratio is low. The authors nevertheless acknowledge that squat pillar behaviour is not correctly simulated by this approach and that this will probably result in pillars being over-designed at depths exceeding 700 m.

A further concern regarding this approach is that the failure criterion is still empirical and the unmodified Equations [10] and [11] may not be appropriate for the pillar material in the Bushveld Complex. Pells42 also expressed some concern about the Hoek and Brown failure criteria and quoted Mostyn and Douglas43 which provided a detailed critique of this failure criterion for intact rock.

'...there are inadequacies in the Hoek-Brown empirical failure criterion as currently proposed for intact rock and, by inference, as extended to rock mass strength. The parameter mi can be misleading, as mi does not appear to be related to rock type. The Hoek-Brown criterion can be generalized by allowing the exponent to vary. This change results in a better model of the experimental data.'

Martin and Maybee22 investigated the strength of hard rock pillars by using both empirical pillar strength formulae and numerical modelling using a Hoek-Brown failure envelope. Figure 18 illustrates a comparison between the pillar stability graph developed by Lunder and Pakalnis3, the Hedley and Grant equation and Phase 2 - Hoek and Brown modelling using different values of GSI. Figure 19 illustrates the same data using the Hoek and Brown brittle parameters. The conclusion reached by Martin and Maybee22 is that twodimensional finite element analyses using conventional Hoek and Brown parameters for typical hard-rock pillars predicted rib pillar failure envelopes that did not agree with empirical pillar failure envelopes. It appears that the conventional Hoek-Brown failure envelopes over-predict the strength of hard rock pillars. This occurs because the failure process is fundamentally controlled by a cohesion loss process in which the frictional strength component is not mobilized. Their recommendation is that Hoek-Brown brittle parameters (mb = 0 and s = 0.11) be used to simulate pillar strength.





Hoek, Kaiser, and Bawden44 summarized the rock mass conditions for which the Hoek-Brown failure criterion can be applied. The criterion is strictly applicable only to intact rock or heavily jointed rock masses that can be considered to be homogenous and isotropic. For cases in which rock mass behaviour is controlled by a single discontinuity or joint set, a criterion that describes the shear strength of joints should rather be used (e.g. see the pillar shown in Figure 20). The implication of this is that for the three case studies of pillar failure in the Bushveld Complex discussed above, explicit simulation of the clay layer will be required. An example was conducted for the authors by Dr John Ryder using the FLAC code, and this is presented below to illustrate the value of modelling.



The pillar composition simulated was the LG6/LG6A 'sandwich' shown in Figure 10. The qualitative effect of a strong pyroxenite layer within a chromitite pillar (with weak contacts, including weak hangingwall and footwall contacts) was therefore modelled in plane strain. A generic model was built to investigate the effect of an inhomogeneous pillar with weak interfaces. Estimated in situ strain-softening parameters were drawn directly from studies carried out in the Bushveld Complex. The hangingwall and footwall were assigned the same properties as the pyroxenite layer, and symmetry was assumed for both the vertical and horizontal centrelines in the following layout in the FLAC finite difference code (Figure 21). The grid size was 0.1 m × 0.1 m. By applying slow displacement loading (velocity 5 × 10-4 mm/step), complete stress-strain and lateral deformation curves could be modelled (Figure 22). Lateral deformations showed no dramatic effects owing to the presence of the 'strong' layer of pyroxenite in the pillar, possibly because the modelled contrasts in strength, Poisson's ratio and dilatancy were not particularly marked. (Note that the horizontal scale in Figure 22b is in millimeters whereas the vertical scale is in metres). Likewise, the presence of a weak interface between the layer and the body of the pillar had virtually no effect, even if the friction angle of interface 1 was set as low as 6º. In contrast, low friction angles on the hangingwall contact (interface 2) had a powerful effect, reducing the peak strength p of the pillar by allowing lateral deformation and reduction of confinement, and reducing also the residual strength (Figure 23).





A further illustration of the value of numerical modelling to investigate pillar strength is given by Esterhuizen45. The UDEC code was used to simulate pillars of different width to height ratios and the effect of jointing on these pillars. The most significant result obtained was that the effect of jointing on pillar strength decreases as the width to height ratio of the pillar increases. This indicated that when designing pillars, it would not be correct to simply determine the rock mass strength and use it in the pillar formulae regardless of the width to height ratio of these pillars. Further work was conducted by Esterhuizen et al.46 to investigate on the effect of large discontinuities on pillar strength. UDEC simulations were used to investigate the effect of these large discontinuities on pillars with width to height ratios of 0.5 to 1.5. A revised empirical pillar strength formula was developed for the slender pillars in the stone mines in the USA. The form of the equation is:

where UCS is the uniaxial compressive strength of the intact rock and LDF is a large discontinuity factor. The authors recommended that a factor of safety of at least 1.8 be used for designs where the pillar strength is assumed to be given by Equation [12]. This will ensure that the pillars remain within the limits of previously successful stone pillar layouts.


The need for monitoring in pillar mining sections

The report of the Government Mining Engineer47, published five years after the Coalbrook disaster, made the following statements regarding instrumentation in the experimental section at the mine:

'There is no evidence that during these operations any systematic inspections were made of the northern portions of No. 10 Section or that any measuring apparatus was installed to ascertain if stability was being affected by the reduction of the abutment support.'

'It appears that no scientific controls were applied to the experiment nor was an attempt made to ascertain what reduction in their support strength was caused by cutting the pillars and increasing their height. The effects of the experiment on the roof and pillars were judged by eye and ear and the results were gauged by the tally board.'

Although monitoring is currently being conducted in some of the Bushveld Complex mines on a limited scale, routine monitoring is probably required in the many bord and pillar sections owing to the uncertainties associated with pillar strength and loads. One example of the type of monitoring that may be valuable is closure monitoring using sensitive instruments. It has been shown that in shallow mining areas, the rate of closure is essentially zero (less than the sensitivity of the instruments) in the stable back areas. Once an area becomes unstable, the rate increases from zero to a very small rate, e.g. 0.3 mm/day. This rate may persist for many days or even weeks48. As this rate is so small, sensitive instruments are required to give early warning of impending instability. In pillar areas, compiling a photographic record of selected pillars over time will also be extremely useful to provide an unbiased opinion whether the pillars are scaling or not (care should be taken that the pillars are appropriately marked and that the successive photographs are taken from the same positions).

As a final comment, in contrast to rock engineering, routine monitoring is commonly used in many of the other engineering disciplines. Civil engineers work with materials that are better characterized and better understood than rock masses, but still conduct routine monitoring of their structures. Structural health monitoring of civil infrastructure systems is a very active field and readers will be able to find a large number of references by searching for this topic on the internet. Many other examples can be quoted from other engineering disciplines, but in general the rule applies that if uncertainty is associated with an engineering design, monitoring should be used to reduce the risk associated with failure of the design.



This paper gives an overview of the difficulties associated with determining the strength of hard-rock pillars. Although a number of pillar design tools are available, pillar collapses still occur. Recent examples of large scale pillar collapses in South Africa were caused by weak partings that traversed the pillars. Currently, two different methods to determine the strength of pillars are used; namely, empirical equations derived from the back analysis of failed and stable cases and numerical modelling using appropriate failure criteria. Both techniques have their limitations and additional work is required to obtain a better understanding of pillar strength.

Empirical methods are popular and easy to use, but care should be exercised that the results are not extrapolated beyond the range of the data which were used to derive them. An example is the Hedley and Grant formula (derived for the Canadian uranium mines) that has been used for many years in the South African platinum and chrome mines (albeit with some adaptation of the K-value) to design pillar layouts in these mines. A careful study of the original publication by Hedley and Grant indicates that this formula was derived for conditions and rock types that are very different to those found in the South African mines. Nevertheless, very few collapses have been reported and in some cases it appears that the Hedley and Grant formula might underestimate pillar strength significantly.

As an alternative, some engineers strongly advocate the use of numerical techniques to determine pillar strength. A close examination unfortunately reveals that these techniques also rely on many assumptions, and extreme care needs to be exercised when using this approach. Pillar load is another unknown not discussed in this paper and care should also be exercised when this is estimated using numerical procedures (some difficulties are outlined in a companion paper in this issue). An area where numerical modelling is invaluable, however, is to study specific pillar failure mechanisms, such as the influence of weak partings on pillar strength. The modelling indicated that a likely failure mechanism in the case studies was that the presence of a parting with a low friction angle that allowed lateral deformation and a reduction of confinement in the pillars. This reduced the peak strength of the pillars considerably.

In conclusion, it appears that neither empirical techniques nor numerical modelling can be used solely to provide a solid basis to predict pillar strength. It is therefore recommended that both these techniques be utilized when addressing pillar design problems in order to obtain the best possible insights. Owing to the uncertainties regarding pillar strength, pillar stress, and loading stiffness, monitoring in trial mining sections (and even in established mining areas) is considered to be an essential tool to assess the stability of pillar layouts in particular geotechnical areas. The need for additional research into pillar strength should also be emphasized strongly as this problem has clearly not yet been solved!



1. ZIPF, R.K. Toward pillar design to prevent collpase in room-and-pillar mines. 108th Annual Exhibit and Meeting, Society for Mining, Metallurgy and Exploration, Denver, 2001.        [ Links ]

2. DISMUKE, S.R., FORSYTH, W.W., and STEWART, S.B.V. The Evolution of Mining Strategy Following the Collapse of the Window Area at the Magmont Mine, Missouri. Proceedings, CIM Annual General Meeting, District 6, Metal Mining, 1994. pp. 3-8.        [ Links ]

3. LUNDER, P.J. and PAKALNIS, R. Determination of of the strength of hard rock mine pillars. Bulletin of the Canadian Institute of Mining and Metallurgy, vol. 68, 1997. pp. 55-67.        [ Links ]

4. ESTERHUIZEN, G.S. An evaluation of the strength of slender pillars. SME Annual Meeting, St. Louis, Missouri., 2006.        [ Links ]

5. MADDEN, B.J. A re-assessment of coal-pillar design. Journal of the South African Institute of Mining and Metallurgy, vol. 91, no. 1, 1991. pp. 27-37.        [ Links ]

6. SALAMON, M.D.G. and WILSON, J.W. An analysis based on a survey of mining dimensions in collieries. Coal Mining Research Controlling Council, No 52/65, 1965.        [ Links ]

7. VAN DER MERWE, N. Beyond Coalbrook: what did we really learn? Journal of the Southern African Institute of Mining and Metallurgy, vol. 106, 2006. pp. 857-868.        [ Links ]

8. BIENIAWSKI, Z.T. and VAN HEERDEN, W.L. The significance of in situ tests on large rock specimens. International Journal of Rock Mechanics and Mining Sciences & Geomechanical Abstracts, vol. 12,1975.        [ Links ]

9. WAGNER, H. Determination of the complete load-deformation characteristics of coal pillars. 3rd ISRM Conference, vol. IIB, 1974. pp. 1076-1081.        [ Links ]

10. SALAMON, M.D.G. and MUNRO, A.H. A study of the strength of coal pillars, Journal of the South African Institute of Mining and Metallurgy, vol. 68, 1967. pp. 56-67.        [ Links ]

11. GREENWALD, H.P., HOWARD, H.C., and HARTMANN, I. Experiments on the strength of small pillars of coal in the Pittsburg Bed. US Bureau of Mines, paper No. 605, 1939.        [ Links ]

12. MADDEN, B.J. The 'squat' pillar design formula and some advantages of its application. The South African Colliery Managers' Association, Circular No 1/88, 1988. pp. 1-25.        [ Links ]

13. WAGNER, H. and MADDEN, B.J. Fifteen years' experience with the design of coal pillars in shallow South African collieries. Design and Performance of Underground Excavations. International Society for Rock Mechanics, Cambridge. 1984.        [ Links ]

14. SALAMON, M.D.G. Unpublished report to Wankie Colliery, 1982.        [ Links ]

15. RYDER, J.A. and JAGER, A.J. A textbook on rock mechanics for tabular hard rock mines. Safety in Mines Research Advisory Committee, Johannesburg, 2002.        [ Links ]

16. BIENIAWSKI, Z.T. A method revisited: coal pillar strength formula based on filed investigations. Workshop on Coal Pillar Mechanics and Design. US Bureau of Mines, IC 9315, 1992.        [ Links ]

17. HUSTRULID, W.A. A review of coal pillar strength formulas. Rock Mechanics, vol. 8, 1976. pp 115-145.        [ Links ]

18. YORK, G. and CANBULAT, I. The scale effect, critical rock mass strength and pillar system design, Journal of the South African Institute of Mining and Metallurgy, vol. 98, 1998.        [ Links ]

19. HILL, F.G. A system of longwall stoping in a deep-level mine, with special reference to its bearing on the pressure burst and ventilation problems. Association of Mine Managers of the Transvaal Papers and Discussions, 1942-1945, vol. 1, 1948. pp. 256-276.        [ Links ]

20. GAY, N.C., ORTLEPP, W.D., RYDER, J.A., STEFFEN, O., and WAGNER, H. The development of rock mechanics in the South African mineral industry, Proceedings of the 12th Congress of the Council of Mining and Metallurgy Institutions, Johannesburg, South Africa, 3-7 May 1982. South African Institute of Mining and Metallurgy, Johannesburg, 1982. pp. 365-376.        [ Links ]

21. KOTZE, T.J. Some observations on the in situ performance of 'hard rock' pillars subject to 'soft' loading conditions. Proceedings of the 8th ISRM Congress, 1995. pp. 1407-1410.        [ Links ]

22. MARTIN, C.D. and MAYBEE, W.G. The strength of hard rock pillars, International Journal of Rock Mechanics and Mining Sciences, vol. 37, 2000. pp. 1239-1246.        [ Links ]

23. HEDLEY, D.G.F. and GRANT, F. Stope-and-pillar design for Elliot Lake Uranium Mines. Bulletin of the Canadian Institute of Mining and Metallurgy, vol. 65, 1972. pp. 37-44.        [ Links ]

24. HEDLEY, D.G.F., ROXBURGH, J.W., and MUPPALANENI, S.N. A case history of rockbursts at Elliot Lake, Szwilski, A.B. and Brawner, C.O. (eds.). 2nd International Conference on Stability in Underground Mining. Society of Mining Engineers, 1984. pp. 210-234.        [ Links ]

25. OZBAY, M.U., RYDER, J.A., and JAGER, A.J. The design of pillar systems as practiced in shallow hard-rock tabular mines in South Africa. Journal of the South African Institute of Mining and Metallurgy, Jan/Feb 1995, pp. 7-18.        [ Links ]

26. WAGNER, H. and SALAMON, M.D.G. Pillar design. Series of lectures given at Sappers Rust, Pretoria, Aug. 1979.        [ Links ]

27. KERSTEN, R.W.O. The design of pillars in the shrinkage stoping of a South African gold mine. Journal of the South African Institute of Mining and Metallurgy, vol. 84, no. 11, 1984. pp. 365-368.        [ Links ]

28. STACEY, T.R. and PAGE, C.H. Practical handbook for underground rock mechanics. Trans Tech Publications, Germany. 1986.        [ Links ]

29. STACEY, T.R. and SWART, A.H. Practical rock engineering practice for Shallow and Opencast Mines, Safety in Mines Research Advisory Committee, Johannesburg, 2001.        [ Links ]

30. LAUBSCHER, D.H. A geomechanics classification system for the rating of rock mass in mine design. Journal of the South African Institute of Mining and Metallurgy, vol. 90, no. 10, 1990. pp. 257-273.        [ Links ]

31. WESSELOO, J. and SWART, A.H. Risk based chromitite pillar design: Part I- Application of locally empirically derived pillar formula, Proceedings SANIRE 2000 Symposium 'Keeping it up in the Bushveld', 2000.        [ Links ]

32. WATSON, B.P., RYDER, J.A., KATAKA, M.O., KUIJPERS, J.S., and LETEANE, F.P. Merensky pillar strength formulae based on back analysis of pillar failures at Impala Platinum, Journal of the Southern African Institute of Mining and Metallurgy, vol. 108, 2008. pp. 449-461.        [ Links ]

33. HOLLAND, C.T. and GADDY, F.L. Some aspects of permanent support of overburden on coal beds. Proceedings of the West Virginia Coal Mining Institute, 1957. pp. 43-65.        [ Links ]

34. RYDER, J.A. and OZBAY, M.U. A methodology for designing pillar layouts for shallow mining. International Symposium on Static and Dynamic Considerations in Rock Engineering, Swaziland, ISRM, 1990.        [ Links ]

35. SPENCER, D. A case study of a pillar system failure at shallow depth in a chrome mine. Hagan, T.O.(ed.). Proceedings of SARES99, 2nd Southern African Rock Engineering Symposium, 1999. pp. 53-59.        [ Links ]

36. NAPIER, J.A.L and MALAN, D.F. The computational analysis of shallow depth tabular mining problems, Journal of the Southern African Institute of Mining and Metallurgy, vol. 107, no. 11, pp. 725-742.        [ Links ]

37. ROBERTS, M.K.C. and CLARK-MOSTERT, V. Is there some commonality with respect to geological structure in the Bushveld Complex and the Great Dyke? The 4th International Platinum Conference, Platinum in transition'Boom or Bust', The Southern African Institute of Mining and Metallurgy, Johannesburg, 2010. pp. 149-156.        [ Links ]

38. JOUGHIN, W.C., SWART, A.H., and WESSELOO, J. Risk based chromitite pillar design: Part II-Non-linear modelling, Proceedings SANIRE 2000 Symposium 'Keeping it up in the Bushveld', 2000.        [ Links ]

39. DAY, A.P. and GODDEN, S.J. The design of panel pillars on Lonmin Platinum's Mines. Proceedings SANIRE 2000 Symposium, Rustenburg, 2000.        [ Links ]

40. HOEK, E. and BROWN, E.T. Underground Excavations in Rock. Institution of Mining and Metallurgy, London, 1980.        [ Links ]

41. PRIEST, S.D. and BROWN, E.T. Probabilistic stability analysis of variable rock slopes. Transaction of the Institute of Mining and Metallurgy, vol. 92, 1983. pp. A1-A12.        [ Links ]

42. PELLS, P.J.N. What happened to the mechanics in rock mechanics and the geology in engineering geology? Stacey, T.R. and Malan, D.F. (eds.) 6th International Symposium on Ground Support in mining and Civil Engineering Construction, Cape Town, April 2008, pp. 1-36.        [ Links ]

43. MOSTYN, G. and DOUGLAS, K. Strength of intact rock and rock masses. GeoEng 2000, Melbourne, 2000.        [ Links ]

44. HOEK, E., KAISER, P.K., and BAWDEN, W.F. Support of underground excavations in hard rock. A.A. Balkema, Rotterdam, 1995.        [ Links ]

45. ESTERHUIZEN, G.S. Jointing effects on pillar strength. Proceedings of the 19th Conference on Ground control in Mining, Morgantown, West Virginia, 2000. pp. 286-290.        [ Links ]

46. ESTERHUIZEN, G.S., DOLINAR, D.R., and ELLENBERGER, J.L. Pillar strength in underground stone mines in the United States, International Journal of Rock Mechanics and Mining Sciences, vol. 48, 2011. pp. 42-50.        [ Links ]

47. GOVERNMENT MINING ENGINEER. Historical and technical review of the Coalbrook disaster. Department of Minerals and Energy, Republic of South Africa, 1965.        [ Links ]

48. MALAN, D.F., PIPER, P.S., POTGIETER, G.M., and DU TOIT, M. Mine-wide panel stability monitoring at Anglo Platinum Union JV, Proceedings of the 1st Hard Rock Safe Safety Conference, Sun City, Sept 2009.        [ Links ]



Paper received Jun. 2010; revised paper received Jul. 2011.



© The Southern African Institute of Mining and Metallurgy, 2011. SA ISSN 0038-223X/3.00 + 0.00.

Creative Commons License All the contents of this journal, except where otherwise noted, is licensed under a Creative Commons Attribution License