Search Documents

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 B. Mishra

This paper presents a rock mechanics investigation of the time-dependent behavior of the rock immediately surrounding the excavation in a hypothetical one-entry coal mine based on the geology of the Illinois coalfields Research in the area of time-dependency is almost nonexistent and, therefore, this research was performed to make inroads into this neglected research area. Numerical analysis of a typical one-entry coal mine in an Illinois basin was performed using the commercial finite difference program FLAC3D. The immediate roof of the coal entry was assumed to be a shale formation. Shale in other formations has shown a tendency to creep with time. Therefore, in this study the authors modeled the shale in the roof using time-dependent characteristics. The input properties were taken from laboratory creep tests performed on shale. Model runs were performed at different stages, mimicking the progressive excavation of the entry. The model results showed that stresses in the roof attained the limit of yielding. Using Mohr-Coulomb failure criteria, potential zones of failure were also predicted. The results showed that the current computational capacities are sufficient to explore this previously uncharted territory of ground control.

Mar 1, 2014

By Tian Tang, Sergio D. Felicelli

"The laser cladding is an important surface treatment for enhancing resistance against corrosion, abrasion, and wear. The focus of the present study is to develop a three dimensional finite element model to analyze the temperature history and thermal stress evolution during multi-pass laser cladding process on H13-H13 tool steel sample. A sequentially coupled thermal-stress analysis methodology was adopted in this study. The thermal model was firstly established to investigate the temperature history of the laser cladding process. Then the temperature history was used as input to mechanical model to calculate the thermal stress evolution during the process and the residual stress after the sample completely cooled down to room temperature. The influences of process parameters including the scanning speed and scanning strategies on temperature history, thermal stress evolution, and residual stress were investigated in details.IntroductionLaser cladding is an advanced material process technique, which is often used to improve mechanical properties or increase corrosion resistance, repair worn out parts, and fabricate metal matrix composites. In this process, the injected powder or wire feedstock material and a thin layer of substrate are melted and consolidated by use of a laser in order to coat part of a substrate or fabricate a near-net shape part.Numerous mathematical models have been established for analyzing the effects of process parameters, microstructure evolution, and thermo-mechanical behavior of laser cladding. Zhou et al. [1] established an analytical model to investigate the laser induction hybrid rapid cladding for Ni-based WC composite coatings. Huang et al. [2] also developed an analytical model to theoretically investigate the interaction between the laser beam and the powder stream of laser cladding process. A simple and realistic model of laser cladding was proposed by Picasso et al. [3]. The cladding layer is very sensitive to cracks due to the accumulation of residual stress, caused by non-uniform temperature changes, which have limited the wide application of laser cladding technique in industry. Recently, fmite element method has been increasingly used to simulate the laser cladding process [ 4-6]."

Jan 1, 2014

By R. Chen

Most tunnels built in Sydney encounter similar ground situations; they penetrate either the Hawkesbury Sandstone or the shales that lie beneath the city. Tunnels in these materials are generally sup-ported by rock bolts in combination with a shotcrete lining. A method of analysis is presented to improve our understanding of the long term behaviour of fiber reinforced shotcrete (FRS) materials and the tunnel support systems. A finite element method SAFEA (Semi-analytical FE Analysis) is developed to model the combined creep behaviour of a shotcrete lining and the Hawkesbury Sandstone. To model the time-dependent creep behaviour, the bulk modulus and shear modulus of the materials are represented by logarithmic functions in time. A Laplace transform is then applied to these functions of time and the solution of the problem found from a finite element analysis. Tablot?s method is then used to invert the transformed values and find the solution in real time. The time dependent creep deflection of the numerical tunnel models reflects the interaction of the shotcrete lining and rock which have different creep parameters.

Jan 1, 2008

By Jun-ichiro Yagi

A mathematical model has been developed for simulating the blast furnace operation with carbon composite agglomerates charging based on multi-fluid theory and reaction kinetics. In this model, the behaviours of carbon composite agglomerates are considered on the basis of previously reported conservation equations and experimental research on the reduction rate of carbon composite agglomerates. A series of calculations are performed to examine the effect of charging carbon composite agglomerates. The model calculation gives two-dimensional distributions of process variables and information on the overall operational parameters. It reveals that in-furnace temperature levels significantly decrease and the reduction of carbon composite agglomerates is retarded with carbon composite agglomerates charging ratio. The furnace performance is remarkably improved with the increase in carbon composite agglomerates rate. The productivity tends to increase while coke rate and total reducing agent rate show decreases at the different degree.

Jan 1, 2006

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 Go-Gi Lee, Jae-Young Jung, Myung-Duk Seo, Hyun-Na Bae, Seon-Hyo Kim

"Bubbles form at the anode during the electrowinning of magnesium and cause hydrodynamic acceleration and electrical field disturbance. A three-dimensional computational model was developed to investigate the effects of the geometry of the electrochemical chamber and the size distribution of bubbles on the fluid flow in the pilot plant. The results are validated by measurements of the surface characteristics of molten salt and the bubbles formed. Fluid flow at the surface was recorded with a high speed camera, and the distribution of bubble size is determined by inspecting the video images frame by frame. The measured size distribution of the bubbles is entered into the computer simulation of a multi-phase flow. A sloped roof is installed above the electrodes, to enhance the flow of the molten salt, thus resulting in higher current efficiency, lower energy consumption, and higher quality of magnesium metal.IntroductionTitanium's excellent properties of low density, high specific strength, high temperature strength, and exceptional resistance to corrosion are highly beneficial to diverse applications including aerospace, automotive, construction, and medical implants. Although the Kroll process[l] and Hunter process have been used for the extraction of titanium in the 20th century, their complicated engineering stages and high production cost, have led to further research into new processes such as FCC Cambridge, Armstrong and TiRO. Despite these intense research efforts, no process has yet been developed that is more economical than the Kroll process, by which most of the global commercial titanium is still being produced."

Jan 1, 2011

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 F. Guevara

The paper presents results of physical and turbulent three-dimensional numerical simulation of the gas flow in the wind box and reaction shaft of the concentrate burner of a flash furnace. Physical simulation work was carried out in a 1:5 linear scale model, where the velocity field was measured using a Pitot tube. The mathematical model consists of the three-dimensional averaged Navier-Stokes equations, turbulent kinetic energy and dissipation of turbulent kinetic energy using the Renormalization Group (RNG) theory for the k-e turbulence model. Computations were performed using FLUENT, based in a finite- volume code. The geometry of the concentrate burner was separated in two sections: wind box and reaction tower. Results obtained at the exit of the wind box were used as boundary conditions for the reaction shaft. The effect of boundary conditions on the numerical results was also studied. Results of physical and numerical simulations show an asymmetry of the velocity field at the exit of wind box. This asymmetry is possibly due to existence of a separated flow in the back of the exit ring. The agreement of physical and numerical results assess the capacity of the turbulence RNG k-e model to reproduce correctly the behaviour of the velocity field at the exit of the concentrate burner.

Jan 1, 1999

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 Peter Griesenbrock

In German underground coal mining mostly arch and some rectangular roadways are used in the development work. Due to its positive results from numerical and physical modeling techniques as well as good underground experiences, the numbers of rectangular roadways has increased significantly. Today the bolt support - independent of its cross-sectional shape - represents an undeniable support element. During the introduction of new bolt techniques a detailed investigation of the rock and support mechanical effectiveness of bolt support in roadways at all operation stages is of basic importance. This effectiveness can be predicted and verified using different simulation techniques. For the planning of rockbolted roadways two supplementary simulation techniques are presented in this paper. By combination of all available planning tools the best result for planning assurance can be achieved. Keywords: numerical modeling, physical modeling, arched roadway, rectangular roadway, roofbolting, advancing longwall.

Jan 1, 2002

By S. Bradshaw, A. Mabentsela

"Physical and numerical modelling methods were followed to identify mechanism(s) for tundish slag entrainment in a bare tundish and one with a flow control device (FCD).The physical and numerical models made use of water and paraffin to model steel and slag respectively. Observations from the physical model showed that the steel-slag interface remains immobile in both cases. Entrained paraffin formed droplets approximately 1 mm in diameter. Results from both models (numerical and physical) showed that in both cases (bare and FCD case), areas of high entrained slag concentration exist near the inlet region. The entrained slag concentration decreases towards the tundish endwalls.Flow patterns and velocities tangential to the steel-slag interface from the numerical model showed that slag entrainment in both the bare tundish and tundish with a FCD possibly takes place via two mechanisms. First, the slag moves across the steel-slag interface via mass transfer; secondly small velocities tangential to the interface at depths greater than 10 mm below the interface carry the already ‘entrained’ slag into the bulk steel phase. These tangential flow patterns are dominant in the inlet region, hence the high concentration of entrained slag in this region. IntroductionThe tundish serves as the last metallurgical vessel through which steel passes before solidifying in the moulds. It is in the tundish that the last traces of non-metallic inclusions should be removed, otherwise the inclusions carry over to the moulds and cause defects in the steel product. In the tundish, the slag is less dense than the melt and thus resides on the top of the melt. The slag provides a sink into which the non-metallic inclusions float to and dissolve. The slag also protects the melt from air and heat loss. However, when slag is entrained in the melt as a result of increased turbulence or shearing at the steel-slag interface, slag can become a source of nonmetallic inclusions.Numerous studies on flow patterns in the tundish have been carried out since the introduction of the tundish. These studies were aimed at improving the flow characteristics in the tundish (Chattopadhyay, Isac, and Guthri, 2010; Cloete, 2014; Jha, Dash, and Kumar, 2001; Jha, Rao, and Dewan, 2008;"

Jan 1, 2017

By R. J. Shaffer

Hydraulic fracturing in rock masses involves complex and coupled processes of fractures propagating in discontinuous media and of fluid flow in discrete channels. The Unconventional Gas Program at LLNL is actively developing numerical models to gain insight in the physics of these coupled processes. This paper describes the current status of the FEFFLAP code, which models fluid-driven discrete fracture propagation in jointed media. We also discuss the results of a preliminary series of physical tests in which hydrofractures were driven across slanted interfaces between dissimilar materials, in blocks loaded in biaxial compression.

Jan 1, 1984

By KeerthiRaaj Subrmanian, Muthukkumaran Kasinathan

"A railway embankment is under construction between Nagapattinam to Thiruthuraipoondi covering 36.5 km. The sub-soil strata comprises of soft clay extending up to 10 m below existing ground level, which is subjected to stability and settlement problems. A trial stretch covering about 5 kilometer was improved using Prefabricated Vertical Drains (PVD) with surcharge was recommended to accelerate the consolidation settlement. The geotechnical characteristics of the subsoil are examined through laboratory tests on representative soil samples. Recent improvement in the finite element analysis has facilitated modeling of PVD in finite element software for computing the settlement magnitude and time rate of consolidation. In this study the theoretically predicted settlement magnitude and time rate of consolidation from Terzaghi.s vertical consolidation, Barron's radial consolidation are compared with the analyzed FEM results for both the cases with and without PVD improved embankment. The comparative analysis between FEM and predicted data hold good for stresses and the associated deformations. The excess pore water pressure was significantly lower and their dissipation was faster in case of embankment with vertical drains. Thus PLAXIS 2D tool shows fairly convincing results in comparison to the theoretical solutions.1. INTRODUCTIONFrom the earlier times civilization has been developed in a fashion either at the banks of river or near to the shore region, because of the easier availability of food and water. This over centuries made them as huge cities, consequently population rise made the increased necessity of highway, railways, port and even more ground for construction. As a stake researchers are pushed to make a low bearing ground into a massive carrier of skyscrapers, embankment, bridges etc. These low bearing ground suffers from total and differential settlement, slope instability. Hence, numerous ground improvement techniques were followed to eradicate these concerns which include preloading, piling, stone column, vertical drains and geosynthetic reinforcement. Among them Prefabricated vertical drain (PVD) are often used to accelerate consolidation of soft subsoil and has been the main soil improvement method used in improving the soft marine clay deposit (Chu et al.(2004); Chai (2001) and Arulrajah et al.(2003)). Literature review of embankment over soft clay deposit reveals that combination of PVD with preloading is widely considered as an economical and less time consuming ground improvement technique."

Jan 1, 2015

By T. H. Yang, K. K. Lee, Z. Z. Liang, C. A. Tang, L. C. Li, Y. Tsui, L. G. Tham

Abstract: A series of numerical tests of hydraulic fracturing in heterogeneous rocks were performed to investigate the influence of pore pressure magnitude and gradient on initiation and propagation of tensile fractures. The numerical tests were run with a Rock Failure Process Analysis code, F-RFPA2D, incorporated with a Flow-Stress-Damage coupling model (FSD). To investigate the influence of pore pressure magnitude on fracture behaviour, “pinch-off” breaking tests, originally employed by Bridgeman to investigate the effective stress law for tensile fracture in a uniform pore pressure field, were numerically modelled. In another numerical test, a double-notched sample, with fluid pressure in one notch while keeping the other one open to the atmosphere, is numerically extended to investigate how the fluid flow direction or the pore pressure gradient will influence the fracture behaviour. The simulation results show that fracture is strongly influenced by both pore pressure magnitude and pore pressure gradient. Keywords: seepage mechanics; heterogeneous; pore pressure; numerical simulation; failure process ZUSAMMENFASSUNG: Um den Einfluss der Wasserdruckgroesse der Oeffnung und der Porisitaet auf den Bruchentstanden und Verbreiten der Tension von Gebirgen zu beurteilen, wurde viele nummerische Pruefungen auf den hydraulischen Bruch in dem Gebirgen aufgenommen. Diese Pruefungen wurden mittels der Analyse- Software F-RFPA2D und Flow-stress-Damage Coupling model ( FSD) durchgefuehrt. Um den Einfluess der Wasserdruckgroesse der Oeffnung auf den Bruch zu beurteilen, wurden auch Bruchpruefung und nummerische Analogexperimenten gemacht, und dadurch wurde die ganze Phase der Bruchentstanden und Verbreiten wiedergezeigt, und die Ergebnissen ergeben sich, dass die Wasserdruckgroesse der Oeffnung und der Porisitaet sehr viel auf den Bruchentstanden und Verbreiten der Tension von Gebirgen beeinfluessen. Schluesselwoerter: Durchlaessigkeit , ungeleichmaessig, Wasserdruck der Oeffnung, nummerische Analog, Zerstoerungspahse Résumé : Une série de simulations numériques a été ménée pour analyser l’influence de la grandeur de la pression de pore et du gradient sur la naissance et la propagation des fractures de tensions par le biais de la simulation numérique des fractures hydroliques dans des rochers hétérogènes . Le logiciel F-RFPA (Rock Failure Process Analysis code ) intégrant FSD ( Flow-Stress-Damage coupling model ) y a été employé. Les résultats des simultations montrent tout le processus de la naissance et de la propagation des fractures sous la pression de pore. Ces résultats nous ont permis de constater que la grandeur de la pression de pore et le gradient exercent une influence importante sur la naissance des fractures dans les rochers et sur leur mode de propagation . Mots de clé : le mécanisme de transfusion ; l’hétéogénéité ; la pression de pore ; la simultation numérique ; le processus destructif

Jan 1, 2003

By S. Peng, Vaughan R. Voller

"Free boundary problems are a common occurrence in materials processing operations. An important example is the filling of material into molds. This paper looks at filling problems associated with polymer processing, in particular Resin Transfer Molding (RTM) and Compression Molding. A general, single governing equation is developed which on numerical solution provides predictions for the position of the polymer front in the mold during filling. Both fixed and deforming finite element solutions are developed. Application to polymer filling test problems reveal that the fixed and deforming grid scheme are accurate and compare closely with each other.IntroductionProblems arising from the analysis of material processing operations often involve the tracking of a moving boundary, on which a balance conditions needs to be met, or the locating of a free surface that defines the solution domain. Material processing examples of so called ""free and moving boundary problems"" include; the tracking of solidification/melting fronts[l], the tracking of the metal/gas interface during the filling of a mold[2] and location and shape of a skull (freeze line) in a melting furnace[3].Computational solution of free and moving boundary problems usually fit into one of two classes, viz.1. Fixed grid scheme: In which the spatial domain is covered by a fixed grid (mesh) of node points and the solution uses an auxiliary variable (e.g. a local liquid fraction) to track the moving boundary or locate the free surface.2. Deforming grid scheme: In which the grid in the spatial domain is continuously deformed (or adapted in a steady state problem) such that a specified set of grid points lie on the moving boundary or the free surface.Details on the basic underlying features in these two classes of solution can be found in the excellent monograph by Crank [4]."

Jan 1, 1994

By W. G. Pariseau, C. H. Schmuck, Fei Duan

The Homestake Mine is located in steeply dipping Precambrian metasediaents, an environment common to a number of world class ore bodies. Development of a pronounced plane of schistosity raises a question concerning the importance of anisotropy to the design of vertical crater retreat stopes. Little is known about the actual influence of elastic anisotropy in such -circumstances; even less is known about strength anisotropy. A cautious step- by-step approach was therefore taken towards a numerical assessment of the role of anisotropy in the analysis of stress and stability. The sequence of analyses begins with elementary test problems for computer program checks and proceeds to the analysis of a study stope area between the 6950 and 7100 Levels of the Homestake Mine. Anisotropy (orthotropy) is important at the Homestake Yine. No simple guidelines for evaluating in advance the need for plastic and strength anisotropy in design analysis appear to be forthcoming, ;o that some physical properties testing is necessary where the question is in doubt

Jan 1, 1984

A number of issues relating to the computational aspects of pillar design are addressed in this paper. The computation of average pillar stress values is important when attempting to establish criteria for pillar design and in the analysis of the stability of tabular pillar layouts. One of the default ?classic? numerical methods that are used to determine pillar stresses is the displacement discontinuity method. In many instances it is not clearly understood that this approach does suffer from some limitations, particularly in relation to the fact that in coarse element simulations, the simulated average pillar stress (APS) can depend on the chosen mesh size. The nature of this error is highlighted in this paper and some strategies are suggested to bound the magnitude of these errors. It is demonstrated as well that the popular linear stiffness approximation to pillar or seam compressibility does appear to allow reasonably accurate estimates of the average pillar stress when either the pillar height is varied or when the seam modulus differs from the host rock modulus. A practical implication of this study is that if the seam modulus is noticeably lower than that of the host rock, such as for coal seams, it is important to use a linear stiffness constitutive model for the pillars rather than a ?rigid? pillar assumption. This added complexity seems unnecessary, however, when simulating hard rock pillars in mines where the seam modulus is very similar to that of the surrounding rock.

Jan 1, 2011