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

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

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 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 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 Yuting Xue, Danqing Gao, Brijes Mishra

Heterogeneity and discontinuity significantly affect rock strength. For accurate stability prediction, intact rock behavior is imperatively included in rock mass behavior. However, past research largely used arbitrary scaling approaches to produce rock strength. This paper considers a stochastic approach in order to determine the strength of a rock. Based on the laboratory data, a MATLAB® with extreme value stochastic model generates a database for each physico-mechanical property. Then FLAC® simulates laboratory-sized rock specimens. The grids developed in the numerical model can in turn develop random material properties in MATLAB®, which researchers then applied to the final FLAC® model. Model runs simulate the approach performed in the laboratory. The results from the model indicate that a stochastic approach produces strengths that are lower than a deterministic approach. Failure modes for each specimen are different, similar to observations in the laboratory. In addition, random density also influences the failure mode, highlighting the importance of stochastic analysis in rocks.

By Yuting Xue, Danqing Gao, Brijes Mishra

Heterogeneity and discontinuity significantly affect rock strength. For accurate stability prediction, intact rock behavior is imperatively included in rock mass behavior. However, past research largely used arbitrary scaling approaches to produce rock strength. This paper considers a stochastic approach in order to determine the strength of a rock. Based on the laboratory data, a MATLAB® with extreme value stochastic model generates a database for each physico-mechanical property. Then FLAC® simulates laboratory-sized rock specimens. The grids developed in the numerical model can in turn develop random material properties in MATLAB®, which researchers then applied to the final FLAC® model. Model runs simulate the approach performed in the laboratory. The results from the model indicate that a stochastic approach produces strengths that are lower than a deterministic approach. Failure modes for each specimen are different, similar to observations in the laboratory. In addition, random density also influences the failure mode, highlighting the importance of stochastic analysis in rocks.

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

In recent years, Computational Fluid Dynamics, CFD, has been commonly utilized in the mining industry to model the fluid flow behavior in underground mine workings. This paper uses CFD modeling to simulate the airflow behavior in underground crosscut regions, where brattice sails are used to direct the airflow into these regions. Brattice sails are cost effective ventilation control devices for temporary or permanent use in underground mining. They can be used to deflect air into the unventilated areas such as crosscut regions. Their design and installation is a fundamental issue for maintaining a sufficient supply of fresh air and achieving effective air circulation and contaminant removal. At the same time, they should have little impact on the mine ventilation system. This paper presents the results of a two-dimensional CFD model, which examines the effects of brattice length on fluid flow behavior in the crosscut regions. The results of this study will assist in understanding the ventilation air behavior and in determining the optimum size of brattice curtain (sails), which provide a highly effective contaminant removal from the unventilated mine areas. This, in turn, helps the mine ventilation designers to meet the mine safety requirements.

Jan 1, 2009

By C. Sun

In this study, numerical simulations of yield pillars employing the finite element method with time-dependent Drucker-Prager material model are conducted. The stress control mechanism of yield pillars on weak floor strata is studied. The development of pillar yield zone, the distribution of pillar stress and the development of failure zone in the weak floor strata are investigated. The study indicates that a properly selected yield pillar size may significantly reduce the peak stress induced in the weak floor strata and effectively control the development of failure zone in the floor.

Jan 1, 1999

By Congling Shi, Marek Borowski, Yue Wang, Jianwei Cheng, Xixi Zhang, Xuan Xu

As the first barrier to prevent accidents caused by explosions in underground coal mines, the integrity and stability of mine seals have great effects on mine workers’ lives and underground production safety. In this paper, the mechanical responding performance of typical mine seals under the impact load are studied by using numerical simulation method. By such simulation studies, the safety of different seals subjected to the explosion load is evaluated. Seals’stability are also investigated under effects caused by various influence factors including thickness of mine seal, cutting depth into surrounding rocks, rock type, and the roof-tofloor convergence under given the explosion load. The research results can provide theoretical basis and reasonable improvements for mine seal’s construction practices for improving its anti-explosion performance.

By Congling Shi, Marek Borowski, Yue Wang, Jianwei Cheng, Xixi Zhang, Xuan Xu

As the first barrier to prevent accidents caused by explosions in underground coal mines, the integrity and stability of mine seals have great effects on mine workers’ lives and underground production safety. In this paper, the mechanical responding performance of typical mine seals under the impact load are studied by using numerical simulation method. By such simulation studies, the safety of different seals subjected to the explosion load is evaluated. Seals’stability are also investigated under effects caused by various influence factors including thickness of mine seal, cutting depth into surrounding rocks, rock type, and the roof-tofloor convergence under given the explosion load. The research results can provide theoretical basis and reasonable improvements for mine seal’s construction practices for improving its anti-explosion performance.