Search Documents

By R. J. Pine, J. M. Rance, J. S. Coggan

ABSTRACT: The safe design of excavations created in a fractured rock mass is of paramount importance in rock engineering. Depending on the nature of the stress regime acting on an excavation, and the associated deformational behaviour of the rock mass surrounding the excavation, a number of potential failure mechanisms exist. These involve a combination of shear failure on existing joints, extension of critically oriented joints and propagation of new fractures through previously intact rock. Conventional numerical codes usually adopt either a continuum or discontinuum approach; this often results in significant over simplification and kinematic limitations particularly when modelling complex rock mass behaviour in the close vicinity of excavations. This paper describes the results from several case examples using a combined finite element-discrete element code, ELFEN that incorporates a crack propagation mode. Through the use of selected rock engineering examples, the authors demonstrate the potential of the code in the realistic simulation of complex rock failure both in underground and surface excavation applications. RÉSUMÉ: La conception sûre des excavations créées dans une masse de roche fracturée est d'importance primordiale dans la géotechnique. Selon la nature du régime d'effort autour d’une excavation, et le déformation associé de la masse de roche autour de l'excavation, il exsite un certain nombre de mécanismes potentiels d'échec. Ceux-ci comportent une combinaison d'échec de cisaillement sur les joints existants, la prolongation des joints d’orientation critique et la propagation de nouvelles ruptures dans la roche précédemment intacte. Les codes numériques conventionnels adoptent habituellement l’approche, soit d’un continuum, soit de discontinuum; ceci a souvent comme conséquence une sursimplification significative et les limitations cinématiques en particulier en modelant le comportement complexe de la masse de roche très proche des excavations. Cet article décrit les résultats de plusieurs exemples de cas en utilisant un code qui s’agit d’une combination d’élément discret et d'élément fini, ELFEN qui incorpore un mode de propagation de fissure. Par l'utilisation des exemples choisis géotechniques, les auteurs démontrent le potentiel du code dans la simulation réaliste de l'échec complexe de roche pour les applications aux excavations soit au fond soit à l’extérieure. ZUSAMMENFASSUNG: Das sichere Design der Aushöhlungen, die in einer zerbrochenen Felsenmasse verursacht werden, ist vom paramount Wert in der Felsentechnik. Abhängig von der Natur des Druckregimes, das auf einer Aushöhlung und dem dazugehörigen deformationalen Verhalten der Felsenmasse umgibt die Aushöhlung fungiert, bestehen eine Anzahl von möglichen Ausfalleinheiten. Diese beziehen eine Kombination des Scherausfalls auf vorhandenen Verbindungen, Verlängerung der kritisch orientierten Verbindungen und Ausbreitung der neuen Brüche durch vorher intakten Felsen. Herkömmliche numerische Codes nehmen normalerweise entweder ein kontinuum oder discontinuum Annäherung an; dieses ergibt häufig bedeutende Übervereinfachung und kinematische Beschränkungen, besonders wenn es kompliziertes Felsenmasseverhalten in der Nähe der Aushöhlungen modelliert. Dieses Papier beschreibt die Resultate von einigen Fallbeispielen mit einem kombinierten begrenzten Element - getrennten Elementcode, ELFEN, das einen Sprungausbreitungmodus enthält. Durch den Gebrauch von vorgewählten Felsentechnikbeispielen, demonstrieren die Autoren das Potential des Codes in der realistischen Simulation des komplizierten Felsenausfalls bei Anwendungen für beide Tagebaue und Untergrundsaushöhlungen.

Jan 1, 2003

One of the most important goals of the experiments is to determine deformation and strength of the rock specimen under applying triaxial load, the status of failure condition is one of the subjects which can be used in soil mechanic and foundation engineering applications. To investigate the effect of confining pressure on the triaxial test, the rock has been considered as bonded-particles method and has been simulated by particle flow code three dimensions (PFC3D). This model, does not have pre-existing microcracks and all the particles have the same size (2 mm diameter), equal to the real specimen shape. To compare the mechanical behavior and failure pattern of the numerical and experimental results at uniaxial and triaxial tests at the same conditions, the Young's modulus, Possion's ratio as well as maximum axial stress have been considered. On the other hand, the microcracks growth and change of failure pattern at the modelling of the uniaxial and triaxial tests with different confining pressures up to failure point were reported. According to the experimental work, the number density of microcracks decreases from starting value and remains almost constant up to the failure beyond that point. The comparison of the numerical and experimental results of maximum axial stress and the Possion's ratio revealed good accordance. The initial number density of microcracks at the model was zero and increased while microcracking. Also, the sudden increment around the maximum stress was observed which is because of unstable growing of microcracks near the maximum stress. In all uniaxial and triaxial tests with different confining stresses, the number density of microcracks during the failure remains almost constant which can be considered as a proof of failure occurring in the model. The experimental results indicate the similar trend as well.

Nov 25, 2010

By Kubik K

The objective of this paper was to determine if recent advances in numerical electromagnetic wave propagation modelling techniques could be used to study the propagation of GPS signals over arbitrary terrain. Specifically the effects of diffraction and multipath were to be modelled accurately for propagation over arbitrary terrain surfaces defined by a digital terrain model. Ideally the model was to be in three dimensions and cover an area of a few square kilometres with a simulation height of a few hundred metres. Simulations were to run on current technology PC's with run times less than a few hours. A relatively new technique that involves a numerical solution to the Parabolic Equations can be used to solve for two-dimensional propagation over any type of terrain. The PE provides a direct solution of Maxwell's wave equations by approximating the Helmholtz scalar wave equation. This technique does not rely on the study of individual ray paths as used in the Geometrical Theory of Diffraction (GTD). It was decided to develop a two dimensional Parabolic Equation code, using the Fourier Split/Step method that allowed propagation over arbitrary terrain (; ; ). The model boundary conditions assumed that the satellite transmission was a plane wave, and that the lower boundary was a perfect reflector. Propagation simulations from the model accurately provide the amplitude and phase of the propagated plane wave at all points within the model domain. The PE model was used to study the effects of diffraction and multipath caused by various types of terrain commonly found in an open cut mining environment. The study of diffraction errors was relatively simple, because we are simply concerned with amplitude of the signal as it diffracts about an obstacle. If the amplitude is above the level receivable by a GPS receiver then the diffraction error can be obtained from the geometry of the terrain obstacle. The study of multipath is more complex, as the computed phase from the PE model is ambiguous and the presence of multipath can only be inferred. Thus the absolute magnitude of the nmltipath error is not currently determinable, though the results of some new work is shown here. This paper shows that diffraction effects can cause significant errors in range. It also shows the effects that multipath has on the received amplitude of the signal and the received phase of the signal. A discussion is made of further work that is presently underway that will allow the direct determination of the multipath error by comparing the ambiguous phase with a reference phase.

Jan 1, 1997

This paper discusses support requirements for the proposed monorail system to be used in decline development. The monorail drilling and loading systems are systems that move on the rail (monorail) installed in the roof of the decline and supported by roof bolts, suspension chains and steel supports. However, due to the weight of the components of the two systems, it is imperative that the force in each roof bolt, suspension chain and steel support capable of suspending the weight of the heaviest component is determined. Numerical models that relate the weight of the monorail drilling and loading components to the required strength in the support system have been developed. Using these developed models, numerical values of the forces in each roof bolt, suspension chain and steel support, required to suspend the weight of the heaviest component of the monorail drilling and loading systems are determined.

Dec 6, 2010

By Andrzej Walentek

This paper presents the results of numerical modeling of the fractured zone around longwall gateroads. The numerical computations have been based on the results of underground measurements of the fractured zone around longwall gateroads, which were performed for various geological and mining conditions and measured with the use of the endoscopic method with a digital camera. Having results of the underground investigations, trials of duplicating them by the Finite Element Method using software PHASE and the Hoek-Brown failure criterion were performed. In order that the fractured zone predicted by the model-based computations closely simulate that measured underground, some selected parameters in the Hoek-Brown failure criterion were modified based on the earlier determined relations used for prediction of the roof-to-floor convergence in longwall gateroads. As a result, the model-based computations produce results very close to those obtained from the underground measurements. Determination of the fractured-rock zone around the mine roads is of great importance to mining practice, especially in planning and designing supports for longwall gateroads. The influence of the front abutment pressure causes the fractured rock zone around the roads to increase and results in an increase of support load. Moreover, proper determination of the fractured rock zone is one of basic elements in the process of planning roof bolting for mine workings.

Jan 1, 2009

By Andrzej Walentek

This paper presents the results of numerical modeling of the fractured zone around longwall gateroads. The numerical computations have been based on the results of underground measurements of the fractured zone around longwall gateroads, which were performed for various geological and mining conditions and measured with the use of the endoscopic method with a digital camera. Having results of the underground investigations, trials of duplicating them by the Finite Element Method using software PHASE and the Hoek-Brown failure criterion were performed. In order that the fractured zone predicted by the model-based computations closely simulate that measured underground, some selected parameters in the Hoek-Brown failure criterion were modified based on the earlier determined relations used for prediction of the roof-to-floor convergence in longwall gateroads. As a result, the model-based computations produce results very close to those obtained from the underground measurements. Determination of the fractured-rock zone around the mine roads is of great importance to mining practice, especially in planning and designing supports for longwall gateroads. The influence of the front abutment pressure causes the fractured rock zone around the roads to increase and results in an increase of support load. Moreover, proper determination of the fractured rock zone is one of basic elements in the process of planning roof bolting for mine workings.

Jan 1, 2009

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

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

By A. T. Iannacchione

Numerical simulation of coal pi1lar loading has traditionally been a difficult task due to the unique and highly variable properties of coal and the inability of numerical procedures to duplicate these properties. This simulation was based upon material properties developed from laboratory tests of coal cubes and examined against coal pillar strength data developed from actual underground measurements. The coal cube measurements indicated the coal followed a nonlinear Mohr- Coulomb failure criterion. The in situ measurements showed that a pillar yield zone developed similar to Wilson&apos;s hypothesis and, at high confinements, the pillar core seemed to followed a pseudo-ductile behavior. The strain-softening material model, based on a 1inear Mohr- Coulomb failure criterion, simulated coal cube and pillar behavior at moderate load conditions. However, at high loads the models were inaccurate, due to an inadequate failure mechanism and/or to changes in the pillar constraints at the coalbed interface.

Jan 1, 1989

By Nan Wang

Based on the commercial computational fluid dynamics software, Fluent, the fluid flow and temperature field in an iron-bath reactor were simulated numerically, and the effects of top-blown and side-blown operating parameters on the mixing time of iron bath, distribution range of oxygen in top space, and post-combustion ratio were investigated. The simulated results showed that the mixing time of iron bath increased initially and then decreased with the increase of inclined degree and immersion depth of side-blown lance. The distribution range of oxygen became wider as pulling the top lance away from the iron bath. Furthermore, post-combustion ratio was reduced given the top lance was too close the bath. The simulation results provide theoretical base for the process optimization of iron-bath reactor and further development of new type iron-bath smelting reduction process.

Jan 1, 2009

By K. S. Rao, R. P. Tiwari, V. K. Shrivastava

ABSRACT: In this paper numerical simulation of jointed block mass behaviour obtained in an experimental study has been presented. Singh et al. (2002) conducted extensive tests on the jointed model materials with several joint configurations. The results of this experimental study were used to validate the strength behaviour of numerical block mass model developed using professionally available software. This software basically contains different kinds of material models. From the results of numerical simulation the existence of strength anisotropy was observed. The strength anisotropy so obtained has been compared with the experimental values reported by Singh et al. (2002). The anisotropic curves in both cases are U-shaped with flat base. They show good agreement in sliding mode of failure zone in which inclination ß (= 90° - ?) of joint set-I lie between 20° to 40° for all cases of interlocking. ZUSAMMENFASSUNG: In dieser numerischen Papiersimulation des verbundenen Blockmasse Verhaltens, das in einer experimentellen Studie erhalten wird, ist dargestellt worden. Singh et al. (2002) leiteten umfangreiche Tests auf den verbundenen vorbildlichen Materialien mit einigen gemeinsamen Konfigurationen. Die Resultate dieser experimentellen Studie wurden verwendet, um das Stärke Verhalten des numerischen Blockmasse Modells zu validieren, das mit professionell vorhandener Software entwickelt wurde. Diese Software enthält im Allgemeinen unterschiedliche Arten der materiellen Modelle. Von den Resultaten der numerischen Simulation wurde das Bestehen von Stärke Anisotrophie beobachtet. Die Stärke Anisotrophie also erreicht ist mit den experimentellen Werten verglichen worden, die von Singh et al. (2002) berichtet werden. Die anisotropen Kurven in beiden Fällen sind mit flacher Unterseite u-förmig. Sie zeigen gute Vereinbarung in gleitender Erscheinungsform des Ausfalls Zone, in der Neigung ß (= 90° - ?) von der Verbindung stellenSie Lüge zwischen 20° auf 40° für alle Fälle vom Ineinandergreifen ein. RESUME: Dans cette simulation numérique de papier du comportement joint de la masse de bloc obtenu en étude expérimentale a été présenté. Singh et autres (2002) ont effectué les essais étendus sur les matériaux modèles joints avec plusieurs configurations communes. Les résultats de cette étude expérimentale ont été employés pour valider le comportement de résistance du modèle numérique de la masse de bloc développé en utilisant le logiciel professionnellement disponible. Ce logiciel contient fondamentalement différents genres de modèles matériels. Des résultats de la simulation numérique on a observé l&apos;existence de l&apos;anisotropie de force. L&apos;anisotropie de force ainsi obtenu a été comparée aux valeurs expérimentales rapportées par Singh et autres (2002). Les courbes anisotropes dans les deux cas sont en U avec la base plate. Elles montrent la bonne concordance dans le mode coulissant de la zone d&apos;échec en lequel la inclination ß (= 90° - ?) du joint placez-Je le mensonge entre 20° à 40° pour tous les cas d&apos;enclencher.

Jan 1, 2003

By Hwasoo Park

Keywords: FC mold, Argon, Magnetic, Solidification Fluid flow plays a significant role for steel surface quality in continuous casting. Although water models are able to model fluid flow, the flow pattern is not of practical interest since the moving solidifying shell cannot be accounted. Computer modeling has become a popular and dependable tool to analyze the complex phenomena occurring in continuous casting. In this numerical simulation we coupled momentum and energy equations considering the solidification of steel along with turbulence ?-e model. The vertical bending mold is equipped with bifurcated nozzle with 15 degrees of port angle. FC mold is supplied with DC induced magnetic fields. Influence of argon and magnetic fields on molten steel flow and heat transfer is investigated by solving momentum equations adding source terms due to multiphase ? magneto hydrodynamics. The results show that magnetic fields reduce the jet velocity from port by lower magnetic core induced field but the influence on meniscus velocity is complex to describe. Magnetic fields increase surface temperature and homogenize its distribution on surface and inside a mold.

Jan 1, 2005

By Y. M. Chen, H. L. Zhang, W. U. Bing, C. Q. Zhou

Multiphase flow in the widely used Vanyukov furnace was numerically studied. An unsteady three-dimensional and three-phase flow model was firstly built using the computational fluid dynamics (CFD) software ANSYS FLUENT®, and then solved with the volume of fluid (VOF) and k - ??model. The results showed that the proposed model could be used to predict the multiphase movement, the slag/air fluctuation, the vortex formation, and effects of structural and operational parameters. By fast Fourier transform (FFT), the dominant frequency of density with time signal was calculated as 0.29 Hz. The analysis of different injection flow rates of enriched air indicated that this variable has a major effect on the mean slag velocity. The peak mean velocity increased from 2.17 to 4.99 m/s while the flow rate of enriched air varied from 70 to 160 m/s. The proposed model provides a method to optimize the furnace structure and operating conditions for the best furnace performance and lowest energy consumption.

Jan 1, 2015

By Toshio Suzuki, Machiko Ode, Won Tae Kim, Seong Gyoon Kim

"The Ostwald ripening and growth and coarsening of the secondary dendrite arm are studied using a phase-field model. The phase-field parameters are obtained at a thin interface limit proposed by Kim et al. In Ostwald ripening analysis, the difference of the concentration due ·to· the curvature effect is shown. The third power of the average particle radius Rav is proportional to time and the ripening constant shows a good agreement with the value by Deppers et al. The particle radii distribution shows the transition behavior to attain a self-similar stationary shape. The growth and coarsening behavior of the secondary dendrite arm is also analyzed using a phase-field model. The estimated arm spacing and its dependence on time are in good agreements with experimental data.IntroductionThe ripening or coarsening process controls the microstructure during the late stage of solidification. The ripening phenomena are recognized firstly by the Ostwald [1]. The first quantitative analysis for the Ostwald ripening in the limit of vanishing solid fraction was provided by Lifshitz, Slyozov and Wagner (LSW) [2,3]. In their study, that the mean radius grows proportional to the cubic root of time and the size distribution of precipitations shows a ""self-similar"" shape are shown. Most of the studies for the Ostwald ripening confirm or predict these two general rules above. In case of finite solid fraction, the general rules are observed but the ripening constant is larger and the self-similar shape is wider than the LSW prediction respectively [4-6]. To analyze· the high solid fraction system, the theories that consider the solute field correlation are proposed [4,7-8]. In the theory, however, since the particle is treated as the point, the effect of the concentration field deviation from the monopole field should be taken into account. Recently the numerical simulation, in which not only diffusive transportation of solute but also convection effect is considered, has been done using a phase-field model [9]."

Jan 1, 2000

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