Search Documents

By Guo Wenbing

"INTRODUCTIONThere is a great amount of coal (about 140 billions tons) left unmined under the surface structures, water bodies, railways (referred to as the “three-body”). Coal mining under “threebody” especially that under surface structures has become a major problem in the mining area. The key problem of mining under surface structures is to control overburden strata and surface movements. At present, the major surface subsidence control methods are strip pillar mining, backfilling mining method, room and pillar mining, grouting of bed separations, harmonic mining, and so on. Wongawilli mining method is a high-efficient shortwall mining method developed based on the room and pillar mining, and named after the first successful trial in Australia “Wongawilli” coal seam. The room and pillar mining equipment such as continuous miner is employed in this method. Its major advantage is flexible face layout. It can be used to mine the small wedge coal blocks and those that are difficult to be mined with retreat longwall mining method. It requires less capital investment, short lead time for production, flexible equipment operation and face move and high productivity. It has been applied in some China’s coal mines.The “Wongawilli strip coal mining technology” is a new coal mining method for mining under surface structures. It combines the strip mining longwall layout and Wongawilli high efficient mining technology. This technology can fully utilize the advantages of both the strip pillar and Wongawilli mining methods, and overcomes the disadvantages such as frequent face move, and low mining efficiency in strip pillar mining, and poor ventilation condition and poor long-term stability problems of coal pillars. In this paper, the Wongawilli strip mining method was designed for Wangtaipu mine and its feasibility was studied using numerical modeling and surface subsidence prediction.Geological and mining conditionsAt present, most of the No.15 coal seam reserves totaling 377.07 million tons in Wangtaipu mine are located under surface structures. In order to mine the coal under the surface structures, the Wongawilli strip mining method was employed. An appropriate panel was selected for trial and the surface movement monitoring lines were set up. The trial area is located in panel XV2214 (East) ( re 1). In this area, the floor elevation of coal seam is +634~646m, the surface elevation is +810~817.5m, and the average depth of coal seam is about 174m.The average panel length is 282m and average panel width is 73.5m. The seam dips 1~3° and is 1.8~2.5m thick with an average of 2.15m. The roof and floor rock characteristics are shown in Table 1."

Jan 1, 2015

By N. Q. Huy

This paper describes the application of finite element analysis to simulate the static and quasi-static pile load tests in different drainage conditions. The possible effect of excess pore water pressure will be studied from the simulations results. The results show that the drainage condition has effect on the bearing capacity of a pile derived from the quasi-static pile load test. The excess of pore pressure reduces the pile toe resistance while it increases the pile shaft resistance. The factors which influence the effect of excess pore pressure are also discussed.

Jan 1, 2006

By Khaled Morsy

The tragic incidents at Sago and Darby coal mines in 2006 pointed out the need for stronger seals to resist the blast effects of violent gob explosions. Mine Safety and Health Administration (MSHA) has published the final rule on Sealing of Abandoned Areas in April 2008. According to the final rule, the mine seals for unmonitored gob areas should withstand an instantaneous pressure of 120 psi. Reinforced concrete (RC) seal is one of the feasible options that can be used to meet the final rule. Concrete is a heterogeneous, cohesive-frictional material and exhibits complex non-linear inelastic behavior under multi-axial state of stress. The increased use of concrete as primary structural material requires the development of sophisticated material models to have an accurate prediction of the material response to a variety of loading situations. The constitutive modeling of concrete has advanced considerably in the past three decades. Choosing a reasonable concrete material model is important in order to obtain an accurate mine seal design. This paper attempts to address the effect of different concrete material models on the response of RC seal subjected to dynamic loads from methane explosion. Three concrete models, from ABAQUS library, were selected to address the effect of non-linear behavior of concrete and damages due to tensile and compressive stresses. In addition, a typical and verified concrete material model was used for conducting 3D finite element models simulating the RC seal by ABAQUS, EXPLICIT ver. 6.6.

Jan 1, 2008

By T. G. Sitharam, Resmi Sebastian

"Wave propagation studies across jointed rocks have significance in the fields of rock mechanics and geophysics. The passage of a seismic wave induces strains in rocks and often, these strains are of low magnitudes except at the vicinity of source. These low magnitude strains are generally in the elastic range and hence elastic constants necessary for the design and analysis of strength of rock mass, can be derived from wave velocities. Measurement of travel time velocities and measurement of resonance frequencies are the two methods usually adopted for obtaining wave propagation velocities.Resonant column apparatus measures wave velocities from the resonant frequencies obtained during testing. Depending on the type vibration applied, torsional or flexural, shear and compression wave velocities are secured. By varying the amplitude of vibrations applied, various strain levels can be introduced and corresponding resonant frequencies and wave velocities can be calculated. The resonant column experiments were performed on a model material, hardened Plaster of Paris (POP) samples, in laboratory. The paper explains the preparation of samples and the working mechanism of resonant column apparatus along with the estimation of wave velocity from resonant frequency. The numerical simulation of resonant column apparatus using distinct element method and its validation with the experimental results are presented in the paper.A distinct element code is usually used to simulate the movement of materials with intersecting discontinuities. Finite rotations and displacements of individual blocks representing the material are permitted and new contacts between the blocks are identified as the calculation progresses. Three dimensional Distinct Element Code (3DEC), the software intended for the numerical simulation of rock engineering projects, can well simulate the behavior of jointed rocks under static and dynamic loading conditions. The numerical simulation of resonant column apparatus has been performed with 3DEC considering its efficiency in modeling jointed rocks.The numerical modeling of the resonant column apparatus with application of torsional and flexural vibrations to the sample is elaborated in the paper. The response of samples to a range of frequencies is obtained and the frequency that induces maximum displacement in the sample is measured as the resonant frequency, for any particular strain level. The displacements measured at various points of the sample during application of vibrations are compared and presented. The numerical modeling of resonant column apparatus with an intact sample (without joint) and a sample with a joint incorporated are explained at a particular strain level. The resonant frequencies obtained from laboratory experiments and numerical simulations are compared and illustrated. It has been found that resonant frequencies obtained using numerical and laboratory experiments are in good comparison and can be extended for various strain levels and sample conditions."

Jan 1, 2015

"A numerical procedure based on Smoothed Particle Hydrodynamics (SPH) is developed for rock blasting phenomena involving multiple intersecting joints. In this procedure, an elasto-plastic damage constitutive model has been applied in the framework of SPH to simulate yielding response of the surrounding rock medium. In this paper, a square shape rock medium with centrally located emulsion explosive in a borehole is numerically blasted with and without the presence of two diagonal rock joints. The model is simulated by assuming all four sides of the rock medium as free surfaces. The paper elaborates the propagation of stress wave, expansion dynamics of gaseous products, and development of fracture, fragmentation of the medium and dissipation of energy due to damage.INTRODUCTIONDamage prediction of rock mass and its mean fragmentation size are of important to design blast parameters for a given explosive. Current blasting practice solely depends on empirical relationships between fragmentation sizes to those blast design parameters while mathematical descriptions of the physical laws governing the process of rock blasting are limited. A constitutive relationship based on properties of explosives and rock mass is of special interest because of its capabilities in describing the key features of rock fracture and fragmentation, resulting in understanding of the physical process involved in rock blasting. It is well known that blast-induced damage and fragmentation of rock mass occur due to detonation induced stress wave and product gas driven fracture propagation. Nevertheless, the fundamental mechanism of the failure process and the behavior of rock characteristics under these two different loading are not well understood in many cases, which are closely related to growth of fractures and interaction of cracks under different loading conditions.Continuum damage mechanics is considered to be appropriate to describe the failure process under such high-amplitude stress wave. Based on this approach, several theoretical models have been developed to analyze the dynamic damage of brittle material such as rock and concrete with microstructure (Chen, 1999; Taylor et al., 1986; Wu et al., 2004; Yazdchi et al., 1996). These models are mostly based on the theory of damage mechanics or fracture mechanics with some statistical treatment to incorporate random distribution of microcracks. In the application of fracture mechanics theory, sometimes it is difficult to estimate the micromechanical parameters such as crack density and fracture toughness in different modes. In this study, damage in rock material is evaluated according to the earlier model developed by Grady and Kipp (1980), which is a fracture model coupled with a material description for stress wave propagation to predict quantitatively fracture and fragmentation under explosive loading conditions."

Jan 1, 2015

By B. Panton, J. Carvalho, E. T. Brown

"A major consideration in the design of high capacity tiedown anchors for dam, bridge and tower foundations is the tensile resistance of the rock mass to pullout, typically as a result of overturning moments or hydrostatic uplift. Rock mass pullout capacity for installed anchors is developed from the tensile strength and fracture propagation properties of intact rock, the orientation and physical properties of the discontinuities and anchor confinement at depth. The typically assumed, but generally conservative, design approach is to calculate anchor pullout capacity using the dead weight of a uniformly shaped inverted “cone” with an assumed initiation point and breakout angle. In some cases, some tensile or shear strength is considered across the area of the assumed pullout cone, however, these strength estimates are based primarily on based on limited historic testing data, empirical relationships for rock mass strength developed for other applications, and engineering judgement.As an alternative to the current foundation anchor design method, this paper presents a new approach to the rock cone pullout problem using Discrete Fracture Networks (DFN) combined with numerical simulations. The simulations presented in the paper investigate the influence of fractures in a synthetic rock mass on ultimate anchor strength. In addition, the global influence on fracture network statistics on rock mass pullout failure is also discussed and it is contended that some of the aleatory and epistemic uncertainty inherent in the design of rock anchors can be reduced by considering fracture statistics in anchor design."

Jan 1, 2015

By A. J. Basu

Rock mass properties testing for determining strength parameters is an integral part of geomechanics. These tests are performed either in a laboratory environment or in a field (in-situ). The authors share their experience in using numerical modelling in simulating the laboratory test procedures. This approach helps studying the rock breakage and failure mechanism more efficiently and provides a high degree of visualization. Moreover, imperfections and discontinuities can be easily incorporated in these models. This capability of providing arbitrary fractures and discontinuities facilitates in building a large knowledgebase of rock strength and failure mechanisms. This knowledgebase may then be used to generate heuristics for determining rock strengths for rock mass with fractures and discontinuities.

Jan 1, 1997

By D. P. Roberts

This investigation aimed to simulate the experiments performed by Nic Gay (Gay, N.C. 1976. Fracture growth around openings in large blocks of rock subjected to uniaxial and biaxial compression. International Journal of Rock Mechanics and Minining Sciences and Geomechical Abstracts, vol. 13. pp. 231-243) on fracture growth around openings in blocks of rock subjected to uniaxial and biaxial compression, using the Elfen discrete element code. The results of the physical experiments indicated a number of trends that were successfully replicated in the numerical simulations, including crushing ahead of the face, the formation of a sub-vertical tensile fracture at centre-span, the formation of groups of fractures originating from the crushed face zone, the formation of large corner-to-corner fractures which precede sample failure, and linking of slot corner and sample corner fractures leading to sample failure. In addition, other features were apparent that were not reported by Gay. It is believed that these took the form of conjugate shear fractures (as reported by Gay in a previous paper) and other features that may have been altered or inhibited by the presence of pre-existing structure in the physical samples. It was found that the samples, though nominally made up of the same rock type, must have had significantly differing properties. The presence of pre-existing structures may have contributed to this discrepancy. Strain measurements from the physical experiments were reasonably well matched, except where the gauge lay directly in the path of advancing fractures.

Jan 1, 2012

By D. P. Roberts

This investigation aimed to simulate the experiments performed by Nic Gay (Gay, N.C.1976. Fracture growth around openings in large blocks of rock subjected to uniaxial and biaxial compression. International journal of Rock Mechanics and Minining Sciences and Geomechical Abstracts, vol. 13. pp. 231-243) on fracture growth around openings in blocks of rock subjected to uniaxial and biaxial compression, using the Elfen discrete element code. The results of the physical experiments indicated a number of trends that were successfully replicated in the numerical simulations, including crushing ahead of the face, the formation of a sub-vertical tensile fracture at centre-span, the formation of groups of fractures originating from the crushed face zone, the formation of large corner-to-corner fractures which precede sample failure, and linking of slot corner and sample corner fractures leading to sample failure. In addition, other features were apparent that were not reported by Gay. It is believed that these took the form of conjugate shear fractures (as reported by Gay in a previous paper) and other features that may have been altered or inhibited by the presence of pre-existing structure in the physical samples. It was found that the samples, though nominally made up of the same rock type, must have had significantly differing properties. The presence of pre-existing structures may have contributed to this discrepancy. Strain measurements from the physical experiments were reasonably well matched, except where the gauge lay directly in the path of advancing fractures.

Jan 1, 2012

By M. Kilfoil

Many minerals processing operations result in crushed rock being mixed with a liquid. Sometimes, the rock particles have to be dried before the next process begins and there may be a considerable distance between the equipment carrying out these processes. Bateman Engineered Technologies (BET) developed, patented and commercialized a method of simultaneous pneumatic conveying and drying. A joint project to assist with optimizing the process is under-way. Part of it is a numerical simulation of the process. The simulation is done using a dedicated program written in MATLAB, (commercial software which is widely used in engineering). It accounts for heat and mass transfer as well as particle transport. The main part of the program consists of a loop containing a series of equations which are evaluated sequentially. During any particular iteration, the properties are changed one at a time. Since the properties are changing relatively slowly, this method gives good results. The outputs of the simulation are consistent with the operational experience of BET’s engineers. The output is very comprehensive and all numerical. Since the program is currently being used for research purposes, flexibility is very important. Thus the output is converted into a spreadsheet file. Variations in the output quantities are displayed graphically using MS Excel. Some simulation could have been done using computational fluid dynamics (CFD) software. The advantages and disadvantages (technical, financial and skills requirements) of the two approaches for this particular type of simulation are compared briefly. Keywords: Pneumatic conveying, Heat and mass transfer, Computer model

Jan 1, 2003

By Bernhard Hiebl

Global trends that are constantly challenging the iron and steel producers are related to worsening of the iron ore grades, higher substitution rates of main iron ores to cheaper raw materials due to economic pressure and increased environmental standards. Changing the raw materials in an operating plant can significantly influence the sintering or pelletizing process. Adjusting all process parameters to regain stable and high productive operation takes considerable time and causes production losses. Furthermore changes in raw materials lead to variations in the process gas flows and emission concentrations.Therefore SVAI developed two models for numerical simulation of sintering process and pelletizing process. Amongst others the models are based on Ergun equation and include chemical reactions of the related process. For already operating plants the simulations shall predict changes of the process parameters when raw materials are exchanged or when operational parameters are varied. Furthermore the simulations give advantageous information for the design of green field plants and allow the simulation of different plant configurations.The objective of this paper is to demonstrate what the latest developments are in these numerical simulations and what Siemens VAI is doing to overcome the remaining challenges in this field.

Aug 2, 2018

By Li. Jiarong, Zuojian Liang, Qingyan Xu, Baicheng Liu, Shizhong Liu, Hailong Yuan

Single crystal solidification of turbine blades is a key technology for the production of reliable and high efficiency gas turbines. The solidification process of single crystal investment castings causes a number of problems that have to be solved such as the increase of casting defects and complexity of production process. Computer modeling and simulation is helpful to predict the defects and improve the process design. In the paper, a mathematical model for three dimensional simulation of solidification and single crystal growth of investment castings was developed based on CA, mechanism of crystal growth and basic transfer equations such as heat, mass transfer equations. Many factors including constitutional undercooling, curvature undercooling and anisotropy, which had vital influences on microstructure evolution, were considered in the model. Temperature, crystal growth interface and structure defects of investment castings of single crystal superalloy were investigated at withdrawal rates of 4.5mm/min and 7mm/min. Simulated results showed that microstructure evolution and structure defects were predicted reasonably and calculated results were coincident with experimental ones. The study indicated that three dimensional solidification simulation technique is a powerful tool for understanding the fundamental of solidification and the formation of structure defects in single crystal investment castings and finally for optimizing process design.

Jan 1, 2004

By Miaoyong Zhu, Sen Luo

"A two-dimensional CAFD model was developed for simulating the solute diffusion-controlled dentritic growth in the low Peclet number regime. Mass conservation and local equilibrium at the solid/liquid interface were solved to determine the growth velocity of interface with consideration of constitutional and curvature undercooling. The modified decentered square growth proposed by Wang et al was adopted to capture the first eight neighboring cells. The comparison between the classic Lipton-Glicksman-Kurz (LGK) analytical model for steady state tip variables and the simulated values showed that the present model was capable of quantitatively determining the dendritic growth, and the simulation tests demonstrated that the effects of the mesh anisotropy in growth kinetics and crystallographic orientation were eliminated. The numerical simulations were performed to the dendritic growth ofFe-0.6%C binary alloy in undercooled melt, and the results showed that the present model could predict the dendrite morphology with artificial growth orientation in acceptable level.IntroductionDuring the solidification of metals and alloys, the dendrite grows from the substance or the undercooling melt and solute redistributes between the dendrite and the interdendritic liquid. The solidification microstructure is close related to the phase and solute evenness and has a close linkage with the properties of final production.In the past few decades, various kinds of deterministic and stochastic models have been developed to predict the evolution of microstructures in solidification, with the fast development of computer power. Among these methods, Phase Field (PF) model [1, 2] has been known as one of the most adequate deterministic models for directly simulating the dendritic growth. However, there is a vital drawback that the PF method requires considerable computation time and is limited to simulate very small domains with a few dendrites [3]. On the other hand, stochastic models, such as Monte Carlo (MC) method [4, 5] and Cellular Automaton (CA) method [6-8], have been usually applied to predict the solidification microstructure evolution. Comparing with the MC method, the dendritic growth kinetics is explicitly integrated in the CA model, while the MC model doesn't. Therefore, the CA method is more reasonable to predict the dendritic growth during the solidification of melt."

Jan 1, 2013

By Min Chen

A numerical simulation of the stress field of the solidified shell in the mold during continuous casting of peritectic steel slab was calculated with commercial software. The results showed that the maximal stress on wide face of the solidified shell was located at the corner and near middle of the slab, and it increased with increasing the width of slab. In addition, the stress was also influenced by operating parameters such as superheat degree, drawing speed and cooling intensity of the mold, and the maximal stress value and its position changed with these parameters. For the slab with section size of 3200mm×150mm, the proper drawing speed was around 1.2m/min, with cooling intensity of 5500L/min while the superheat degree was 15°C to 25°C in order to prevent longitudinal crack occurring during peritectic steel continuous casting.

Jan 1, 2009

By Huimin Lu, Lian Zhou, Qiang Li

"A three-dimensional mathematical model is developed for the regenerative rotary hearth furnace with 40 meters medium diameter. The influence of combustion switchover and outlet pressure on the fluctuation and distribution of temperature and pressure is studied with numerical simulation method. The results showed that the temperature field and flow field is suitable for requirement of industrial production; the fluctuation of the temperature field and pressure is acceptable; the influence of outlet pressure on the furnace pressure is great.IntroduceRotary hearth furnace which is originally utilized as heat treating and reheating industrial furnace, now is also used for carbon based DRl production. With the need of market and the concern of resource recycling and environmental issues, the DRl production from ore or from waste material by rotary hearth furnace attracts more and more attentions from steel industry [1,2].When rotary hearth furnace works, the feed pellets are charged on the hearth and rotate within an annular furnace chamber with the process of reduction. For the sake of rapid reduction, usually the depth of the pellet is one or two layers, accounting for less than 90% of the volume of furnace. As result, the exhaust gases carry away a lot of heat, so that energy utilization efficient is low [3]. Currently, in order to reduce the consumption of energy and improve the production efficiency, many methods have been adopted. The regenerative combustion technology is one of the most promising."

Jan 1, 2013

By Dong-dong Liu, Jian-ping Niu, Qi Liu, Shu-sen Cheng

A numerical model was developed to simulate the three-dimensional heat transfer in the belly and lower shaft region of the blast furnace wall consisting of copper steel composite stave, packing and shell. The temperature distribution of copper steel composite stave without the slag layer or refractory lining was investigated and then compared with the copper stave. The results showed that copper steel composite stave can satisfy the requirement of overheatingfree and had an analogous heat transfer behavior as copper stave but more stable cooling intensity. This model had been validated by the good agreement between the computed and measured temperature in the thermal test.

Jan 1, 2015

By Dengfu Chen, Huabiao Chen, Yunwei Huang, Mujun Long, Youguang Ma, Lintao Gui

In the current work, a three-dimensional model coupling turbulent flow, heat and solute transport was developed to investigate the solute transport and redistribution in the mold of slab continuous casting where the Reynolds number is very high and the molten steel flows strongly. The conservation equations of momentum, energy and species for a multicomponent system which includes solute element C, Si, Mn, P, S were solved with the commercial software ANSYS Fluent. The parameters used in the model are based on the actual continuous casting. The fluid flow, temperature distribution and solute elements distribution were analyzed. Results showed that negative segregation occurred near the strand surface, and the species concentration reached a peak value at the solidification front during solidification. The segregation degree of solute element S is higher than the other solute elements because of its lower equilibrium partition coefficient. The species concentration in the liquid pool is homogeneous because of the turbulent flow effect.

Jan 1, 2015

By S Servant

In the mine of this study, cemented hydraulic backfill is placed under the form of a slurry at a percentage of solid by weight varying from 50 to 65 per cent. This slurry settles and cures in the stope. During this process an important amount of water must be drained from the stope. Field observations in the case studied here show that the geometry of the drains system has a significant impact on the drainage of the stopes. The mode of filling produces a backfill showing an evident stratification typical of a sedimentation process. In situ observations on this material show that this stratification translates into a very strong heterogeneity in the physical properties of the cemented backfill. This heterogeneity should affect the hydraulic conductivity of the backfill in the stope. In order to study this aspect, finite element simulations were performed. This paper presents the results of these simulations, where the impact of this stratification on hydraulic conductivity is related to the performance of the backfill drainage system for various geometry. This preliminary study is aimed at providing a better understanding of the role of the drains in the dewatering process of the backfilled stope in the case of this mine.

Jan 1, 1998

By Servant S

In the mine of this study, cemented hydraulic backfill is placed under the form of a slurry at a percentage of solid by weight varying from 50 to 65 per cent. This slurry settles and cures in the stope. During this process an important amount of water must be drained from the scope. Field observations in the case studied here show that the geometry of the drains system has a significant impact on the drainage of the stopes. The mode of filling produces a backfill showing an evident stratification typical of a sedimentation process. In situ observations on this material show that this stratification translates into a very strong heterogeneity in the physical properties of the cemented backfill. This heterogeneity should affect the hydraulic conductivity of the backfill in the scope. In order to study this aspect, finite element simulations were performed. This paper presents the results of these simulations, where the impact of this stratification on hydraulic conductivity is related to the performance of the backfill drainage system for various geometry. This preliminary study is aimed at providing a better understanding of the role of the drains in the dewatering process of the backfilled stope in the case of this Mine.

Jan 1, 1998

By Kumar Ranjan, Kumari A. V. Asha, G. E. Sreedhar, Banerjee R. Raja, M. Narasimha, R. Shivakumar

In this paper, numerical simulation of magnetite, medium segregation and coal partitioning have been made in a 350 mm diameter cyclone. An attempt is made to assess the nature of impairment in efficiency in relation to the increasing amount of near gravity materials in coal. The CFD model approaches the Algebraic Mixture (ASM) model with the granular option and Reynolds stress model (RSM) to resolve the air core and turbulent mixing of the medium and coal particles. The partition characteristics of the DMC were modelled using both Lagrangian & Eulerian particle tracking for particles ranging in size from 0.5 to 8 mm and in density from 1200 to 2000 kg.m-3. Lagrangian frame work based Discrete Phase model used to simulate the coal particle trajectories inside the DMC and attempted to understand the residence time of different size and density coal particles. It was observed that coal particles having density near to separation-density (having fifty percent of probability to under flow and overflow) are exhibiting increased residence time than any density particles. As expected the smaller size coal particles shows higher residence time than the coarse coal particles. CFD simulations on effect of NGM mass fraction are simulated and analyzed using ASM model including for coal and magnetite. Coal particles with high NGM content show significant effect in misplacement of coal particles towards wrong products at all relative densities. The residence time of the particles increases because of its increased interaction with near dense particles. The effect of NGM fraction on overall cyclone performance as well as detailed flow assessment in terms of residence time distribution, NGM particle concentration profiles with respect to LZVV is presented.

Jan 1, 2018