Search Documents

By M. E. Duncan Fama

The aim of the paper was to model numerically some aspects of the extraction of coal pillars by split and fendering. A two-dimensional elasto-plastic plane strain model with strain softening and interface modelling was used. The model was able to simulate the improvement of conditions experienced in lifting from the fender, when the roof adjacent to the fender has caved. In contrast, when a stook was left, providing support to a span of uncaved roof, the model predicted the fender would suffer considerably more deformation and strain.

Jan 1, 1992

By Ian Masters, Saul B. Davies, David T. Gethin

The purpose of this model is to simulate and analyze the energy flows within a rotary aluminum recycling furnace, using the global properties of its main components, the gas, metal and refractory walls. The model will simulate a series of operations such as loading, heating and stirring which will constitute the processing of a batch. Each main component of the furnace is treated as a one-dimensional heat transfer medium, the thickness and area of each medium being calculated from the actual geometry of the furnace. The physical properties are taken as the weighted average of those for the various layers forming the medium. Incorporated into the model are elements describing the radiative exchange within the furnace, melting of the charge, and forced convection in the melt due to the rotation of the furnace. This approach has produced a model with a reasonable computing time that can be used to run a whole series of simulations for the prediction and comparison of the furnace performance under different operating conditions.

Jan 1, 2000

By K Koskinen, S Ronkainen

The flow field induced in an Outokumpu TankCell and by an OK rotor-stator mixing device in a laboratory and a process size flotation cell was modelled numerically using computational fluid dynamics (CFD). The primary target was to discover and reproduce the various interactions between the flow field and the flotation cell by numerical means. As power consumption depends sensitively on flow field, it can be regarded as a top-level derived quantity (or measure/metric) indicating the validity of the numerically obtained flow field. Having properly validated numerical means to produce the flow field the influence of various design parameters can be captured with certain accuracy. In addition, separate gas and solid phase models for bubbles and solids, respectively, were implemented. The comparisons between velocities and turbulence quantities measured using the laser doppler velocimetry (LDV) technique and CFD modelling were performed. The grid generation was performed with the commercial Gambit 1.3 and the subsequent integration and solving of the resulting algebraic equations with the commercial CFD-code Fluent 6.22. The computationally economic multiple reference frame method was included to model the moving rotor blades. Gas and solid phase models were modified mixture models.

Jan 1, 2005

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

A mathematical model for spontaneous combustion of coal with moisture included is presented in this paper. The one-dimensional unsteady state model consists of conservation equations for oxygen, water vapour and inherent moisture of coal and energy for both gaseous and solid phases. A first order Arrhenius reaction rate for oxidation under both pore diffusion and chemically controlled reaction regime is considered. The rate of evaporation or condensation is also considered as a function of temperature of coal, water content of coal and gas streams. The equation of the model is solved numerically by the finite difference technique. Influences of different parameters on the process of spontaneous combustion can be examined by using this model. Also the model can be used to simulate full-scale storage conditions.

Jan 1, 2005

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 Mohamed Sassi, I. Kolo, R. L. Sousa, R. K. Abu Al-Rub, Manhal Sirat

"Fractures induced by folds are pivotal in hydrocarbon exploration, ground water transport and harnessing of geothermal energy. This is because of the need to predict and understand fracture propagation in reservoirs which have folded rock formations. Despite its prominence, there is a lack of reliable numerical models for simulating fractures resulting from rock folding. This is partially due to the difficulty involved in fracture modelling. As a contribution to fold-fracture modelling, this work applies a coupled plasticity-damage model to rock fracturing and anticlinal folds. Parametric studies on a single folded layer give insight into the fold formation and behaviour. The anisotropic continuum damage model which considers both tension and compression is formulated using the power damage evolution law. For ease in formulation, strain equivalence hypothesis is adopted whereby the strain is the same for both damaged and undamaged configurations. Lubliner plasticity yield criterion is adopted for plastic deformation. The model is coded in Abaqus user subroutine UMAT and is applicable to quasi-brittle materials.INTRODUCTIONThe understanding of fractures in the earth’s crust helps in the prediction, evaluation and characterization of fractured reservoirs, whether they are petroleum, geothermal or underground water reservoirs. This is partly because reservoir permeability and stability are dependent on fracture properties (Nelson, 2001). To achieve this, various approaches have been developed to analyze reservoir rock fractures. These could be grouped into (a) discrete approaches – where the rock mass is represented by a finite number of well-defined components and (b) continuum approaches – where the mathematical assumption of an infinitesimal element is made (Jing, 2003). Ivanova (Ivanova, 1998) comprehensively studied fracture systems in nature and highlighted folds as one of the major geologic settings that produce fractures. Naturally occurring folds in reservoirs sometimes compound fracture analysis. It might not be clear whether the fractures formed before, during or after rock folding. Fractures formed during the folding process are expected to show stress orientations related to folding and are thus of primary concern in rock folding simulation. It could be argued that there are no numerical models that as"

Jan 1, 2015

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

By D. J. Reddish

The paper describes a case study where numerical modelling has been applied to predict potential risks associated with inflows of water in active mine workings from an overlying sandstone aquifer. The objective of the modelling was to predict the development of strata fractures and strains as the longwall panels were extracted and investigate their relationship to an sandstone aquifer. The numerical modelling study also included the investigation of the effect of mining on a fault that existed close to the start position of one of the panels. The potential effects of mining on the fault were undertaken by parametric studies adjusting properties for the fault and the relative position of extractions, which allowed the mine planner to position the panel layout accordingly.

May 1, 2003

By W. Keilich

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

Jan 1, 2006

By C. V. Deutsch, O. Babak, J. G. Manchuk

High levels of geological heterogeneity occurring over relatively short distances pose a challenge to the unconventional open-pit techniques used in oil sands mining. There is strong motivation to mine selectively and minimize dilution and loss; however, the large equipment and high mining rate used to meet production targets tend to reduce selectivity. A numerical model is developed to account for geological complexity, mining equipment, regulatory criteria, and other operational parameters. This framework is used for more accurate forecasts, optimization of equipment, and other operational decisions. An example demonstrates the approach in light of Alberta Energy Regulator criteria.

Jan 1, 2014

This paper discusses support requirements for the proposed monorail system to be used in decline development. The monorail drilling and loading systems are systems that move on the rail (monorail) installed in the roof of the decline and supported by roof bolts, suspension chains and steel supports. However, due to the weight of the components of the two systems, it is imperative that the force in each roof bolt, suspension chain and steel support capable of suspending the weight of the heaviest component is determined. Numerical models that relate the weight of the monorail drilling and loading components to the required strength in the support system have been developed. Using these developed models, numerical values of the forces in each roof bolt, suspension chain and steel support, required to suspend the weight of the heaviest component of the monorail drilling and loading systems are determined.

Dec 6, 2010

By S. P. Lin

The numerical models of nonlinear piles for deep foundations are studied in this paper. The purpose of these is to clarify the mechanical behaviour of deep foundation nonlinear piles subjected to axial and lateral loadings. Based on the site data of p-y curve as well as t-z curve of foundation soils, it is possible to evaluate pile performance under above mentioned loads by means of combined FEM/FDM models. Numerical examples are presented in tall building and offshore jacket for illustrating the technique, it show high efficiency of these models. Some information to be interested by design engineers are given, especially for those who pay very concerns on the interaction behaviours of structure and its pile foundation.

Jan 1, 2010

By T. Saksala

In this paper, the embedded discontinuity approach is applied to finite element modeling of rock fracture. A rate-dependent constitutive model based on the embedded displacement discontinuity theory is developed to describe the mode I, mode II and mixed mode fracture of rock in tension and compression. The bulk material is described as linear elastic until reaching the elastic limit. Beyond this limit, a ratedependent exponential softening law governs the evolution of the embedded displacement jump. The present approach incorporates the rock heterogeneity by random description of the mineral texture of rock. Moreover, the inherent initial microcrack populations of natural rocks are accounted for as randomly oriented embedded discontinuities. Numerical examples demonstrate the model behavior in uniaxial compression and tension. The effect of loading rate and confining pressure is tested as well in 2D numerical simulations. These simulations show that the model captures the main features of rock in confined compression and uniaxial tension. The developed method has the computational efficiency of continuum plasticity models. However, it has an important advantage of accounting for the orientation of introduced microcracks.

Jan 1, 2015

By Andrzej Walentek

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

Jan 1, 2009

By Andrzej Walentek

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

Jan 1, 2009

By J. A. Jankowski

The situation in the raw material market has brought into consideration the commercial exploitation of the mineral deposits on the deep sea bed and in particular the manganese nodules. They occur generally in depths of 4000-5000 m, partially covered or buried in bottom sediments. The deposits and their potential recovery have been intensively explored and researched during the past twenty years (Bischoff 1979, Halbach 1988, Padan 1990). Although the present market situation limits the rentability of deep sea mining on an industrial scale, the mining technology is being further developed in many countries. The stagnation phase in the deep mining issue allows sufficient time for conducting research on the precautionary protection of the deep sea environment and for consultation in the development of technologies with the least possible impact on it. The subject has gained in importance since autumn 1994, when the United Nations International Law of the Sea was signed. This also regulates human activities in international waters from the environmental point of view. As fields tests have shown in the past, even with the most modern technologies, mining of the manganese nodules will inevitably be accompanied by a resuspension of some bottom sediments. The generated sediment plumes will be further transported with the ocean currents and in the worst case form stable turbidity layers in different depths (Ozturgut 1981a, Baturin 1991). One of the problems of the environmental impact assessment is the estimation of the plume persistence before it eventually dilutes to approximately the ambient concentration level or settles at the bottom (Lavelle 1981a, Lavelle 1981b, Lavelle 1982, Thiel 1991a, Thiel 1991b). The aim of the research is to develop a numerical sediment transport model in connection with potential manganese nodule mining in a reference area located in the southeastern equatorial Pacific Ocean off Peru (Peru Basin, Discol Experimental Area). In this region, a German mining claim is registered and interdisciplinary field works of the TUSCH (German: TiefseeUmweltSCHutz, deep sea environmental protection) research group are carried out in order to obtain the required scientific background for evaluating of the deep sea mining impacts and to provide the necessary supporting data for the modelling (Thiel 1991a, Thiel 1991h).

Jan 1, 1995

By J. H. Curran

In this paper, a numerical model for determining the stress and displacement fields around a crack of arbitrary geometry in a linear elastic medium is described. Accuracy and computational speed were the main considerations in selecting the numerical method and Propagation criteria used to model the growing crack. Some results regarding the stresses and displacements around a two-dimensional crack of arbitrary geometry are discussed.

Jan 1, 1984

Following a short review of selected literature, a back analysis of stress measurements of stress regeneration in a closed mini-longwall is shown. The agreement is good, lending confidence to the parameters used in the proposed numerical model. Foundation yielding was modelled in the absence of total stope closure. A ?Prandtl? wedge was formed, coinciding with the yield point of the pillar/foundation system. Some theoretical work is used as a starting point for the formulation of a new design methodology for strike or dip stabilizing pillar/foundation systems in deep level goldmines. The existing guideline for the stability of foundations (APS< 2.5UCS) is critically examined. Finally, the advantages of the new approach, and the work required to substantiate the proposed new design methodology, are outlined.

Jan 1, 1998

By D. F. Malan

Considerable practical knowledge and effective use of empirical evidence is required in the design of tunnel support systems. In many instances, it is unclear how individual support elements interact with fractured rock and which elements of the support system are critical in controlling and containing the movement of material proximate to the excavation. The present paper outlines a displacement discontinuity approach that allows the coupling of surface constraints, in the form of liner or mesh material, to an explicit representation of time-dependent rock discontinuity creep movements. The method is illustrated for a simple square-shaped tunnel profile in which the fracture zone is mobilized by an imposed stress field. A crucial aspect of the model is the computational treatment of the coupling between the liner material and there presentation of the fractured rock mass. This allows the response of the liner material to be assessed quantitatively in terms of the possible reduction of movement in the rock mass. The importance of installing the liner timeously after the excavation is formed is illustrated in a specific example. A number of computational difficulties in simulating liner-rock interaction mechanisms are identified.

Jan 1, 2008