Search Documents

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'existence de l'anisotropie de force. L'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'échec en lequel la inclination ß (= 90° - ?) du joint placez-Je le mensonge entre 20° à 40° pour tous les cas d'enclencher.

Jan 1, 2003

By Yuting Xue, Danqing Gao, Brijes Mishra

Heterogeneity and discontinuity significantly affect rockmass strength. Past studies have largely used arbitrary scaling approaches to produce rockmass strength from laboratory rock tests. This paper proposes a stochastic approach to produce rockmass strength from laboratory rock tests. Based on the laboratory data, MATLAB with an extreme value stochastic model generated a database for each physico-mechanical property of the rock. Then finite-difference software FLAC2D simulated laboratory-sized rock specimens. Random material properties in MATLAB were generated using the grids developed in the numerical model. The random properties were then applied to the final FLAC2D model. Model runs simulated the approach performed in the laboratory. The results from the model indicate that a stochastic approach produced compressive strengths that were lower than those from the deterministic approach. Failure modes for each specimen were different from those observed in the laboratory. In addition, random density also influenced 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 D. F. Malan, J. A. L. Napier

A number of shallow coal or hard rock mines employ pillar mining systems as a strategy for roof failure control. In certain platinum mine layouts, pillars are designed to ‘crush’ in a stable manner as they become loaded in the panel back area. The correct sizing of pillars demands some knowledge of the pillar strength and the overall layout stress distribution. It is particularly important to understand the impact of the layout geometry on the effective regional ‘stiffness’ of the rock mass around each pillar. An important design strategy is to model relatively detailed layout configurations which include a precise representation of the local pillar layout geometry and to analyse multiple mining scenarios and extraction sequences to select optimal pillar sizes and barrier pillar spacing. Although computational solution techniques are now impressive in terms of run time efficiency, a major difficulty is often encountered in assigning suitable material properties to the pillars and in devising an effective material description of the layered rock strata overlying the mine excavations. This paper outlines an efficient numerical strategy that can be used to assess large-scale pillar layout performance while retaining the ability to modify individual pillar constitutive behaviour. The proposed method is applied to selected layouts to compare estimated average pillar stress values against values determined by detailed modelling and against observed behaviour.

May 3, 2023

By R. L. Lewis

The White Pine Mine is the second largest active underground copper mine in the United States, covering an area of approximately 34 km2 in the Upper Peninsula of Michigan. A numerical flow model of the White Pine Mine was constructed to evaluate engineering controls that would reduce the rate of mine water discharge after closure. Of the passive controls evaluated, the placement of low-flow bulkheads within the mine workings was predicted to offer the greatest benefit, reducing the rate of discharge by as much as 35 percent.

Jan 1, 1996

By Laurentiu Nastac

"An advanced casting simulation-based optimization approach was applied in this study to assist in the improvement of the mold design of aerospace components made of A3 56 alloy. By using this approach, mold filling and solidification related defects (including macro-shrinkage and shrinkage porosity) were significantly minimized and hence it helped in cost reduction, performance enhancement and quality assurance of complex A3 56 cast parts.An experimental validation and detailed calibration procedures of the models for prediction of macro-shrinkage and shrinkage porosities were performed using A3 5 6 plates cast in fiiran-silica sand molds using the Prometal's Rapid Casting Technology (RCT) - mold printing technology. Correlations between Niyama values and the pore percentage were also developed. Thus, the severity level of shrinkage porosity can be determined via the Niyama criterion. Predictions were then compared with the macro-shrinkage and porosity measurements in plates of various plate thicknesses and in other commercial A3 56 casting components.IntroductionThis paper describes the process modeling approach used in the development of the mold design technology used to produce assist in the improvement of the mold design of aerospace components made of A356 alloy. Nova Flow&Solid™ software [3] was applied in this study. This casting simulation software was previously validated for several alloys including Ti, steel and superalloys [1-4]. The objective of this paper is to calibrate and validate the simulation software for A356 alloy sand mold castings. For this purpose, an experimental validation and detailed calibration procedures of the models for prediction of macro-shrinkage and microshrinkage (e.g., shrinkage porosities) were performed using A356 plates cast in fiiran-silica sand molds using the Prometal Rapid Casting Technology (RCT) mold printing technology. The details of the RCT technology are shown in Fig. 1 [5]. Prometal's RCT technology is used for rapid production of sand molds and cores directly from CAD files. The RCT technology relies upon ink jet printer technology to selectively bind the grains of sand to create three-dimensional sand molds and cores. Mold filling and solidification modeling as well as stress analysis are typically performed to properly design the casting and the rig system of each prototype component."

Jan 1, 2012

By Cheng Li, Yanlei Wang, Lingrui Kong, Jianxin Tang, Fang Yuan

Coal mined-out areas have a significant impact on mining safety and affect the pressure distributions of neighboring coal faces. The PFC3D software is used to study the goaf mechanical characteristics, and the core loading test verifies the rationality of the numerical simulation method. The results from the numerical modeling showed that the void ratio, particle size ratio, and particle size distribution of the rock particle sample (RPS) had a significantly remarkable influence on the compression deformation characteristics of the rock block. The parallel-bond tensile strength and parallel-bond cohesion had a notable impact on the compression deformation traits of the rock block only when the axial stress level was lower than 12 MPa. By contrast, when the axial stress exceeded 5 MPa, the effective modulus significantly influenced the compression deformation features. However, the stiffness ratio, loading rate, and parallel-bond friction angle had almost no effect on the compression deformation characteristics. Furthermore, when the stress level was less than 10 MPa, the stress–strain curves from the rock core samples’ numerical simulation and loading tests were approximately linear. In addition, the Pearson correlation coefficients of the stress–strain curves suggest that the goaf zone’s core numerical simulation results had a robust correlation with the loading test results. The findings of this study can serve as a reference for the field application of PFC3D software and safety assessments in goaves.

Jan 18, 2022

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 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 K. Hashiba, Y. Z. Liang, K. Fukui, M. Koizumi

"Efficiency and penetration rate during percussive rock drilling significantly affect the speed of tunneling and excavation. Many researchers and engineers have continuously engaged in the advancement of rock drilling performance. Development and modification of the rock drills based on test data are timeconsuming and costly. Hence it is essential to make effective use of numerical simulation. However, most of the numerical models proposed in previous studies were too simple in comparison with the actual percussive rock drilling in which many factors are involved. The objective of our research project is to construct a new realistic numerical model of continuous percussive drilling, which is based on the previous Okubo-Nishimatsu’s model. In this study, the parts of the previous model involved in stress wave propagation and attenuation were modified. The two models were prepared for HD712 hydraulic rock drill mounted on crawler drills and HD210 hydraulic rock drill mounted on drilling jumbos. Stress wave propagation was computed with the discretized one-dimensional wave equation. The calculated results were compared to the stress waves measured on a rod during percussive drilling. In both models for HD712 and HD210, the measured stress waves can be reproduced with the one-dimensional model of the same acoustic impedance as the actual piston, shank rod and rod combined with the CI+spring model as a rod joint. The calculated results indicate that the spring constant in the model of a rod joint needs to be set appropriately depending on drilling conditions. INTRODUCTIONEfficiency and penetration rate during percussive rock drilling significantly affect the speed of tunneling and excavation. Many researchers and engineers have continuously engaged in the advancement of rock drilling performance. Components related to the percussive drilling are a drill body, a piston, a shank rod, rods, rod joints, a bit and rock. The piston in the drill body is reciprocated by hydraulic pressure and collides with the shank rod a few thousand times per minute. Stress waves generated by the collision propagate and attenuate in the rods and rod joints, and then reach the bit. The bit penetrates and crushes rock. During percussive drilling, the components mutually affect each other and self-induced oscillation occurs. The number of blows and the rate of penetration are determined by the practical drilling conditions such as hydraulic pressure, thrust force, number of rods and mechanical properties of rock. Development and modification of the rock drills based on test data are time-consuming and costly. Hence it is essential to make effective use of numerical simulation. However, most of the numerical models proposed in previous studies were too simple in comparison with the actual percussive rock drilling in which many factors are involved."

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