Search Documents

By C De-Vitry

"The large Yandi channel iron deposit (CID) is located in the north-west of Western Australia and is currently mined via open pit methods. Drilling at Yandi is widely spaced and current estimates are mostly linear averages of the surrounding local samples. Linear estimates such as ordinary kriging (OK) based on block sizes that are significantly smaller than the dimensions of the drilling grid are less variable than the unknown reality (Armstrong and Champigny, 1989).Conditional simulation (CS) provides a series of plausible realisations that represent the actual grade variability at a small scale (eg 1 m + 1 m + 1 m). These simulations are very useful for a number of purposes, eg: optimising grade control, stockpiling and blending selection of optimal drill hole spacing assessing resource risk on contained tonnage and metal.At Yandi the major constituents of the ore (Fe, P2O5, SiO2, Al2O3 and LOI) are correlated to varying degrees. Independent CS of each assay variable will not reproduce these correlations. Reproducing the correlations (and the shape of scatterplots) is important if the simulations are used for the abovementioned purposes. This paper examines and compares various methods of CS that aim to reproduce the correlations between the assay variables.The CS methods evaluated in this paper are conditional co-simulation (CCS) using the linear model of co-regionalisation and two other, alternative multivariate simulation methods, which use min/max autocorrelation factors (MAF) or the stepwise conditional transform (SCT) to de-correlate assay variables before performing independent CS. It was found that CCS best reproduced the characteristics (histograms, variogram, scatterplots, etc) of the declustered drill hole data (Figure 1). Whilst simulation using MAF generally produced acceptable results, the method did not adequately reproduce the variance of the assay variables. The univariate statistics, histograms and experimental variograms were acceptable for simulations generated from the SCT. Bivariate statistics, cross-variograms and scatterplot reproduction for the SCT were unacceptable. The latter method produced inferior results because there were insufficient samples to obtain a reliable transform of the third assay variable in the transformation sequence.None of the abovementioned simulation methods provides good reproduction of all the scatterplots. This is due to the presence of clay pods, which cause some scatterplots to display two trends. These clay pods contain higher percentages of SiO2 and Al2O3 than the surrounding CID. The clay pods are generally smaller than the drill-grid dimensions. Thus their locations are typically unknown and they cannot be defined by manual interpretation, eg via wireframes. The locations of the clay pods were simulated via indicator CS. CCS was then performed separately for the clay pods and the CID. This significantly improved the reproduction of the correlation coefficients, the scatterplots and visually the simulation was more geologically realistic. Simulating the location of clay pods followed by CCS of assay variables is the recommended approach for simulation at the Yandi CID.CITATION:De-Vitry, C, 2013. Simulation of correlated assay variables - a case study from the Yandi Channel iron deposit, in Proceedings Iron Ore 2013, pp 137-148 (The Australasian Institute of Mining and Metallurgy: Melbourne)."

Aug 12, 2013

By H S. Mitri, A Sainoki

Destress blasting techniques have been widely employed in underground mines for the purpose of preconditioning highly stressed rock mass. Appropriate applications of the technique successfully contribute to increasing mine safety (eg mitigating the risks for rock bursts). This paper presents a simulation technique to estimate the extent of damage zones induced by destress blasting. A 3D numerical model with a single blasthole is constructed with FLAC3D code, whereby static analysis is carried out in which horizontal and vertical stresses are applied to the model in order to simulate stress conditions in a deep underground mine. Afterwards, dynamic analysis is performed whilst applying time-varying blast pressure along the blasthole. The blast pressure profile takes into consideration an instantaneous increase and gradual decrease in blast pressure as well as the propagation of detonation along the blasthole. In this paper, constitutive models allowing for shear and tensile strengths of rock mass that are dependent upon strain rate are newly implemented into FLAC3D code with C++ programming language, whereby blast-induced damage is assessed by the extent of tensile and shear failure zones. The results obtained from the dynamic analysis show significant decrease in the extent of yielding zones induced by destress blasting with increasing the magnitude of in situ stress. They indicate that considering in situ stress state is indispensable in order to propose proper blast design for each level in underground mines. Based on these results, the relationship between the extent of the damage zone and depth at which destress blasting is performed is established. Furthermore, stress dissipation induced by destress blasting is examined with static analyses, in which the modulus of elasticity is decreased according to a rock fragmentation factor, a, only for the area where tensile and shear failures takes place during the dynamic analyses. A model parametrical study with respect to a is undertaken. The relationship between the rock fragmentation factor and the stress dissipation in the yielding zones is developed from the parametrical study. The developed relationship can be useful when destress blasting is designed for massive rock masses.CITATION:Sainoki, A and Mitri, H S, 2015. Simulation of damage zones induced by destress blasting, in Proceedings 11th International Symposium on Rock Fragmentation by Blasting, pp 245–256 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Aug 24, 2015

By R. Noorani, M. Cai

"Excavation of deep underground openings in highly stressed grounds can generate new fractures near the openings. As a result, rock mass bulking may occur, leading to large wall deformation. A better understanding of brittle rock failure and the associated fracturing process at the laboratory scale is an essential step towards evaluating dilation behavior of rocks. A grain-based model that can capture dilation behavior of rocks under various loading conditions at the laboratory scale can be a useful approach to study rock mass bulking at the field scale. The aim of the paper is to study the dilation behavior of rocks at the laboratory scale employing numerical simulation. A modified grain-based model is used to simulate the macroscopic response of crystalline rocks by considering micro-fractures at the grain scale. Grains generated with the Voronoi tessellation are allowed to fail in tension or in shear using a strain-softening model. Damage is represented by either grain boundary fracture or grain failure, or the combination of the two. As stress increases, the process of damage, fracture initiation and coalescence in rocks can be captured.The grain-based model was calibrated using test data of the Lac du Bonnet granite and a parametric study was conducted. Numerical simulation results showed that the proposed modeling approach is able to capture most of the mechanical properties, fracture pattern, and dilation behavior in both unconfined and confined loading conditions. In an unconfined compression test, interactions of microscopic fractures may lead to macroscopic shear or axial splitting failure depending on the chosen micromechanical parameters. Numerical experimental results show that the dilation behavior of the model is strongly related to the failure mechanism of rocks. Tensile-failure-dominated failure mode generates large dilation due to the formation of new cracks that are sub-parallel to the loading direction. On the other hand, shear-failure-dominated failure mode produces small volume increase. One significant finding from this study is that in order to capture large dilation using grain-based UDEC models, intragranular grain failure should be considered."

Jan 1, 2015

By P. K. Chatterjee

The overall success of many strip coal mining operations depends primarily upon the efficient use of draglines to remove overburden. These machines require enormous capital investment and unless used efficiently the profitability of a venture is adversely affected. It is important that the dragline works continuously and it is imperative that the walking path and other facets of its operation be predetermined exactly so that the desired production rate is achieved with a minimum of non-productive time. Details of dragline operations which are of particular interest are walking path, prediction of best cut and positioning to maximize production, the angle of swing which produces the swiftest operation, optimum fill factor allowable for the type of material and conditions of ground, the maximum depth from which it can excavate and operating level with respect to the elevation of the high wall. The aims of the present paper are: 1. To establish a cut-off point in the cost-depth ratio. 2. Predict the most economic method of stripping in a given area using graphical plots of walking path implemented with information necessary to ensure optimum production. The optimum operating efficiency of the dragline in relation to the terrain and other conditions imposed on its operation will be evaluated. The model is based upon operating data which has been collected from site observations of a large strip coal mining operation in Central Queensland, Australia.

Jan 1, 1977

By A. Gakwaya, P. Vachon-Joannette, P. Gagnon, J. Lévesque

"Predictive numerical methods are essential tools for reduced development cost and improved high performance designs of modern aircrafts. Their use leads to savings in terms of design cycle time and experimental trials. Composite materials subjected to impact loadings demonstrate various damage modes such as matrix cracking, fiber breakage, delamination and, ultimately, fabric fracture. Composite laminates are extensively used in the aerospace industry, and their damage modeling is complex due to their inhomogeneity and anisotropic behavior. Damage evolution affects both in-service properties and performance of these composites. In this work, the drop-weight impact test has been simulated. Using a material model that accounts for both delamination and intralaminar damage, laminates response is analysed using numerical simulations based on 3-D finite-element models of the commercial code Abaqus/Explicit. Experimental and simulation results are compared, in terms of energy absorption and contact force. The obtained numerical results show good agreement with experimental data.INTRODUCTIONAs primary aircraft structures are increasingly being made from advanced fiber-reinforced polymer matrix composites and are in general exposed to a variety of impact scenarios (bird, hail, debris etc.), an understanding of the impact event and the performance of composite materials during impact is important in the development of new structures. This work is focussed on low velocity impacts that typically occur when foreign objects hit the aircraft during maintenance or other on-the-ground operations. These impacts can cause non visible damage that affect the structural integrity of the composite components. While for metallic structures, damaging impact events typically leave noticeable dents due to plastic deformation, in contrast, for high strength carbon fiber reinforced composites, they can deform elastically to strain levels exceeding 1% and the internal damage can often occur without the formation of a visible dent. The identification and effect of impact damage on the performance of the structure is thus of great concern to the aircraft industry as the durability and damage tolerance of these composite structures are strongly influenced by the size of detectable internal damage."

Jan 1, 2016

By A. Ockewitz, A. Trondl, D. -Z. Sun

"Numerical simulation of extrusion processes offers an efficient method for tool design and optimization of geometries and properties of extruded parts. However, until now there are few systematic investigations about influences of material model and different parameters e.g. for friction and thermal processes on calculated local and global responses. Since extrusion processes are essentially influenced by friction between billet and extrusion tools, new friction tests were performed in this work to develop a reliable friction model. Compression and torsion tests were performed at different temperatures and velocities to calibrate the material model for extrusion simulations. Extrusion processes for circular rods and U-profiles were conducted at TU Berlin with backward and forward extrusion for different extrusion ratios, billet temperatures, product velocities and cooling conditions. This work focuses on numerical simulations of the extrusion tests with the FE code HyperXtrude. Based on the calculated distributions of local strain, strain rate and temperature the recrystallized volume fraction and grain sizes were calculated with empirical models and compared with the experimental results from optical micrographs and EBSD scans. The calculated total and friction forces and temperatures were compared with the experimental results and the deviations between experiment and simulation were analyzed.INTRODUCTION Extruded aluminum components are widely used in vehicle constructions due to light weight and manufacturing consideration. Local mechanical properties in the components are inhomogeneous as a consequence of the spatial distribution of microstructure which is influenced by numerous factors like material composition, temperature, tool design, friction, extrusion ratio and speed during the extrusion process. For the optimization of the extrusion process with regard to all these factors numerical simulation can be an essential tool. Based on results of finite element simulations of the extrusion process the microstructure in the extruded components can be predicted and then the resulting mechanical properties can be numerically determined instead of by costly trial and error methods. The aim of the work presented in this paper is the development of a numerical method for the simulation of extrusion processes with modeling of microstructure."

Jan 1, 2018

By S. Malinov, J. J. McKeown, W. Sha, S. McShane

"A model is developed for prediction of fatigue stress life S-N diagrams for Ti-6Al-4V alloy at different conditions, using Artificial Neural Network (ANN). Different factors which have influence on the S-N diagram are involved as input parameters of the NN, namely microstructure, environment, texture, test/work temperature, surface treatment and stress amplitude. As the output the S-N fatigue diagram is simulated and plotted.The model is based on standard multilayer feedforward neural network. The NN is trained with a comprehensive data set collected from the literature, using Levenberg-Marquardt training algorithm in combination with Bayesian regularisation. A very good performance of the trained Neural Network is achieved. Good correspondence between ANN"" prediction and experimental S-N diagrams is observed. Using the model S-N diagrams for different conditions are predicted and analysed. Graphical User Interface is developed for use of the model. We believe that the model will be important for practical applications in solving various problems on fatigue behaviour of the titanium alloy.IntroductionTitanium and its alloys have an increasing range of applications due to their very good combination of mechanical properties, corrosion and erosion resistance. They are currently finding widespread applications in many industries due to their desirable and versatile properties. Most important among the alloys is Ti-6AI-4Y. Jet aircraft manufacturers are the principal 'consumers of Ti-6Al-4V and it is often specified for critical parts, the failure of which could result in the loss of an entire system. A generic class of Ti-6AI-4V has been developed over the past years and it has become the best understood of all the titanium alloys. The advantages of titanium alloys include: (i) low densities, giving very attractive strength to weight ratios allowing lighter and stronger structures; (ii) superior corrosion and erosion resistance in many environments; (iii) high temperature capability. Although the titanium materials are still considered as expensive materials, for many applications the cost of titanium alloys can be justified on the basis of their desirable properties."

Jan 1, 2001

By M. Rahman, P. C. Kapur

The fixed-time locked-cycle grinding test is employed to simulate the industrial closed grinding circuit on the laboratory scale. A mathematical algorithm for 'cyclewise' computer simulation of the locked-cycle test has been derived and used to study an idealised laboratory-scale model of a closed ball-mill circuit. The agreement between simulated and experimental results is reasonably good, except when the grinding times are relatively long

Dec 1, 1971

By MáAndrusiewicz, M M. Hilden, M S. Powell, P D. Kay, B Foggiatto

Circuits with unconventional and potentially more complex designs are predicted to become more common in the future as plant designers and operators look to process ores more efficiently. Ideas have been presented on how circuit design may be changed to utilise todayÆs technology to address tomorrowÆs challenges (Powell and Bye, 2009; Powell, Benzer and Mainza, 2011a). High pressure grinding roll (HPGR) installations in various configurations; application of preconcentration plants eg sorting, magnetic or gravity separation; and preweakening devices eg microwave processing; are all examples of technologies having the potential to deliver a significant reduction in comminution energy for ores amenable to these processes. Not only will future novel circuits contain equipment for which models are not as well developed (if at all), many will likely require the ability to respond to variations in feed, and interact with geometallurgical models requiring process adjustments when the ore characteristics change.In order to confidently simulate more complex circuits and accommodate variations in feed type, a significant advancement in modelling and simulation capability is required. This includes dealing with multiple components in the feed, accommodating different competencies of ore, utilising physical separation opportunities and tracking the grade reporting along the different potential processing routes. Developing the means to simulate complex circuits has been a key priority of a collaborative group of comminution researchers working to solve these challenges.Simulations are presented utilising a new configurable simulation interface called the model developers kit (MDK). The MDK is designed to accommodate multi-components and incorporate new models under development. These will be used to demonstrate the potential of flexible circuit design, based on measured ore properties and equipment performance.CITATION:Foggiatto, B, Hilden, M M, Powell, M S, Kay, P D and Andrusiewicz, M, 2012. Simulation of flexible circuits, in Proceedings 11th AusIMM Mill OperatorsÆ Conference, pp 291-298 (The Australasian Institute of Mining and Metallurgy: Melbourne).

Oct 29, 2012

By Sylvain Duchesne, Claude Bazin, Robert Grant

Grinding simulators are now widely available to assist in the selection of circuit operating conditions or flowsheet designs that yield a given product size distribution. These simulators do not inform of the impact of the grinding product size distribution on subsequent flotation stages. This paper presents a method to predict the effect of fineness of grind on a flotation bank operation. Application of the method is illustrated with data from two Canadian concentrators.

Jan 1, 1995

By J Arola, B Kettunen, W M. Bond

Separation of fine and coarse particles as separate streams prior to flotation, and their differential treatments, could substantially increase recovery and reduce collector consumption. A series of simulations was carried out to quantify the benefits from plant scale implementation of such a treatment at the Ispat Inland flotation circuit (Minorca mine). The study involved detailed plant sampling and sample analysis to generate the data needed for simulation studies. Samples were taken from all the streams from flotation feed to final concentrate, as well as cell lips and cell pulp. Sampling was repeated for the second ore blend that was processed by the plant during different periods of a year. Sample analysis involved sieve analysis and size by size chemical analysis for total iron, magnetic iron and silica. Raw data were mass balanced, and flow rates within the circuit were calculated. Mass balanced data were used to calculate model parameters representing the existing conditions. Usim Pac was used for both mass balancing and simulation studies. Available pilot and bench scale data were utilised to define the modified flow sheet and expected performance. The simulated flow sheet included classification of flotation circuit feed by hydrocyclones. While the coarse fraction followed the existing flow sheet, the fine fraction was sent to a magnetic separator. Magnetic and rougher flotation concentrates were combined to form the final concentrate. Two simulations for each blend were carried out. The first one assumed that there would be no change in size by size kinetic behaviours of both gangue and magnetite due to this modification. This was named as the worst case scenario, in which particles had a longer residence time due to reduced flow rates within the flotation cells. The second scenario anticipated that there would be an increase in recoveries of gangue and magnetite going into froth launder due to more efficient use of collector in absence of fines. Results showed that pre-classification of flotation feed could improve the flotation recovery by 0.3 - 0.5 per cent, with the major benefit being significant reduction (40 - 60 per cent) in the amount of collector used.

Jan 1, 2005

By R. P. King

The two-phase distributed-parameter kinetic model for flotation has been wed as the basis for a digital simulation of flotation plants of arbitrary configuration. Ores that contain any number of minerals can be handled. The effects of incomplete liberation and particle-size distribution are specifically allowed for. The simulator allows the investigation of variations in throughput and in the slurry density of the feed, and of variations in stage volumes and aeration rates. It can handle cell behavior even under conditions of hindered flotation. The requirements for data preparation are completely documented so that the computer program is suitable for general-purpose use. The program is written in Fortran. The output from the simulator consists of the complete material flows in the plant, as well as the grade, recovery, and particle-size distribution in each stream.

Jan 1, 1976

By J. B. Rubinstein, V. J. Badenicov, S. T. Hall

A simulation technique of the flotation process and models for the solution of various problems are developed. A modular model of column flotation is constructed taking into account the kinetic relationship of the process and the hydrodynamic regime in the cell. Based on this, a technique of calculating the column parameters based on laboratory tests is developed (scaling-up), and applied for calcetest results is developed (scaling-up), and applied for calculating the design of commercial-scale columns.

Jan 1, 1996

By A. M. Engelbrecht

Paper written on project work carried out in partial fulfilment of B Eng. (Chemical) Membrane processes are vital to a diverse range of technical and biological systems, which also feature in mineral processing applications such as the treatment of mine effluent and extraction of metals. In this study, the use of cellular automata to model a semi-permeable membrane that separates two compartments is considered. Cellular automata are simple models that enable the analyst to represent well-understood local effects, such as the chemical or steric interaction between molecules themselves and between molecules or particles and their environment, which can lead to very complex behaviour on a macroscopic scale. Various physical phenomena, including reverse osmosis, were simulated. Initial results for the effect of temperature, membrane structure and chemical properties of the solutes on the performance of the membrane appear promising. Keywords: cellular automata, dynamic modelling, filtration, membrane fouling, osmosis.

Jan 1, 2007

By M. W. Hildyard, J. A. L. Napier

The improved design of support systems and stoping strategies when changes in depth, stratigraphy, or mine layouts are considered must rely on the designer's fundamental understanding of the deformation mechanisms adjacent to underground openings so that he can extrapolate and predict expected mining conditions. This paper presents a numerical approach to the modeling of the growth and interaction of mining-induced fractures near the edges of deep tabular excavations in brittle rock. Fractures are represented by displacement-discontinuity elements that are allowed to propagate according to prescribed growth rules. This procedure is justified by the reproduction of several features observed in laboratory simulations of fracturing around openings in blocks of rock subjected to uniaxial and biaxial compression. An attempt is then made to explain some of the attributes of the fracture zone surrounding tabular stopes. Particular questions of interest are the sequence in which fractures form and are mobilized in front of the mining face, the effect of parting planes parallel to the plane of the excavation, and the distribution of stress parallel to the hangingwall and footwall.

Jan 1, 1992

By Alan P. Prosser

The heap leaching of gold ores by cyanide solution was simulated by computer program. The program uses only basic m1neralogical and chemical data; no process testwork data, such as column leach results, are required. The simulation can be used to forecast the consequences of various changes before and during the testwork, design and operation stages of a particular ore-process system. Also, data about the state of different parts of the heap at various times help create an understanding of what is happening. The effects of mineralogical and process variables on the overall performance of the ?heap leaching process are illustrated. Comparisons are made to column leach tests and production data. The potential of the simulation strategy is discussed for making better use of operating experience and for suggesting specific research objectives. It has become clear that more attention should be paid to other metals that dissolve in the solution.

Jan 1, 1989

A universal grinding circuit model suitable for an arbitrary structure of flowsheet has been suggested. In this model classifiers and mills are described in a most general form. What puts the model in a class by itself is the automatic mass balancing in all streams of the circuit. In order to provide this a special classifier model is adopted. A previously published continuous steady state process model is used as a grinding model. Non-linearity of classifier model leads to a resultant circuit model which is equivalent to a non-linear integral equation. It complicates somewhat calculation algorithm in comparison to that of linear model, but provides a physically meaningful solution.

Jan 1, 1993

By B. Sonmez, B. Aksanl, S. Ersaym, I. F. Erkal, S. L. Ergun

The comminution circuit at the plant has experienced severe problems arising from a feed containing much larger amounts of clay and having a higher moisture content than expected. Three solutions were proposed: replacement of the crushing circuit by tooth roll-type crushers; using washing drums; and introducing primary autogenous/semi-autogenous millling instead of crushing. Computer simulations studied the effects of the first two solutions on the circuit capacity. Laboratory and industrial scale ball mill batch grinding tests were carried out. The parameters of a kinetic model were determined. The ball mill circuit was sampled and operating model parameters were calculated. The ability of batch grinding parameters to represent the operating conditions was studied. Tests failed to find acceptable model parameters, but parameters from an industrial size batch mill were found to be suitable for simulation studies. The simulations indicated that a continuous feed of -12 mm provided by tooth roll crushers would maintain the desired capacity. Increasing the feed size to -25 mm would reduce the capacity by 8-15%. The use of washing drums could also solve the capacity problem, but small fluctuations in water flow rates might result in problems of excess water in the grinding circuit.

Jan 4, 1993

By L. G. Austin, P. T. Luckie, R. R. Klimpel, R. S. C. Rogers

1. INTRODUCTION The Bond method(1.2) for sizing tumbling ball mills has been used successfully for many years. The primary purpose of the Bond calculation is to predict the mill size and mill power to give a desired capacity (Q tph) from a wet overflow ball mill in normal closed circuit operated with a circulation ratio of 2.5. The input to the calculation is the make-up feed size characterized by the 80%-passing size in micrometers (xG), the desired circuit product size characterized by the 80%-passing size (xQ), and an empirical grindability number (the standard Bond Work Index) determined by a standard test. A series of empirical correcting factors(2) are used to convert the results to open-circuit or dry grinding in grate discharge mills. The Bond method has two major engineering advantages: (i) it is very simple; (ii) experience has shown that it works for many circumstances, as a reasonable approximation, because it is based on direct conversion of test data to results measured on industrial mills. There are four main disadvantages of the Bond method. First, there will be a range of error for any specific mill and material and set of operating conditions because the method does not explicitly include several variables which are obviously important: (i) recycle ratio and classifier efficiency (ii) mixture of ball sizes in the mill (iii) variations of residence time distributions with mill geometry and slurry density (iv) lifter design (v) slurry rheology (vi) different underfilling or overfilling of the mill as flow rate is changed. Second, the Bond method explicitly assumes that specific grinding energy E is not a function of ball load, whereas industrial practice and laboratory tests(3) show that the specific grinding energy (to go from a specified feed to a specified product) is less for lower ball loadings than the ball loading for maximum mill capacity. Third, the method uses only the 80% passing size of circuit feed and product to characterize the size distributions, whereas mill capacity in general must depend on the shape of the feed size distribution and the product size distribution. The prime example of this is the use of the reverse closed circuit which removes some fraction of the make-up feed before it is

Jan 1, 1982

By Yusuf R

Many heap and dump leaching processes involve competitive, complementary and sequential reactions between multiple minerals and reagents. A comprehensive model has been developed to deal with these complex chemical interactions which is based on diffusion control of the reaction rates. All of the mineralogical and chemical data are obtained from representative samples of the ore or the chemical literature. These properties have been incorporated in a general way, making the model applicable to a range of ore-process systems. The simulation is used to predict the absolute reaction rates and can be used before any leach tests are carried out on the ore itself. The applicability of the model and some tests of its validity are presented.

Jan 1, 1984