Search Documents

By Joshua R. Omer

A numerical model is developed for predicting the load transfer and settlement behaviour of vertically loaded, bored cast in-situ piles. Based on analysis of full-scale pile tests, closed-form equations are formulated to uncouple the shaft and base transfer resistance components of the pile/soil system. Appropriate parameters are introduced to take into account: (i) the shaft resistance distribution in unlined and partially lined piles, (ii) the effects of ?soft toe? (voided toe pile) behaviour on base resistance development. For convenience, an interactive computer program has been developed in Matlab? with a user friendly interface system. The program can be applied to predict the complete load-settlement curve and axial force distribution in bored cast in-situ piles if reliable soil test data are available. To demonstrate its validity, the proposed method is applied to 4 test piles, which are also analysed using two previously published methods: (a) elastic solution method and (b) layer discretisation method. Methods (a) and (b) above are based on load transfer (t-z) and finite element models respectively. It is shown that the proposed method yields axial force distributions and load-settlement curves that are in good agreement with the observed pile test data.

Jan 1, 2010

By P Lindsay

A numerical model has been developed to simulate the effectiveness and costs of a number of conceptual remedial measures aimed at reducing acidity and metal concentrations in the Ngakawau River and its tributaries near the Stockton Coal Mine, situated on the West Coast of New Zealand. The Stockton Coal Mine has been in operation over the past 100 years. During that time, seepage from waste rock dumps has contributed to elevated acidity and metal concentrations in the Ngakawau River and some of its tributaries. The yearly mean rainfall at the mine site is 6000 mm. The model, developed using GoldSim software, simulates rainfall, run-off and stream discharge for 26 tributaries. It then calculates acidity and metal concentrations in the mine impacted tributaries and the dilution downstream from the mine. The model has been calibrated to measured flow and acidity and metal concentrations in the Ngakawau River and some of its tributaries. Acidity and metal loads are removed from the system in accordance with specified efficiencies and capacities of the respective conceptual remedial measures. The user can compare specific remedial measures, or a combination of conceptual remedial measures. The model predicts the acidity and metal concentrations at a number of monitoring stations for any remedial strategy selected by the user. It also calculates the capital and operational costs for that particular remedial strategy. This paper presents the development of the numerical model and describes how the end-user can use it to compare the efficiency and costs of various remedial strategies by running different æwhat ifÆ scenarios.

Jan 1, 2006

By S. K. Sharma, ANKUSH GALAV, G. S. P. Singh

Providing a suitable protective water barrier pillar (PWBP) is common to reduce inundation hazards in underground coal mines. The imperative factors influencing its performance include the water head acting on the pillar, cover depth, pillar width, strength properties, and permeability characteristics. The mechanical failure of such pillar is a stress-controlled phenomenon, whereas the hydraulic failure is a strain-based phenomenon. A finite-difference numerical modeling approach was developed to study the hydro-mechanical coupled behavior of protective water barrier pillars. The mechanical stability was evaluated in terms of the percentage of failed (ZoF) and intact zones. The influence of the strain-controlled weakening on the permeability of the flow medium was studied through the coupling of the mechanical and hydraulic effects. The coupled steady-state model was used to estimate the outflow rate and its hydraulic stability. The adequacy of the protective pillar was also investigated by assessing mechanical stability and capability to resist hydraulic pressure against the maximum expected water head. A seepage rate-based classification system has also been proposed to evaluate the seepage potential and assess the hydraulic stability. The model has been validated for two case studies at the cover depth of 136–189.5 m and the existing pillar width of 16–42 m against the water head of 25–141 m.

Aug 18, 2022

By S.K. SHARMA, ANKUSH GALAV, G.S.P. SINGH

This paper proposes a hydromechanically coupled numerical modeling approach to quantify the influence of various parameters on the mechanical and hydraulic performance of protective water barrier pillars (PWBP). The results are validated for two case studies having pillar widths of 16 and 42 m at cover depths of 136 and 190 m and water heads of 25 and 141 m. The outcome can be significantly helpful in the rational design of these pillars in underground coal mines.

Mar 1, 2023

By Gang Mei, Nengxiong Xu, Yanlin Zhang

After open-pit mining finishes, an iron mine with steeply inclined ore-bodies will generally transform to underground mining. Non-pillar sublevel caving method, with horizontal slice mining from up to down, is commonly used for mining deep ore-bodies under the open-pit bottom. As horizontal slice mining proceeds, a gradually enlarging goaf will be formed under the open-pit bottom, the rocks surrounding the handing and heading sides of the goaf will collapse continuously, and then these results in rock-mass failure with large deformation and displacement in the open-pit slope. Based on FLAC3D, a modeling approach to predict pit slope deformation and failure induced by transitioning open-pit to underground mining is proposed in this paper. This method mainly consists of three procedures of modeling, i.e., mining modeling, landslide modeling, and filling modeling, of which landslide modeling is the core of the proposed approach. In the procedures of modeling with FLAC3D, iteration non-convergence is taken as a local failure criterion of pit slope. When iteration calculation cannot converge, the sliding modeling will start by deleting some elements with large displacements, and finish until convergence. Before modeling landslide, several displacement contours of rock-masses in the pit slope are searched out, and each contour is an approximate possible sliding surface of the open-pit slope. From these contours, the most proper one is selected as the slope sliding surface by using the iteration calculation convergence and minimum element deletion as the judgment condition. Landslide modeling is then realized by removing the elements outside the slope sliding surface. Finally, to demonstrate the effectiveness, the proposed approach is used to simulate the process of open-pit to underground mining at Yanqianshan Iron Mine, and to predict the deformation and failure characteristics of open-pit slope during underground mining. The case study results show that: the prediction of the slope stability determined by using our method is consistent with that calculated by using the limit equilibrium methods.

Jan 1, 2015

By Janusz Cygankiewicz, Marian Branny, Józef Wactawik

The paper presents mathematical model of concealed fire. We consider a porous medium (goaf, coal fissures) through which air flows. The set of equations is constituted by transient, two dimensional mass, momentum and energy balance equations with sources describing oxygen loss and heat generation caused by coal oxidation. Paper presents several examples illustrating the phenomenon of spontaneous heating in goaf areas behind longwall faces.

Jan 1, 1997

By Ch. Arvanitidis

This paper presents the design solution for an actual case of settlement tanks, for a wastewater treatment plant under construction, actually founded on swelling clay (the nature of which had not been appreciated beforehand) and suffering considerable damage as a result. Two points of interest prompted the use of a numerical model as an alternative approach to the traditional method of analysis: First, the desire to compare the in situ measured deformation and displacement of the tanks with the figures calculated by the numerical model. Second the need to check the validity of the suggested solution as a remedy, i.e. the replacement of a certain thickness of clay by graded material.

Jan 1, 2005

By Patrick Linzer

The majority of existing models describing (methane) gas transport within, and liberation from, coal structures treat these as continua wherein two main gas flow modes (diffusive and Darcy permeable) prevail, and use finite difference techniques to generate solutions to the resulting multivariate parabolic partial differential equation(s). Isotropy and homogeneity in the physical and geometric coal structure properties are often assumed, and the existence of fundamental random fluctuations in these proper- ties is not catered for. The model proposed herein does not suffer these inadequacies, and its efficiency is enhanced by the replacement of the diffusion equation by an accurate approximation. Dynamic mining activity, gas contributions made by the geology surrounding the mining seam, goafing of overlying strata and the effects of induced strain are all taken into account.

Jan 1, 1995

By Dale S. Preece

The computer program. DMC Qistinct Motion &de), which w&s developed for simulating the rock motion associated with blasting, has been used to study the intluence of row delay timing On rock mOtiOn. The munerical simulations correspond with field observations in that vmy short delays (< 5Oms) and vety long delays(> 3oOms) produce a lower percent-cast thao a medium delay (100 to ZOO ms). The DMC predicted relatio&lp between row delay timing and percent-cast is more complex than expected with a dip in the auve where the optimum timing might be expected. More study is required to gain a full undasranding of this phenomenoo.

Jan 1, 1994

By M. D. Voegele, C. Fairhurst

INTRODUCTION A computer study of excavation support loads in tunnels and other excavations in discontinuous rock masses was undertaken with the distinct element method introduced by Cundall (1971). In the distinct element method the rock mass is modeled as a composition of individual blocks of rock defined by the joint planes. Normal and shear deformation laws are ascribed to the joint surfaces; the blocks of rock themselves are considered to be rigid. The particular version of the program utilized in this study ran a mini computer and used inter, active techniques for the input of geometrical data such as joint orientation and spacing and subsequent display of rock mass deformation and other data. Two major aspects of the behavior of excavations in jointed rock masses are examined in this paper: i) the behavior of unsupported excavations in rock masses with varying joint configuration, joint properties and in situ stress fields; and ii) the modification of the behavior of similar excavations when support is provided. In the latter case the excavation behavior was examined through the use of ground reaction curves. THE DISTINCT ELEMENT METHOD The Distinct Element method introduced by Cundall (1971) is a computer based analysis that simulates the behavior of a system of discrete, semi-rigid rock blocks. Block interactions are governed by realistic friction and stiffness laws; each block may undergo unlimited displacement and rotation while progressive failure is modeled.

Jan 1, 1982

By H. Gercek, M. Genis

In this paper, the effects of dynamic loads on the geometry of the failure zone occurring around underground openings, excavated in deep and massive rock masses, have been investigated. Following an elasto-plastic stress analysis of a two-dimensional underground opening under static conditions, dynamic loads have been applied. In the static and dynamic analyses, FLAC2D computer program has been utilized. Maximum ground acceleration value is selected as an input parameter. An index determining the relative size of the yield zone around the opening has been used for comparison of instabilities. In diesem Bericht wurden die Auswirkungen von dynamischen Lasten, auf die Geometrie der Bruchzone in der Umgebung von Untertagehohlräumen, die in tiefen und massiven Gebirge aufgefahren wurden, untersucht. Nach einer elasto-plastischen Spannungsanalyse eines zwei-dimensionalen Untertagehohlraumes unter statischen Beanspruchungen, wurden dynamische Lasten aufgebracht. Für die statischen und dynamischen Analysen wurde das FLAC2D Computer Programm verwendet. Die maximale Bodenbeschleunigung wurde als Eingabe Parameter gewählt. Zum Vergleichen der Unstabilitäten wurde ein Index verwendet, der die relative Grösse der Bruchzone in der Umgebung des Hohlraumes bestimmt. Dans ce papier, l’influence des chargements dynamique sur la géométrie de la zone en rupture observée autour des cavités profondes excavées dans des massifs rocheux ont été étudiées. Des chargements dynamiques ont été appliqués, dans le cas d’une modélisation bi-dimensionnelle de l’excavation d’une cavité profonde, sur la base d’une analyse élastoplastiques des contraintes en conditions statiques. Le code de calcul FLAC2D a été utilisé pour les analyses statiques et dynamiques. La valeur maximale de l’accélération du sol est choisie comme paramètre d’entrée. Un indice déterminant la taille relative de la zone sollicitée autour de la cavité a été utilisé pour comparer les différentes instabilités.

Jan 1, 2003

By Saba Gharedash

In this study the results of a series of numerical simulations which are conducted for one of the tunnels of the Mashhad metro line 2 project are presented to evaluate the role of tunnel face support pressure and soil strength parameters on deformation of tunnel lining and face deformation. To achieve this goal, comparative numerical simulations of a mechanized tunnel advance in soft non- homogeneous soil on top of the ground water table and beneath the building are conducted and sensitivity analysis is performed. Development of a step by step tunnel construction process is modeled using a three dimensional finite difference model (FLAC3D) taking into account all the relevant components of shield tunneling. In order to understand the mechanism of the effects of a variety of support pressures on face deformation and lining stresses, numerical simulations are conducted. It is found that the use of different face supports causes a considerable influence on the forces of tunnel lining. In addition, it is revealed that building weight help to change the stress paths and the stress distributions in the tunnel lining.

Aug 1, 2013

By Edson Charikinya, Pierre-Henri Koch

"The extraction of mineral values from ore requires liberation followed by separation steps. Liberation is achieved by size reduction operations which are energy inefficient processes typically dominating the energy consumption in a mineral concentrator. As the grade of ore reserves declines, future viability of mineral operations will be determined by energy costs of comminution. The application of high power microwave energy to secondary crusher products has been suggested as a possible commercially viable thermal treatment method for reducing comminution energy and improving value mineral liberation. Recent studies have shown that microwave pre-treatment of coarse sphalerite ore particles (> 5mm) at specific microwave heating energies (1-3 kWh/t), induces microfractures and creates new crack surfaces. This suggests that subsequent crushing of these microwave treated particles could yield enhanced liberation. However, limited studies have been carried out investigating the mode of breakage and the extent of enhanced liberation in that case. The objective of this study is to develop numerical methods for quantifying the extent of enhanced liberation and mode of breakage in crushed microwave treated and untreated particles. Sphalerite ore particles representing small (-5+4.75) mm, medium (-16+9.5) mm, and large (-25+19) mm HPGR and cone crushed particles were microwave treated at specific energies between 1-3 kWh/t. Cracks in the ore particles before and after microwave treatment were analysed with QEMSCAN and numerical models of the measured particles (before and after microwave treatment) were developed in MATLAB. The propagation of random and non-random cracks was investigated by simulating comminution of the modelled treated and untreated particles. Results of this study demonstrate that the breakage mechanism and the liberation of valuable minerals from gangue in microwave treated and untreated particles is significantly different. This study contributes to the development of numerical tools to quantify crack propagation in heterogeneous ore particles."

Jan 1, 2018

By M. A. Sobhi, L. Li

With the continuous reduction of mineral resources at shallow depth, more and more mines have to go deeper beneath the earth's surface to reach the ore zones. This increase in mine depth leads to a natural stress increase, which also raises the stresses around mine openings. This may result in ground instability, including the risk of rock bursts or exceedingly large opening closure. Stope backfilling is a well known solution to improve the ground stability condition. It also has the advantage of reducing the mining wastes disposed on surface, thus helping to reduce the environment impact from the industry. However, most available solutions for backfill design were developed for an isolated stope without taking into account the influence of mine depth, except for a few numerical analyses performed recently by Falaknaz and coworkers for the case of multiple backfilled stopes laying in the same level. In this paper, a numerical study is carried out to evaluate the stresses in backfilled stopes overlying a sill mat. The mine depth and the excavation of the underlying stope below the sill mat are both taken into consideration. Compared to the case of a single isolated backfilled stope, the obtained numerical results show that the stress magnitudes in the overlying backfill are considerably increased due to the excavation of the underlying stope. In general, the stresses increase with mine depth and backfill stiffness, but decrease with an increase in the surrounding rock mass stiffness. These results indicate that the barricade and sill mat backfill design should take into account the excavation of the underlying stope and the mine depth.

Jan 1, 2015

By Hyu-Soung Shin

This study is concerned with the characteristics of the behavior of steel-concrete composite drilled shafts, a material that is used for large bridge foundations, and especially for weak submerged ground conditions. Depending on the loading or boundary conditions considered, the behavior of foundation will vary. The bearing capacity of steel-concrete composite drilled shaft is influenced by initiation and propagation of fractures along the interface between steel casing and filled concrete. The interface effect is taken into account with the following two simplified cases: (1) detached model with free interface (no shear stiffness of interface) between concrete core and steel casing, (2) confined model by an upper steel cap with free interface. In the confined model, it is intended to have an unified behavior of composite shaft between concrete core and steel casing but there is interface between them. Intensive numerical studies have been undertaken to investigate detailed behavior of the composite shaft. In this study, 3-D finite element analysis is employed. Ground conditions typical of those encountered for the design of drilled shaft under water are considered affecting the behavior. Also, an attempt is made to generate possible combinations of the loads which have been considered in fields. Detailed characteristics of the stress and displacement distributions are evaluated to understand the mechanisms of composite shafts.

Jan 1, 2007

By M. Ghamgosar, Erarslan N.

Rock cutting includes many unknown parameters which result in different consequences of the different cutting conditions. However, based on the current suggested models for rock cutting, there are some critical parameters that control the fracture initiation and propagation in the rock domain. In general, there are two major cutting tools are used in rock cutting: drag bits and indenters. Different mechanical attacks of cutters result in different fracture models. For example drag bits break the rock parallel to the rock surface, which is called ‘undercutting’ while indenters cut rocks by pressing normal to the rock surface. Both tools cause rock failure by inducing tensile cracks in the rock. However, drag bits induce more tensile cracks compared with the indenters. The full benefits of improvement in rock cutting efficiency will arise from a clear understanding of rock fracture mechanics principles. Therefore, fracture mechanics principles will be applied to the rock cutting as the rock fragmentation process is based on the rock fracturing. The Oscillating Disc Cutting (ODC) has more advantages than drag bits and disc cutters. The forces transmitted back through the cutting head are much lower with ODC. Although ODC attacks rock in much the same way as drag bits, a combination of the oscillating action of the cutter and the undercutting mode produces the lower cutting forces and causes rock fatigue. Fatigue phenomenon in rock has advantages for the rock breaking by applying less effort and energy due to the cyclic loading. The aim of this study is using the Extended Finite Elements Analyses (XEFM) numerical modelling for simulating the ODC cutting with different cyclic loading conditions. An accurate prediction of numerical investigation requires detail laboratory experiments. A series of laboratory tests were performed to obtain the fracture toughness values of our Monsonite rock samples under various cyclic amplitude and frequency values. The numerical results showed different types fatigue damage mechanisms due to the various cyclic loading conditions. Thus, modelling and understanding the crack propagation behaviour depending on different amplitude and frequencies are believed very important for the dynamic rock cutting researches and rock cutting machine manufacturers.

Jan 1, 2015

By Juhyung Lee

Recently, large diameter drilled shafts are often used for foundations of large-scale marine bridges with high capacity and required to resist lateral loads. Steel casings are essential in supporting the hole-diving concrete installation of the drill shafts constructed in the sea. The casings, however, are usually not considered as the structural elements and their load bearing capacity is often overlooked in the design of drilled shafts. As the diameter of drilled shafts becomes larger, there are higher demands for more efficient and economical design of large-scaled drill shafts in the sea. Steel-concrete composite columns are popular for the superstructure of bridges, and the outer steel attached to the shaft increases the shaft resistance through confining the concrete. In this study, the lateral resistance of steel-concrete composite drilled shafts was evaluated quantitatively based on the numerical analysis. When steel casings are used as structural elements such as the composite columns. Ultimate lateral bearing capacity of composite drilled shafts with various diameters was numerically calculated through the 3D finite element analysis. For that, elasto-plastic model with the perfect plasticity is involved to capture the ultimate load. A commercial FEA program, MIDAS-GTS, is used in this study. Real field conditions of the west coast, Korea were examined to set up the ground conditions and pile lengths required for this parametric studies. Detailed characteristics of the stress and displacement distributions are evaluated for the better understanding on the mechanisms of the composite shaft behavior.

Jan 1, 2008

By F. Mackay, S. A. B. da Fontoura

"The theoretical solution for a viscoelastic well was presented by Gnirk and Johnson in 1964. From this an exercise was executed to validate the numerical simulation of a cemented salt well. The exercise comprised the geomechanical behavior through time of a cemented salt well, this was achieved using a finite element commercial software. This simulation showed the accuracy of the theoretical solution for viscoelastic materials such as salt rock. The theoretical solution presented elastic and viscoelastic analytic equations for the tangential and radial stresses. The equations for the elastic material (casing and cement) and for the viscoelastic material (salt rock) are presented here: Where equations 1a & 1b are the radial and tangential stresses of the elastic material; equations 2a & 2b are the radial and tangential stresses of the viscoelastic material; and a, r and b are the geometrical parameters explained in Figure 1. These equations are function of Pb, a stress time dependent factor that comes from the Burgers model of viscoelasticity. This is well explained in Gnirk and Johnson (1964) and for the sake of simplicity is not disclosed here. Po is the stress at section A-A’ dependent of the depth h."

Jan 1, 2015

By Emile Scheepers, Chris Aldrich, Jacques J. Eksteen

The pneumatic injection of CaC2 powder into molten iron has become a very reliable and effective way to carry out desulphurisation. The key variables that influence the kinetics of this process are the condition and composition of the top slag, the oxygen activity in the hot metal and the injection conditions. Notable injection parameters include carrier gas flowrate, reagent flowrate, lance depth and initial sulphur concentration. A one-dimensional, steady-state model for momentum, heat- and mass transfer in rising gasliquid- powder plumes was developed for conditions relevant to the CaC2 refining process at Saldanha Steel?, South Africa. Combined with a model that predicts the contribution of the topslag to the process, the overall rate of desulphurisation as a function of time was determined, thereby quantitatively demonstrating the influence of the above-mentioned injection parameters on the kinetics of the desulphurisation process.

Jan 1, 2003

By Whiten WJ

A method that allows data to be mass balanced and screen models to be fitted in a single step is described. The inclusion of several variables to describe the feed data in the fitting procedure allows a mass balance to be conducted at the same time as the fitting of the screen model. Analysis of Groote Eylandt Mining Company screening data has demonstrated that this method is ideally suited to applications where complete stream data are not available. Also, certain random fluctuations in screen model parameters arising from the fitting of consecutive mass balances and screening models are significantly reduced. A cubic spline function with fixed knot positions corresponding to particle sizes and variable (fitted) cumulative percent passing values is used to describe the feed size distribution. When the spline function is interpolated at the experimental feed particle sizes a fitted feed size distribution is obtained. This feed size distribution is then used as input to the screen model calculations from which product size distributions are generated. Since both feed and product size distribution data are used to minimise the residual errors between the observed and fitted stream size distributions, a mass balance is Conducted simultaneously with the screen model parameter fitting.

Jan 1, 1984