Search Documents

By J. H. Curran

In this paper, a numerical model for determining the stress and displacement fields around a crack of arbitrary geometry in a linear elastic medium is described. Accuracy and computational speed were the main considerations in selecting the numerical method and Propagation criteria used to model the growing crack. Some results regarding the stresses and displacements around a two-dimensional crack of arbitrary geometry are discussed.

Jan 1, 1984

Following a short review of selected literature, a back analysis of stress measurements of stress regeneration in a closed mini-longwall is shown. The agreement is good, lending confidence to the parameters used in the proposed numerical model. Foundation yielding was modelled in the absence of total stope closure. A ?Prandtl? wedge was formed, coinciding with the yield point of the pillar/foundation system. Some theoretical work is used as a starting point for the formulation of a new design methodology for strike or dip stabilizing pillar/foundation systems in deep level goldmines. The existing guideline for the stability of foundations (APS< 2.5UCS) is critically examined. Finally, the advantages of the new approach, and the work required to substantiate the proposed new design methodology, are outlined.

Jan 1, 1998

By D. F. Malan

Considerable practical knowledge and effective use of empirical evidence is required in the design of tunnel support systems. In many instances, it is unclear how individual support elements interact with fractured rock and which elements of the support system are critical in controlling and containing the movement of material proximate to the excavation. The present paper outlines a displacement discontinuity approach that allows the coupling of surface constraints, in the form of liner or mesh material, to an explicit representation of time-dependent rock discontinuity creep movements. The method is illustrated for a simple square-shaped tunnel profile in which the fracture zone is mobilized by an imposed stress field. A crucial aspect of the model is the computational treatment of the coupling between the liner material and there presentation of the fractured rock mass. This allows the response of the liner material to be assessed quantitatively in terms of the possible reduction of movement in the rock mass. The importance of installing the liner timeously after the excavation is formed is illustrated in a specific example. A number of computational difficulties in simulating liner-rock interaction mechanisms are identified.

Jan 1, 2008

By J W. Lee, E Baafi, N Aziz

Ground subsidence due to mining is a common dilemma confronting the underground coal mining industry. The effects of underground longwall mining on river valleys have come under scrutiny, especially when mining is located underneath catchment areas. As public awareness of mining impacts on the environment increases, there is the need to develop damage mitigation strategies. Results of numerical modelling show that river valley response to underground longwall mining can be simulated and evaluated.

Jan 1, 2005

By B BussiFre, J Molson, M Aubertin, G W. Wils

Waste rock piles containing reactive sulfide minerals can serve as long-term sources of acid mine drainage (AMD). Within these piles, many physical, geochemical and biological processes can contribute to the production of AMD. The flow of water, for example, depends on internal pile structure, which depends on the dumping methods used to construct the pile. In this paper, the internal structure and grain size distribution of typical waste rock piles are first reviewed with respect to their effect on internal water flow. A numerical model is then applied to simulate unsaturated flow within a waste pile assuming two primary waste materials, different atmospheric conditions and several geometric configurations. The simulations suggest that water movement can be controlled by introducing fine-grained material layers within the coarser host material. However, these layers should not be horizontal because they can lead to localised high-saturation regions as well as preferential flow favouring the production of AMD. A shallow slope towards the pile boundary may be introduced to force drainage to the outside while maintaining low saturation within the interior, thus reducing the production of AMD.

Jan 1, 2003

By Horst Peter Wölfel

The goal of biodynamic models is to simulate the vibration behaviour of the human body. In combination with experimental studies biodynamical models can be a powerful tool for the analysis of the effects of vibration exposure on health [1] and comfort. This paper gives an overview of the state of the art of biodynamic whole-body vibration models of humans, addressing both numerical models and hardware dummies. Method Two approaches are distinguished, the phenomenological and the anatomical, as illustrated in Figure 1. [ ] Phenomenological models aim to reproduce the vibration behaviour of humans with respect to particular physical quantities, chiefly the driving-point impedance at the interface to the seat, and partly with respect to other transfer functions. Discrete systems of masses, springs, and dampers with several degrees of freedom whose topology and parameters are determined by structure- and parameter identification methods are used in the sense that the functions derived from measurements are reproduced as well as possible. This paper provides an evaluation of this methodology and defines its range of application as well as its limits.

Jan 6, 2006

By Henry Hu, Alfred Yu

"A numerical model has been developed to simulate heat and mass transfer phenomena occurring during a novel magnesium reduction process-the vertical retort technology. The model was based on the control - volume finite difference approach. Simulations were run to determine the effect of various processing parameters on magnesium reduction cycle time. The model predicted temperature distributions, heating curves, recovery ratio of magnesium, and total process times. The predictions were used to optimize the magnesium reduction process including dimensions of the retort, geometries of the reactant compounds, and productivities. The results indicate that the implementation of the optimized processing parameters significantly decreases the reduction time, and consequently increases the production rate as well as lowers the energy consumption.IntroductionMagnesium consumption has been growing considerably in the past decade [1]. Of the total magnesium production, 30% is produced by the thermal reduction (Pidgeon) process, and it is anticipated to be still the main technology for producing magnesium in future [2]. One of the advantages of the Pidgeon process is that it is capable of producing high purity magnesium demanded by both the automotive and aerospace industries. However, the inefficiency of heat and mass transfer taking place in the horizontal Pidgeon reduction process results in low productivity, high energy consumption, labor-intensive operation, large floor space and high capital investment. This limits the production scale of a horizontal Pidgeon reduction process [3].Vertical reduction processTo overcome the shortcomings of the horizontal magnesium reduction process, a novel vertical magnesium reduction process [3] for large-scale production has been developed, in which heat and mass transfer in the retort is enhanced. The new process shortens reduction time, reduces energy consumption and increases productivity. The vertical reduction system is comprised of a vertically-placed retort, a heating furnace and a vacuum system. A section of the vertical retort is illustrated in Figure 1. Since the retort is set vertically, instead of horizontally, the entire procedure of loading the reactants material in the new process is different from the one for the horizontal reduction process."

Jan 1, 2001

Floor heave events in underground mining often happen suddenly and cause risk to workers and equipment. Reducing the risk of floor heave must be addressed in mining layout design. To improve that design one must understand the significant parameters that are correlated to the occurrence of floor heave. To accomplish that goal, researchers from the National Institute for Occupational Safety and Health (NIOSH) completed a numerical study using FLAC3D, taking advantage of its strain softening and ubiquitous joint capabilities to better simulate the floor heave and caving mechanisms involved at an actual mine site.

Jan 7, 2020

By S. P. Goplen, M. Henriksen, T. F. Barron

In the past decade, numerous large scale drilling projects have made use of reverse circulation air lift methods. Reverse air lift pumping can be established by injecting gas into a column of liquid; the subsequently decreased density of the gas-liquid mixture allows atmospheric pressure outside the lift conduit to establish flow. The rate of flow can be controlled by varying the gas to liquid ratio. A typical reverse circulation air lift system is depicted in [Fig. 1]. Allen 1 has described a number of reverse air lift drilling projects. Although cited examples are generally related to large diameter holes, the method has seen use in other applications such as sand dredging, chemical pumping, water well pumping. and continuous coring. One particular benefit obtained with enclosed reverse circulation drilling is elimination of lost circulation. The required air/water ratio for pumping at a given submergence and lift has been established experimentally by Gibbs.2 His data, however, do not reflect sensitivity of flow rate to air/water ratio. Further, in Gibbs&apos; work, air apparently is injec-

Jan 1, 1980

By R. I. Jeffery, J. W. Summers

An account is given of the verification and calibration of a numerical modelling technique that has been developed for the prediction of subsurface strata deformation following the extraction of coal by longwall mining methods.

Jan 4, 1992

By S. P. Girrens

A coupled finite element geomechanical- hydrology code is currently under development for application to the problem of predicting groundwater disturbances associated with mine subsidence. The structural-fluid coupling is addressed by calculating the subsided mine geometry, with emphasis placed on determining the strata disturbance and locating damaged regions, for input into a hydrology code, which determines localized volume flow rates and aquifer fluctuations. Benefits from coupling will be best realized when field measurements, an additional aspect of the study concurrent with analytical investigations, indicating the relationship between increasing rock strain and increasing permeability are incorporated into hydraulic material descriptions. Hydrologic and struc- tural calculations are presented to demonstrate computational capabilities applicable to mine subsidence

Jan 1, 1982

By G. D. Schena, M Zanin

A numerical approach is presented for predicting the liberation of phases during comminution, with direct calculation of a liberation spectrum giving the fraction of the particulate system in each narrow size-composition class. The numerical model uses a Poisson flat process for stochastic generation of the non-convex three-dimensional Poisson polyhedral textures and the system of fractures that give rise to convex fragments. With this numerical approach, there is no need to assume random uniform isotopic fracture (lack of selective liberation); only computational geometry primitives are exploited for the derivation. The agreement with relationships obtained previously under the random uniform isotopic fracture hypothesis demonstrates that the approach is capable of deriving full liberation spectra directly. Both volumetric and superficial composition can be calculated

Jun 21, 1905

By W. M. Ma, Daniel W. H. Su, K. Centofanti, Yi Luo, W. L. Zhong, Syd S. Peng

NTRODUCTION Surface subsidence monitoring is the most effective means for investigating the nature of, and developing the techniques for pre¬dicting the surface movement and deformation process. However, because of imperfection in survey equipment and field practices, the surface movement data collected from field survey are not error-free. Also due to lack of proper numerical techniques for computing surface deformation from the measured surface move¬ment data, a relatively large portion of the subsidence data col¬lected could not provide information other than surface movement. In a paper by Luo and Peng (1992), the sources and magnitude of errors in the measured surface movements using an electronic total station in subsidence survey are discussed. If the errors in the measured surface movement components are left untreated before calculation for deformation indices (i.e., slope, strain and curva¬ture), the errors will be carried over to the deformation indices and enlarged if improper numerical procedures are used in the calcu¬lation. It should be noticed that in protecting structures affected by ground subsidence, accurate surface deformation indices are much more important than surface movement. In order to improve the accuracy of the collected surface movement data and the calcu¬lated deformation indices, a number of numerical procedures are also proposed in the same paper. The proposed techniques include: a data smoothing technique for data pre-processing; and a number of numerical methods and/or procedures for calculating surface deformation indices from the measured surface movement data. DATA SMOOTHING Because of various factors, the measured surface movement data always contain errors, small or large. In order to better utilize the collected data, a data smoothing method is proposed to correct the errors in the measured movement data. The method uses the same concept as the influence function method (e.g., Knothe, 1957) on which many existing subsidence prediction methods are based. The data smoothing method makes necessary corrections on the measured movement for the point of interest according to not only the measured movement at this point but also those at the surrounding points. The contribution of a surrounding point in making the correction for the point of interest is inversely propor¬tional, in exponential form, to the distance between them. The new magnitude of a movement component at the point of interest (Zi*) is calculated by

Jan 1, 1992

By Hongbo Yang, Tongbo Zhang, Ting Huang, Jingshe Li, Fangfang Song

"Liquid core length is a very important parameter of the process of continuous casting, which directly reflects the cooling status in secondary cooling zone and solidification rate of the slab. In this paper, the billet caster in some steel plant in China was chosen as the prototype of simulating. A numerical simulating model for calculating the liquid core length of continuous casting billet was built by ProCAST software, which was verified reasonable and accurate by pin-shooting technique. Based on the model, the relationships between casting speed, casting temperature and liquid core length of casting billet were found, which was very practically useful for the workers on sit to improve casting billet quality. IntroductionContinuous casting process developed very rapidly because of the high metal yield, energy saving, and low-cost 111• And the quality requirements of continuous casting billet for customers became higher and higher with the fast development of economic121• Internal defect usually is the major quality problem for continuous casting billet[3•61 , which is influenced by the solidification process of casting billet, especially the solidification process of liquid core17•81 • So it is very important to study on the relationship between casting parameters and the liquid core length which means the distance between the crystallizer meniscus and the casting billet solidification end-point. In practical production process, the workers usually control the quality of casting billet according to the liquid core length."

Jan 1, 2013

By S. B. Johnson

Keywords: Spray, Aluminum, Simulation A stochastic, droplet-resolved model has been developed to describe the behavior of a binary aluminum alloy spray during the spray-rolling process. In this process, a molten aluminum alloy is atomized and the resulting spray is deposited on the rolls of a twin-roll caster to produce aluminum strip. The one-way coupled spray model allows the prediction of spray characteristics such as enthalpy and solid fraction, and their distribution between the nozzle and the deposition surface. This paper outlines the model development and compares the predicted spray dynamics to PDI measurements performed in a controlled configuration. Predicted and measured droplet velocity and size distributions are presented for two points along the spray centerline along with predicted spray averaged specific enthalpy and solid fraction curves.

Jan 1, 2005

By Huang Fenglei, Jin Qiankun

A new linear charge method has been applied in coyote blasting for highway excavation in which presplitting boreholes are set along the designed rock slope with some empty chambers at the end of linear charges whose axes are perpendicular to the axis of the highway. Three-dimensional nonlinear finite element analysis for this method has been conducted using LS-DYNA3D. The results indicate that the blasting damage to the rock slope is reduced greatly using the new linear charge method in coyote blasting.

Jan 1, 1998

By H. Huang, J. C. Wu, C. H. Xie, Y. Q. Xue, Z. H. Zhang

Two- and three-dimensional groundwater flow models for dewatering in an opencast mine are presented. The models focus on solving special problems using numerical simulations of the opencast mine. These special problems include the drying up of the aquifer, the movement of the model boundary and the prediction of the residual head. These models are used to describe the decline of groumdwater levels in the Yuanbaoshan Opencast Coal Mine and its neighboring area during the drainage period. The simulated values agree very well with the field data. There have been several successful examples of numerical simulation for predicting the mine dewatering rate in an opencast mine. However, The hydrogeological conditions for the Yuanbaoshan Opencast Coal Mine are rather complex, and the total average flow rate is very high at 4.676 x 105 m3/d (1.65 x 107 cu ft/day). In this paper, a two-dimensional groundwater flow model in a large area is developed to solve the problems mentioned above by simulating regional drawdowns of the groundwater level during the drainage period, and a three-dimensional flow model is developed to simulate the residual head near the open-pit bench.

Jan 1, 1998

By F. Liu, H. A. Meier, D. M. Valiveti, P. A. Gordon

Hydraulic fracture is a frequently used technique in the Oil and Gas industry to create high permeability pathways for extraction of hydrocarbons. Numerical simulation of this method poses challenges such as representation of 3D fractures, propagation of randomly oriented fractures in 3D media, enabling fluid leak-off from fracture to poro-elastic medium. We have developed a numerical model using the Generalized Finite Element Method (GFEM) to simulate fluid driven fracture in a poro-elastic medium that addresses these challenges. In the GFEM framework, special enrichment functions are used to augment the solution space in the displacement and pore-pressure fields in order to accurately capture the response in the vicinity of the fracture. Adaptive refinement of the finite element mesh is also employed to facilitate accurate and efficient solution of the model. Our model fully couples fluid flow in the fracture, represented with the lubrication approximation, to the geo-mechanical behavior of the bulk poroelastic rock, comprised of an elastic rock skeleton saturated with fluid. Fracturing fluid exerts traction on the fracture surface and is also allowed to leak-off into the poroelastic medium using a transfer function approach. Stress intensity factors (SIFs) are evaluated based on the Contour Integral Method, accounting for the fluid pressure on the fracture surface. The SIFs are used to govern a full 3-D evolution of the fracture. Verification of the model is performed for known behaviors for the poroelastic response and fluid flow in the fracture.

Jan 1, 2015

By Kathy S. Kalenchuk

Damages have been observed in shafts, backfill infrastructure and the slick line at the Leeville Mine. Shaft damage (shearing and spalling of the concrete liner) is attributed to mine-induced rock mass deformations concentrated where regional faults intersect the shafts. Aggregate bin, mixer station and conveyor drive damage (shotcrete shearing, heave of concrete floor slabs and bending of a vertical steel beams) is attributed to rock mass squeezing induced by elevated magnitudes of the vertically orientated major principal stresses. Separation of the slick line has occurred do to extension of the stressshadowed rock mass immediately above the mined out areas. Three-dimensional numerical models have been developed and calibrated to simulated rock mass deformations with magnitudes similar to those believed to have produced the observed damages. By stepping through the mine sequence, models reproduce deformations at the same mining stage as when damage has occurred. Calibrated models have been used to forward project further damage over the remaining mine life.

Jun 1, 2012

By Y Wada, Y Ogata, T Nishida

"Shaped charge is used for blast demolition of steel construction in general. It is important tocontrol movement of cutting part with shaped charge for the blast demolition of steel construction.Movement of cutting part was studied in the small scale experiment with steel bar. A numericalsimulation using the Discontinuous Deformation Analysis (DDA) was carried out. The DDAmethod was applied to the full scale blast demolition of steel construction. Results of thesenumerical simulation was compared with the blasting demolition work of steel tower."

Jan 1, 1999