Search Documents

By K. Fotoohi

This paper presents a basic, and yet practical, three-dimensional finite element model that was developed by the authors for stress/stability analysis of mine structures. While the code is linear elastic at this stage, it has some practical features that make it an attractive tool for mine design. Such features include: simulation of mine backfill support, mining induced strain energy calculation, powerful out-of-core frontal solver and failure analysis by Mohr-Coulomb or Hoek-Brown failure criteria. Model development criteria are discussed, and the potential application of the model to mining is demonstrated by design examples.

Jan 1, 1998

By Wilson Blake

A mathematical model based on the finite element method of stress analysis has been used to describe the behavior of the western wall of the Kimbley Pit as its slope was steepened from 45 degrees to 57 degrees. Theoretical stress changes because of slope modification are slight except at the toe of the modified slope. Magnitudes of stress created are small when compared with the strength of a competent rock. Resulting displacement patterns point out the need for determining the in situ state of stress at a site prior to excavation or modification. A comparison of model results with field measurements shows that, in general, the model has described the measured behavior of the west wall of the Kimbley Pit remarkably well.

Jan 1, 1968

By I. S. Laitinen

Outokumpu Technology has recently introduced new technology in both copper electrorefining and electrowinning. This technology includes, among other things, a new permanent cathode (Outokumpu Permanent Cathode) and new busbar technology (Outokumpu Double Contact System). The new technology has already been tested in commercial scale cells. To get a better understanding of the busbar system a multiphysical simulation model of the electrowinning and electrorefining cells was developed. A mathematical model of electrical conduction and heat transfer in a copper electrolysis cell was developed from partial differential equations. The model was expanded to a cell group including seven electrolysis cells connected electrically in series. The model was solved numerically using the finite element method (FEM). Using the model, different intercell busbar systems were simulated. Typical electrical disturbances, like loose contacts between cathodes and the main busbar and short circuits between electrodes, were included in the simulation. The results of the simulation runs with conclusions are presented in this paper.

Jan 1, 2007

By Shad M. Sargand, Glenn A. Hazen, A M. M Jayasuriya

Various types of waves originate from a blast source to propagate through rock/soil media away from the blast. Input vibration to a structure depends on the soil-structure interaction which varies with the dynamic properties at the intersection as well as the stress and strain levels developed from the vibration of the structure-ground system. Hence it is important to study such effects on the response of a structure due to blast loads and this knowledge can be used in the design and control of blasts that are carried out close to buildings, highways and other structures.

Jan 1, 1994

By B. Q. Li, H. C. de Groh, X. Lu

"A fully 3-D, transient finite element model is developed to represent the magnetic damping effects on complex fluid flow, heat transfer and electromagnetic field distributions in a metallic melt. The model development is based on our in-house finite element code for the fluid flow, heat transfer and electromagnetic field calculations. The numerical model is tested against numerical and experimental results for water as reported in literature. Various numerical simulations are carried out for the melt convection and temperature distribution with and without the presence of a transverse magnetic field. Numerical results show that magnetic damping can be effectively applied to stabilize melt flow, reduce turbulence and flow levels in the melt and over a certain threshold value a higher magnetic field resulted in a greater reduction in velocity. Also, for the study of melt flow instability, a long enough running time is needed to ensure the final fluid flow recirculation pattern. Moreover, numerical results suggest that there seems to exist a threshold value of applied magnetic field, above which magnetic damping becomes possible and below which the convection in the melt is actually enhanced.1. INTRODUCTIONThe imposition of a DC magnetic field on a moving, electrically conducting fluid, such as molten metal or semiconductor melt, produces an opposing Lorentz force that results in a reduction of flow velocity in the liquid. This effect has been explored in both the semiconductor and metals industries to control natural convection and/or surface tension induced Marangoni convection as well as thermally induced flow instability during solidification processing of these melts [1-6]. A mathematical model can be of great use in developing both fundamental understanding of the fluid-magnetic field interactions in these magnetically-assisted processing systems and basic guidelines for process design, control and optimization.This paper presents a 3-D mathematical model for magnetically damped fluid flow and heat transfer in a melt cavity, which is useful for interpreting results for unidirectional solidification processes at slow solidification rate, and some preliminary results concerning the magnetic damping effects on the flow instability in the melt. The model development is based on the finite element solution of the Navier-Stokes equations, energy equations and equations for induced electric field. The model is tested against reported results and is then applied to study the fluid flow and thermal distribution in the metallic melt with and without an external magnetic field."

Jan 1, 1999

By A. S. Usmani

The flow associated with the pouring of molten metal into moulds has been modelled on the basis of certain simplifying assumptions. A finite element based Navier-Stokes equation solver has been used to analyse the flow. The velocities obtained from the Navier-Stokes equation solver are used to advect a pseudo-concentration function for modelling the fluid front. A simple problem with a known analytical solution is solved first to test the model. Subsequently, several more examples with varied geometrical configurations are modelled, including one with curved flow boundaries.

Jan 1, 1992

By D Frantzos

Pressure time histories were recorded for low density ammonium nitrate/fuel oil, detonated in long heavy walled steel cannons of various bores. These were shell used in a finite element model of a horizontally layered limestone rock mass to predict the crack propagation limits in presplit blasting for a range of borehole diameters. Apart from showing very good agreement with field results, the model clearly demonstrated the strong dependence of the results on the pressure-time curve. The important elements were the peak pressure, the rise time to it and its duration. Control of these characteristics offers the possible future production of more ideally suited pressure-time profiles for the optimization of crack propagation distances and borehole spacings for various ground conditions. Also the results to date provide the basis for investigating the characteristics of the explosives presently used in presplit blasting, and finding ways to modify them with the purpose of optimizing the field results.

Jan 1, 1987

By Francois E. Heuze

Excavation of large rock caverns is a stepwise process. Whether one uses a heading and bench method or a multiple drift system, the use of excavation steps enables one to introduce the reinforcement, such as bolts or shotcrete, at the earliest possible stage. Knowing the spatial distribution and properties of rock materials and geological discontinuities, the cavern geometry and its reinforcement must be optimized. The finite element method of analysis is a powerful tool for such optimization. However, the models are only as good as their ability to account realistically for the physical attributes of the prototype and for the excavation sequence. They must include the following: - ability to model an incremental construction process. Thus requires the use of "Restart" procedures where displacements and stresses from a previous equilibrium are saved for the next step. - ability to model non-linear rock material behavior; in particular, the post-peak behavior of the rock which retains a residual strength when con-fined. The stability of many under-ground structures rests with the "broken" rock carrying significant loads. - ability to account for the non-linear behavior of geological discontinuities. These are two-fold. First, up to peak, non smooth rock joints will dilate. In a confined situation, dilation will increase joint strength by a consider-able amount (Obert, et al, Brady, 1971, Heuze,1978a), Then, beyond the peak, the joint also retains residual strength.

Jan 1, 1978

By Daniel W. H. Su

Surface subsidence induced by multiple-panel coal extraction was calculated with finite element stress analysis. The use of nonlinear material behavior and GAP elements, which provide a realistic representation of shearing along existing planes of weakness, allows the finite element program to accurately reproduce observed subsidence profiles. A reduction factor of 1/6 is applied to the measured rock strength and modulus for use in the finite element model. The models satisfactorily predict a significant change in maximum subsidence and subsidence profile associated with an increase from subcritical to near critical panel width. Modeled shearing along planes of weakness is in good agreement with available time domain reflectometry field observations and is demon¬strated to have a very profound effect on the maximum subsidence and the shape of the subsidence profile. When both material nonlinearity and planes of weakness are incorporated, the finite element model appears to possess great potential for accurately predicting surface subsidence and subsurface strata deformation in environments where little subsidence information is available.

Jan 1, 1990

By Daniel W. H. Su

The ability to predict the effects of longwall mining on topography and surface structures is important for any coal company in making permit applications and anticipating potential mining problems. The sophisticated finite element model described and evaluated in this paper is based upon five years of underground and surface observations and evolutionary development of modeling techniques and attributes. The model provides a very powerful tool to address subsidence and other ground control questions. The model can be used to calculate post-mining stress and strain conditions at any horizon between the mine and the ground surface. This holds the promise of assisting in the prediction of mining-related hydrological effects.

Jan 1, 1986

By Youliang He

Twin-roll casting is an integrated casting and rolling process that involves the continuous solidification of a melt and the plastic deformation of the solidified metal between two rotating rolls. Computer modeling is a valuable tool to help engineers understand the metallurgical process during twin-roll casting. However, commercially available software can only model either the fluid flow/heat transfer of the melt or the plastic deformation of the solid, but not both. In this study, an object-oriented finite element program was developed to simulate the steady-state fluid flow, solidification and plastic deformation of an A7,31 magnesium alloy in a horizontal twin-roll caster. The melt, the mush zone and the solid were treated as a generalized Newtonian flow in an Eulerian reference frame and different constitutive equations were employed to describe the material behaviours. The effect of operational parameters on the solidification process was discussed based on the simulation results.

Jan 1, 2010

By S. Behera, K. Dey

This paper presents a methodology for finite element (FE) modelling to analyse the total deformation of rockmass under blast-load, with a focus on predicting backbreak that may occur during surface blasting operations. Backbreak can result from improper use of explosive energy, leading to damage beyond the excavation zone and jeopardizing operational safety and productivity. The paper highlights the impacts of backbreak on bench stability, excavator operation, drilling operation, flyrock generation, and improper fragmentation and boulder formation in surface blasting. However, prior estimation of extent of rock breakage zone can provide an idea to formulate the blast design to control the backbreak. While vibration-based prediction models are commonly used to estimate backbreak, they require physical blasting experimentation. To overcome this limitation, the paper describes the development of a FE model using ANSYS Explicit dynamics, which simulates rock blasting and calculates the total deformation of the rockmass. The rockmass deformation is modelled using Drucker-Prager (D-P) strength model with the Jones-Wilkins-Lee (JWL) equation of state for the explosives. The accuracy of the FE model is evaluated through physical experiments, and the results are compared with the vibration models. The FE model demonstrates close proximity to physical observations as well as vibration model result in predicting backbreak, thereby offering a promising alternative to vibration-based prediction models. The primary objective of the study is to predict the response of the insitu rockmass to explosive blast loading and determine the effective rock breakage zone.

Jan 1, 2024

The paper describes a practical methodology based on numerical modelling for the prediction of progressive failure of underground mine stopes with cable bolt support. A case study of the 40-761 stope at Kidd Creek Mine, Timmins, Ontario is presented. The stope was mined by blasthole open stoping in 1986 as a lead stope within a mining block between 36-1 and 38-2 levels. Initially designed to peak at 38-2 level, the stope back caved upwards approximately 40m towards 36-1 level. Cable bolts were installed above the stope from 36-1 level after caving had begun. The stope eventually stopped to cave when the falling rock mass encountered a sufficiently high cable density. The paper provides details of the geomechanical data and modelling results of the case study. Numerical modelling is carried out using a finite element package known as e-z tools which has, among other things, a special cable/rockbolt simulation capability.

Jan 1, 1994

By P. K. Currie, J. V. Walters

Solution mining of magnesium salts (Carnallite and Bischofite) is taking place near Veendam, in the northeast Netherlands. The magnesium salts are used for the production of magnesia sinter, used as a base material for the manufacture of refractory bricks for furnaces. This paper gives, first, an introduction to the operational procedure followed in Veendam, with the emphasis on some of the unique aspects, and then discusses rock-mechanical aspects of the cavity behaviour. The creep behaviour of magnesium salts is described, and Finite Element calculations on cavity stability and potential surface subsidence are discussed.

Jan 1, 1985

By I. Mc Dougall

Finite element models were generated of two submerged arc furnaces of different geometries. A 48MW circular furnace and a 68MW 6-in-line rectangular furnace were studied. The electrodes, raw material, slag and molten metal were included in the model. ANSYS/Multiphysics was used to predict the current density distribution in the electrodes, raw material, slag and molten metal as a result of the three phase AC current. The effect of the electrode immersion on the current path was studied. The ?proximity effect? caused a pronounced asymmetry in the electrode currents in the circular furnace, whilst the electrode currents in the rectangular furnace exhibited radial symmetry. The current distribution in the melt in both furnaces indicated that the bulk of the current flows between phases. The current density is greatest between electrodes. Results were compared to qualitatively similar published work, operating experience and results of furnace dig-outs.

Jan 1, 2007

By F. Fournier Dit Chabert, F. Tancret, J. François Castagné, R. Le Gall

During the repair of turbine blades by a laser cladding process, it is observed that interdendritic cracking occurs in the heat affected zone (HAZ). Auger electron spectroscopy, performed on sample broken at high temperature, reveals interfacial sulfur segregation, which is known to be one of the main causes of microcraking in the HAZ during welding processes. Moreover, thermal gradients generated by the process are responsible for important thermomechanical stresses. In the present work, the finite element code Comsol Multiphysics is used to simulate the laser cladding process in order to determine both the thermal field and, using an elastoplastic law, the stress field inside the component. The originality is that the temperature history calculated in the thermomechanical model is extracted then injected in a segregation model in order to calculate the local segregation behavior. This model can be eventually used to predict the risk of cracking as a function of process parameters and alloy chemistry.

Jan 1, 2006

By Richard K. Komarek, Antonios Zavaliangos, Roman T. Dec

Current industrial compacting and briquetting practices are largely based on trial-and¬error techniques. This paper introduces a two dimensional finite element model of compacting/briquetting process in the roller press. Compared with other modeling techniques used in the past, finite element method offers the most versatile approach because it incorporates adequate information about powder behavior, geometry and frictional conditions and offers predictions for the density distribution, velocity field, stress distribution, roll force, roll torque and other important process variables. More over, when the model includes elasticity, it offers predictions on the internal stress state in the briquet after its exit from the pressing zone. This enables to perform realistic computer experiments which are helpful for establishing the pressing system design, forming cavity geometry and process conditions necessary to produce compacted or briquetted product with the desired properties.

Jan 1, 2004

By C. Arson

"Optimizing hydraulic fracture injection parameters in order to maximize hydrocarbon extraction and to avoid extensive borehole spalling is still an open issue. In this paper, the Differential Stress Induced Damage (DSID) model is employed to simulate the anisotropic damage distribution around pressurized boreholes drilled in fractured shale. The loading path is purely mechanical: excavation and fluid injection are simulated by stress relaxation followed by pressurization. We investigate the effect of pre-existing fracture orientation, in-situ stress and injection pressure on the propagation of anisotropic damage during pressurization. Finite Element results indicate that (1) damage initiates and develops when the ratio of two in-situ stress component deviates from 1, and the magnitude of damage is proportional to this stress deviation; (2) For a given in-situ stress field, damage initiates for a certain threshold of injection pressure, and the damage zone extends as the injection pressure increases; (3) A pre-existing fracture only affects the damage zone locally; the orientation of that fracture does not affect the damaged zone, which is controlled by in situ stress and injection pressure.INTRODUCTIONHydraulic fracturing is used to increase shale permeability and to enhance hydrocarbon extraction. In order to maximize the extraction ratio and to avoid wellbore spalling and clogging during fluid injection, the parameters related to injection for each particular site need to be optimized. Therefore, it is necessary to clearly understand how shale responses upon excavation and pressurization.It is well known that breakouts grow in the direction of minimum compression and that plastic deformation around wellbores tends to arrest breakouts (Zoback et al., 1985). Laboratory acoustic emission tests revealed that fractures propagate along the direction of maximum compression stress, and fractures (breakouts) are surrounded by the Excavation Damaged Zone (EDZ), in which rock properties are weakened by the presence of micro cracks (Eberhardt et al., 1999; Martino et al., 2004). The anisotropy from orientated foliations and cracks may influence the development and propagation of fractures and breakouts (Everitt et al., 2004). Sonic televiewer logs proved that during hydraulic fracturing, fluid injection further alters the stress distribution (Evans et al., 2005). Shale natural fabric strongly affects the evolution of the EDZ around boreholes."

Jan 1, 2015

By K. Pietruszka

The article presents the concept of the application of the Finite Element Method (FEM) and ABAQUS software in the global monitoring of changes in the rock mass and on the surface of the area, as a result of backfilling the workings. Waste materials applied for stabilizing old mining workings can have a strong influence on the rock mass and surface. According to the regulations, this process should be carefully monitored, especially in terms of the safety of surface structures and rock mass stability. The article presents the concept of an analytic and graphic presentation of the results of the measurements of old workings, applying the Finite Element Method and ABAQUS software in model studies. The main purpose of such interpretation of the geodetic observations results is the application of methods providing wide possibilities to present the geometry of workings, their changes and the influence of the rock mass and surface.

Jan 1, 2008

By Victor E. Saouma

Chip formation in rock under a line load and in front of a drag bit cutter is numerically investigated. Analysis is accomplished by a special purpose interactive graphics finite element code, SICRAP, written for the simulation of mixed mode crack propagation under 1inear elastic fracture mechanics assumptions. The first study provides some interesting qualitative results, and in the second study, correlation with experimental tests on chip formation by drag bit cutter in Berea sandstone is found to be very satisfactory. This indicates that elastic analysis, coupled with fracture mechanics, is capable of modeling rock cutting.

Jan 1, 1984