Search Documents

By Wen-Chieh Cheng

Two bored piles with base grouting are loaded axially and the results are used to compare with numerical analyses. One compressive pile with a diameter of 2 meters is drilled through clay deposit and into bedrock to a depth of 13.7 meters. One tensile pile with a diameter of 1.5 meters is drilled through the same clay deposit and into bedrock to a depth of 17 meters. Both piles are installed with monitoring system such as rebar strain meters, concrete strain gauges and extensometers. The monitored data during pile load test are used to obtain average frictional stress (t) a displacement (z) curve, pile tip bearing (q) as well as pile head load (p) and a displacement (s) curve. FLAC3D is used to simulate the pile loading and base grouting. The soil and bedrock are modeled with Mohr-Coulomb. The parameters of the constitutive law are established from soil biaxial tests and rock bore hole deformation tests. The pile is assumed as an elastic body. The interface element is used to simulate the pile/soil interface and pile/rock interface. The numerical results are compared with the pile load test in terms of t-z curve, q-z curve and p-s curve. Finally, the pile base grouting is simulated and the uplift is analyzed. The result from this study is promising and the numerical modeling is successful.

Jan 1, 2006

Roof falls in coal mines typically occur in or around intersections and this has not changed in the last two decades. This research develops an improved scientific understanding of stress distribution and associated failure behavior around 4-way and 3-way coal mine intersections. Three-dimensional numerical analyses were performed to determine factors that influence intersection stability. The analyses used Hoek-Brown failure criterion. Failure zones around 3-way and 4-way intersections were developed. Intersection span and horizontal stress have a major influence on intersection stability. For the 4-way intersection, pillar corners across the intersection fail first and lead to progressive failure of immediate roof and floor layers. The mechanism of failure is similar for the 3-way entry but the shape and extension of failed zones differ. Coal ribs mostly fail due to tensile stress, while roof and floor strata fail due to shear stresses. Rib corners fail due to a combination of shear and tensile stresses. In addition to stress distribution; safety factor contours analysis was performed to assess stability of intersections.

Feb 27, 2013

Two miners were fatally injured when a pillar bump occurred during retreat mining in a southern West Virginia coal mine. The mine was operating in the Eagle seam with overmining in the No. 2 Gas and Powellton seams. A coal bump is defined as a sudden and violent failure of coal caused by the release of stored strain energy in the pillar. While significant strides have been made by academia, industry, and regulatory agencies to better understand bump conditions and mitigation techniques, coal bumps represent a long standing, highly site-specific engineering problem in which the exact failure mechanism(s) is not clearly understood.

Jan 7, 2020

By Yazen Khasawneh, Yao Zhang

Upgrading and maintenance programs were proposed at a power generating plant in Indiana. To perform construction activities associated with the proposed upgrading and maintenance, a LR11350 crawler crane with a dead load of 1,500 ton was required in the vicinity of existing plant structures. The contact pressure from the crane tracks for the heaviest lift, maximum reach and angle (worst loading condition) were estimated by the rigging engineer to be on the order of 84 psi. Several concerns were raised regarding the total and differential settlements while the crane is in operation and the induced stresses on existing underground utilities due to crane loads. Three-dimensional finite element simulations were conducted using hardening soil model for the subsurface soils and elastic model for the other structural elements (underground utilities and crane mats) to optimize the configuration of the proposed crane foundation. It was found 1-foot thick engineered steel mats underlain by a 1.5-foot thick crushed stone pad will provide enough support to the crane. It will also control the settlement to the allowable limit and maintain the stresses in the underground utilities within tolerance. INTRODUCTION A LR11350 heavy lift crawler crane was proposed as part of an upgrade program at an existing power plant in the state of Indiana. The proposed crane will be sitting on an open ground in the vicinity of existing plant facilities, with multiple underground utilities going beneath the proposed footprint of the new crane pad. The existing utilities are of diverse kinds, including sanitary sewer corrugated metal pipes, electrical conduits and variable steel water/oil drain pipes. The safety and functionality of these underground utilities are critical to maintain the proper operation of the power plant. Therefore, the design of the new crane pad foundation system is of essential to control the ground settlement and the imposed stresses to the existing utilities (Davie et al., 1998, Scott 2017). To ensure all utility locations and depths are identified, the proposed footprint was thoroughly surveyed. Information in regard to each identified utility was provided by the owner. Subsurface explorations were conducted using geotechnical soil borings and soundings to evaluate the ground conditions and generate soil profiles. Numerical simulation and analysis utilizing finite element method were performed to design the crane pad foundation system and evaluate the responses of the utilities. This paper summarizes project information, geotechnical conditions and then focuses on the finite element analysis in assisting the design of crane pad and the utilities evaluation.

Jan 1, 2019

By Eugie Kabwe

Rock bolts are components of tunnel support structures and their efficacy in reinforcement requires the appropriate estimation of the bolt pre-tension. However, cases of ineffective rock bolt application due to loosening are common which leads to the potential failure of the support structures, thus negatively affecting the tunnel stability. The effective employment of pre-tensioned bolt within support structures requires an appropriate installation tension as a result of sufficient bolt head/nut pre-load. Henceforth, this paper presents a load indicator washer (LIW) capable of withstanding assigned static loading as well as provides an appropriate pre-tension. The different LIW breadths are numerically analysed in the finite element code (ABAQUS) to estimate deformation when a designated 60-kN pre-load is applied. The numerical simulation employs a bi-linear isotropic hardening constitutive law governed by the von Mises yield criterion. The simulation also considers the nominal temperature of the nutLIW-rock bolt connection. Therefore, a long-term relaxation effect on bolt pre-tension as a result of stress relaxation is neglected. The compressive force required to achieve the LIW-predetermined deformation is the bolt pre-tension. It is drawn from the simulation results that the 1.5-, 2.0-, 3.0-, 3.5-, and 4.0-mm-thick LIWs are completely flattened when subjected to 28-, 42.5-, 65-, 76-, and 88-kN pre-load respectively.

By Eugie Kabwe

Rock bolts are components of tunnel support structures and their efficacy in reinforcement requires the appropriate estimation of the bolt pre-tension. However, cases of ineffective rock bolt application due to loosening are common which leads to the potential failure of the support structures, thus negatively affecting the tunnel stability. The effective employment of pre-tensioned bolt within support structures requires an appropriate installation tension as a result of sufficient bolt head/nut pre-load. Henceforth, this paper presents a load indicator washer (LIW) capable of withstanding assigned static loading as well as provides an appropriate pre-tension. The different LIW breadths are numerically analysed in the finite element code (ABAQUS) to estimate deformation when a designated 60-kN pre-load is applied. The numerical simulation employs a bi-linear isotropic hardening constitutive law governed by the von Mises yield criterion. The simulation also considers the nominal temperature of the nutLIW-rock bolt connection. Therefore, a long-term relaxation effect on bolt pre-tension as a result of stress relaxation is neglected. The compressive force required to achieve the LIW-predetermined deformation is the bolt pre-tension. It is drawn from the simulation results that the 1.5-, 2.0-, 3.0-, 3.5-, and 4.0-mm-thick LIWs are completely flattened when subjected to 28-, 42.5-, 65-, 76-, and 88-kN pre-load respectively.

By J. Stacho, Frankovská. J.

"Abstract The use of continuous flight auger (CFA) piles in Slovakia has increased significantly in recent decades. High-rise building of administrative centre Lakeside Park in Bratislava is based on piled-raft foundation. The technology of CFA has been used due to the advantages of this technology in comparison with driven and bored piles in local geotechnical conditions. Piles have been designed in accordance with Eurocode 7. Resistance of piles was calculated by using an analytical calculation model for bored piles according to Slovak technical standard and CFA technology has been taken into account by coefficient equal to 1.1. Settlement of piles was calculated by using a calculation model which is being used in Czech Republic and Slovakia. Two well instrumented static load tests were carried out to verify the pile design. The results of static load tests showed that the design was conservative. The paper presents analysis of numerical modelling of CFA pile and static load tests. Load-settlement curves and load distribution curves over pile shafts are compared. Software Plaxis 2011 has been used for numerical modelling. Hardening soil model allows considering over consolidation of Tertiary Silts and water pressure in deep sandy lenses. The article also presents the parametric study where the influences of pile-soil interface reduction factor, over consolidation ratio and material models of soils are analysed. The results of numerical modelling showed that the most appropriate results have been obtained when the interface reduction factor equal 0.9 and real behaviour of groundwater were considered. The conclusion of the paper confirms that mechanical behaviour of the ground can be better understood and determined by using FEM than analytical methods.1. IntroductionPiles are among the most frequently used elements of deep foundations. Piles transfer loads down to a suitable bearing stratum when the soil immediately below a construction is unsuitable. Piles are also used to reduce settlement of structure. At present, the requirements of limit settlement of structures are more stringent. Examples of these structures are commercial and administrative centres that form a complex of different tall buildings. High-rise buildings cause greater contact stresses than connected low-rise buildings. Pile foundations in such cases should be placed under the high-rise buildings to achieve an appropriate settlement.Resistance of piles can be determined from static load test results, in-situ tests and empirical experiences or calculations. Static load tests give the most reliable results, but these tests are also the most expensive. This is the reason for constant development of computational methods. Especially numerical methods give more accurate results by using more sophisticated material models e. g. hardening soil or hypoplastic model."

Jan 1, 2014

By P. C. Sui

Numerical investigation and parametric simulations of the coupled electrical conduction and mass diffusion in the GDL of a PEMFC are reported. The coupled transport is solved using a 2D, "under-the-rib" computational domain for the cathode GDL. The current density on the GDL/catalyst layer interface, which constitutes one of the boundary conditions for the GDL domain and reflects the activation overpotential in the catalyst layer and the ohmic loss in the membrane, is solved iteratively using a novel numerical algorithm. A parametric study is performed to investigate the effects on current density distribution of various operating conditions such as oxygen concentration and membrane resistance, and of design factors such as GDL geometry, anisotropic transport properties, and deformation under the land area. The results show that the current density distribution under the land area can be dominated by either electron transport or mass transport, depending on the operating regime.

Jan 1, 2005

By Baokuan Li, Fang Wang, Meilong Shan

"As the advantage of uniform bus loading, the electroslag remelting (ESR) furnaces with three - phase installation by delta or star connection are mainly used to produce the large or super large ingot. The electrical efficiency is significantly improved as the short circuiting through the molten slag. However, little is found about fiindamental research on the electromagnetic field of the ESR system with three -phase installation. A three-dimensional (3D) finite element model was developed to simulate the current density, magnetic field, electromagnetic force and Joule heating for a system of electrode, slag and ingot in the ESR system with three -phase installation.IntroductionThe metallurgical performance of the ingot is the fundament and precondition to obtain the large forging of high quality, which has a critical role in the producing of forging. The electroslag process (ESR) to product the metal materials of high homogenization and clarity become more and more popular[1].The electroslag process is a special smelting technology of producing small ingots, which originate in 1950[2].At present the three, four or six electrode is used in large scale electroslag furnace. It has been found that the deeper molten pool and the wider of the two-phase zone with only one large diameter of electrode to product large ingot effect the solidification microstructure of ingot[3]. However, the molten pool becomes shallow relatively with several electrodes producing the large ingot, which is benefit to the solidification microstructure of ingot. Three electrodes electroslag remelting furnace take the three phase equilibrium power supply to making the load and electric power network balance, avoiding the distortion phenomenon induced by single phase power supply and saving electric power. It also can uniformize the distribution of heat source."

Jan 1, 2012

The N3500 orebody is the northern continuation of the 3500 orebody at Xstrata Copper's Enterprise Mine and occurs at a depth of 1.3 to 1.8 km. At the nearby S3500 orebody, extreme deformation was observed in access drives where there was a high number of bedding planes per metre combined with either graphite or talc infill. The observed failure mechanism was largely controlled by the geometry of the inclined bedding submitted to high horizontal stresses. Due to the anisotropic strength and deformation behaviour of the shale rock mass, traditional isotropic methods of geomechanical property estimation based upon GSI are not appropriate. A numerical modelling approach, based upon the distinct element codes UDEC (2D) and 3DEC (3D), has been developed to accurately simulate the failure mechanisms and deformation observed in the initial infrastructure surrounding the N3500 orebody.

Nov 25, 2010

By Hwa-Soo Park, J. K. Yoon, Hoseok Nam

In a thin slab casting process, control of fluid flow in mold is particularly difficult due to the high casting speed and large aspect ratio of the mold and the understanding of the fluid flow motion as well as heat transfer in mold is regarded as the most important information because the melt delivery is directly related to the problems of shell growth nonuniformity, surface turbulence and mold powder entrapment. A mathematical model has been developed for the coupled analysis of fluid flow, heat transfer and solidification in funnel type mold using finite volume method based on a body fitted coordinate and the characteristics of the transport phenomena in the mold of thin slab caster were analyzed by numerical simulation. As a result of simulation, the basic flow pattern is characterized as four large recirculations and two small eddies near the narrow faces of mold, which is somewhat different from that of parallel type mold or other conventional slab caster. The heat transfer model accurately predicted the mold positions susceptible to cracking as a result of thermal stress and the predicted shell thickness shows good agreements with the measured shell thickness on breakout products. We also investigated the influences of some operation parameters such as casting speed, nozzle type, submerging depth and mold material on fluid flow, heat transfer and distribution of solidifying shell and various information for the determination of optimal operating conditions were obtained.

Jan 1, 2000

By Khamis Y. Haramy

Quality assurance and quality control of drilled shafts has become a concern due to difficulties in accurately locating construction anomalies and determining load bearing capacity of defected drilled shafts. Various non-destructive evaluation techniques have been developed to estimate the integrity of the concrete in drilled shafts. While these techniques have been widely accepted, variables and unknowns can affect the measurement results. Results are typically difficult to interpret, leading to unnecessary construction delays and possible litigation over shaft integrity. In addition, influences of surrounding ground materials, stress states under different load conditions, and crack and residual stress development during concrete curing further complicate determination of shaft performance. This paper focuses on the evaluation of load bearing capacity of defected drilled shafts under various conditions by numerical analysis and modeling.

Jan 1, 2007

By V. A. Shea-Albin

Sandstone escarpments surrounding western coal-bearing plateaus serve as critical habitat for animal species. Effects of mining subsidence on escarpment stability are not well documented; therefore, regulatory agencies may limit extraction of coal reserves underlying escarpments. Because millions of tons of this coal exist on Federal Lands, associated public revenues from Federal royalties, and valuable coal reserves, risk being lost. The U. S. Bureau of Mines is using numerical modeling techniques to analyze effects of subsidence on escarpment stability. Finite element models that incorporate bedding plane slip, gob formation, and movement along joints simulate escarpment subsidence observed in the field.

Jan 1, 1994

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 Alexander B. Reeb, Richard J. Varuso, George M. Filz, Jehu B. Johnson

"T-walls, which are concrete floodwalls in an inverted “T” shape, are an integral component of the Greater New Orleans Hurricane and Storm Damage Risk Reduction System (HSDRRS). In order to resist large lateral flood loads, T-walls are typically supported on long batter piles that extend through the weak and compressible clays and peats commonly encountered in southern Louisiana. Since the piles are battered, settlement of the compressible foundation soils can result in loads acting normal to the piles and induce pile bending moments. Therefore, the design of pile-supported T-walls is a complex soil-structure interaction problem that requires detailed and careful numerical analysis to accurately estimate the stresses and bending moments in the piles. For this paper, numerical analyses were performed for three cross sections that each included a two-pile T wall system that was installed within a 15-ft tall new levee embankment. The cross sections were identical except that the elevation of the T-wall base was varied. This paper discusses the analyses, and it compares the influence of settlements and flood loading on the pile stresses and bending moments. In addition, this paper includes general recommendations for the analysis and design of pile-supported T-walls in practice.INTRODUCTION AND BACKGROUNDThe recently completed Greater New Orleans Hurricane and Storm Damage Risk Reduction System (HSDRRS) includes 350 miles of levees and floodwalls and is designed to defend against a 100-year storm surge event. The floodwalls predominantly consist of I-walls and T-walls. I-walls are comprised of a sheet pile wall curtain capped with a monolithic concrete wall above grade. T-walls are similar to I-walls, but they include a wide concrete base at grade, which forms an inverted “T” shape, and they are generally supported by long batter piles that provide additional lateral and rotational stability for resisting flood loads.During Hurricane Katrina, several I-walls failed, resulting in large breaches and devastating flooding (Duncan et al. 2008, IPET 2009). T-walls on the other hand performed quite well, including T-walls that were overtopped and eroded on the protected side, yet still did not collapse. As part of the HSDRSS, more conservative designs and height restrictions were required for I-walls (USACE 2011), and a number of I-wall sections were replaced with T-walls. However, the design of pile-supported T-walls is complex and requires careful consideration of both flood loading and settlements."

Jan 1, 2017

By Zhou Haizuo, Diao Yu, Zheng Gang

"The conventional design method for the embankment supported by DM columns usually assumes the simultaneous failure of all the columns. However, the progressive failure may occur in practice, especially considering that many researches have justified the critical load of DM columns is bending moment. In this study, a finite difference method was adopted to simulate the embankment on soft clay improved by DM columns. An increasing surcharge was applied on the embankment surface to investigate the mechanism of the progressive failure of columns beneath embankment. To capture the brittle bending failure of the DM columns, the model employed a cut-off method, in which the column was separated on failed section as long as the ultimate bending failure was reached. Furthermore, high-strength DM columns were used to enhance the global stability of embankment in this study. The numerical results showed that the ultimate bending failure was first reached at the column beneath embankment toe, and then a continuous collapse extended to the column beneath the embankment shoulder. It is found that the reinforcement could be obviously effective when the reinforced area of high-strength DM columns was extended to a certain area.INTRODUCTIONDeep mixed (DM) columns have been commonly used as a soil-improvement technique to increase stability and to reduce settlement of embankment. Research on failure modes of DM columns owing to embankment fill and surcharge has been the subject of interest for a number of investigators (Miyake 1991; Kitazume 2000; Han 2005; 2006).The potential failure modes of DM column include bending, inclination, shearing, rotating, soil flowing and a combination of these modes (Kitazume and Terashi 1991; Broms 1999). Many researchers have found that bending failure of DM columns was more critical than the shear failure. A centrifuge test was performed by Kitazume (2000) and the results showed that the columns failed in bending before shear failure. In addition, several numerical computations illustrated that design method using the shear strength overestimated the stability of embankment over DM columns (Han and Huang 2005; 2006; Navin 2006; Liu 2007; Wang 1996). In practice, an equivalent shear strength method based on bending failure could better reflect the stability of embankment (Zheng 2010; China Code JGJ79-2011)."

Jan 1, 2015

By Dao Hong Quang

This paper describes the numerical analysis of some geotechnical factors influencing the application of Longwall Top Coal Caving (LTCC) method. The study was aimed at comparing the differences in stress distribution in the chain pillars of the face ends and caving mechanism of the top coal in LTCC method between the flat and inclined thick coal seams. The outcomes from this analysis will benefit mine planners in designing the location of roadway development in the next longwall panel. In addition, results from the study will extend the understanding of caving mechanism in top coal section when applying the LTCC method (with the direction of the longwall face is down dip, retreating along strike) in inclined coal seams. The results demonstrate a clear effect of dipping factor on cavability of the top coal and on the stress distribution in the chain pillar of the face ends.

Jan 1, 2008

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 Qingwen Shi, Brijes Mishra

The abstract discusses the use of numerical simulations to investigate roof fall in underground mines. It highlights the advantages of the discrete element method (DEM) in simulating large strain problems and fracturing processes, while acknowledging its computational limitations for large-scale problems. The study presents a 3D numerical model coupling the finite difference method (FDM) and DEM to simulate laminated roof failure in a coal mine. The model represents the laminated roof with bonded particles and discontinuities, while the surrounding rock is represented by continuum zones. The research analyzes the fracturing of the laminated roof as the entry advances, demonstrating the coupled method's ability to reproduce laminated roof failure.

Mar 17, 2022

By B. Mishra, Q. Shi

Numerical simulations have been used to investigate roof fall in the past decades. Comparing to the continuum based methods, the discrete element method has been well recognized for its excellent performance in simulating large strain problems and fracturing processes in the rock mass. However, the discrete element method has its limitations. When simulating a large-scale problem, a large number of particles are required, which consumes significant computational resources. In the present study, a 3D numerical model coupled with the finite-difference method (FDM) and the discreteelement method (DEM) was created to simulate the laminated roof failure in an eastern coal mine in the US. In the FDM-DEM coupled model, the laminated roof of the entry is represented by an assembly of bonded particles using PFC3D. The laminated roof is simulated by adding parallel discontinuities. The discontinuities are given certain stiffness and strength in the model. The ribs, floor, and far-field surrounding rock are represented by continuum zones using FLAC3D. The strength of the discontinuities is varied sequentially and the fracturing of the laminated roof was analyzed with respect to entry advance. The results showed that the laminated roof fractures along the discontinuities while an intact roof fractures in a dome shape. Increasing the strength of discontinuity could enlarge the tensile zone inside the laminated roof and hence change the crack distribution. In addition, the crack quantity in the laminated roof behind the advancing face decreases with the increment of discontinuity strength while this trend is not occurring ahead of the advancing face. The simulation demonstrated that coupled method using PFC3D/FLAC3D can reproduce the laminated roof failure in underground coal mines.

Mar 2, 2022