Search Documents

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 Adriana Amaro Diacenco

This work presents the finite element modeling of Euler-Bernoulli beams constituted of a metallic material, discussing theoretical and numerical aspects. The numerical analysis uses Hermitian interpolation functions and finite elements with two degrees of freedom per node. A numerical simulation will be done for uniform loading considering two kinds beams (steel and aluminum). The paper is organized as follows: introductory comments and numerical simulations to evaluate the static and dynamical behavior of the metallic beams.

Aug 7, 2018

By Giovanni B. Barla, Stefan H. Boshkov

The block caving method is examined in this paper on the basis of experimental results and observations in the field, and through the use of numerical modeling by the Finite Element Method. The Grace Mine, developed by a panel caving method, served as a reference case. However, the results obtained and the simulation techniques which have been developed may find useful application in more general conditions as encountered in other mining operations by block caving. Following a brief description of the Grace Mine and of the mining method, the problems associated with the behavior of underground drifts at the production levels, during undercutting and caving, are addressed. The most important results obtained through monitoring of strains and loads in the supports of a runway are discussed, while giving relevance to the sequence of mining operations in the immediate vicinity. The most important geomechanical data, as provided for subsequent modeling, are described with reference to: (a) rock mass characterization and (b) in situ state of stress. Then, the numerical modeling procedures, developed by the Finite Element Method in order to predict the caving phenomena and the behavior of underground openings during undercutting, are discussed. Finally, an attempt is made to compare the pre- dictions of the numerical analyses with, the measurements carried out at the instrumented site. Suggestions as to how to improve the results both qualitatively and quantitatively are also given.

Jan 1, 1984

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

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 Y. Fautrelle, F. Felten, O. Metais, Y. Du Terrail

"We deal with the prediction of electromagnetically-driven turbulent flows by means of a Large-Eddy- Simulation method (LES). The model is applied in the case of a liquid metal pool submitted to a polyphase linear electromagnetic stirrer. We investigate two cases: (i) the frequency of the magnetic field is equal to 50 Hz; in that case the effects of the pulsating part of the Lorentz forces may be neglected; (ii) the frequency of the magnetic field is low (f= 1 Hz), then the oscillating part of the electromagnetic forces is taken into account. The LES predictions agree well with the mean velocity measurements, as does the standard k-a model. However, as for the turbulent kinetic energy predictions, there is a large discrepancy between the two models. When the oscillating part of the Lorentz forces is taken into account, the computations show that the fluid flow is sensitive to the unsteady part of the forces. We find that the mean velocity is not affected by the fluctuating component of the force. As for the turbulence parameters, the presence of the pulsating part leads to a significant reduction of the turbulent kinetic energy, whilst the turbulence length scale decreases.IntroductionElectromagnetic stirring is commonly used in various metallurgical processes. In continuous casting for example, it may improve the structure of the solidified metal products. Fluid flows generated by electromagnetic stirrers have been extensively investigated (1 ). As for numerical predictions, many works have been published on the subject. In most of the previous works, onepoint closure models have been used. The most used model is the well-known k-i:: model. Those models yield fairly good results as for the mean velocity prediction. However, it has two main drawbacks:(i) the turbulence is not predicted accurately,(ii) it is not well fitted for unstationary flows as well as transitional turbulent motions, where theReynolds number is moderate."

Jan 1, 1999

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

In order to investigate the design parameters of a future longwall in a thick coal seam a numerical model of gas emission from the goaf area was developed. A Longwall goaf was modelled as a collection of different size coal blocks. For each coal block the partial differential equation of gas flow was independently solved. In developing the flow equations, both pressure-driven and diffusion driven gas flow regimes were considered. The total gas emission rate was calculated by integrating the gas release from each individual coal block over the whole goaf area. The main output of the computer program are gas flowrate from goaf and residual content of caved coal in the goaf as a function of distance from the face. The model was applied to a longwall face driven in a thick seam where large volumes of coal in the roof were to be left behind the face. The gas contained in the coal was mainly CO, with a high desorption rate. Simulations were under- taken to quantify the effect of size of the caved coal on the flowrate from the goaf. The flow properties of coal, such as diffusivity, permeability and coal swelling due to gas adsorption, were directly measured on coal samples from the seam.

Jan 1, 1996

By P. Ruz

A melt transfer ladle was 3D modeled by means of a mathematical fluid dynamics and heat transfer program. Internal temperature profiles were calculated as a time function for a normal ladle and for a pre-heated refractory lined ladle. It was shown that in a normal ladle the molten metal cools down quickly forming accretions. By using an inner lining there are practically no accretions build up within a typical time of operation. Therefore, by applying a refractory lining to a melt transfer ladle and keeping it at a moderate to high temperature it is possible to reduce cold dope circulation in the smelter.

Jan 1, 1999

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 R. L. McNearny, K. A. Barker

A study was conducted to compare physical and numerical models of the block-caving mining method. Results from the numerical analysis closely simulated the blocky behavior of the physical model, including bridging and interlocking of the simulated orebody during draw. The draw-down patterns and the rate of draw developed during numerical modeling also closely simulated the physical models. An unexpected result from the numerical modeling was the development of a zone of downward deflection which preceded the development of the zone of material movement.

Jan 1, 1999

By R. L. McNearny

A comparison study was conducted between an actual physical model and a numerical model based on the physical model of the block cave mining method. Results from the numerical analyses closely simulated the blocky behavior of the physical models, including bridging and interlocking of the simulated orebody during draw. The draw-down patterns and the rate of draw developed during numerical modeling also closely simulated the physical models. An unexpected result from the numerical modeling was the development of a zone of downward deflection which preceded the development of the zone of material movement.

Jan 1, 1997

By Salah Badr

Longwall mine layouts with their entries, chain pillars, face support systems, advancing mining faces and compacting gobs represent a geomechanically complex mining operation. Geomechanical aspects are further complicated if mining occurs at depths exceeding 350 m, as stresses induced adjacent to the excavations exceed the unconfined strength of the coal and result in fracturing of the sidewalls of entries and chain pillars during development. Traditional methods of geomechanical analysis, whether empirical or numerical, are unable to represent such behavior satisfactorily. Achieving some meaningful results requires explicit representation of the three-dimensional boundary conditions and implementation of appropriate constitutive models for the rock mass. This paper describes a numerical modeling methodology for deep longwall mines operation using the FLAC3D finite difference de, which incorporates a strain softening constitutive law for modeling failure of coal. An important aspect of the modeling method includes accounting for stress relief provided to the pillars and entries due to the increasing load applied to the compacting gob.

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 J S. Kuszmaul

The economic development of modified in situ oil shale retorting will benefit from the ability to design a blasting scheme that creates a rubble bed of uniform permeability. Preparing such a design depends upon successfully predicting how a given explosive will break rock. Numerical models are used to make predictions of the extent of damage caused by an explosion.

Jan 1, 1985