Search Documents

By Paul Hughes

This paper describes a focus of work presently being conducted at the Rock Mechanics Research Group at the University of British Columbia. The underhand method under consolidated fill ensures a high recovery under an engineered back that is comprised of cemented rock fill and/or cemented paste fill. This method of mining has been employed successfully in mines throughout North America as a method of mitigating exposure to the operator having to work under a seismically active ?back?. This paper reviews design methodology in the placement and analysis of sill mats with reference to site observation coupled with on going numerical and analytically derived solutions.

Jan 1, 2006

By Ehsan Khajeh

In this study, the permeabilities of Al-15.5wt%Cu and Al-19.5wt%Cu have been determined by solving the full Navier-Stokes equations on 3D geometries of interdendritic channels obtained by X-ray microtomography (XMT). The samples for XMT were sectioned from different positions in directionally solidified cylinders in order to obtain a variety of microstructures for each composition. The flow has been solved using a 2nd-order accurate Finite Volume Method (FVM) approach. The effect of employing uniform or unstructured mesh on the calculated permeability has been studied. The Marching-cube triangulation method was necessary to extract accurate surfaces for unstructured volume meshing. Calculated values of permeability for this range of solid fractions show partial agreement with previous experiments. Observed deviations are attributed to the difficulties of experimental permeability measurements.

Jan 1, 2009

By S. Finsterle, R. C. Trautz

Numerical modeling studies were performed to investigate water seepage into underground openings excavated in unsaturated fractured rock. Water seepage is an important factor affecting the performance of a nuclear waste repository, such as the one proposed at Yucca Mountain, NV. The amount of water dripping into an underground opening depends on the connectivity and permeability of the fracture network, as well as on the capillarity of individual fractures intersecting the opening. A high-resolution numerical model of an unsaturated fracture network was developed, and the appropriateness of seepage predictions using a fracture-continuum model was examined. Model calibrations to liquid-release tests performed at Yucca Mountain are also discussed.

Jan 1, 2002

By S. Finsterle, R. C. Trautz

We performed numerical modeling studies to investigate water seepage into underground openings excavated in unsaturated fractured rock. Water seepage is a key factor affecting the performance of a nuclear waste repository such as the one proposed at Yucca Mountain, Nevada. The amount of water dripping into an underground opening depends on the connectivity and permeability of the fracture network as well as on the capillarity of individual fractures intersecting the opening. We developed a high-resolution numerical model of an unsaturated fracture network and examined the appropriateness of seepage predictions using a fracture-continuum model. Model calibrations to liquid-release tests performed at Yucca Mountain are also discussed.

Jan 1, 2000

By S. Stefanizzi

In the present paper, numerical simulations based on a finite/discrete element approach are used to study deformation and failure process during standard rock mechanics. FEM/DEM modeling has been used as available with the ELFEN code, which explicitly model the transition from continuum to discontinuum. The transition is simulated through a crack nucleation and propagation that obeys to the Griffith?s failure criterion. If the fracture criterion within the intact rock (represented by FEM) is met, a new crack (represented by DEM) is initiated. Remeshing allows the fracture process through the FEM mesh to be tracked and visualized, thus contact properties can be assigned to pre-existing fractures and newly generated fractures. In this study, compressive tests of rock samples are simulated. The applicability of ELFEN to properly and efficiently simulate typical laboratory rock mechanics tests, with reference to the rock fracture process, is described in detail. In addition, an indirect validation of the extensional strain criterion proposed by Stacey (1981) is given.

May 1, 2009

By Benjamin Mirabile

"Current and emerging practices to ventilate longwall gobs are prompting increased design and analysis of standing support requirements related to longwall mining. Specifically, additional design methodologies are being employed to evaluate standing supports in the #2 entry of a 3-entry longwall system.Numerical modeling will be used to evaluate the performance of standing support in the #2 entry between gobs of longwalls in the Herrin No. 6 seam. The performance of standing support in the #2 entry will be evaluated based on the following: (1) The effectiveness of standing support in maintaining an open airway in the #2 entry; (2) The interaction of standing support in the #2 entry with standing support in the longwall tailgate; and (3) Loading conditions of longwall pillars. Initial numerical modeling activities are focused on evaluating the roof to floor convergence in the #2 Entry and using the ground reaction curve concept for support design.Numerical modeling indicates that, like standing support in tailgates, supports in the #2 entry require some type of yield mechanism to withstand initial uncontrolled convergence. Numerical modeling also suggests that support density can be lower than that used in a tailgate support application. The numerical models used for initial investigation of support performance in the #2 entry are limited in their ability to evaluate support performance. The magnitude of uncontrollable convergence can be estimated with LaModel software. Near-seam conditions have not been investigated in this research since the LaModel software does not have the ability to evaluate inelastic rock behavior.IntroductionDuring the past three decades, many researchers have investigated standing support performance in longwall tailgate entries. Very little attention was paid to the middle entry of a threeentry longwall gateroad. As longwalls retreated, no personnel would travel the #2 entry, and the entry was considered part of the longwall gob. In response to recent mine disasters, regulatory agencies and operators have devoted more attention to these entries. Prior to the Upper Big Branch mine disaster, little or no secondary support was applied in these entries. Unlike secondary support in longwall tailgates, there is little empirical knowledge or quantified research specifically applied to secondary support in the #2 entry. This paper outlines initial steps used to develop a design methodology for secondary support in the #2 entry."

Jan 1, 2015

By Sylvain Gouttebroze, Mohammed M'Hamdi

In a previous work, a small scale Bridgman furnace has been used to study silicon bi-grain crystallization at different cooling rate. This work expands the analysis by studying the mechanical interaction between the crucible and ingot during the solidification and cooling. The thermal model is based on a 2D-axisymetric heat-transfer model. The flux histories are then transferred to the ingot-crucible 3D-model. Anisotropic Elastic and Crystal Plasticity model are used to model the silicon deformation. Four different assumptions are applied to model the mechanical contact crucible-ingot and three grain misorientations are considered. The results show the strong impact of the alumina crucible contraction on the stresses and deformations in the silicon ingot.

Jan 1, 2015

By Ximeng Li

The entries that are subjected to dynamic pressure are much more difficult to maintain due to repeated disturbance. A case study of an entry stress distribution was performed using the FLAC3D program. By analyzing the distribution of stress field around the entries, this paper provides guidelines for entry support design in the lower (#16) seam subjected to the repeated disturbance from longwall mining in the immediate overlying (#17) seam, 72ft (22 m) above. The results showed that there was stress superposition inside the barrier pillar in the lower seam, 15ft (5m) to 45ft (15m) deep, during development. The influence zone of the front abutment pressure in the upper (#16) seam was about 100ft (30m) outby the faceline. The stress-relieved zone was from 0 to 100ft (30m) inby the face. The area further inby the face (i.e., 230 ft (70 m) inby the face) was within the gob compaction zone. Along the face advancing direction, the affected zones of the entries in the lower (#17) seam that were subjected to the dynamic pressure from longwall mining of the immediate overlying seam could be divided into four stages: pillar development stress stage, front abutment pressure stage, unstable gob stage, and stable gob stage. Those entries within an area from 100ft (30m) outby to 100ft (30m) to 230ft (70m) inby the faceline should be reinforced with special supports. The numerical modelling matched the monitoring results well.

Jan 1, 2013

By Young-Dong Cho, Yonghun Jong, Dong-See Shim, Seokwon Jeon, Chung-In Lee

The most important operation in the tunnel blasting procedure is to create an opening in the face in order to develop more free surface in the rock mass. This operation is highly dependent on the cut holes. In general, the burn-cut or V-cut is commonly adopted to blast the cut section. But, these traditional methods often raise problems. There are instances (a) that the actual advance is considerably less than the designed advance depending on the rock mass, (b) that they cannot sometimes enlarge the free surface efficiently and (c) that the blast vibration resulted from the cut blast exceeds the allowable vibration level. In order to resolve the problems, many different kinds of the cut blast can be considered.

Jan 1, 2005

By R. L. McNearny, K. A. Barker

Four numerical models written with the "PFC2D " program were compared against four similar physical models. The numerical models were constrained numerically to closely follow the caving behavior of the physical models. Results from the numerical analyses closely simulated both the granular and blocky behaviors of the physical models, including bridging and interlocking of the simulated orebody during draw. The drawdown patterns and the rate of draw developed during numerical modeling also closely simulated the physical models. The same assumptions and governing parameters used in analyzing the first four models were applied to an additional model. The fifth model employed a completely different, uneven, draw sequence than the other four models, allowing the material to flow in a more granular manner with a demonstrated instability of the undercut sections between each drawpoint opening. Drawdown patterns and behavior of the numerical material of the first four models progressed at uneven and inconsistent rates, similar to those seen in the physical models. By inspection, the overall shape and flow lines of both the numerical and physical models were close in the area removed and in the flow characteristics. The fifth model proved to be the most unstable of all the models tested, with poor correlation to the physical models. The uneven draw sequence, coupled with an additional frictionless wall condition and an added drawpoint, caused the material to flow in a very granular manner. Information generated by this study can be used as a first approach in determining the proper draw rate, drawpoint spacing and undercut thickness for a block-caving operation. Numerical models are often difficult to predict and interpret. However, this study should increase the confidence level for use of the PFC2D program in the future.

Jan 1, 2000

By A. S. Tawadrous

The stress field around a detonating charge in a borehole was calculated using AUTODYN for a variety of in-situ stress conditions. Calculations were conducted for zero in-situ stress, hydrostatic conditions at 100MPa, uniaxial stress conditions of 100MPa, and bi-axial state of stress with amplitudes of 200MPa and 100MPa along the axes. Two different explosives, ANFO and emulsion, were used for the calculations. It was shown that in-situ stress fields had an effect on the tensile and compressional damage zones around the blasthole, calculated on the basis of the maximum amplitude of the stress pulses in the rock after the detonation of the explosive in the blasthole. The extent of compressional damage increased around the borehole, whereas, the extent of the tensile damage decreased. The damage zones are of elliptical nature with their axes parallel to the direction of the principal stresses. The difference between the major and minor axis of the elliptical damage was not very significant in the case of the compressional damage; however it was significant in the case of the tensile damage.

Jan 1, 2007

By Khaled Morsy

Longwall mining is one of the most widely practiced underground coal mining method. The behavior of the longwall gob is very critical in the understanding of the complex ground response to longwall mining. In the limited data available, the values of gob moduli used in numerical modeling ranged from 1,000 psi to over 300,000 psi. Such a wide variation in moduli would greatly affect the results of the numerical analysis. A better understanding of the behavior of the gob will allow numerical models to be more accurate for simulating longwall mining conditions. In this study, a gob model based on the Terzaghi's model was incorporated into the comprehensive finite element package, ABAQUS, to simulate the loading behavior of longwali gob material. A case study was used to verify the proposed gob model.

Jan 1, 2002

By Carlos Schvezov, Eliana Agaliotis, Alicia E. Ares, Mario Rosenberger

"The interaction between a solidification front and a spherical particle was modeled and simulated in order to study the phenomenon of ""pushing"". Finite element methods were employed in axisymetric domains. The drag force is calculated from the calculated fluid flow field and the repulsion force is the Lifshitz-Van der Waal force. The thermal field is decoupled from the force field. The model is applied to a metallic matrix containing particles with similar, lower or higher thermal conductivities than the matrix, resulting in planar, convex and concave interface shapes, respectively.The results show that the critical velocities for pushing for concave and convex interfaces are about one order of magnitude larger and one order of magnitude lower than the critical velocity for flat interfaces, respectively. The simulated results were in very good agreement with reported experimental results.IntroductionThe solidification interface of materials containing spherical particles in some cases can push the particle ahead the interface. This phenomenon generally occurs during the solidification of different materials such as, metallic alloys containing native particles, metal matrix composites containing particles or fibers as reinforcements and crystallization of semiconductor, optical, optoelectronic or biological materials. In all cases the interaction strongly determines the distribution of particles and the resulting physical and physic-chemical properties of the resulting material. The phenomenon has been investigated for many years both experimentally and by mathematical modeling [1-21]. The process is complex due to the variety of phenomena involved during pushing. In all cases, the physical problem depends on the properties, nature, and morphology of the interacting media; the particle, melt and solid; and external fields such as gravity, thermal and electromagnetic fields."

Jan 1, 2012

By Andre Zingano

It is known that the backfilling confinement increases the pillar strength and also reduces the rib lateral displacement. With backfilling, it is possible to increase the room-and-pillar mining recovery. This paper aims to study the behavior of coal pillars confinement using rock fill as backfilling material. Laboratory tests were carried out to determine the behavior of backfilling materials. Constitutive models were developed for the tested backfilling materials. Numerical simulations for the effect of backfilling materials on the ground stability in room-and-pillar coal mine entry were performed. The backfill method is a good strategy to (i) improve the ground stability in room-and-pillar coal mines for long term, especially for weak coal seams, (ii) increase the coal mining recovery and (iii) reduce the quantity of waste material on waste piles at the surface.

Jan 1, 2011

By S. Guo

In this paper, a model for the resin anchor bolt system is developed in which the smooth bar is modeled by using one-dimensional truss elements and the resin annulus is modeled by using axisymmetric elements. These two submodels are combined together and then incorporated into the' existing computer code ROOFBOLT. After this development, all of the components in the resin anchor bolt system including smooth bar, bearing plate and resin annulus can be correctly modeled.

Jan 1, 1995

By J. C. Stankus, S. Guo, S. S. Peng

In this paper, a model is developed for a resin anchor bolt system in which the smooth bar is modeled by using one-dimensional truss elements, and the resin annulus is modeled by using axisymmetric elements. These two submodels are combined and then incorporated into an existing computer code called Roofbolt. After this development, all of the components in the resin anchor bolt system, including smooth bar, bearing plate and resin annulus, can be correctly modeled.

Jan 1, 1997

By Riqing Ye

Laser Engineered Net-shaping (LENS TM) is an emerging manufacturing technique that ensures significant reduction of process time between initial design and final components. The fabrication of fully dense parts with appropriate properties using LENS TM process requires an in-depth understanding of the entire thermal behavior of process. In this paper, the thermal behavior during LENSTM processing of 316 stainless steel was studied both numerically and experimentally. Temperature distribution and gradient in the fabricated 316 stainless steel part were obtained by Finite Element Method (FEM) simulation. The numerical results are in good agreement with the experimental observations. The numerical method may be used to optirnize process parameters and predict the thermal response of LENSTM fabricated components. Key words: Laser Engineered Net-shaping, Thermal Behavior, Finite Element Analysis

Jan 1, 2003

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 Matthew E. Fowler, Jean Habimana, Roy F. Cook

The West Side CSO Tunnel, Pump Station & Pipeline Project is currently beginning construction along the downtown waterfront of Portland, Oregon. This4.3 meter diameter, 5,500 meter-long soft ground tunnel will be driven along the westbank of the Willamette River adjacent to buildings and utilities and several major pile-supported bridges. High groundwater pressure, extended lengths of gravel to cobble-size alluvium, and the need to control ground loss and surface settlements were key considerations in the selection of slurry tunnel boring machines for the project This paper presents an overview of the numerical modeling that was performed during design to estimate the strains and additional loads in the pile-supported structures and to assess the efficiency of proposed mitigation measures.

Jan 1, 2003