Search Documents

By Essie Esterhuizen, Michael M. Murphy, Chris Mark

"The design of underground coal mines requires a clear understanding of the overburden response, the loading of pillars, the loading of the gob, the pillar failure process, and the ultimate load carried by partly or fully yielding pillars. Very few highquality stress measurements of yielding pillars and gob loading have been made in full extraction mining. Well-calibrated numerical models can assist in providing a better understanding of the load and failure processes, provided the coal, the overburden, and the gob are all modeled with sufficient realism. A program of numerical model calibration and validation was carried out using FLAC3D.1 The models were calibrated against observed and measured performance of coal pillars and the overburden in operating mines to provide a basic set of input parameters that can be used to provide a realistic first estimate of expected ground response and pillar loading. Input parameters for modeling coal pillar response were based on data from triaxial testing on coal samples, combined with both matching the depth of failure in the coal ribs to observations as well as matching the peak pillar resistance to an empirical equation. The models were calibrated against strong roof and floor case histories in which the pillar strength is governed by failure and yielding of the coal within the pillar and the surrounding strata only had a limited impact on pillar strength. Input parameters for the overburden were determined from a large database of laboratory tests and model calibration against maximum subsidence and subsidence curvature. Further overburden calibration was carried out by matching stresses in the mining horizon to field measurements. Three examples of the application of the calibrated dataset and modeling methodology to field measurements are presented. The results show that a reasopable estimate of the in-seam stress distribution and overburden response can be obtained for both strong and weak overburden scenarios at various depths of cover.INTRODUCTIONThe planning and design of coal mine excavations requires reliable estimates of the expected strength and loading of the mine structures to achieve global stability. Empirical methods are widely used to estimate the strength and loading of coal pillars and have been incorporated into pillar design procedures, such as Analysis of Longwall Pillar Stability (ALPS) (Mark, 1987) and Analysis of Retreat Mining Pillar Stability (ARMPS) (Mark and Chase, 1997), that are widely used in the United States. ~umerical models are finding increasing application as a tool for underground mine design because of their. versatility and the ever increasing computational power available to mine designers. A prerequisite for the application of numerical models is the calibration of the models against observed rock mass response (Hoek et al., 1990; Skiles and Stricklin, 2009). The calibration process may include comparison of model results to measured stress and deformation and modifying the input parameters in a systematic manner to achieve a satisfactory agreement between the model results and measurements. Once models have been calibrated, they can be applied to evaluating similar mining layouts in similar geological conditions."

Jan 1, 2010

By K. M. Mohamed, I. B. Tulu

"Failures of coal ribs are major hazards in underground coal mines. For the period of 2010 to 2014, rib falls were the leading cause of groundfall fatalities and accounted for 35% of all groundfall-related fatalities. Currently, the National Institute for Occupational Safety and Health (NIOSH) is conducting a research project to develop an engineering-based design methodology for rib control in coal mines to reduce rib fall accidents and fatalities in underground coal mines. Numerical models can be used as part of the design methodology provided the models are appropriately calibrated and verified. The paper presents a coal rib model to simulate the mechanism of the coal deformation and fracturing that can lead to stress-driven rib falls. The proposed coal rib model was calibrated against published data from an instrumented rib site at the West Cliff mine in Australia. Fracturing of the coal rib was defined by critical values of shear strain of the coal material and rib displacement. A critical plastic shear strain of 2.8% and a sudden increase of rib displacement of 30 to 45 mm were sufficient to develop a fracture plane in the rib model. Accordingly, a discontinuous rib deformation was achieved and it compared favorably against field monitoring data. INTRODUCTION Coal ribs are the vertical walls of coal formed when entries are developed in a coal bed. Failures of coal ribs are a major hazard in underground coal mines. Furthermore, failing ribs can contribute indirectly to roof and floor instability by increasing opening widths across intersections and entryways. Over the last 18 years, there has been a continual occurrence of injuries and fatalities in underground coal mines due to rib falls. For the period of 2010 to 2014, rib falls (excluding longwall face and bursts) were the leading cause of ground fall fatalities and accounted for 35% of all groundfall-related fatalities [1]. Despite the threat of rib falls, the amount of rib falls and rib support research conducted in the U.S. is limited compared with what has been done in pillar and roof support design. Currently, the National Institute for Occupational Safety and Health (NIOSH) is conducting a research project to develop an engineering-based design methodology for rib control in coal mines to minimize the number of rib fall accidents and fatalities in underground coal mines."

Jan 1, 2016

By Michael Murphy, Khaled M. Mohamed, Berk Tulu

"Failures of coal ribs are major hazards in underground coal mines. For the period of 2010 to 2014, rib falls were the leading cause of groundfall fatalities and accounted for 35% of all groundfall-related fatalities. Currently, the National Institute for Occupational Safety and Health (NIOSH) is conducting a research project to develop an engineering-based design methodology for rib control in underground coal mines to reduce rib fall accidents and fatalities. Numerical models can be used as part of the design methodology provided the models are appropriately calibrated and verified.This paper presents a coal rib model to simulate the mechanism of coal deformation and fracturing that can lead to stress-driven rib falls. The proposed coal rib model was calibrated against published data from an instrumented rib site at the West Cliff Mine in Australia. Fracturing of the coal rib was defined by critical VALUES (of plastic shear strain of the coal material and rib displacement. A critical plastic shear strain of 2.8% and a sudden increase of rib displacement of 30 to 45 mm were sufficient to develop a fracture plane in the rib model. Accordingly, a discontinuous rib deformation was achieved and it compared favorably against field monitoring data.INTRODUCTIONCoal ribs are the vertical walls of coal formed when entries are developed in a coal bed. Failures of coal ribs are a major hazard in underground coal mines. Furthermore, failing ribs can contribute indirectly to roof and floor instability by increasing opening widths across intersections and entryways. Over the last 18 years, there has been a continual occurrence of injuries and fatalities in underground coal mines due to rib falls. For the period of 2010 to 2014, rib falls (excluding longwall face failure and bursts) were the leading cause of ground fall fatalities and accounted for 35% of all groundfall-related fatalities (MSHA, 2015)."

Jan 1, 2016

By Morgan M. Sears, Gamal Rashed, Jim Winfield, Brent Slaker

"Over the past decade, ground falls have resulted in five fatalities, representing over 40% of all fatalities occurring in underground stone mines in the United States. In conjunction with recent high-profile pillar collapses, the need for further ground control research in the nation’s stone mines is evident. Ground conditions are growing more difficult. Therefore, this research is focused on deep cover, steeply dipping, and multi-level underground nonmetal mines. The objective of this paper is to introduce a geomechanical study of pillar development in a steeply dipping limestone deposit. To accomplish this, a numerical model was developed using input from modified laboratory test data and in situ stress measurements provided by the mine to simulate the development of a pillar under these conditions. The model was then verified against field observations and stress mapping. The calibrated model was then used to calculate the estimated stresses in the pillar and immediate roof as development progresses around an instrumented pillar.INTRODUCTIONGround falls represent a significant hazard in the nation’s underground stone mines. In response to recent large-scale pillar collapses, as well as stakeholder interactions, researchers from the National Institute for Occupational Safety and Health (NIOSH) are currently conducting detailed investigations into the complex loading conditions at mines operating in challenging operational environments. These include deep-cover, steeply dipping, and multi-level mines where there is the potential for an increased risk of ground failure. The Pleasant Gap Mine, shown in Figure 1, operates in the Valentine Formation, which is typically light grey, extremely fine-grained limestone, and approximately 21.3 m (70 ft.) thick. The Centre Hall Formation makes up the immediate and main roof members and is approximately 9.1 m (30 ft.) thick, consisting of medium-dark-gray limestone (Murphy et al., 2018). Situated in a syncline, the top of the Valentine Formation dips at approximately 19° with a strike of approximately N54E in the area of interest (see Figure 2)."

Jan 1, 2018

By C. Caceres

Fibre reinforced shotcrete in combination with other types of support systems is becoming the leading support method for underground excavations, especially those present in low quality rock masses. Sudden or rapid ground support failure can expose workers to safety hazards, damage to equipment and others resulting in production interruptions that will impact the cost of operation. Knowing the peak and post-peak strength characteristics of the support system over time are of major importance. Round Determinate Panel Testing of fiber reinforce shotcrete is the primary test used by the mining industry to determine strength characteristics; however, the tests are costly and significant manpower is required in a mine environment to build, store, handle and test the panels. For these reasons numerical modeling can result in an effective method to produce realistic results. This paper examines the strength characteristics Round Determinate Panel tests and its ultimate application in underground excavations, using Itasca?s three dimensional particle flow code (PFC3D).

Jan 1, 2010

By Xinghai Lei, Nan Wang, Xiangyang Jin, Yuqi Ren, Yang Li

The upward mining method, as an effective mining method, is widely employed in close-multiple coal seams to extend the service life of coal mines and improve the recovery rate of residual coal resources. However, the upward mining feasibility of the upper coal seam is related to the height of caved zone of interlayer strata after lower coal seam extraction. Therefore, it is extremely important to determine the height of the caved zone after the lower coal seam. The No. 9 coal seam (upper coal seam) and the No. 12- 1 coal seam (lower coal seam) in the Qianjiaying coal mine are selected as the research object. The characteristics of the interlayer strata structure field, fracture field, and subsidence field are studied using UDEC numerical modeling, in which the height of the caved zone is gotten after the No. 12- 1 coal seam extraction. Furthermore, the height of the caved zone is verified by onsite detection of GPR and BTVI. In the numerical modeling, the height of the caved zone is 19.6 m. And the height of caved zone is 14.7 m and 15.4–21.2 m, respectively, by using the onsite detection of GPR and BTVI. The heights of the caved zone onsite detection results are consistent with the numerical modeling result. The No. 9 coal seam is always located above the caved zone after the No. 12- 1 coal seam extraction. Consequently, the upward mining of the No. 9 coal seam is feasible.

Jun 28, 2023

By M. Niinae

Electrokinetic soil remediation by electroosmosis and ionic migration is restricted to soluble sub-stances. At a higher pH, the soil particles sorb more heavy metals, such as cadmium, than at a lower pH and the precipitation reactions of heavy metals are promoted. Both mechanisms make heavy metals immobile, rendering clean up more difficult and electrokinetic extraction inefficient. There-fore, the acidification of soils is very important to increase the clean up efficiencies of heavy metals, such as cadmium by electrokinetic method from contaminated soils. In the present study, acidification of soils and cadmium transport are predicted with a numerical method in which electrical flow is coupled with hydraulic flow, and the sorption process in the soils is induced in the model. The effects of some factors such as electrode configuration, electric field strength and electrode spacing for saturated soils are investigated using this numerical model. The results obtained this study are summarized as follows: The two-dimensional electrode configurations containing cathode surrounded by three or more anodes rapidly acidify the soils and remove cadmium from soils compared with one-dimensional electrode configuration. However, two-dimensional electrode configurations develop spots of inactive electric field between cathodes. A larger electrode spacing increases the processing time required, and it was also found that a larger electric field strength reduced the processing time required.

Jan 1, 2014

By Ed Van Eeckhout

As part of a U.S. Bureau of Mines-sponsored project on the utilization of "stiff" backfill to minimize potential rock burst hazards, a finite element study was undertaken to predict stresses and displacements around a stope mined up from the 5400 level of the Sunshine mine, Kellogg, Idaho. This stope was to be backfilled using a high-modulus fill. In preparation for that configuration, certain instrumentation was installed near the 5400 level to monitor stress changes. Due to untenable economics, only 2 cuts of the stope were ever made. This paper presents the results of the instrument readings taken during the first cut and compares them with the numerical predictions, as well as presents the modeling results of closure and fill stresses in the stope if it had been mined.

Jan 1, 1984

By Ed Van Eeckhout

As part of a U.S. Bureau of Mines-sponsored project on the utilization of "stiff' backflll to minimize potential rock burst hazards, a finite element study was undertaken to predict stresses and displacements around a stope mined up from the 5400 level of the Sunshine mine, Kellogg, Idaho. This stope was to be backfilled using a high-modulus fill. In preparation for that configuration, certain instrumentation was installed near the 5400 level to monitor stress changes. Due to untenable economics, only 2 cuts of the stope were ever made. This paper presents the results of the instrument readings taken during the first cut and compares them with the numerical predictions, as well as presents the modeling results of closure and fill stresses in the stope if it had been mined.

Jan 1, 1984

By Thomas Novak, Sampurna Arya, Joseph Sottile

Dust generation and its dispersion in a longwall mine pose significant health and safety issues. The concern is prominent at both the working face, where coal excavation exposes workers to high dust concentrations, and also in the return, where float dust accumulates. To address this issue, researchers at theUniversity of Kentucky designed and fabricated a flooded-bed dust scrubber integrated into a full-scale physical model of a longwall shearer. This mock-up was installed in the longwall test gallery at the Pittsburgh Research Laboratory (PRL) of the National Institute for Occupational Safety and Health (NIOSH). A full-factorial experiment was designed to measure the effectiveness of the scrubber on dust capture. Test results indicated dust reduction of up to 57% in the return airways. This paper presents the development of a computational fluid dynamics (CFD) model and its validation by comparing the simulation results against the experimental results. The CFD simulations show good agreement with experimental results with differences between + 6.9 and − 9.7%. The validated CFD model was used further to predict the scrubber performance under additional operating conditions, resulting in a dust reduction of up to 66.5%. With this level of success, the scrubber concept can be a game changer for the mining industry, and it can be utilized to improve health and safety in longwall mines.

May 26, 2020

By Eleir Mundim Bortoleto

A finite element modeling approach was developed to study adhesion phenomena during the unidirectional contact between a two-dimensional plane-strain square and a flat slab. Surfaces were metallic or ceramic, different pairs of materials were analyzed and their adhesion was considered by means of a FORTRAN subroutine (DLOAD), which was connected with commercial finite element code Abaqus and provided surface attractive forces based on the Lennard-Jones interatomic potential. Adhesive loads were considered both during approximation and separation of the surfaces. During the separation step, material transfer between surfaces, due to adhesion, was modeled by means of damage initiation and propagation at the flat slab. The parameters considered in the simulations were normal load, chemical affinity and system size and the different conditions were analyzed by comparison of the interaction forces during approach and withdrawal. This work also presents detailed descriptions of: (i) the evolution of energy dissipation due to adhesion hysteresis, (ii) the formation-growth-breakage of adhesive junctions and (iii) the evolution of stress distribution during the process. Results have indicated that chemical affinity between the bodies in contact is more relevant for adhesion than the applied load. Besides, smaller asperities are less prone to promote adhesive wear, but they may be subjected to high stress concentration effects.

Jul 31, 2018

By N. Manopulo, J. Maier, L. Tong, C. Becker, P. Hora

The contribution presents a new finite element approach in which the common updated Lagrangian formulation is replaced by a spatial flow formulation applying a special arbitrary Lagrangian-Eulerian for-mulation (ALE). The description separates spatially the process in a die inside zone and in a profile zone. For the description of the inside zone, where the filling process occurs, a special Euler time independent mesh is applied. By this formulation the mesh has to be generated only once, which simplifies the generation of spatially optimized meshes. For mapping of the material flow during the filling process the method uses a flow surface traction method as applied in casting simulations. A further strength of the new method is the computational efficiency which allows the modelling of extrusion process with the highest complexity with a response time of 12 hours, which is important for the industrial optimization process.

Jan 1, 2015

By Adriana Amaro Diacenco

This work presents the finite element modeling of Euler-Bernoulli beams constituted of a metallic material, discussing theoretical and numerical aspects. The numerical analysis uses Hermitian interpolation functions and finite elements with two degrees of freedom per node. A numerical simulation will be done for uniform loading considering two kinds beams (steel and aluminum). The paper is organized as follows: introductory comments and numerical simulations to evaluate the static and dynamical behavior of the metallic beams.

Aug 7, 2018

By N. Golchinfar

Tunnels in deep mines are subjected to high stresses. Fracture and failure of the rock mass around an opening can occur when mining-induced stresses are high. In burst-prone grounds, mining-induced seismicity can cause additional dynamic loading which can further increase the stress around the tunnels. When the dynamic stress reaches the rock mass strength, fracturing can occur and rockburst may happen (depending on loading stiffness around the failed rock mass). Rock support installed in the tunnel must be capable of dissipating dynamic energy and holding the failed rocks. Hence, it is important to estimate the depth of failure for rock support design. This paper focuses on how to model depth of failure under dynamic loading. In the modeling, two scenarios of rock failure are considered. In the first scenario, rock failure occurs under static loading and subsequent dynamic loading further increases the depth of failure. In the second scenario, no failure occurs under static loading but in the subsequent dynamic loading, failure around the tunnel is created.

Aug 1, 2013

By M. N. Vu, G. Armand, A. Noiret, D. Seyedi

Extensive experimental observations have been performed around drifts and shafts at the Meuse Haute-Marne (North-Eastern France) Underground Research Laboratory (URL) to assess the extent and pattern of the induced fracture networks in Callovo-Oxfordian claystone. Two types of fractures can be distinguished, namely; shear (mode II) and extension unloading (mode I) fractures. The main goal of the present paper is to provide insights on the extent of damaged zone around the drifts by numerical simulations. A two-step approach is considered. First, 2D elastic calculations are performed for simulating the excavation of a drift following the direction of the major in situ horizontal stress. The evolution of the stress and strain tensors around drifts are calculated by applying an elliptic convergence function on the drift wall. Two failure criteria are then verified for the induced damage. An equivalent strain and an elastic energy based damage criteria are considered. The obtained results show that the employed simplified approach provides rather good results for the extent of the damaged zone for drifts parallel to major horizontal stress. However, post-processing elastic results leads to overestimating the vertical extent of the damage zone and underestimating its lateral extent, where the in situ stress state is quasi-isotropic in a section of drift. Nevertheless, the obtained results show that an energy based damage criterion could be considered as a potential way for reproducing the extent of the damaged zones if a convergence function can be defined for the drift.

Jan 1, 2015

By Tu Bao

HP-Pile foundations are the most favorable type of foundations for bridges because of their constructability, versatility with regard to subsoil conditions, and settlement controllability. For many years, bridge engineers have puzzled over the most economical way to orient HP-Piles: Should the strong axis be oriented in the longitudinal direction or in the transverse direction? The accurate design of foundations may considerably reduce construction costs. In order to answer this question, 3-D finite element models will be set up to analyze the effects of pile orientation on bridge behavior. Because seismic loading usually controls the structural design in seismic zones, it is important to investigate the effects of pile orientation under such loading. HP-Piles are designed conservatively due to the uncertainty of the dynamic soil-structure interaction during earthquakes. For the bridge substructure, the soil?s stiffness plays an important role in load distribution as the soil and piles act as a combined system to resist the seismic loading. The dynamic soil-structure interaction mechanisms of bridges in both the longitudinal and the transverse directions will be investigated through the use of 3-D finite element models. The numerical models will realistically simulate the soil-pile-abutment-wingwall-superstructure system subjected to seismic loads. Nonlinear soil springs will be employed to model the non-elastic soil behavior. The analysis will provide theoretical support to the design practice. KEYWORDS: HP driven piles, finite element, pile orientation, soil-structure interaction

By Andrew Rietz, Yu Bao, Steven Halewski

"HP-Pile foundations are the most favorable type of foundations for bridges because of their constructability, versatility with regard to subsoil conditions, and settlement controllability. For many years, bridge engineers have puzzled over the most economical way to orient HP-Piles: Should the strong axis be oriented in the longitudinal direction or in the transverse direction? The accurate design of foundations may considerably reduce construction costs. In order to answer this question, 3-D finite element models will be set up to analyze the effects of pile orientation on bridge behavior. Because seismic loading usually controls the structural design in seismic zones, it is important to investigate the effects of pile orientation under such loading. HP-Piles are designed conservatively due to the uncertainty of the dynamic soil-structure interaction during earthquakes. For the bridge substructure, the soil’s stiffness plays an important role in load distribution as the soil and piles act as a combined system to resist the seismic loading. The dynamic soil-structure interaction mechanisms of bridges in both the longitudinal and the transverse directions will be investigated through the use of 3-D finite element models. The numerical models will realistically simulate the soil-pileabutment- wingwall-superstructure system subjected to seismic loads. Nonlinear soil springs will be employed to model the non-elastic soil behavior. The analysis will provide theoretical support to the design practice.IntroductionAs sustainability has asserted itself in the public’s eye, bridge professionals have rushed to develop criteria for sustainable bridges. The three main categories that define and distinguish sustainable bridges are: lower energy input, increased durability, and simplified deconstruction. By optimizing the structural members, the engineer is able to reduce the material quantities for the bridge, effectively reducing the energy input for the bridge (Hunt; Whittemore). A deeper understanding of soilfoundation- structure mechanism will facilitate sustainable bridge design by lowering energy input and increasing the service life of bridges (Rietz et al.)."

Jan 1, 2017

By Y. Fautrelle, F. Felten, O. Metais, Y. Du Terrail

"We deal with the prediction of electromagnetically-driven turbulent flows by means of a Large-Eddy- Simulation method (LES). The model is applied in the case of a liquid metal pool submitted to a polyphase linear electromagnetic stirrer. We investigate two cases: (i) the frequency of the magnetic field is equal to 50 Hz; in that case the effects of the pulsating part of the Lorentz forces may be neglected; (ii) the frequency of the magnetic field is low (f= 1 Hz), then the oscillating part of the electromagnetic forces is taken into account. The LES predictions agree well with the mean velocity measurements, as does the standard k-a model. However, as for the turbulent kinetic energy predictions, there is a large discrepancy between the two models. When the oscillating part of the Lorentz forces is taken into account, the computations show that the fluid flow is sensitive to the unsteady part of the forces. We find that the mean velocity is not affected by the fluctuating component of the force. As for the turbulence parameters, the presence of the pulsating part leads to a significant reduction of the turbulent kinetic energy, whilst the turbulence length scale decreases.IntroductionElectromagnetic stirring is commonly used in various metallurgical processes. In continuous casting for example, it may improve the structure of the solidified metal products. Fluid flows generated by electromagnetic stirrers have been extensively investigated (1 ). As for numerical predictions, many works have been published on the subject. In most of the previous works, onepoint closure models have been used. The most used model is the well-known k-i:: model. Those models yield fairly good results as for the mean velocity prediction. However, it has two main drawbacks:(i) the turbulence is not predicted accurately,(ii) it is not well fitted for unstationary flows as well as transitional turbulent motions, where theReynolds number is moderate."

Jan 1, 1999

By PDF

Flotation columns are widely used in mineral processing and coal preparation industry today. Mathematical modeling of flotation columns has conventionally been done using axial dispersion or tank-in-series models. However, for design and scale up purposes these models are not adequate enough because they do not provide detailed information of the flow fields within the column. This paper presents a Two-Dimensional, Two-Phase. Fluid-dynamic Model and its numerical solutions for some cases encountered in column flotation. The MAC (Marker and Cell) method was employed to obtain the numerical solutions, and based on this method a computer code COLUMNFLOW was developed. The simulated flow field reveals the presence of a circulation flow structure in flotation columns, which could be one of the main reasons for mixing in column flotation. The effects of some parameters such as superficial gas and liquid velocities on the characteristics of the flow structure in flotation columns were also tested with the numerical model.

Jan 1, 1994

By Amalia Giannakou

When designing structures in seismically active areas, foundations of critical structures are typically located away from known faults. Indeed many seismic codes prohibit construction in the immediate vicinity of seismically active faults establishing minimum setbacks from the faults. However, for long structures such as bridges, tunnels and pipelines, a fault maybe unavoidable, and fault rupture risk impossible to preclude. This paper presents a case study of a bridge foundation design that incorporates the consequences of fault deformations passing through the foundation. A series of site surveys were conducted to characterize the fault environment, and fault displacement hazard analyses were conducted to develop design displacement demands. Subsequently, numerical analyses were performed to evaluate the fault rupture-induced demands on a pier foundation in terms of displacements and rotations. Initially, analyses were performed that model the free field fault rupture propagation to the ground surface of the fault to ensure that the basic site characteristics are captured. Subsequently, the foundation elements were introduced in the model and the fault rupture is simulated again. The interaction and influence of the foundation on the fault rupture propagation is calculated, and the consequent demands and distortions of the foundation block provide a basis for structural engineers to design the bridge.