Search Documents

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.

By R. L. Lewis

The White Pine Mine is the second largest active underground copper mine in the United States, covering an area of approximately 34 km2 in the Upper Peninsula of Michigan. A numerical flow model of the White Pine Mine was constructed to evaluate engineering controls that would reduce the rate of mine water discharge after closure. Of the passive controls evaluated, the placement of low-flow bulkheads within the mine workings was predicted to offer the greatest benefit, reducing the rate of discharge by as much as 35 percent.

Jan 1, 1996

By Cheng Li, Yanlei Wang, Lingrui Kong, Jianxin Tang, Fang Yuan

Coal mined-out areas have a significant impact on mining safety and affect the pressure distributions of neighboring coal faces. The PFC3D software is used to study the goaf mechanical characteristics, and the core loading test verifies the rationality of the numerical simulation method. The results from the numerical modeling showed that the void ratio, particle size ratio, and particle size distribution of the rock particle sample (RPS) had a significantly remarkable influence on the compression deformation characteristics of the rock block. The parallel-bond tensile strength and parallel-bond cohesion had a notable impact on the compression deformation traits of the rock block only when the axial stress level was lower than 12 MPa. By contrast, when the axial stress exceeded 5 MPa, the effective modulus significantly influenced the compression deformation features. However, the stiffness ratio, loading rate, and parallel-bond friction angle had almost no effect on the compression deformation characteristics. Furthermore, when the stress level was less than 10 MPa, the stress–strain curves from the rock core samples’ numerical simulation and loading tests were approximately linear. In addition, the Pearson correlation coefficients of the stress–strain curves suggest that the goaf zone’s core numerical simulation results had a robust correlation with the loading test results. The findings of this study can serve as a reference for the field application of PFC3D software and safety assessments in goaves.

Jan 18, 2022

By A. J. Basu

Rock mass properties testing for determining strength parameters is an integral part of geomechanics. These tests are performed either in a laboratory environment or in a field (in-situ). The authors share their experience in using numerical modelling in simulating the laboratory test procedures. This approach helps studying the rock breakage and failure mechanism more efficiently and provides a high degree of visualization. Moreover, imperfections and discontinuities can be easily incorporated in these models. This capability of providing arbitrary fractures and discontinuities facilitates in building a large knowledgebase of rock strength and failure mechanisms. This knowledgebase may then be used to generate heuristics for determining rock strengths for rock mass with fractures and discontinuities.

Jan 1, 1997

By Yuting Xue, Brijes Mishra

The relaxation test on sandstone specimens showed typical behavior in pre-failure region and stepwise behavior in the postfailure region. Overall, more significant stress relaxation occurred within the failed specimens than the intact ones. Numerical simulations were conducted with pre-defined failure plane and with Voronoi tessellation to visualize the relaxation behavior. The model with pre-defined failure plane showed the key role of failure plane and asperity in simulating the step-wise post-failure relaxation behavior. The inhomogeneous stress distribution within failed specimen and the viscous deformation of intact rocks induced high stress concentration at the asperities. The observed step-wise relaxation initiated from the failure of asperity. Furthermore, the simulations with Voronoi tessellation showed time-dependent fracture development during relaxation in post-failure region. The presence of fractures completely changed the stress distribution. Stress concentration occurred at the front area of fractures and at the interacting areas between fractures. During relaxation, fractures still developed with time and the sudden significant increase in the “damage” coincided with the acceleration of stress relaxation leading to step-wise relaxation. Finally, the results showed the possibility of using residual strength as the long-term strength of failed rock.

Jun 25, 2020

In recent years, Computational Fluid Dynamics, CFD, has been commonly utilized in the mining industry to model the fluid flow behavior in underground mine workings. This paper uses CFD modeling to simulate the airflow behavior in underground crosscut regions, where brattice sails are used to direct the airflow into these regions. Brattice sails are cost effective ventilation control devices for temporary or permanent use in underground mining. They can be used to deflect air into the unventilated areas such as crosscut regions. Their design and installation is a fundamental issue for maintaining a sufficient supply of fresh air and achieving effective air circulation and contaminant removal. At the same time, they should have little impact on the mine ventilation system. This paper presents the results of a two-dimensional CFD model, which examines the effects of brattice length on fluid flow behavior in the crosscut regions. The results of this study will assist in understanding the ventilation air behavior and in determining the optimum size of brattice curtain (sails), which provide a highly effective contaminant removal from the unventilated mine areas. This, in turn, helps the mine ventilation designers to meet the mine safety requirements.

Jan 1, 2009

By C. Sun

In this study, numerical simulations of yield pillars employing the finite-element method with a time-dependent Drucker-Prager material model are conducted. The stress¬control mechanism of yield pillars on weak floor strata is studied. The development of a pillar yield zone, the distribution of pillar stress and the development of a failure zone in the weak floor strata are investigated. The study indicates that a properly selected yield pillar size may significantly reduce the peak stress induced in the weak floor strata and may effectively control the development of a failure zone in the floor.

Jan 1, 2000

By B. Pesic

Northwest Alloys, a wholly owned subsidiary of Alcoa, produces magnesium from dolomite by reduction with ferro silicon and aluminum. This process produces alkaline by-products in quantities that allow their consideration for use as commercial by-products. Two such products are dicalcium silicate slag, Di-Cal, and magnesium oxide condensate, MgO. The alkaline properties of these materials offer acid neutralizing capability comparable to the traditionally used lime, therefore, warranting a study for their use in applications where the neutralization of acidic waters or process solutions is required, such as treatment of mine waters. In this ongoing study, three mine waters (Bunker Hill Mine, Berkeley Pit, and Blackbird Mine) were selected for treatment with Di-Cal and MgO. Lime was used as the background control. The solid precipitates were subsequently contacted with pure water at four different pH conditions for stability testing of prepared metal precipitates (expanded TCLP). It was found that neutralization with Di-Cal results in precipitates which are of equivalent to superior stability than those from lime treatment.

Jan 1, 1999