Search Documents

By Pierre Roy, Claude Gosselin, Maryam Sadeghi, Michel Bourassa, Christian Rochefort, Claude Bazin, Gautier Mahieu, François Lavoie, Marilène Renaud, David Cataford

"Spirals are gravity concentrators used for the concentration of iron oxide minerals from ores. This paper presents the mineral size recovery curves observed for spirals operating in two Canadian iron ore processing plants. The size recovery curves are presented for rougher, cleaner and re-cleaner spirals. The partition curves observed for the spirals are significantly different than those observed for typical gravity concentration devices such as hydrocyclones or hydraulic classifiers. The size recovery curves show the complex behaviour of the minerals as they flow down a spiral and are analysed with a view toward process modelling and simulation. Initial simulation results show that gains in the iron recovery can be achieved by increasing the recovery of coarse iron oxide particles in the rougher spirals.INTRODUCTION Spiral concentrators are gravity based separation devices used for the preparation of coal, iron and different heavy minerals ores (Wills, 1992). The classification of the particles in spirals is based on the combination of gravity, centrifuge, drag, lift and Bagnold forces (Atasoy and Spottiswood, 1995; Burt, 1984; Bouchard, 2001). Several papers are dedicated to the analysis of the forces acting on the particles flowing down a spiral (Das and al., 2007; Sivamohan and Forssberg, 1985; Holland-Batt and Holtham, 1991; Holtham , 1992). Most of these papers aim at the development of a general model for the simulation of the classification of the particles for coal ore processing. Few papers deal with the operation of spirals for processing iron ores (Hyma and Meech, 1989; Srivastava, 2001; Dallaire et al., 1978) and fewer papers present results for spirals used as cleaners following a rougher stage. Since the classification in spiral classifiers is dependent on the particle size and specific gravity a partition or performance curve commonly applied to classifiers (Mai Manga et al., 2012) is used here to analyse spiral operation. Few papers were found to describe the operation of spirals in terms of partition curves (Atasoy and Spottiswood, 1995; Holland-Batt and Holtham, 1991; Richards et al., 2000). Most of the investigations presented the application of a model to simulate the position of coal particles of various sizes and specific gravities along the trough of a spiral (Kapur et al., 1996; Holland-Batt and Holtham, 1991). In addition most of the available literature on spirals deals with one stage of concentration by spirals (Dallaire et al., 1978; Das and al., 2007; Sivamohan and Forssberg, 1985) while industrial circuits usually consist of rougher and cleaner spirals."

Jan 1, 2014

By Hossein Aminian, Daniel Hodouin, Claude Bazin, Claude Jacob

"Solvent extraction followed by electrowinning (SXEW) is an economical option for the processing of low grade and oxidized copper ore. Phenomenological models were developed to simulate the copper solvent extraction and the electrowinning processes. The two models were linked together to make a SX/EW simulator. The simulator is used to predict the operation of the mines Gaspe pilot plant put in operation in 1996. Results confirm the potential of the simulation that would subsequently be used for student training, process optimization and to assess the performances of control strategies.INTRODUCTIONSolvent extraction followed by electrowinning (SXEW) has become in the last two decades a key process for the recovery of copper from solution obtained by leaching low grade copper ore. Several operations have been commissioned in North and South Americas in the last decade (Arbiter, 1994). In Canada, Gibraltar Mine (Scott, 1991) operates a SX plant followed by electrowinning (EW) to produce copper cathode from leached low-grade ore dumps. From 1996 to 1997 Mines Gaspe operated a pilot plant to test the economical feasibility of leaching an oxidized copper ore, followed by solvent extraction and electrowinning of copper. Despite copper SX/EW processes have the reputation of being low cost operations, plant engineers should always look at way to make the plant more efficient. Process models provide an economic platform to learn, control and ultimately optimize the operation of a plant.Several models have been independently proposed to simulate the SX process (Cadman and Hsu, 1970; Medina Gomez, 1980; Wilkinson and Ingham, 1983; Garcia et al., 1988; Aminian et al., 1997) and the EW process (Bryson, 1981; Barton and Scott, 1992; Barton and Scott, 1993). However few papers, if any, deal with the integration of the two process models into a single simulator of a SX-EW plant. This paper describes the development of an integrated simulator for the SX-EW process and shows some steady-state simulation runs of an actual SX-EW pilot plant. The paper consists of four sections. The first section shows the general layout and principles of operation of a SX-EW plant.The second section describes the SX model while the third section deals with the process model used for the copper electrowinning process. The last section shows how the two process models are linked to simulate the operation of the Gaspe pilot plant."

Jan 1, 1999

By B. K. Mishra

The motion of the charge inside a grinding mill is of considerable practical importance. The dynamics of grinding mills have been the subject of study by a number of workers because it can accurately predict: i) the power drawn by the mill, ii) the distribution of forces inside the mill, and hence the wear rate of balls as well as the liners, and iii) mass transport of milled product through the porous ball charge. Some of the earlier analyses of the motion of a charge of particulate material in a tumbling mill were reported by White (1904), Davis (1919), and Dyer (1945) and have been reviewed by Rose and Sullivan (1958). In these works, frictional force at the wall of the mill and between the colliding balls, and the interaction of different ball trajectories, was not considered. Moreover, it was assumed that the profile of the tumbling charge is independent of the degree of filling, which is generally not true. However, all this work led to the identification of so-called "equilibrium surface" which is the plane that separates the rotating charge into two regions. The inner region is static; here there is no motion relative to the mill shell. The other region is a shear region where either the balls flow down the surface of the charge (cascading motion), or are projected out of the surface and follow approximately parabolic paths (cataracting motion).

Jan 1, 1990

By Tomohide Takeyama, Yuki Yoshida, Masaki Kitazume

"Abstract A series of centrifuge model tests was conducted to investigate the failure pattern and stability of cement treated clay ground having an unsupported vertical slope. The model ground was subjected to monotonically increasing centrifugal acceleration until it failed due to the enhanced self-weight. The test results were compared with the numerical analyses in which the treated soil was simulated by the Smoothed Particle Hydrodynamics (SPH) method. In the simulation, the Dracker- Prager model is employed as the constitutive model. In order to simulate the brittle behavior of the cement treated soil, it is assumed that the unconfined strength is gradually reduced to a certain value which corresponds to the residual strength due to accumulated plastic strain. It was found that the simulation by using SPH method can express the brittle failure observed in the centrifuge model tests.IntroductionBecause of economical and environmental concerns it has become more difficult in recent years to obtain enough soil to construct reclaimed land for airports, electric power plants, manufacture plants and so on. On the other hand, it has become difficult to find and construct disposal sites for waste subsoil or dredged soil. These circumstances have encourage geotechnical engineers at present to use these poor quality materials as filling materials more than they were used one or two decades ago. Several soil admixture methods have been developed to improve their properties and to manufacture high quality fill materials.The deformation and strength characteristics of cement treated soils have been investigated in Japan, which have revealed that the treated soil is characterized by a brittle behavior with quite high compressive strength and quite small strain at failure, and the treated soil with relatively high strength shows tensile failure instead of shear failure. Regarding to evaluate the stability of structures adjacent to or buried in the treated ground, Kitazume (1998) and Kitazume et al. (2003) investigated the earth pressure of cement treated ground by the centrifuge model tests. And they emphasized that the failure pattern of the cement treated ground was much different from that of ordinary ground. After the Hyogoken- Nambu earthquake in 1995, many design codes for the seismic stability were revised in Japan, so that seismic stability must be evaluated more precisely. The dynamic earth pressure of ordinary grounds is usually calculated by a pseudo-static limit equilibrium method based on the Mononobe- Okabe’s method (Okabe, 1924; Mononobe and Matsuo, 1929). This calculation method is also applied to the cement treated ground, although its applicability has not been thoroughly clarified."

Jan 1, 2014

By Wilhelm Brandstätter, Johannes Steiner, Alois Triessnig, Christian Redl

"The results of numerical simulations of bubble detachment from an orifice in a horizontal plate submerged in liquid are presented. The bubbles are created by forcing gas through this orifice into stagnant liquid. The simulations are performed with a two-dimensional incompressible Lattice Boltzmann (LB) scheme. In previous works the same problem was investigated by two-phase simulations in which surface tension effects were considered by an interaction potential formulation between the phases. In contrast, in the present work a volume of fluid (VOF) model consisting of a carrier fluid and a passively advected scalar function for the calculation of gravitational and surface tension forces, respectively, is used. The aims of this study are to examine the applicability of the one-phase LB method to the numerical analysis of bubble formation and to verify the theoretical and experimental predictions of the dependencies of the departure diameter on gravitational forces, surface forces and wettability effects. To this purpose the simulations were performed at different values of gravitational acceleration g, surface tension coefficient s and contact angle ß. Nonlinear fits to the data yield functional dependencies which are in good agreement with the predictions. IntroductionBubble formation behavior is an extensively studied subject, especially in chemical engineering, power industries and metallurgical processing. Knowledge of the underlying physics of bubble motion is, for example, required for the design of chemical reaction towers and agitation tanks. Another important application is the argon gas injection through the wall of the tundish nozzle into liquid steel stream to deter clogging with solid inclusions. Argon bubble motion and its effect on flow depends greatly on different parameters, e.g. gas flow rate and size of the bubbles. Thus, if these bubbles are not managed properly, they can become encapsulated upon solidification and result in undesirable impurities in the material [1, 2]."

Jan 1, 2004

Carbon anode cells were studied by two-dimensional computer simulation to evaluate their suitability for alumina electrolysis at high superficial current densities. The most efficient anode arrangement was needed to overcome the electrochemical disadvantages of alumina's low solubility in chlorides. Cells with packed-bed anodes might provide practical superficial current densities greater than 0.67 A/cm2 when the granule size was less than 1 cm, but in the process of electrolysis the granules within the active zone would be consumed and become smaller, affecting alumina concentration transport. For zones of very small granules a concentration-dependent model should be considered; for more complete consideration a time-dependent model can be introduced to describe the variation of granule size and its influence on mass transfer and polarisation

Jun 19, 1905

By Donald Leroux, Colin Hardie

"New or improved mathematical models of mineral processing machines are published almost on a monthly basis. Often, these models can be calibrated from bench-scale test results and programmed using spreadsheet software to simulate different designs and operating conditions in open-circuit configuration. However, since mineral processing plants frequently operate comminution, separation and/or extraction machines in closed-circuit configurations, open-circuit simulations often have limited applicability.The simulation of closed-circuit configurations necessitates the programming of an iterative loop to bring the circuit to equilibrium. This can be done by using spreadsheet software but often requires advanced programming skills. As an alternative, LIMN® software, which simplifies the programming of closed-circuit configurations within spreadsheets can be used. This software was specifically designed to turn a Microsoft EXCEL® spreadsheet into a closed-circuit modeling and simulation platform with minimal programming.During the past few years, the Mineral Processing Group at the Noranda Technology Centre, has used this software for numerous projects which include the simulation of multi-stage cyclone circuits, closed-circuit regrinding and locked-cycle flotation tests. Using these projects as examples, the flexibility and the easy to use nature of the LIMN® software is demonstrated.INTRODUCTIONMineral processing engineers routinely have to make predictions regarding comminution and separation processes. Predicting the appropriate dimensions of equipment and the process response to changes in operating mode are probably the most common situations that they face. In all cases these predictions are expected to be rapid and accurate.Predictions made by engineers are frequently based on mathematical models. Many first principles, empirical, semi-empirical, steady state and dynamic models of mineral processes can be found in the literature. In fact, new or improved mineral processing models are continuously being proposed and in spite of this variety, process engineers often have to either develop their own or tailor existing models to suit their particular needs."

Jan 1, 2003

By J. P. Gu, C. Beckermann

"Melt convection and macrosegregation in casting of a large steel ingot is numerically simulated. The simulation is based on a previously developed model for multicomponent steel solidification with melt convection, and involves the solution of fully coupled conservation equations for the transport phenomena in the liquid, mush, and solid. Heat transfer in the mold and insulation materials, as well as the formation of a shrinkage cavity at the top are taken into account. The numerical results show the evolution of the temperature, melt velocity, and species concentration fields during solidification. The predicted variation of the macrosegregation of carbon and sulfur along the vertical centerline is compared with measurements from an industrial steel ingot that was sectioned and analyzed. Although generally good agreement is obtained, the neglect of sedimentation of free equiaxed grains prevents the prediction of the zone of negative macrosegregation observed in the lower part of the ingot. It is also shown that the inclusion of the shrinkage cavity at the top is important in obtaining good agreement between the predictions and measurements.IntroductionMacrosegregation in large steel ingots is one of the most well-known and now classical problems in the field of solidification and casting [1]. Commonly found macrosegregation patterns in steel ingots are a positively segregated zone at the top, negative segregation near the bottom in the equiaxed zone, inverse segregation near the ingot surface, V-segregates along the centerline, and A-segregates in the columnar zone [2, 3]. It is well established that the macroscopic transport of species leading to macrosegregation can occur by a variety of mechanisms, including thermally and solutally driven natural convection of the melt in the mushy zone, flow due to solidification contraction, and the sedimentation of free equiaxed crystals. Fredriksson and coworkers [ 4,5] provide recent overviews of the solidification transport phenomena and macrosegregation in ingot casting."

Jan 1, 1999

Conveyor belt haulage offers the advan tages of continuity of operation, minimum labour requirements, and the ability to handle large volumes of material over extended periods. However, breakdown of critical con- veyors in the system can cause undue and substantial losses in production output at the working faces, unless adequate surge bin capacity is incorporated in the haulage system. A conveyor belt simulation programme, code name 'Beltsim' was developed by Bucklen et al (1968) at the Virginia Polytechnic Institute in the U.S.A. and provides a quantitative measure of the performance of a conveyor belt network design. A limitation still exists with the computer programme and involves the inability of the programme to predict the necessary surge bin capacity to overcome conveyor belt downtime. The computer programme was modified to overcome this limitation, to simulate a con- veyor belt breakdown, and to provide a quantitative prediction of the surge bin capacity necessary to permit a modern mining system to operate effectively and continuously.

Jan 1, 1977

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 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 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

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 Ji Hoe Kwon

Hindered-settling column process has been used for the separation of mineral particles for many years. In this study, a model based on DEM (Discrete Element Method) was used to simulate the movements of particles in a settling column separation process. DEM model is capable of describing the motion of particles and calculating the impact forces between particle and particle, and particle and wall contact under various operating conditions. Drag force, rotation induced Magnus lift force, contact force and the gravitational force acting on individual particles were considered for establishing the mathematic models. The overall motion of the particles in water is satisfactorily simulated by the DEM, and the result shows good agreements with the result of lap-scale experiment.

Jan 1, 2014

By Harushige Kusumi, Yuzuru Ashida, Toshiki Watanabe, Toshifumi Matsuoka

ABSTRACT: The technique of various kinds of numerical analysis offered remarkable results to explication of many rock mechanical phenomena. These techniques regarded a rock stratum as a continuum body and solved the governing equation of continuum mechanics by using finite differ-ence method, finite element method and/or others. Using numerical analysis based on continuum mechanics, however, the simulations of destruction phenomena caused by elastic waves and on the other hand, the elastic wave propagation which generated by fracturing events are quite difficult. In this paper, we express a rock stratum by a granular body model and develop a simulation tech-nique can simulate both wave phenomena and destruction phenomena at once. The DEM (discrete element method) simulation results are compared with finite difference solutions of elastic wave equation. Furthermore, we simulate Hopkinson effect by this proposed method. Hopkinson effect is famous as a destruction phenomenon caused by elastic wave propagation. ABRISS: In dieser Studie drücken wird eine Felsschicht durch ein Granularkörpermodell aus und entwickeln eine Simulationstechnik, die sowohl Wellenphänomene als auch Zerstörungsphänomene simulieren kann. Die Ergebnisse der Simulation mit DEM (Diskrete-Elemente-Methode) werden mit denen der Finit-Differenz-Lösung von Gleichungen elastischer Wellen verglichen. Desweiteren simulierten wir den Hopkinson-Effekt mit dieser vorgeschlagenen Methode. Der Hopkinson-Effekt ist als Zerstörungsphänomen berühmt, verursacht durch Fortpflanzung elastischer Wellen. RESUME: Dans cette étude, nous exprimons un lit de roches par un modèle à corps granulaire et développons une technique, permettant de simuler d'un coup le phénomène des ondes et le phé-nomène de destruction. Les résultats de la simulation DEM (méthode d'éléments discrets) sont comparés aux solutions à différence finie de l'équation d'ondes élastiques. En outre, nous simulons l'effet d'Hopkinson par la méthode proposée. L'effet d'Hopkinson est connu comme phénomène de destruction, causé par la propagation d'ondes élastiques.

Jan 1, 2003

By Alexander Siefert, Steffen Pankoke

In the development of commercial vehicles as well as of passenger cars, the effects of vehicle vibrations on operating safety, health and comfort can only be predicted by numerical simulation when appropriate occupant models are available. Such models must be based on human anatomy and have dynamic properties of real humans in order to achieve realistic results. Since human dynamic behavior depend on posture and percentile, the occupant model needs to be adjustable to these parameters with respect to geometry and dynamic properties [1,2]. Dynamic Human Finite-Element-Model CASIMIR CASIMIR is a non-linear, dynamic finite-element-model of the human body. It consists of a dynamic model of the upper torso with head, neck, shoulders and arms as well as of a dynamic model of the lower extremity with pelvis and legs. The most important part is the lumbar area with dynamic non-linear models of the lumbar spine and of back and abdominal musculature. The frequency-dependent characteristics of the intervertebral discs and the effects of muscle activation and non-linear frequency-dependent muscle properties are included. In the latest stage of development, CASIMIR has been equipped with a compliant model of the body surface in the contact areas to the seat. This results in a very realistic transmission of static and vibrational forces into the human body, see fig. 1. Intense model verification and validation has been performed in all stages of model development, starting with validation of small components like intervertebral disc, ending with validation of whole-body-vibrations using measurements of the dynamic mass / mechanical impedance [4]. For an in-detail examination of stresses in the vertebral bodies and discs, a non-linear submodel of the lumbar spine with an increased number of degrees of freedom can be coupled to the whole-body-model, enabling the researcher to examine local effects of vibrations and single shocks on the lumbar materials. [ ] Since it is well known that human dynamic behavior is significantly affected by anthropometric data and posture, CASIMIR can be individualized to the anthropometric status of single individuals or to the mean values for specific percentile groups. Furthermore, posture can be adjusted to the seating conditions applicable to a specific vehicle. Posture modification capabilities include the variation of the lumbar lordosis [3].

Jan 6, 2006

By Kenji Matsuzaki

A mathematical simulation model of the Imperial Smelting Furnace operation has been developed for purpose of improving the furnace operation. The model consists of top, shaft and bottom regions, Analysis in the shaft region is based on differential heat and materials balances taking account of reaction kinetics and heat transfer rates, and those in the top and the bottom regions are based on heat and materials balances under a steady-state condition, Longitudinal distributions of process variables such as temperatures of gas, sinter and coke, gas composition and fractional conversion under a practical operating condition were calculated, With this model, the furnace performance under various operating conditions also can be shown reasonably predicted,

Jan 1, 2003