Search Documents

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

By Tu Bao

HP-Pile foundations are the most favorable type of foundations for bridges because of their constructability, versatility with regard to subsoil conditions, and settlement controllability. For many years, bridge engineers have puzzled over the most economical way to orient HP-Piles: Should the strong axis be oriented in the longitudinal direction or in the transverse direction? The accurate design of foundations may considerably reduce construction costs. In order to answer this question, 3-D finite element models will be set up to analyze the effects of pile orientation on bridge behavior. Because seismic loading usually controls the structural design in seismic zones, it is important to investigate the effects of pile orientation under such loading. HP-Piles are designed conservatively due to the uncertainty of the dynamic soil-structure interaction during earthquakes. For the bridge substructure, the soil?s stiffness plays an important role in load distribution as the soil and piles act as a combined system to resist the seismic loading. The dynamic soil-structure interaction mechanisms of bridges in both the longitudinal and the transverse directions will be investigated through the use of 3-D finite element models. The numerical models will realistically simulate the soil-pile-abutment-wingwall-superstructure system subjected to seismic loads. Nonlinear soil springs will be employed to model the non-elastic soil behavior. The analysis will provide theoretical support to the design practice. KEYWORDS: HP driven piles, finite element, pile orientation, soil-structure interaction

By Andrew Rietz, Yu Bao, Steven Halewski

"HP-Pile foundations are the most favorable type of foundations for bridges because of their constructability, versatility with regard to subsoil conditions, and settlement controllability. For many years, bridge engineers have puzzled over the most economical way to orient HP-Piles: Should the strong axis be oriented in the longitudinal direction or in the transverse direction? The accurate design of foundations may considerably reduce construction costs. In order to answer this question, 3-D finite element models will be set up to analyze the effects of pile orientation on bridge behavior. Because seismic loading usually controls the structural design in seismic zones, it is important to investigate the effects of pile orientation under such loading. HP-Piles are designed conservatively due to the uncertainty of the dynamic soil-structure interaction during earthquakes. For the bridge substructure, the soil’s stiffness plays an important role in load distribution as the soil and piles act as a combined system to resist the seismic loading. The dynamic soil-structure interaction mechanisms of bridges in both the longitudinal and the transverse directions will be investigated through the use of 3-D finite element models. The numerical models will realistically simulate the soil-pileabutment- wingwall-superstructure system subjected to seismic loads. Nonlinear soil springs will be employed to model the non-elastic soil behavior. The analysis will provide theoretical support to the design practice.IntroductionAs sustainability has asserted itself in the public’s eye, bridge professionals have rushed to develop criteria for sustainable bridges. The three main categories that define and distinguish sustainable bridges are: lower energy input, increased durability, and simplified deconstruction. By optimizing the structural members, the engineer is able to reduce the material quantities for the bridge, effectively reducing the energy input for the bridge (Hunt; Whittemore). A deeper understanding of soilfoundation- structure mechanism will facilitate sustainable bridge design by lowering energy input and increasing the service life of bridges (Rietz et al.)."

Jan 1, 2017

By PDF

Flotation columns are widely used in mineral processing and coal preparation industry today. Mathematical modeling of flotation columns has conventionally been done using axial dispersion or tank-in-series models. However, for design and scale up purposes these models are not adequate enough because they do not provide detailed information of the flow fields within the column. This paper presents a Two-Dimensional, Two-Phase. Fluid-dynamic Model and its numerical solutions for some cases encountered in column flotation. The MAC (Marker and Cell) method was employed to obtain the numerical solutions, and based on this method a computer code COLUMNFLOW was developed. The simulated flow field reveals the presence of a circulation flow structure in flotation columns, which could be one of the main reasons for mixing in column flotation. The effects of some parameters such as superficial gas and liquid velocities on the characteristics of the flow structure in flotation columns were also tested with the numerical model.

Jan 1, 1994

By Amalia Giannakou

When designing structures in seismically active areas, foundations of critical structures are typically located away from known faults. Indeed many seismic codes prohibit construction in the immediate vicinity of seismically active faults establishing minimum setbacks from the faults. However, for long structures such as bridges, tunnels and pipelines, a fault maybe unavoidable, and fault rupture risk impossible to preclude. This paper presents a case study of a bridge foundation design that incorporates the consequences of fault deformations passing through the foundation. A series of site surveys were conducted to characterize the fault environment, and fault displacement hazard analyses were conducted to develop design displacement demands. Subsequently, numerical analyses were performed to evaluate the fault rupture-induced demands on a pier foundation in terms of displacements and rotations. Initially, analyses were performed that model the free field fault rupture propagation to the ground surface of the fault to ensure that the basic site characteristics are captured. Subsequently, the foundation elements were introduced in the model and the fault rupture is simulated again. The interaction and influence of the foundation on the fault rupture propagation is calculated, and the consequent demands and distortions of the foundation block provide a basis for structural engineers to design the bridge.

By Orestis Zarzouras, Amalia Giannakou, Jacob Chacko

"When designing structures in seismically active areas, foundations of critical structures are typically located away from known faults. Indeed many seismic codes prohibit construction in the immediate vicinity of seismically active faults establishing minimum setbacks from the faults. However, for long structures such as bridges, tunnels and pipelines, a fault maybe unavoidable, and fault rupture risk impossible to preclude. This paper presents a case study of a bridge foundation design that incorporates the consequences of fault deformations passing through the foundation. A series of site surveys were conducted to characterize the fault environment, and fault displacement hazard analyses were conducted to develop design displacement demands. Subsequently, numerical analyses were performed to evaluate the fault rupture-induced demands on a pier foundation in terms of displacements and rotations. Initially, analyses were performed that model the free field fault rupture propagation to the ground surface of the fault to ensure that the basic site characteristics are captured. Subsequently, the foundation elements were introduced in the model and the fault rupture is simulated again. The interaction and influence of the foundation on the fault rupture propagation is calculated, and the consequent demands and distortions of the foundation block provide a basis for structural engineers to design the bridge.IntroductionThe Izmit Bay bridge, a 3-km-long suspension bridge that crosses Izmit Bay in Turkey, is planned to be constructed in one of the most seismically active places in the world. The south approach structures for the Izmit Bay Bridge along the western side of the Hersek peninsula bring the bridge down from on the order of El. 60 meters at the South Anchorage to an elevated embankment approximately 1.3 km farther south (Figure 1). The south approach structures are located within a zone of secondary deformation around the primary trace of the North Anatolian Fault. The concept design team has developed a concept design for a “viaduct” structure that consists of 10 spans with lengths varying between 136 m and 100 m and two back spans one of 125 m, attached to the main bridge and the other of 74 m attached to the south embankment."

Jan 1, 2017

By A. Pellegrino, F. Collin, C. Guiducci, R. Charlier

A numerical approach for modeling coupled Hydro-Mechanical fracture behavior in 2D is presented in this paper. This approach accounts for the strong dependence of the flow characteristics on fracture apertures and offers an anisotropic description of the flow propagation in the fracture. An interface finite element and a constitutive model are developed . In order to validate the applicability of the model to general geomechanic problems, 2D simulations were performed. Results show the good application of the mathematical model and the fundamental hydraulic role played by the fractures. In diesem Beitrag wird ein Stoffgesetz vorgeschlagen um hydromechanisches Kluftverhalten numerisch zu Modellieren. Dieser Vorgang nimmt die starke Abhängigkeit der Fliesseigenschaften von der Kluftöffnung mit in Anbetracht. Ein Interfaceelement und ein Stoffgesetz wurden entwickelt. Zum validieren der Anwendbarkeit dieses Modells wurden 2D Simulationen durchgeführt. Ergebnisse haben gezeigt dass das mathematische Model anwendbar ist und dass Klüfte eine fundamentale hydraulische Rolle spielen. Une approche numérique pour la modélisation du comportement couplé hydro-mécanique des fractures en 2D est présentée dans ce papier. Cette approche prend en compte la forte dépendance des caractéristiques du flux de fluide sur les ouvertures des fractures et présente une description anisotrope de la propagation du flux dans la fracture. Un élément fini d'interface et un modèle constitutif sont développés. Dans le but de valider l'applicabilité du modèle à des problèmes généraux de géomécanique, des simulations 2D ont été réalisées. Les résultats montrent la bonne application du modèle mathématique et le rôle hydraulique fondamental joué par les fractures.

Jan 1, 2003

By M. Z. Emad, M. U. Khan

"Coal mines opting for longwall mining or its variations rely on the main and tail gates for production. Stability of these entries or crosscuts is crucial in coal mining operations. Roof or strata failure can be fatal, may lead to property loss and delay in production. Failure of the gates depends upon many factors including the mining method, panel dimensions, pillar width, mine depth, caving rate, dimensions of gates and roof strength properties. In Pakistan coal is mined through a variation of longwall mining method with almost no mechanization while using timber as a primary support. Timber support as chock is very useful in coal mines especially in overhanging face. Similarly gate entries and coal faces are supported by chock supports throughout their length. In this paper the stability of the longwall gate is assessed for a typical coal mine in Pakistan by using finite element code. Numerical modelling considerations for gate entries are discussed. A model parametric study is performed for key factors. The results showed that the existing mining practice is not suitable and there is a need to review the mining approach. INTRODUCTION Main and tail gate entries leading to longwall mining panel play a crucial role in production cycle. Failure of any of these may lead to loss of life, property and production delays. The stability of gate entries depends primarily on the nature of loading and stresses incurred. It is relatively easier to assess the ground response to production for common longwall mining panel. In thin-seam coal mines the gate entries are generally supported throughout their length by timber made chock supports. The tailgates receive huge deformations and stresses during the mining and backfilling cycle. The chock supports are required to take care of rock falls leading to blockage of entries causing production delays and loss of capital. A number of supports are available for supporting the gate entries (Hoek and Brown 1980; Brady and Brown 1985). The support system can be selected based on load-deformation characteristics including the yielding behavior, pre- and post-yielding behavior (Barczak, 2003). The principal of support selection is that the support characteristics must match the ground response in order to maintain the gate entries (Mucho, et al., 1999). It has been observed that thin-seam mining (<1 m) is very challenging and it is necessary to develop an understanding of ground mechanics before the selection of an efficient support system Dolinar (2003). The chock supports have been tested for load-deformation characteristics in the laboratory (Barczak 2003). But the ground and support interaction, response and loads must be assessed for individual site. This can be accomplished by in-situ testing and observations at few locations. The generated data is generally insufficient for complete understanding. Numerical models can combine the data generated in the laboratory and on-site to fulfill the requirements."

Jan 1, 2017

By Shahé Shnorhokian

A study was conducted on the seismic activities and associated rockburst events that took place at an underground Canadian metal mine over a period of five years. In the first stage, a back analysis of the microseismic database led to the identification of clusters of major events. The locations of these clusters were identified with respect to the local geology and the planned mining sequence. In the second stage, a numerical model was constructed using the finite difference code FLAC3D to simulate the in-situ stress, mine geometry, local geology, and stope sequence. Mining-induced stress fields were compared against two rockburst events that took place during this period, as well as against the various geologic units present. It was shown that by combining microseismic analysis and numerical modeling, it was possible to provide a basis for evaluating the influence of the in-situ stress regimes and stope sequencing on mining-induced seismic response.

May 1, 2011

By Ehsan Khajeh

In this study, the permeabilities of Al-15.5wt%Cu and Al-19.5wt%Cu have been determined by solving the full Navier-Stokes equations on 3D geometries of interdendritic channels obtained by X-ray microtomography (XMT). The samples for XMT were sectioned from different positions in directionally solidified cylinders in order to obtain a variety of microstructures for each composition. The flow has been solved using a 2nd-order accurate Finite Volume Method (FVM) approach. The effect of employing uniform or unstructured mesh on the calculated permeability has been studied. The Marching-cube triangulation method was necessary to extract accurate surfaces for unstructured volume meshing. Calculated values of permeability for this range of solid fractions show partial agreement with previous experiments. Observed deviations are attributed to the difficulties of experimental permeability measurements.

Jan 1, 2009

By S. Finsterle, R. C. Trautz

We performed numerical modeling studies to investigate water seepage into underground openings excavated in unsaturated fractured rock. Water seepage is a key factor affecting the performance of a nuclear waste repository such as the one proposed at Yucca Mountain, Nevada. The amount of water dripping into an underground opening depends on the connectivity and permeability of the fracture network as well as on the capillarity of individual fractures intersecting the opening. We developed a high-resolution numerical model of an unsaturated fracture network and examined the appropriateness of seepage predictions using a fracture-continuum model. Model calibrations to liquid-release tests performed at Yucca Mountain are also discussed.

Jan 1, 2000

By S. Stefanizzi

In the present paper, numerical simulations based on a finite/discrete element approach are used to study deformation and failure process during standard rock mechanics. FEM/DEM modeling has been used as available with the ELFEN code, which explicitly model the transition from continuum to discontinuum. The transition is simulated through a crack nucleation and propagation that obeys to the Griffith?s failure criterion. If the fracture criterion within the intact rock (represented by FEM) is met, a new crack (represented by DEM) is initiated. Remeshing allows the fracture process through the FEM mesh to be tracked and visualized, thus contact properties can be assigned to pre-existing fractures and newly generated fractures. In this study, compressive tests of rock samples are simulated. The applicability of ELFEN to properly and efficiently simulate typical laboratory rock mechanics tests, with reference to the rock fracture process, is described in detail. In addition, an indirect validation of the extensional strain criterion proposed by Stacey (1981) is given.

May 1, 2009

By Ximeng Li

The entries that are subjected to dynamic pressure are much more difficult to maintain due to repeated disturbance. A case study of an entry stress distribution was performed using the FLAC3D program. By analyzing the distribution of stress field around the entries, this paper provides guidelines for entry support design in the lower (#16) seam subjected to the repeated disturbance from longwall mining in the immediate overlying (#17) seam, 72ft (22 m) above. The results showed that there was stress superposition inside the barrier pillar in the lower seam, 15ft (5m) to 45ft (15m) deep, during development. The influence zone of the front abutment pressure in the upper (#16) seam was about 100ft (30m) outby the faceline. The stress-relieved zone was from 0 to 100ft (30m) inby the face. The area further inby the face (i.e., 230 ft (70 m) inby the face) was within the gob compaction zone. Along the face advancing direction, the affected zones of the entries in the lower (#17) seam that were subjected to the dynamic pressure from longwall mining of the immediate overlying seam could be divided into four stages: pillar development stress stage, front abutment pressure stage, unstable gob stage, and stable gob stage. Those entries within an area from 100ft (30m) outby to 100ft (30m) to 230ft (70m) inby the faceline should be reinforced with special supports. The numerical modelling matched the monitoring results well.

Jan 1, 2013

By Young-Dong Cho, Yonghun Jong, Dong-See Shim, Seokwon Jeon, Chung-In Lee

The most important operation in the tunnel blasting procedure is to create an opening in the face in order to develop more free surface in the rock mass. This operation is highly dependent on the cut holes. In general, the burn-cut or V-cut is commonly adopted to blast the cut section. But, these traditional methods often raise problems. There are instances (a) that the actual advance is considerably less than the designed advance depending on the rock mass, (b) that they cannot sometimes enlarge the free surface efficiently and (c) that the blast vibration resulted from the cut blast exceeds the allowable vibration level. In order to resolve the problems, many different kinds of the cut blast can be considered.

Jan 1, 2005

By Khaled Morsy

Longwall mining is one of the most widely practiced underground coal mining method. The behavior of the longwall gob is very critical in the understanding of the complex ground response to longwall mining. In the limited data available, the values of gob moduli used in numerical modeling ranged from 1,000 psi to over 300,000 psi. Such a wide variation in moduli would greatly affect the results of the numerical analysis. A better understanding of the behavior of the gob will allow numerical models to be more accurate for simulating longwall mining conditions. In this study, a gob model based on the Terzaghi&apos;s model was incorporated into the comprehensive finite element package, ABAQUS, to simulate the loading behavior of longwali gob material. A case study was used to verify the proposed gob model.

Jan 1, 2002

By Riqing Ye

Laser Engineered Net-shaping (LENS TM) is an emerging manufacturing technique that ensures significant reduction of process time between initial design and final components. The fabrication of fully dense parts with appropriate properties using LENS TM process requires an in-depth understanding of the entire thermal behavior of process. In this paper, the thermal behavior during LENSTM processing of 316 stainless steel was studied both numerically and experimentally. Temperature distribution and gradient in the fabricated 316 stainless steel part were obtained by Finite Element Method (FEM) simulation. The numerical results are in good agreement with the experimental observations. The numerical method may be used to optirnize process parameters and predict the thermal response of LENSTM fabricated components. Key words: Laser Engineered Net-shaping, Thermal Behavior, Finite Element Analysis

Jan 1, 2003

By R. Karl Zipf, Francois E. Heuze

Stress analysis programs such as MULSIM/NL, LAMODEL, MinSim 2000, and EXAMINETAB are used in the mining industry to analyze stresses and displacements in coal mines, platinum mines, gold reefs, and tabular-type deposits. These relatively simple numerical models can efficiently simulate yielding and failure of a rock mass near a mine opening and subsequent stress transfer. The main input parameters for these models are a family of nonlinear stress-strain curves for the in-seam material. For many applications of these models, such as in coal mining, extensive field measurements and observations exist from which to estimate these stress-strain curves and calibrate the numerical models; however, such measurements are lacking for other types of mines. A three-step method is presented to determine nonlinear stress-strain curves for boundary-element (BE) programs used in many mining applications. The method requires a suite of laboratory-scale strength tests at various confining pressures. Note that this analysis uses laboratory-measured strength and modulus properties of alluvium. The dependency of the mechanical properties of alluvium on scale has not been established, although more tests are being planned. However, it is believed that the relative comparisons performed are valuable regardless of scale effects. First, the FLAC2D computer program is used to model the laboratory tests and determine cohesion (c) and friction angle () for a Mohr-Coulomb material model. Second, another FLAC2D model uses the c and values to calculate vertical stress and strain around a single opening in the rock mass. This model calculates the stress-strain path of points at various distances from the opening boundary. Finally, based on these stress-strain paths, the stress-strain curves needed for a BE analysis are derived. This three-step method was demonstrated in a large, flat-lying underground test facility in very weak rock. The BE analyses agreed well with observations of failure at the test facility and provided a basis for evaluating the behavior of alternative layouts.

By R. Karl Zipf

Stress analysis programs such as MULSIM/NL, LAMODEL, MinSim 2000, and EXAMINE TAB are used in the mining industry to analyze stresses and displacements in coal mines, platinum mines, gold reefs, and tabular-type deposits. These relatively simple numerical models can efficiently simulate yielding and failure of a rock mass near a mine opening and subsequent stress transfer. The main input parameters for these models are a family of nonlinear stress-strain curves for the in-seam material. For many applications of these models, such as in coal mining, extensive field measurements and observations exist from which to estimate these stress-strain curves and calibrate the numerical models; however, such measurements are lacking for other types of mines. A three-step method is presented to determine nonlinear stress¬strain curves for boundary-element (BE) programs used in many mining applications. The method requires a suite of laboratory¬scale strength tests at various confining pressures. Note that this analysis uses laboratory-measured strength and modulus properties of alluvium. The dependency of the mechanical properties of alluvium on scale has not been established, although more tests are being planned. However, it is believed that the relative comparisons performed are valuable regardless of scale effects. First, the FLAC2D computer program is used to model the laboratory tests and determine cohesion (c) and friction angle (0) for a Mohr-Coulomb material model. Second, another FLAC2D model uses the c and 0 values to calculate vertical stress and strain around a single opening in the rock mass. This model calculates the stress-strain path of points at various distances from the opening boundary. Finally, based on these stress-strain paths, the stress-strain curves needed for a BE analysis are derived. This three-step method was demonstrated in a large, flat-lying underground test facility in very weak rock. The BE analyses agreed well with observations of failure at the test facility and provided a basis for evaluating the behavior of alternative layouts.

Jan 1, 2003

By R. K. Zipf

A method is presented for creating realistic numerical models for practical coal mine ground control. The method includes procedures to collect the necessary mechanical input parameters from a geologic core log, procedures to set up a model and procedures to interpret calculation results. The input parameters come from a detailed geologic core log and extensive point load tests of estimate rock layer strength. A suite of material property input parameters is proposed which allow the user to go from core log to numerical model inputs. Rock bolt anchorage properties are also linked to the material properties of each geologic layer in the model. Following this procedure leads to very realistic calculations of the rock failure process and rock support system behavior. These calculations in turn enable realistic comparison of the effectiveness of alternative rock support systems.

By S. S. SHINDE, P. C. JIKAR, N. B. DHOKEY

Mill scale is an oxide waste product of the steel mills, which contains around 70% iron and other allied impurities. In the present context, a dynamic prototype counter current reactor (CCR) has been used for experimentation and simulation of particle bed by using discrete element method (DEM). Mill scale was pulverized to ASTM mesh range 140 (106 μm)/ + 270 (53 μm), and then it was oxidized in CCR to make single oxide Fe2O3 phase by maintaining required thermodynamic conditions such as pure oxygen at 1100 °C. Then this oxidized powder was subjected to reduction reaction in CCR having gas mixture ( H2:N2) ratio of 1:4 at 875 ± 5 °C for filling degree (fd) and reactor revolution (ω) ranging from 5.38 to 16.14% and 2 to 5 rpm, respectively. A creeping fluid flow condition (Re ˂ 1) and natural heat convection (Gr/Re2 ˃˃ 1) has been maintained while undergoing solid–gas reaction. Besides particle bed behavior has been quantified by Froude number (Fr) for the optimum operating window. The simulation was validated by the experimental results of reactor reaction product by X-ray diffraction (XRD) and scanning electron microscope (SEM) images. Prevailing heat transfer mechanism and reduction reaction mechanism has been established for solid–gas counter current conditions.

Nov 22, 2021