Search Documents

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 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 Juan Andres Jarufe, Carlos Cifuentes

Numerical modelling has been extensively used to estimate plastic displacement in faults and rock mass surrounding mine excavations. Particularly, for fault slip seismicity, modelled plastic displacement is usually utilized to estimate seismic potential, however these methods assume that all the displacement caused during an specific simulated excavation step, is released as a single event at each step. This representation corresponds to a conservative estimation since rock mass displacement does not occur as one single event but as an exponential distribution of many small displacements and a few large ones. When displacement is expressed as seismic magnitude, this exponential distribution is addressed as the Gutenberg Richter (GR) relationship. This paper present a method to decompose the results from numerical models into a distribution of seismic events The reconstruction of the GR relationship from modelling displacement permits not only a more realistic approximation of the expected maximum magnitudes but also provides a tool to develop a probabilistic estimation of the maximum possible magnitudes.

Apr 26, 2022

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 W. Keilich

A methodology of subsidence prediction using the Distinct Element code UDEC has been developed as an alternative for subsidence modelling in the Southern Coalfield, New South Wales, Australia. The models have been validated by comparison with empirical results and observed caving behaviour. At this stage, modelling capability is limited to flat lying terrain. It is planned to apply the methodology to areas of high topographical relief to investigate the mechanics of valley closure.

Jan 1, 2006

By Andrzej Walentek

This paper presents the results of numerical modeling of the fractured zone around longwall gateroads. The numerical computations have been based on the results of underground measurements of the fractured zone around longwall gateroads, which were performed for various geological and mining conditions and measured with the use of the endoscopic method with a digital camera. Having results of the underground investigations, trials of duplicating them by the Finite Element Method using software PHASE and the Hoek-Brown failure criterion were performed. In order that the fractured zone predicted by the model-based computations closely simulate that measured underground, some selected parameters in the Hoek-Brown failure criterion were modified based on the earlier determined relations used for prediction of the roof-to-floor convergence in longwall gateroads. As a result, the model-based computations produce results very close to those obtained from the underground measurements. Determination of the fractured-rock zone around the mine roads is of great importance to mining practice, especially in planning and designing supports for longwall gateroads. The influence of the front abutment pressure causes the fractured rock zone around the roads to increase and results in an increase of support load. Moreover, proper determination of the fractured rock zone is one of basic elements in the process of planning roof bolting for mine workings.

Jan 1, 2009

Floor heave events in underground mining often happen suddenly and cause risk to workers and equipment. Reducing the risk of floor heave must be addressed in mining layout design. To improve that design one must understand the significant parameters that are correlated to the occurrence of floor heave. To accomplish that goal, researchers from the National Institute for Occupational Safety and Health (NIOSH) completed a numerical study using FLAC3D, taking advantage of its strain softening and ubiquitous joint capabilities to better simulate the floor heave and caving mechanisms involved at an actual mine site.

Jan 7, 2020

By S. P. Girrens

A coupled finite element geomechanical- hydrology code is currently under development for application to the problem of predicting groundwater disturbances associated with mine subsidence. The structural-fluid coupling is addressed by calculating the subsided mine geometry, with emphasis placed on determining the strata disturbance and locating damaged regions, for input into a hydrology code, which determines localized volume flow rates and aquifer fluctuations. Benefits from coupling will be best realized when field measurements, an additional aspect of the study concurrent with analytical investigations, indicating the relationship between increasing rock strain and increasing permeability are incorporated into hydraulic material descriptions. Hydrologic and struc- tural calculations are presented to demonstrate computational capabilities applicable to mine subsidence

Jan 1, 1982

By W. M. Ma, Daniel W. H. Su, K. Centofanti, Yi Luo, W. L. Zhong, Syd S. Peng

NTRODUCTION Surface subsidence monitoring is the most effective means for investigating the nature of, and developing the techniques for pre¬dicting the surface movement and deformation process. However, because of imperfection in survey equipment and field practices, the surface movement data collected from field survey are not error-free. Also due to lack of proper numerical techniques for computing surface deformation from the measured surface move¬ment data, a relatively large portion of the subsidence data col¬lected could not provide information other than surface movement. In a paper by Luo and Peng (1992), the sources and magnitude of errors in the measured surface movements using an electronic total station in subsidence survey are discussed. If the errors in the measured surface movement components are left untreated before calculation for deformation indices (i.e., slope, strain and curva¬ture), the errors will be carried over to the deformation indices and enlarged if improper numerical procedures are used in the calcu¬lation. It should be noticed that in protecting structures affected by ground subsidence, accurate surface deformation indices are much more important than surface movement. In order to improve the accuracy of the collected surface movement data and the calcu¬lated deformation indices, a number of numerical procedures are also proposed in the same paper. The proposed techniques include: a data smoothing technique for data pre-processing; and a number of numerical methods and/or procedures for calculating surface deformation indices from the measured surface movement data. DATA SMOOTHING Because of various factors, the measured surface movement data always contain errors, small or large. In order to better utilize the collected data, a data smoothing method is proposed to correct the errors in the measured movement data. The method uses the same concept as the influence function method (e.g., Knothe, 1957) on which many existing subsidence prediction methods are based. The data smoothing method makes necessary corrections on the measured movement for the point of interest according to not only the measured movement at this point but also those at the surrounding points. The contribution of a surrounding point in making the correction for the point of interest is inversely propor¬tional, in exponential form, to the distance between them. The new magnitude of a movement component at the point of interest (Zi*) is calculated by

Jan 1, 1992

By H. Huang, J. C. Wu, C. H. Xie, Y. Q. Xue, Z. H. Zhang

Two- and three-dimensional groundwater flow models for dewatering in an opencast mine are presented. The models focus on solving special problems using numerical simulations of the opencast mine. These special problems include the drying up of the aquifer, the movement of the model boundary and the prediction of the residual head. These models are used to describe the decline of groumdwater levels in the Yuanbaoshan Opencast Coal Mine and its neighboring area during the drainage period. The simulated values agree very well with the field data. There have been several successful examples of numerical simulation for predicting the mine dewatering rate in an opencast mine. However, The hydrogeological conditions for the Yuanbaoshan Opencast Coal Mine are rather complex, and the total average flow rate is very high at 4.676 x 105 m3/d (1.65 x 107 cu ft/day). In this paper, a two-dimensional groundwater flow model in a large area is developed to solve the problems mentioned above by simulating regional drawdowns of the groundwater level during the drainage period, and a three-dimensional flow model is developed to simulate the residual head near the open-pit bench.

Jan 1, 1998

By I. R. McGowan, J. Tilk, J. I. Hershman, R. L. Lewis

The White Pine Mine is a large underground copper mine located in the Upper Peninsula of Michigan. A numerical ground water flow model of the White Pine Mine was constructed to evaluate potential engineering controls that would reduce the rate of mine-water discharge after closure. Of the passive controls evaluated, the placement of bulkheads within the mine workings was predicted to offer the greatest benefit, reducing the rate of discharge by as much as 35%.

Jan 1, 1998

By R. L. Lewis

The White Pine Mine is a large underground copper mine located in the Upper Peninsula of Michigan. A numerical ground water flow model of the White Pine Mine was constructed to evaluate potential engineering controls that would reduce the rate of mine-water discharge after closure. Of the passive controls evaluated, the placement of bulkheads within the mine workings was predicted to offer the greatest benefit, reducing the rate of discharge by as much as 35%.

Jan 1, 1997

By B. K. Mishra, R. K. Rajamani

This paper focuses on the use of a numerical tool known as the discrete element method (DEM) to study the motion of ball charge in ball mills. DEM is employed to simulate the motion of individual balls. While doing so, the effect of liner design on the grinding performance of the mill is studied. Two important, factors, such as the configuration of the liner and the geometry of the lifting surfaces, are correlated with the overall mill performance. An attempt is made to show the effectiveness of the lifter bars in transmitting power from the mill shell to the interior of the mill, with particular emphasis on large diameter mills.

Jan 1, 1994

By B. K. Mishra

This paper focuses on the use of a numerical tool known as the discrete element method (DEM) to study the motion of ball charge in ball mills (14-18 ft). DEM is employed to simulate the motion of individual balls. While doing so, the effect of liner design on the grinding performance of the mill is studied. Two important factors such as the configuration of the liner and the geometry of the lifting surfaces are correlated with the overall mill performance. An attempt is made to show the effectiveness of the lifter bars in transmitting power from the mill shell to the interior of the mill with particular emphasis to large diameter mills.

Jan 1, 1992

By Murali M. Gadde, Syd S. Peng, C. T. Holland

Cutter roof failure is one of the most common ground control problems affecting the safety and economy of an underground coal mine operation. Several attempts were made in the past to understand the processes responsible for cutter development and to devise mitigating measures. Yet, many cutter failures still happen and their advance prediction remains as elusive as before. One of the reasons for the impasse in predicting cutters lies in the limitations of the current numerical modeling technologies. The development and propagation of cutters is a progressive phenomenon and continuum models can only approximate this process. Since the discontinuum modeling is severely restricted by the lack of reliable input data, alternatives must be sought to better understand cutter failures using the continuum models. This paper proposes one approach using strain-softening material model to simulate the cutter initiation and propagation process. An attempt has also been made to provide some insights into the reasons for the change in cutter location at different parts of a development layout.

Jan 1, 2005

By Menglong Bian, Xiaotong Dong

Thermal damage is an urgent problem faced by deep mining. High-temperature rock mass, hot water and electromechanical equipment are the sources of the thermal damage. The increase in air temperature in mine tunnel is the result of fluid–solid coupling heat transfer from various heat sources. This paper takes a roadway with hot water and a mine car as the research object. Based on the Fluent, Static Structural and Steady-State Thermal modules in ANSYS, the thermal-fluid–structure coupling analysis in the roadway is carried out. It is found that the hot water has a significant effect on the thermal environment of the mine. Under the experimental parameters, the maximum heat flux density and equivalent elastic strain at the interface between the hot water and the rock are 9789.3 W/m2 and 0.00054532 m/m, respectively. The heat flux distribution in the rock area near the mine car does not change much, and the equivalent elastic strain value in this area is high. After the air is heated in the hot water area, the temperature curve changes gradually and then rises steadily. By changing the inlet air temperature parameter, it is found that the lower the inlet air temperature, the better the cooling effect, and the change of the inlet air temperature can change the air temperature field in the roadway as a whole. By changing the parameter of inlet air speed, it is found that higher wind speed brings better cooling effect, but with the increase in wind speed, the cooling effect of roadway thermal environment is very limited. When the wind speed increases from 5 m/s to 10 m/s, the air temperature of most measuring points drops by more than 2℃; when the wind speed increases from 15 m/s to 20 m/s, the air temperature of most measuring points drops below 0.7℃. After the wind speed reaches a certain value, the cooling effect of ventilation is very limited.

May 5, 2022

By Yuting Xue, Danqing Gao, Brijes Mishra

Heterogeneity and discontinuity significantly affect rock strength. For accurate stability prediction, intact rock behavior is imperatively included in rock mass behavior. However, past research largely used arbitrary scaling approaches to produce rock strength. This paper considers a stochastic approach in order to determine the strength of a rock. Based on the laboratory data, a MATLAB® with extreme value stochastic model generates a database for each physico-mechanical property. Then FLAC® simulates laboratory-sized rock specimens. The grids developed in the numerical model can in turn develop random material properties in MATLAB®, which researchers then applied to the final FLAC® model. Model runs simulate the approach performed in the laboratory. The results from the model indicate that a stochastic approach produces strengths that are lower than a deterministic approach. Failure modes for each specimen are different, similar to observations in the laboratory. In addition, random density also influences the failure mode, highlighting the importance of stochastic analysis in rocks.

By Yuting Xue, Danqing Gao, Brijes Mishra

Heterogeneity and discontinuity significantly affect rock strength. For accurate stability prediction, intact rock behavior is imperatively included in rock mass behavior. However, past research largely used arbitrary scaling approaches to produce rock strength. This paper considers a stochastic approach in order to determine the strength of a rock. Based on the laboratory data, a MATLAB® with extreme value stochastic model generates a database for each physico-mechanical property. Then FLAC® simulates laboratory-sized rock specimens. The grids developed in the numerical model can in turn develop random material properties in MATLAB®, which researchers then applied to the final FLAC® model. Model runs simulate the approach performed in the laboratory. The results from the model indicate that a stochastic approach produces strengths that are lower than a deterministic approach. Failure modes for each specimen are different, similar to observations in the laboratory. In addition, random density also influences the failure mode, highlighting the importance of stochastic analysis in rocks.

By Yuting Xue, Danqing Gao, Brijes Mishra

Heterogeneity and discontinuity significantly affect rockmass strength. Past studies have largely used arbitrary scaling approaches to produce rockmass strength from laboratory rock tests. This paper proposes a stochastic approach to produce rockmass strength from laboratory rock tests. Based on the laboratory data, MATLAB with an extreme value stochastic model generated a database for each physico-mechanical property of the rock. Then finite-difference software FLAC2D simulated laboratory-sized rock specimens. Random material properties in MATLAB were generated using the grids developed in the numerical model. The random properties were then applied to the final FLAC2D model. Model runs simulated the approach performed in the laboratory. The results from the model indicate that a stochastic approach produced compressive strengths that were lower than those from the deterministic approach. Failure modes for each specimen were different from those observed in the laboratory. In addition, random density also influenced the failure mode, highlighting the importance of stochastic analysis in rocks.