Search Documents

By R. Karl Zipf, Francois E. Heuze

Stress analysis programs such as MULSIM/NL, LAMODEL, MinSim 2000, and EXAMINETAB are used in the mining industry to analyze stresses and displacements in coal mines, platinum mines, gold reefs, and tabular-type deposits. These relatively simple numerical models can efficiently simulate yielding and failure of a rock mass near a mine opening and subsequent stress transfer. The main input parameters for these models are a family of nonlinear stress-strain curves for the in-seam material. For many applications of these models, such as in coal mining, extensive field measurements and observations exist from which to estimate these stress-strain curves and calibrate the numerical models; however, such measurements are lacking for other types of mines. A three-step method is presented to determine nonlinear stress-strain curves for boundary-element (BE) programs used in many mining applications. The method requires a suite of laboratory-scale strength tests at various confining pressures. Note that this analysis uses laboratory-measured strength and modulus properties of alluvium. The dependency of the mechanical properties of alluvium on scale has not been established, although more tests are being planned. However, it is believed that the relative comparisons performed are valuable regardless of scale effects. First, the FLAC2D computer program is used to model the laboratory tests and determine cohesion (c) and friction angle () for a Mohr-Coulomb material model. Second, another FLAC2D model uses the c and values to calculate vertical stress and strain around a single opening in the rock mass. This model calculates the stress-strain path of points at various distances from the opening boundary. Finally, based on these stress-strain paths, the stress-strain curves needed for a BE analysis are derived. This three-step method was demonstrated in a large, flat-lying underground test facility in very weak rock. The BE analyses agreed well with observations of failure at the test facility and provided a basis for evaluating the behavior of alternative layouts.

By R. Karl Zipf

Stress analysis programs such as MULSIM/NL, LAMODEL, MinSim 2000, and EXAMINE TAB are used in the mining industry to analyze stresses and displacements in coal mines, platinum mines, gold reefs, and tabular-type deposits. These relatively simple numerical models can efficiently simulate yielding and failure of a rock mass near a mine opening and subsequent stress transfer. The main input parameters for these models are a family of nonlinear stress-strain curves for the in-seam material. For many applications of these models, such as in coal mining, extensive field measurements and observations exist from which to estimate these stress-strain curves and calibrate the numerical models; however, such measurements are lacking for other types of mines. A three-step method is presented to determine nonlinear stress¬strain curves for boundary-element (BE) programs used in many mining applications. The method requires a suite of laboratory¬scale strength tests at various confining pressures. Note that this analysis uses laboratory-measured strength and modulus properties of alluvium. The dependency of the mechanical properties of alluvium on scale has not been established, although more tests are being planned. However, it is believed that the relative comparisons performed are valuable regardless of scale effects. First, the FLAC2D computer program is used to model the laboratory tests and determine cohesion (c) and friction angle (0) for a Mohr-Coulomb material model. Second, another FLAC2D model uses the c and 0 values to calculate vertical stress and strain around a single opening in the rock mass. This model calculates the stress-strain path of points at various distances from the opening boundary. Finally, based on these stress-strain paths, the stress-strain curves needed for a BE analysis are derived. This three-step method was demonstrated in a large, flat-lying underground test facility in very weak rock. The BE analyses agreed well with observations of failure at the test facility and provided a basis for evaluating the behavior of alternative layouts.

Jan 1, 2003

By R. K. Zipf

A method is presented for creating realistic numerical models for practical coal mine ground control. The method includes procedures to collect the necessary mechanical input parameters from a geologic core log, procedures to set up a model and procedures to interpret calculation results. The input parameters come from a detailed geologic core log and extensive point load tests of estimate rock layer strength. A suite of material property input parameters is proposed which allow the user to go from core log to numerical model inputs. Rock bolt anchorage properties are also linked to the material properties of each geologic layer in the model. Following this procedure leads to very realistic calculations of the rock failure process and rock support system behavior. These calculations in turn enable realistic comparison of the effectiveness of alternative rock support systems.

By S. S. SHINDE, P. C. JIKAR, N. B. DHOKEY

Mill scale is an oxide waste product of the steel mills, which contains around 70% iron and other allied impurities. In the present context, a dynamic prototype counter current reactor (CCR) has been used for experimentation and simulation of particle bed by using discrete element method (DEM). Mill scale was pulverized to ASTM mesh range 140 (106 μm)/ + 270 (53 μm), and then it was oxidized in CCR to make single oxide Fe2O3 phase by maintaining required thermodynamic conditions such as pure oxygen at 1100 °C. Then this oxidized powder was subjected to reduction reaction in CCR having gas mixture ( H2:N2) ratio of 1:4 at 875 ± 5 °C for filling degree (fd) and reactor revolution (ω) ranging from 5.38 to 16.14% and 2 to 5 rpm, respectively. A creeping fluid flow condition (Re ˂ 1) and natural heat convection (Gr/Re2 ˃˃ 1) has been maintained while undergoing solid–gas reaction. Besides particle bed behavior has been quantified by Froude number (Fr) for the optimum operating window. The simulation was validated by the experimental results of reactor reaction product by X-ray diffraction (XRD) and scanning electron microscope (SEM) images. Prevailing heat transfer mechanism and reduction reaction mechanism has been established for solid–gas counter current conditions.

Nov 22, 2021

By S. S. SHINDE, P. C. JIKAR, N. B. DHOKEY

Mill scale is an oxide waste product of steel mills, containing around 70 percent iron and other allied impurities. In the present context, a dynamic prototype counter-current reactor (CCR) was used for the experimentation and simulation of a particle bed using the discrete element method (DEM). Mill scale was pulverized to ASTM mesh range 140 (106 μm)/+ 270 (53 μm) and oxidized in a CCR to make a single oxide Fe2O3 phase by maintaining the required thermodynamic conditions, such as pure oxygen at 1,100 °C. This oxidized powder was then subjected to a reduction reaction in the CCR, with gas mixture (H2:N2) ratio of 1:4 at 875±5 °C for filling degree, fd , and reactor revolution, ω, ranging from 5.38 to 16.14 percent and 2 to 5 rpm, respectively. A creeping fluid flow condition (Re < 1) and natural heat convection (Gr/Re2 >> 1) were maintained during the solid–gas reaction, and the particle-bed behavior was quantified by Froude number for the optimum operating window. The simulation was validated with the experimental results of reactor reaction product by X-ray diffraction (XRD) and scanning electron microscope (SEM) images.

By A. P. E. Dirige

Extensive field study has been performed to assess the capabilities of thin sprayon liners (TSLs) and conventional spray-on support systems for preventing rock and support material damage that often results due to mining-induced rockbursts. Support performance was studied using field scale explosive detonation trials to simulate dynamic failure effects that are known to develop during typical rockburst events. Multiple seismic and high speed photographic monitoring techniques were used to provide detailed information concerning rock motion, surface fracturing, ejected fragment motion, and support liner survivability characteristics. Results of this study have demonstrated that TSL?s and variant layer combinations may be as effective as or better than conventional support materials for mitigating rockburst or like damage in highly stressed mine environments. In this study, field results have been verified using numerical modeling procedures to better understand the support behaviour of TSLs when subjected to highly stressed mine environments and mining-induced rockbursts, a factor that is essential when designing any rock support system. A series of numerical modeling assessments were conducted using FLAC3D modeling of a half circular TSL-lined tunnel, influenced by anisotropic stresses and simulated mining-induced rockbursts, to investigate the support potential of TSLs under prototype conditions. The effect of a mine rockburst was simulated using the three-dimensional dynamic analysis option of FLAC3D. Results of TSL-lined tunnels in highly stressed and mine rockburst conditions indicated exceptional effectiveness of TSLs for suppressing rock deterioration resulting from stresses and rockbursts. The rock support potential of TSLs on tunnel surfaces was compared with that of thin shotcrete linings, and results indicated that thin spray-on lining products currently available may be equally as effective in support as shotcrete materials.

May 1, 2009

By S. K. Sahoo, S. K. Sharma, U. K. Singh, G. S. P. Singh

This research focuses on the strata control issues in a Bord and Pillar depillaring working under the influence of dead load condition of voluminous and fragmented softcover of overburden dump material. It reports the modeling of a mine working through the FLAC 2D software under the Indian geo-mining conditions. The modeling results of an actual mine working under softcover were compared with that under the intact overburden condition to assess the influence of the soft overburden over the depillaring working. The modeling work involved progressive caving of the strata and the cyclic goaf filling following the modelsimulated occurrence of main fall and periodic caving. Mohr-Coulomb Failure criteria along with a simplified strain-softening (MCSS) model was used to study the failure and caving mechanism. Salamon’s compaction model was used to simulate compaction and resultant stress recovery in the goaf material. Comparison of the model findings was done in terms of load on supports, abutment stress, and convergence, apart from the stress redistribution and failure mechanism. The load transfer mechanism of the soft overburden was also studied for an overall understanding of the strata behavior in such workings.

Jun 16, 2020

By SK SAHOO, UK SINGH, SK SHARMA, GSP SINGH

Surface mining accounts for most coal production in Indian geo-mining conditions compared with underground methods, whose implementations are in deeper coal reserves and for lower stripping ratio in a congested operating space. The bord-and-pillar method is the preferred underground method in India for reasons of flexibility and investment cost. While previous researchers conducted field investigations on intact overburden conditions, this study focuses on the strata control issues in a bord-and-pillar depillaring working under the influence of the dead-load condition of voluminous and fragmented softcover of overburden dump material. The modeling results helped in ascertaining the criticality of strata behavior and rib-pillar stability in actual mining conditions when compared with the conventional intact overburden.

By Andrew Theobald, Peter Dare-Bryan, Trevor Byers

The use of advanced modelling techniques and the timing flexibility of electronic initiation systems allows for unprecedented control over blast outcomes as shown in several complex blasting problems experienced by the open pit operation Xstrata Ernest Henry Mine (EHM) in north-west Queensland, Australia.

Jan 1, 2010

By Ashok Kumar, Rajendra Singh, Petr Waclawik

High-grade coal is locked in the shaft protective pillars at 900 m below surface in the Czech part of Upper Silesian Coal Basin. To optimise recovery of coal from the larger size shaft protective pillar, room and pillar (called bord and pillar in India) is trialled to further develop these on smaller pillars. This development followed a constant gallery width of 5.2 m and 3.5 m height. The monitored pillar size in the first developed panel varied from 860 m2 to 1225 m2 of area. The stability of these pillars (i.e., natural supports) at this depth of cover is observed to be a vital issue to control ground movement and safety of the surface structures. A simple calculation of safety factor revealed that application of an empirical method for estimating size and shape of the pillar is inappropriate due to complex geo-mining conditions of the site. Accordingly, geo-technical instruments are installed in the panel for monitoring the stability of the diamond shaped pillar at different stages of the development. Extensive monitoring of the pillar behaviour by these instruments provided an important set of data for such optimum recovery. An established process of simulation-based pillar strength estimation is often used by the mining industry, but calibration is generally done with the help of empirical formulation. However, for the deep-seated pillar, the strength suggested by the existing empirical formulations (except the CMRI formula) is found to be much different from the instrument observations. Using a study conducted at Czech hard coal mines, an attempt is made to develop a strain-softening-based numerical modelling approach to assess pillar performance at deeper cover. The recorded value and nature of the pillar spalling in the field along with the CMRI formula are used to calibrate the developed simulation approach. An assessment of the stability of the developed pillars using simulated models did not indicate any collapse of the natural supports. This paper briefly presents difficulties in pillar design at great depth of cover and field and simulation results of pillar strength estimation in Czech hard coal mines.

Jan 1, 2018

By M. W. Grabinsky, M. S. Diederichs, E. Hoek

Provided that they are used correctly, numerical models can be of considerable assistance to mining engineers in designing underground excavations. Guidance is given to the engineer or geologist who is not a specialist in numerical analysis but who is interested in using numerical modelling tools or in employing someone to analyse the stresses, displacements and failure in the rock surrounding underground excavations. A description of the theoretical background to a number of these numerical models and discussion of the advantages and disadvantages of each are followed by guidelines on the selection of appropriate models for each stage in the development of a mine, from exploration to the completion of mining

Jan 1, 1991

By A. J. Basu

Geotechnical design plays a key role in mine design. I A proper geotechnical design provides a safe working environment and facilitates efficient mining operation. There are mainly three design approaches I) judgmental, based on observational and experiential facts, 2) empirical, based on empirical design formulae, and 3) analytical, incorporating numerical analysis. Most design applications are based on the design methods one and two as mentioned earlier. However, complexities arise in deep underground mining operations and deep open pit mines. These complexities and difficult geological environment introduce more uncertainties and simple empirical formulae would not be adequate as a design tool. These situations require numerical modelling and analysis of the design problems including interpretation of the results obtained from the analysis. There are several tools available for performing the numerical modelling and analysis. A brief review of the available tools was addressed earlier (Basu and Kuruppu, 1995). The authors selected ANSYS (a general purpose finite element analysis tool) for use in mining geomechanics and geotechnical design. This paper includes several examples demonstrating the applicability of this software in designing structures applicable to mining geomechnics. The authors hope that the design examples will demonstrate the use of the numerical modelling and analysis application in excavation design.

Jan 1, 1996

By O. K. Mahabadi

The scope of this study is to simulate the behaviour of a layered rock sample under standard laboratory Brazilian Disc test using an innovative combined finite-discrete element method (FEM/DEM) research code. The influence of layering and the direction of loading on the samples behaviour are studied through various layering and loading configurations. A parametric study including the effect of loading rate on the tensile strength of samples has been carried out. This paper demonstrates the suitability of the numerical approach to explicitly model rock deformation and failure.

May 1, 2009

By Lok M. Sharma, Deep K. Khatri, Bradfield. Brett

"A shored enclosure was required to access a control structure upstream below the dam of a deep lake. The shoring involves the interaction of the structure, soil, and water and the loads imposed include the hydrostatic forces of the water, soil pressures as well as dynamic forces due to currents and waves. Deep excavation support systems have traditionally been designed using classical methods that involve direct computations of earth pressures. However, traditional analysis does not take into account soil-structure interaction. Numerical modelling offers better means to analyze soil-structure interaction. This paper presents a 2D finite element (2D-FE) analysis of a deep excavation support structure that was required to remove and replace a control structure at the upstream toe of a dam. The sheet-piling shoring was designed to withstand the lateral pressures from 14 m deep free water and 5.5 m deep excavation below the free water to reach the required excavation level to replace the gate structure. Elastic-plastic analyses were performed to predict the stresses and deformations of the complicated shoring structure and excavation bottom. The soil parameters, procedures used to develop the 2D-FE model, and the results are discussed in this paper. The shoring system is currently under construction.INTRODUCTIONA deep shoring system enclosure is required to keep out water and soil and to provide a support for unbalanced loadings due to dewatering and/or, excavation in order to permit construction of the permanent structures in the dry. A braced shoring system consists of a single wall sheet piling which is driven into the ground to form a “box” around the excavation site and is braced on the inside. After bracing the sheet piling, seal concrete is placed into the bottom to prevent water from seeping into the open excavation from underneath the sheet piling, and the water is pumped out to provide the dry work area. The braced shoring involves the interaction of the structure, soil, and water and the loads imposed include the hydrostatic forces of the water, soil pressures as well as dynamic forces due to currents and waves. Deep excavation support systems have traditionally been designed using classical methods that involve direct computations of earth pressures, (active, passive, or apparent pressures). There have been some uncertainties regarding the actual distribution of earth pressures due to influence of the support wall flexibility, as these traditional methods do not take into account soil-structure interaction. Numerical modelling offers better means to analyze soilstructure interaction."

Jan 1, 2017

By Medhurst TP, Trueman R

The aim of the paper was to model numerically some aspects of the extraction of coal pillars by split and fendering. A two-dimensional elasto-plastic plane strain model with strain softening and interface mod- elling was used. The model was able to simulate the improvement of conditions experienced in lifting from the fender, when the roof adjacent to the fender has caved. In contrast, when a stook was left, providing support to a span of uncaved roof, the model pre- dicted the fender would suffer considerably more de- formation and strain.

Jan 1, 1992

By Ian Masters, Saul B. Davies, David T. Gethin

The purpose of this model is to simulate and analyze the energy flows within a rotary aluminum recycling furnace, using the global properties of its main components, the gas, metal and refractory walls. The model will simulate a series of operations such as loading, heating and stirring which will constitute the processing of a batch. Each main component of the furnace is treated as a one-dimensional heat transfer medium, the thickness and area of each medium being calculated from the actual geometry of the furnace. The physical properties are taken as the weighted average of those for the various layers forming the medium. Incorporated into the model are elements describing the radiative exchange within the furnace, melting of the charge, and forced convection in the melt due to the rotation of the furnace. This approach has produced a model with a reasonable computing time that can be used to run a whole series of simulations for the prediction and comparison of the furnace performance under different operating conditions.

Jan 1, 2000

By K Koskinen, S Ronkainen

The flow field induced in an Outokumpu TankCell and by an OK rotor-stator mixing device in a laboratory and a process size flotation cell was modelled numerically using computational fluid dynamics (CFD). The primary target was to discover and reproduce the various interactions between the flow field and the flotation cell by numerical means. As power consumption depends sensitively on flow field, it can be regarded as a top-level derived quantity (or measure/metric) indicating the validity of the numerically obtained flow field. Having properly validated numerical means to produce the flow field the influence of various design parameters can be captured with certain accuracy. In addition, separate gas and solid phase models for bubbles and solids, respectively, were implemented. The comparisons between velocities and turbulence quantities measured using the laser doppler velocimetry (LDV) technique and CFD modelling were performed. The grid generation was performed with the commercial Gambit 1.3 and the subsequent integration and solving of the resulting algebraic equations with the commercial CFD-code Fluent 6.22. The computationally economic multiple reference frame method was included to model the moving rotor blades. Gas and solid phase models were modified mixture models.

Jan 1, 2005

By R. J. Pine, J. M. Rance, J. S. Coggan

ABSTRACT: The safe design of excavations created in a fractured rock mass is of paramount importance in rock engineering. Depending on the nature of the stress regime acting on an excavation, and the associated deformational behaviour of the rock mass surrounding the excavation, a number of potential failure mechanisms exist. These involve a combination of shear failure on existing joints, extension of critically oriented joints and propagation of new fractures through previously intact rock. Conventional numerical codes usually adopt either a continuum or discontinuum approach; this often results in significant over simplification and kinematic limitations particularly when modelling complex rock mass behaviour in the close vicinity of excavations. This paper describes the results from several case examples using a combined finite element-discrete element code, ELFEN that incorporates a crack propagation mode. Through the use of selected rock engineering examples, the authors demonstrate the potential of the code in the realistic simulation of complex rock failure both in underground and surface excavation applications. RÉSUMÉ: La conception sûre des excavations créées dans une masse de roche fracturée est d&apos;importance primordiale dans la géotechnique. Selon la nature du régime d&apos;effort autour d’une excavation, et le déformation associé de la masse de roche autour de l&apos;excavation, il exsite un certain nombre de mécanismes potentiels d&apos;échec. Ceux-ci comportent une combinaison d&apos;échec de cisaillement sur les joints existants, la prolongation des joints d’orientation critique et la propagation de nouvelles ruptures dans la roche précédemment intacte. Les codes numériques conventionnels adoptent habituellement l’approche, soit d’un continuum, soit de discontinuum; ceci a souvent comme conséquence une sursimplification significative et les limitations cinématiques en particulier en modelant le comportement complexe de la masse de roche très proche des excavations. Cet article décrit les résultats de plusieurs exemples de cas en utilisant un code qui s’agit d’une combination d’élément discret et d&apos;élément fini, ELFEN qui incorpore un mode de propagation de fissure. Par l&apos;utilisation des exemples choisis géotechniques, les auteurs démontrent le potentiel du code dans la simulation réaliste de l&apos;échec complexe de roche pour les applications aux excavations soit au fond soit à l’extérieure. ZUSAMMENFASSUNG: Das sichere Design der Aushöhlungen, die in einer zerbrochenen Felsenmasse verursacht werden, ist vom paramount Wert in der Felsentechnik. Abhängig von der Natur des Druckregimes, das auf einer Aushöhlung und dem dazugehörigen deformationalen Verhalten der Felsenmasse umgibt die Aushöhlung fungiert, bestehen eine Anzahl von möglichen Ausfalleinheiten. Diese beziehen eine Kombination des Scherausfalls auf vorhandenen Verbindungen, Verlängerung der kritisch orientierten Verbindungen und Ausbreitung der neuen Brüche durch vorher intakten Felsen. Herkömmliche numerische Codes nehmen normalerweise entweder ein kontinuum oder discontinuum Annäherung an; dieses ergibt häufig bedeutende Übervereinfachung und kinematische Beschränkungen, besonders wenn es kompliziertes Felsenmasseverhalten in der Nähe der Aushöhlungen modelliert. Dieses Papier beschreibt die Resultate von einigen Fallbeispielen mit einem kombinierten begrenzten Element - getrennten Elementcode, ELFEN, das einen Sprungausbreitungmodus enthält. Durch den Gebrauch von vorgewählten Felsentechnikbeispielen, demonstrieren die Autoren das Potential des Codes in der realistischen Simulation des komplizierten Felsenausfalls bei Anwendungen für beide Tagebaue und Untergrundsaushöhlungen.

Jan 1, 2003

A mathematical model for spontaneous combustion of coal with moisture included is presented in this paper. The one-dimensional unsteady state model consists of conservation equations for oxygen, water vapour and inherent moisture of coal and energy for both gaseous and solid phases. A first order Arrhenius reaction rate for oxidation under both pore diffusion and chemically controlled reaction regime is considered. The rate of evaporation or condensation is also considered as a function of temperature of coal, water content of coal and gas streams. The equation of the model is solved numerically by the finite difference technique. Influences of different parameters on the process of spontaneous combustion can be examined by using this model. Also the model can be used to simulate full-scale storage conditions.

Jan 1, 2005

One of the most important goals of the experiments is to determine deformation and strength of the rock specimen under applying triaxial load, the status of failure condition is one of the subjects which can be used in soil mechanic and foundation engineering applications. To investigate the effect of confining pressure on the triaxial test, the rock has been considered as bonded-particles method and has been simulated by particle flow code three dimensions (PFC3D). This model, does not have pre-existing microcracks and all the particles have the same size (2 mm diameter), equal to the real specimen shape. To compare the mechanical behavior and failure pattern of the numerical and experimental results at uniaxial and triaxial tests at the same conditions, the Young&apos;s modulus, Possion&apos;s ratio as well as maximum axial stress have been considered. On the other hand, the microcracks growth and change of failure pattern at the modelling of the uniaxial and triaxial tests with different confining pressures up to failure point were reported. According to the experimental work, the number density of microcracks decreases from starting value and remains almost constant up to the failure beyond that point. The comparison of the numerical and experimental results of maximum axial stress and the Possion&apos;s ratio revealed good accordance. The initial number density of microcracks at the model was zero and increased while microcracking. Also, the sudden increment around the maximum stress was observed which is because of unstable growing of microcracks near the maximum stress. In all uniaxial and triaxial tests with different confining stresses, the number density of microcracks during the failure remains almost constant which can be considered as a proof of failure occurring in the model. The experimental results indicate the similar trend as well.

Nov 25, 2010