Search Documents

By F. Kavicka

An original application of ANSYS simulated the forming of the temperature field of a massive casting from ductile cast-iron during the application various methods of its cooling using chills. The numerical model managed to optimize more than one method of cooling but, in addition to that, provided serious results for the successive model of structural and chemical heterogeneity, and so it also contributes to influencing the pouring structure. The file containing the acquired results from both models brings new and, simultaneously, remarkable findings of causal relationships between the structural and chemical heterogeneity and the local solidification time in any point of the casting. The determined relations therefore enable the prediction of the face density of the spheroids of graphite in dependence on the local solidification time. The calculated temperature field of 500x500x1000 mm casting of ductile cast-iron with various methods of cooling has successfully been compared with temperatures obtained experimentally. Acknowledgments. This analysis was conducted using a program devised within the framework of the GA CR projects (No.106/04/1334, No. 106/04/1006 and No. 106/04/094), of the COSTAPOMAT-OC526.10, MPO/FI-IM/021, EUREKA No.2716 COOP and KONTAKT 14/ 2005/06.

Jan 1, 2006

By Axel Studeny

Geomechanical-numerical methods are an established tool for planning of underground excavations worldwide. By describing the interaction between heterogeneous layered strata and support elements it is possible to determine their suitability. The German hard coal mining industry has more than twenty years experience in investigation and using these parameters to carry out numerical calculations. Today, nearly all excavations (shafts, pit bottoms areas, bunkers, coal faces, development drivages, gate roads etc.) are calculated by geomechanical-numerical tools in combination with empirical planning methods. A high accuracy of planning results, which is a special requirement in the German coal industry with its deep-level longwalls, bedded and weak surrounding strata and multiple seam mining layouts, depends on the extensive calibration of the numerical models, which is achieved by: 1) Taking a large number of underground measurements with empirical evaluation. 2) Carrying out physical modeling. 3) Using the characteristic features of the installed support elements. This paper compares the planning approach used for multiple and single seam mining and presents the advantages of this planning tool. The main focus is on the feasibility of rockbolting as the sole means of roadway support for multiple seam mining operations, though the paper will also look at how the costs can be improved by using a low bolting density in single seam layouts.

Jan 1, 2007

By Giovanni B. Barla, Stefan H. Boshkov

The block caving method is examined in this paper on the basis of experimental results and observations in the field, and through the use of numerical modeling by the Finite Element Method. The Grace Mine, developed by a panel caving method, served as a reference case. However, the results obtained and the simulation techniques which have been developed may find useful application in more general conditions as encountered in other mining operations by block caving. Following a brief description of the Grace Mine and of the mining method, the problems associated with the behavior of underground drifts at the production levels, during undercutting and caving, are addressed. The most important results obtained through monitoring of strains and loads in the supports of a runway are discussed, while giving relevance to the sequence of mining operations in the immediate vicinity. The most important geomechanical data, as provided for subsequent modeling, are described with reference to: (a) rock mass characterization and (b) in situ state of stress. Then, the numerical modeling procedures, developed by the Finite Element Method in order to predict the caving phenomena and the behavior of underground openings during undercutting, are discussed. Finally, an attempt is made to compare the pre- dictions of the numerical analyses with, the measurements carried out at the instrumented site. Suggestions as to how to improve the results both qualitatively and quantitatively are also given.

Jan 1, 1984

By A. Pellegrino, F. Collin, C. Guiducci, R. Charlier

A numerical approach for modeling coupled Hydro-Mechanical fracture behavior in 2D is presented in this paper. This approach accounts for the strong dependence of the flow characteristics on fracture apertures and offers an anisotropic description of the flow propagation in the fracture. An interface finite element and a constitutive model are developed . In order to validate the applicability of the model to general geomechanic problems, 2D simulations were performed. Results show the good application of the mathematical model and the fundamental hydraulic role played by the fractures. In diesem Beitrag wird ein Stoffgesetz vorgeschlagen um hydromechanisches Kluftverhalten numerisch zu Modellieren. Dieser Vorgang nimmt die starke Abhängigkeit der Fliesseigenschaften von der Kluftöffnung mit in Anbetracht. Ein Interfaceelement und ein Stoffgesetz wurden entwickelt. Zum validieren der Anwendbarkeit dieses Modells wurden 2D Simulationen durchgeführt. Ergebnisse haben gezeigt dass das mathematische Model anwendbar ist und dass Klüfte eine fundamentale hydraulische Rolle spielen. Une approche numérique pour la modélisation du comportement couplé hydro-mécanique des fractures en 2D est présentée dans ce papier. Cette approche prend en compte la forte dépendance des caractéristiques du flux de fluide sur les ouvertures des fractures et présente une description anisotrope de la propagation du flux dans la fracture. Un élément fini d'interface et un modèle constitutif sont développés. Dans le but de valider l'applicabilité du modèle à des problèmes généraux de géomécanique, des simulations 2D ont été réalisées. Les résultats montrent la bonne application du modèle mathématique et le rôle hydraulique fondamental joué par les fractures.

Jan 1, 2003

By Shahé Shnorhokian

A study was conducted on the seismic activities and associated rockburst events that took place at an underground Canadian metal mine over a period of five years. In the first stage, a back analysis of the microseismic database led to the identification of clusters of major events. The locations of these clusters were identified with respect to the local geology and the planned mining sequence. In the second stage, a numerical model was constructed using the finite difference code FLAC3D to simulate the in-situ stress, mine geometry, local geology, and stope sequence. Mining-induced stress fields were compared against two rockburst events that took place during this period, as well as against the various geologic units present. It was shown that by combining microseismic analysis and numerical modeling, it was possible to provide a basis for evaluating the influence of the in-situ stress regimes and stope sequencing on mining-induced seismic response.

May 1, 2011

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 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 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 R. Karl Zipf, Francois E. Heuze

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

By R. Karl Zipf

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

Jan 1, 2003

By R. K. Zipf

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

By A. P. E. Dirige

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

May 1, 2009

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

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

Jan 1, 1991

By A. J. Basu

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

Jan 1, 1996

By O. K. Mahabadi

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

May 1, 2009

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

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

Jan 1, 2017

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

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

Jan 1, 2003

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

Nov 25, 2010

By Kubik K

The objective of this paper was to determine if recent advances in numerical electromagnetic wave propagation modelling techniques could be used to study the propagation of GPS signals over arbitrary terrain. Specifically the effects of diffraction and multipath were to be modelled accurately for propagation over arbitrary terrain surfaces defined by a digital terrain model. Ideally the model was to be in three dimensions and cover an area of a few square kilometres with a simulation height of a few hundred metres. Simulations were to run on current technology PC's with run times less than a few hours. A relatively new technique that involves a numerical solution to the Parabolic Equations can be used to solve for two-dimensional propagation over any type of terrain. The PE provides a direct solution of Maxwell's wave equations by approximating the Helmholtz scalar wave equation. This technique does not rely on the study of individual ray paths as used in the Geometrical Theory of Diffraction (GTD). It was decided to develop a two dimensional Parabolic Equation code, using the Fourier Split/Step method that allowed propagation over arbitrary terrain (; ; ). The model boundary conditions assumed that the satellite transmission was a plane wave, and that the lower boundary was a perfect reflector. Propagation simulations from the model accurately provide the amplitude and phase of the propagated plane wave at all points within the model domain. The PE model was used to study the effects of diffraction and multipath caused by various types of terrain commonly found in an open cut mining environment. The study of diffraction errors was relatively simple, because we are simply concerned with amplitude of the signal as it diffracts about an obstacle. If the amplitude is above the level receivable by a GPS receiver then the diffraction error can be obtained from the geometry of the terrain obstacle. The study of multipath is more complex, as the computed phase from the PE model is ambiguous and the presence of multipath can only be inferred. Thus the absolute magnitude of the nmltipath error is not currently determinable, though the results of some new work is shown here. This paper shows that diffraction effects can cause significant errors in range. It also shows the effects that multipath has on the received amplitude of the signal and the received phase of the signal. A discussion is made of further work that is presently underway that will allow the direct determination of the multipath error by comparing the ambiguous phase with a reference phase.

Jan 1, 1997