Search Documents

By Shinji Muraishi, Shinji Kumai, Yasutaka Kedo, Takashi Kambe

"Magnetic pulse forming (MPF) is one of the extreme high-speed forming processes. A metal sheet which is accelerated by electromagnetic force is collided onto a mold and fine and sharp details of the product can be obtained. The deformation mechanism of the MPFed metal cannot be understood only by the microstructure observation of the deformed metal. The deformation behavior of the metal should be investigated also by considering the local microstructure change with strain, strain rate and temperature change during deformation. Numerical analysis is an effective method to achieve this requirement. In this study, pure Al and Cu sheets were deformed by MPF. Impact velocity of the metal sheet was calculated by using ANSYS Emag-Mechanical. Based on the obtained impact velocity, the deformation behavior of the metal sheet was reproduced, and local strain, strain rate, and temperature were analyzed by using ANSYS AUTODYN. The relationship between the local microstructure change of Al and Cu sheets and the deformation parameters was examined.INTRODUCTIONMagnetic pulse forming (MPF) is one of the high-speed forming processes using electromagnetic force. Fine and sharp details of the product can be obtained, and springback after deformation is extremely reduced in comparison to the conventional forming processes(Psyk, Risch, Kinsey, Tekkaya, & Kleiner, 2011). Figure 1 shows a schematic diagram of MPF. A metal sheet is placed over the coil in the discharge circuit. The discharge current runs through the coil, and magnetic fieldis producedaround the coil. An eddy current is inducedin the metal sheet. Interaction between magnetic field and eddy current creates electromagnetic force. The metal sheet accelerated by electromagnetic forcecollides ontothe mold. Forming is completed in a few tens of micro-seconds.MPF is suitable for metals with a high electrical conductivity like Al and Cu.The MPFed metal may have characteristic microstructure because of high strain rate deformation."

Jan 1, 2018

By A. Manouchehrian, M. Cai

"Rockburst, which often occurs around underground openings in highly stressed hard rocks, is a sudden and violent unstable rock failure and is associated with energy release. Rockburst is very dangerous to workers and mining equipment, and it has been a major concern for deep underground construction and mining. Before excavation, rocks are loaded under a true-triaxial condition. Upon excavation, the rocks at the excavation boundaries are loaded in the maximum principal stress direction and unloaded in the minimum principal stress direction and this provides a necessary condition for unstable rock failure to occur. Understanding rock behavior under true-triaxial unloading condition is critical for developing measures to control unstable rock failure in underground construction.In this article, an explicit FEM tool is used to simulate rock behavior under true-triaxial unloading condition, with a focus on analyzing stable and unstable failure and the released energy during failure. Firstly, a uniaxial compression test was simulated. The results indicated that the adopted numerical tool was able to simulate stable and unstable failure by varying the Loading System Stiffness (LSS). In this research, Loading System Reaction Intensity (LSRI) was proposed as an indicator to identify stable and unstable failure. Secondly, unstable rock failure under true-triaxial unloading condition was simulated. A rock specimen was loaded in three directions; the lateral pressures (s3 and s2) were kept constant and the vertical stress (s1) was increased to a pre-defined stress. The specimen was suddenly unloaded in the minimum principal stress direction (s3) and loading in the s1 direction was continued until failure occurred. Different LSS values were considered and the results showed that the LSS influenced the failure mode; a soft loading system resulted in more violent failure and released more energy than that in a stiff loading system. The influence of specimen height to width ratio on rock failure was also investigated. It was found that failure in taller specimens was more violent than that in shorter specimens. Insights gained from the simulation are useful for better understanding of the influencing factors on the released kinetic energy during unstable rock failure and to find means to control rockburst damage in engineering practice."

Jan 1, 2015

By N. E. Mkadmi

Backfilling increases the stability of underground stopes and may help reduce the environmental footprint from the mining operation. The barricades holding the backfill in place must be designed to prevent failure based on the total and effective stress states, which depend on fill properties, stope geometry, and filling rate. This paper presents results from a numerical study aimed at analyzing the stress states associated with the response of paste backfill following stope filling, with an emphasis on the short term behavior. The simulations illustrate how stresses evolve with backfill drainage and filling rate in a vertical stope.

May 1, 2011

By M. Gimpel, R. Galler, T. Antretter, H. Kargl, M. Mikl-Resch, R. Tichy, W. Ecker, G. Pittino

"A computational method based on a continuum-mechanical approach to simulate rock cutting is developed and evaluated. The aim is to obtain the local loading acting at the cutting pick during operation. The explicit FEM is chosen to model the pick-rock interaction. The presented method takes advantage of (i) an advanced material model (ii) a robust contact algorithm and (iii) fast calculations due to parallelization. A damage-plasticity law is chosen and modified with user-subroutines. A user integrated element deletion routine evades critical mesh distortions. Infinite elements at the model boundaries avoid unrealistic deformation wave reflection.INTRODUCTIONExcavation in mining and tunneling is a costly and complex process. Special equipment such as the roadheader (fig. 1) is employed to cut the rock in the excavation process but high material strengths limit their usability (Tockner 1999). Nowadays a challenge is to extend their field of application to harder and more abrasive rock materials. The forces on the pick (fig. 2) are of great importance and determine wear and especially the efficiency of the machine. Figure 2 shows the principle of cutting and the predominant forces. The variation of pick geometry, cutting depth or attack angle will result in an alteration of these forces and determine the efficiency of the process (Tockner 1999). The mechanisms of rock cutting follow the principle of indentation (Kou et al. 1999). Directly under the indenter a hydrostatic stress state evolves causing the material to fail only due to crushing and fragmentation mechanisms. In the close surroundings large strains emerge and tensile cracks are formed but in a further distance the material behaves elastically."

Jan 1, 2015

By M. Gimpel, R. Galler, T. Antretter, H. Kargl, R. Tichy, W. Ecker, M. J. Mikl-Resch, G. Pittino

"A computational method based on a continuum-mechanical approach to simulate rock cutting is developed and evaluated. The aim is to obtain the local loading acting at the cutting pick during operation. The explicit FEM is chosen to model the pick-rock interaction. The presented method takes advantage of (i) an advanced material model (ii) a robust contact algorithm and (iii) fast calculations due to parallelization. A damage-plasticity law is chosen and modified with user-subroutines. A user integrated element deletion routine evades critical mesh distortions. Infinite elements at the model boundaries avoid unrealistic deformation wave reflection.INTRODUCTIONExcavation in mining and tunneling is a costly and complex process. Special equipment such as the roadheader (Figure 1) is employed to cut the rock in the excavation process but high material strengths limit their usability (Tockner 1999). Nowadays a challenge is to extend their field of application to harder and more abrasive rock materials. The forces on the pick (Figure 2) are of great importance and determine wear and especially the efficiency of the machine. Figure 2 shows the principle of cutting and the predominant forces. The variation of pick geometry, cutting depth or attack angle will result in an alteration of these forces and determine the efficiency of the process (Tockner 1999). The mechanisms of rock cutting follow the principle of indentation (Kou et al. 1999). Directly under the indenter a hydrostatic stress state evolves causing the material to fail only due to crushing and fragmentation mechanisms. In the close surroundings large strains emerge and tensile cracks are formed but in a further distance the material behaves elastically."

Jan 1, 2015

By F. S. Su, P. Z. Pan, X. T. Feng, F. Yan, W. W. Ji

"This paper aims to simulate the spalling failure behavior in hard rock based on a self-developed numerical model, i.e. a rock discontinuous cellular automaton (RDCA). A crack initiation and propagation criterion is proposed based on the plastic strain, which is used to link the continuity and discontinuity during the failure process. The discontinuities of rock, such as cracks and other internal boundaries, are represented by discontinuous displacement functions, i.e. tip enrichment and step enrichment. The partition of unity is used to improve integral precision of the elements intersected internal boundaries, which is traced by a level set method. Arbitrary crack geometries can be conveniently represented and it does not require any remeshing in the modeling of crack initiation, propagation and coalescence. Numerical modeling of excavation in hard rock using RDCA is conducted. The spalling failure process of hard rock induced by excavation is well simulated. The depth and shape of hard rock spalling failure are good agreement with in situ observation.INTRODUCTIONThe excavation in rock mass is more and more popular. At the same time, the disasters, such as rock burst, rock spalling and falling etc, induced by excavation are quite prominent. Numerical modeling is an important way to help us understand the failure mechanism of rocks. Many numerical methods have been used to simulate the excavation of rock mass. For example, the finite element method (FEM) (Karakus and Fowell, 2003), the discrete element method (DEM)(Cundall and Strack, 1979, Labra et al., 2008), Finite difference method (FDM)(Zhu et al., 2009), Element Free Galerkin Method (Belytschko et al., 1994), Discontinuous displacement analysis (DDA) (Shi, 1988) and Boundary element method (BEM) (Katsikadelis, 2002, Chen and Hong, 1992) are commonly found in the stability analysis of excavation in rocks. In the modeling of crack propagation or excavation, most of these methods are tedious since the internal boundaries have to be consistent with the element boundaries. And in many cases, the remesh has to be conducted. In recent years, people are interested in particular types of FEMs, such as the extended finite-element method (XFEM), the particle discretization scheme finite-element method(PDS-FEM(Hori et al., 2005, Oguni et al., 2009)), the combined finite-discrete element method (FEM/DEM) (Munjiza and Wiley, 2004, Cai, 2008) and spectral-element method(Gharti et al., 2012) etc."

Jan 1, 2015

By P. Carlone, A. Astarita, F. Rubino

Cold gas dynamic spray is an innovative technology allowing the deposition of thin metallic layers on a bulk or sheet substrate. Titanium on aluminum deposition is very attractive to the aeronautic industry, due to the enhanced corrosion resistance and wear properties of titanium coating as well as its improved compatibility with CFRP, preserving the reduced cost of aluminum substrate. One main drawback, however, is related to the micro-porosity of the deposed layer, negatively affecting the corrosion barrier performance. In this paper, the effect of selective laser post-treatment on pure grade 2 titanium coatings on AA2020-T3 sheets was experimentally and numerically investigated. Morphological features, microstructure and chemical composition of the treated zone were assessed by means of optical microscopy, scanning electron microscopy and energy dispersive X-ray spectrometry. Microhardness measurements were also carried out to evaluate the mechanical properties of the treated coating. A numerical model of the laser treatment, based on a finite volume scheme, was implemented and solved to simulate the process and discuss the experimental outcomes. Obtained results highlighted the key role played by heat input and dimensional features on the effectiveness of the treatment.

Jan 1, 2015

By S. Bhattacharjee

"Rock failure process is generally recorded or measured using conventional sensors like strain gauges, seismic, micro-seismic, acoustic, electromagnetic and other techniques. The above sensors are either to be pasted on the surface of the rock sample or placed very near to it to capture the response of the sample with the increment of load. This requires precise tools to cut or polish the locations where sensors are to be placed. Moreover, deformation or acoustic or micro-seismic activities are recorded only at the locations where sensors are pasted.Recent development in Digital Image Correlation (DIC) technique has shown promise to measure deformation / strain at every location of the surface of a sample. In DIC, two images (reference and deformed) are taken within a small time interval of the same object are compared to estimate the change in reflectance as a function of displacement, velocity or rate of strain at each pixel while the object is undergoing through a loading process. In this study, extended digital image correlation (X-DIC) technique has been implemented using 4-noded quadrilateral elements with a discontinuity surface in the deformed image. The continuous and enriched displacements are estimated at each node.The X-DIC method is verified using numerically generated images. The study has been further extended to evaluate the performance of X-DIC method on deformation analysis of experimental samples under uni-axial loading condition. The results of this study show that X-DIC method is an effective tool for determination of displacements of the entire surface of a rock sample having a discrete crack or discontinuity in a non-contact manner. This method demonstrates the applicability to investigate rock failure mechanism for the entire sample in a laboratory set up as well as monitoring of pillar failure mechanism in underground mines."

Jan 1, 2015

By Mahmoud Shakeri1

Various types of thin-walled structures have been frequently used as energy absorbers to mitigate the adverse effects of impact and hence protect the vehicle structure under consideration when subjected to different loading conditions. Compression tubular structures that buckle in a progressive manner, can be used to advantage, providing inexpensive and versatile, but nevertheless weight and stroke efficient energy absorbers. In this paper quasi-static nonlinear finite element simulations are performed to study the effects of number of corners of the cross-section, wall thickness, plastic hardening rate of material, foam-filler and load angle on deformation mode and mean crushing force of tube. Crush strength increases as the number of corners of the cross-section increases, though it almost saturates for the number of corners beyond 11. Higher plastic hardening rate stabilizes the collapse pattern. Foam-filled tubes seem to have higher critical load angles in oblique loading and display a more stable and gradual decrease in mean load over the transition region. Substantial decrease in the initial peak load and dynamic load at any a given deflection was observed as the load angle increases. The investigation of bitubal arrangements and Multi-cell tubes reveals that these may be preferable to monotubal specimens due to the presence of the inner profiles. In addition, a comparison is conducted between the numerical and experimental results of the mean forces of some special angle elements. Good agreement is obtained.

Oct 1, 2011

By D. D. Tannant

An experimental and numerical study was conducted on the geomechanical behaviour of oil sands during penetration tests in which a steel wedge with an apex angle of either 20 degrees or 30 degrees was forced into compacted oil sands. The boundary conditions simulated in the laboratory tests are similar to shovel penetration into oil sands. The objective of the laboratory tests was to measure the force required to push a steel wedge into compacted oil sands. A wedge penetration test differs from conventional triaxial tests in that much more shear is generated. The laboratory data was needed to improve calibration of micro-mechanical input parameters for a two-dimensional discrete element model of oil sands constructed with Particle Flow Code (PFC). Initial model calibration was performed with triaxial test data and then a PFC model of the wedge penetration test was constructed. The measured force-displacement behaviour of the steel wedge was compared with PFC modelling results. It was found that different combinations of PFC particle and bond parameters could equally well replicate the laboratory triaxial compression data for bitumen-rich oil sands. When the wedge model was run, it was found that the resulting force required to penetrate the wedge into the oil sands was about four to six times higher than that measured. This clearly indicates that further research is required to improve the predictive capability of PFC models of oil sands.

Jan 1, 2002

By B. Y. Park

This paper develops a series of three-dimensional simulations for the Bayou Choctaw Strategic Petroleum Reserve. The U.S. Department of Energy plans to leach two new caverns and convert of one of the existing caverns within the Bayou Choctaw salt dome to expand its petroleum reserve storage capacity. An existing finite element mesh from previous analyses is modified by adding two caverns and increasing the size of existing Cavern 102. The structural integrity of the three expansion caverns and the interaction between all the caverns in the dome are investigated. The impacts of the expansion on underground creep closure, surface subsidence, infrastructure stability, and well integrity are quantified. The results show that from a structural viewpoint, the locations of the two newly proposed expansion caverns are acceptable, and all three expansion caverns can be safely constructed and operated.

May 1, 2009

By W. Swidzinski

"Zelazny Most" post-flotation copper mine tailings pond is the only place to of storeing wastes for KGHM POLSKA MI EDZ S.A. The flotation tailings is are the main waste stream in copper mining in Poland. Sufficient storage capacity for tailings is crucial for , thus uninterrupted production of copper is crucially related to the possibility of storing wastes in the pond. The pond covers a Its over 14 km2 arca, almost 50 m high dams making this dame this construction one of the largest in the world. The risk assessment of potential damage due to un stable ility loss of dams or environmental contamination, as well as safely constructing damsincrease of dams within the design criteriaprocess, require quick access to a large amount of the monitoring data , which in vast amount have to be monitored, that needs to be analyzed and processed and reported, every day. This In the paper describes the numerical information system "SyZeM" , used for designated to the "Zelazny Most" pond, baseding on the Oracle database and Geographical Information System , to manage the graphical and descriptive data in coherent way., is presented. The main structure of the system, as well as the methodology of storing, processing and reporting the geological, geotechnical, hydrogeological, geophysical, chemical and survey data are shown. The aspects of Optimizing optimal and risk mitigation and safe management of copper wastes and risk mitigation with the help of "SyZeM" are discussed.

Jan 1, 2007

By Jim Hazzard

A numerical investigation is presented to examine fluid flow in fractured rock slopes. In mining geomechanics, different numerical modelling approaches are used to estimate fluid pressures in open pit mines. The most common technique is the Equivalent Porous Medium (EPM) approach in which the fractured rock is assumed to behave as a porous continuum. Depending on the geometry and properties of the problem, this may not always be the best approach. The paper presents a parametric study to examine the effect of fracture spacing, fracture and matrix diffusivity and excavation rate on the fluid behaviour in rock slopes. Based on the model results, guidelines are provided as to when EPM models can be used and when different modelling approaches (e.g., discrete fracture flow) may be necessary.

Jun 1, 2012

By M. B. Dusseault, A. Pirayehgar

"There is concern that hydraulic fracturing (HF) near an existing critically stressed fault could trigger slip and generate a seismic event of significance. To this end, a mathematical discontinuum model representing a naturally fractured rock (NFR) mass containing a fault is created, and fault slip response to fluid injection is studied. UDECTM is used as a Discrete Element Method (DEM) approach to define the seismic moment, and also to assess the influence of injection and slip on further shear slip propagation. In most of the simulations, slip and dilation take place along a number of oriented fracture surfaces, so it is reasonable to assume that in a real case, microseimic emissions would accompany each slip “event”. The moment magnitude scale (Mw) is used to measure the size of the microseismic event in terms of the energy released. Seismic moment is a measure of the total energy (work) released during a seismic event and is used to measure its size. From another aspect, the injection energy is calculated as the total energy required for achieving the hydraulic fracturing; this is viewed as energy to open fractures against the ambient stress field and much of it is potential stored energy that can be released and create discontinuities (work of fracture), deformations (W=F·d), heat and radiated seismic energy. A small amount of the total released energy is transformed to radiated seismic energy and another small amount results in tensile fractures. Fracture energy is used to calculate the work required to create fracture surfaces. Comparing the energy release by fracturing or seismic events indicates their contribution to the total amount of released energy or total work. The amount of energy associated with the stick-slip event is small compared with the overall work done to generate the distortions (i.e. increases in aperture by forcing them open) in the rock mass, and much of the energy in slip is lost as heat from the plastic deformation, rather than as radiated seismic energy.Microseismic monitoring shows the spatial distribution and magnitude of seismicity associated with slip of bedding planes and natural and incipient fractures. Work calculations have been undertaken to look at the various sources of energy storage (elastic strain energy) or dissipation (slip), and likely the energy losses from viscous dissipation could be addressed as well, helping to clarify HF processes in NFRs."

Jan 1, 2015

By U. Caliskan, S. Miskovic

Gravity separation is a technique that depends on the difference in density of minerals for their separation. Although this method is one of the oldest separation techniques, it is still preferred in the industry because of its low capital and operating cost while it is environmental-friendly. In this paper, we present a new numerical methodology that can be used to investigate dynamics of dense slurry flows in gravity separation systems. A simple channel with a set of riffles is selected as a test platform in order to observe the effect of the riffle shape on segregation profile of particles with different densities. It is believed that the riffle shape would affect the slurry flow and possibly cause one class of particles to segregate within the riffle zone. The turbulence, as a disturbance effect within the riffle zone, causes heavy particles to displace the light, that is, the heavy particles will be trapped in riffles, while the lighter particles are carried by water to the end of the channel, which is an example to uninterrupted suspension. MP-PIC model in OpenFOAM is used to study the effect of riffle shape and slurry percentage of solids in this simple gravity separator on particle segregation profile of four discrete mineral density classes, which are Quartz/Clay (SG 2.6), Pyrite (SG 4.2), Chalcopyrite (SG 5.1), and Gold (SG 19.32). Gold particles are found to be readily captured and retained in the riffles because of their high density, while particles with the intermediate densities are more influenced by the flow. The developed model framework has been proven robust and capable of providing physical results even for the flows with a very high percentage of solids, characteristic for mineral processing operations. INTRODUCTION Gravity-based ore concentration techniques make use of differences in density of minerals within an ore to enable separation of valuables from the gangue. As some of the oldest techniques used for ore beneficiation, they are known to offer a range of benefits such as low specific energy consumption and low capital and operating costs (Burt, 1999; Falconer, 2003; Holland-Batt, 1998), however, they suffer from low selectivity, poor separation efficiency, and low capacity. Gravity separation techniques are common practices in the mining industry and are employed widely for the precious mineral concentration, iron ore recovery and coal preparation (Kroll-Rabotin, Bourgeois, & Climent, 2013; Laplante & Spiller, 2002; Williams, 1997). For example, gravity separation systems are highly applicable for precious metal enrichment operations by being included in primary grinding circuits, shaking table tailing separation, pre-concentration of leach feed and flotation concentrate enrichment for both flash and conventional flotation purposes (Burt, 1999; Falconer, 2003; Holland-Batt, 1998; Laplante & Spiller, 2002; McAlister, 1992). Gravity separation can be employed in the copper processing circuits by being integrated into the pre-flotation circuit to capture the gravity recoverable gold (GRG) and any free copper aiming to enhance the gold and copper recovery of the plant. When native copper and GRGs are present in the deposit, this integration not only facilitates plant overall performance via removal of the valuables early in the circuit, but it also reduces downstream operational expenditures due to less retention time of coarser copper particles in the system.

Jan 1, 2019

By S. Pietruszczak, G. Su, T. S. Nguyen, E. Haghighat

"In this paper, the mechanical properties of Tournemire argillite are modeled using a mathematical framework which employs the microstructure tensor approach (Pietruszczak & Mroz, 2000). The experimental data incorporates the results of a series of triaxial tests (Abdi & Evgin, 2013) conducted at different orientation of bedding planes relative to the direction of loading. First, a plasticity framework is presented which involves the specification of the conditions at failure as well as the description of the irreversible deformation process. Later, a systematic procedure is outlined for identification of material functions/parameters involved. The paper is concluded by presenting the results of numerical simulations of the triaxial tests. It is demonstrated that the model can adequately reproduce the primary effects of anisotropy in the stable deformation regime as well as the onset of fracture.INTRODUCTIONThe understanding of the mechanical behaviour of argillaceous rocks is important in many types of geotechnical projects, including petroleum extraction, carbon dioxide sequestration as well as deep geological disposal of radioactive waste. In all these cases, the primary concern is the onset and propagation of damage due to excavation, transport of pore fluids and/or elevated temperature. The Canadian Nuclear Safety Commission (CNSC), Canada’s nuclear regulator, will be responsible for the licensing of future deep geological repositories. In order to evaluate the safety of deep geological repositories, CNSC staff keep abreast of scientific developments by collaborating with national and international partners on experimental and theoretical research on the long-term performance of those engineered and natural barriers. Through that collaboration, the CNSC has access to data resulting from experiments performed at different Underground Research Laboratories (URL). Such an URL is situated in Tournemire, France, in an argillaceous rock formation, the Tournemire argillite. The Tournemire argillite is characterized by the presence of closely spaced bedding planes and exhibits a strong anisotropic behaviour, so that its strength as well as deformation properties are directionally dependent."

Jan 1, 2015

By H. Zniber El Mouhabbis

Ore dilution or stope overbreak increases the cost of production and ultimately reduces the profitability of a mining operation. The adverse economic impact of ore dilution is due to the added costs associated with the mucking, haulage, crushing, hoisting, milling, and treatment of waste rock or low grade ore in the hangingwall and footwall. One factor that could have a significant influence on stope overbreak is the undercutting of the stope into the walls. Undercutting stope walls is often necessary, in particular in narrow vein mines. Moreover, equipment size has a role in drift sizing where stope undercutting may become a must in order to fit the mobile equipment while respecting legal clearances.

May 1, 2009

By N. Manopulo, J. Maier, L. Tong, C. Becker, P. Hora

The contribution presents a new finite element approach in which the common updated Lagrangian formulation is replaced by a spatial flow formulation applying a special arbitrary Lagrangian-Eulerian for-mulation (ALE). The description separates spatially the process in a die inside zone and in a profile zone. For the description of the inside zone, where the filling process occurs, a special Euler time independent mesh is applied. By this formulation the mesh has to be generated only once, which simplifies the generation of spatially optimized meshes. For mapping of the material flow during the filling process the method uses a flow surface traction method as applied in casting simulations. A further strength of the new method is the computational efficiency which allows the modelling of extrusion process with the highest complexity with a response time of 12 hours, which is important for the industrial optimization process.

Jan 1, 2015

By N. Golchinfar

Tunnels in deep mines are subjected to high stresses. Fracture and failure of the rock mass around an opening can occur when mining-induced stresses are high. In burst-prone grounds, mining-induced seismicity can cause additional dynamic loading which can further increase the stress around the tunnels. When the dynamic stress reaches the rock mass strength, fracturing can occur and rockburst may happen (depending on loading stiffness around the failed rock mass). Rock support installed in the tunnel must be capable of dissipating dynamic energy and holding the failed rocks. Hence, it is important to estimate the depth of failure for rock support design. This paper focuses on how to model depth of failure under dynamic loading. In the modeling, two scenarios of rock failure are considered. In the first scenario, rock failure occurs under static loading and subsequent dynamic loading further increases the depth of failure. In the second scenario, no failure occurs under static loading but in the subsequent dynamic loading, failure around the tunnel is created.

Aug 1, 2013

By M. N. Vu, G. Armand, A. Noiret, D. Seyedi

Extensive experimental observations have been performed around drifts and shafts at the Meuse Haute-Marne (North-Eastern France) Underground Research Laboratory (URL) to assess the extent and pattern of the induced fracture networks in Callovo-Oxfordian claystone. Two types of fractures can be distinguished, namely; shear (mode II) and extension unloading (mode I) fractures. The main goal of the present paper is to provide insights on the extent of damaged zone around the drifts by numerical simulations. A two-step approach is considered. First, 2D elastic calculations are performed for simulating the excavation of a drift following the direction of the major in situ horizontal stress. The evolution of the stress and strain tensors around drifts are calculated by applying an elliptic convergence function on the drift wall. Two failure criteria are then verified for the induced damage. An equivalent strain and an elastic energy based damage criteria are considered. The obtained results show that the employed simplified approach provides rather good results for the extent of the damaged zone for drifts parallel to major horizontal stress. However, post-processing elastic results leads to overestimating the vertical extent of the damage zone and underestimating its lateral extent, where the in situ stress state is quasi-isotropic in a section of drift. Nevertheless, the obtained results show that an energy based damage criterion could be considered as a potential way for reproducing the extent of the damaged zones if a convergence function can be defined for the drift.

Jan 1, 2015