Search Documents

By D. Magoon

The metals industry has enjoyed a resurgence in prices and demand for its commodities during the global surge in economic activity that began in 2005. This success has brought with it an unprecedented demand for professionals in all sectors of our economy. Nowhere is that more true than in the mining and metallurgy sectors. However, this welcome improvement follows on the heels of almost two decades of challenging markets for the mining and metallurgy. Our sector's protracted economic slump created a paucity of demand for professionals that had several consequences. Our industry came to be viewed as a low ?tech, uninteresting place to work. Low demand for professionals encouraged students to seek other career options. Reduced demand for graduates, coupled with diminished program funding by our sector, resulted in dramatic shrinkage in university mining and metallurgy programs. This paper explores the lessons to be learned from our recent economic history and discusses some strategies for addressing our image issue. Attracting and retaining good employees is also strongly influenced by a company's reputation and record as a good employer. This reality is very well illustrated by the findings of an international survey undertaken by Towers Perrin in 2005 of 86,000 employees of mid-sized and large companies in 16 countries. The key findings related to attracting and retaining employees are presented. These findings are also discussed in the context of their specific applicability to the mining and metallurgy industry.

Jan 1, 2007

By Gordon. I. M., J. Blackshaw

"The wet part of the Hollinger mill is housed in two buildings, the grinding plant and the treatment plant. In the grinding plant there are four 7xl5 ball mills which grind the ore to 60 per cent through 200 mesh. There is also one 7xl5 ball mill which grinds table concentrate to 9 5 per cent through 200 mesh.The primary classifier overflow is pumped to the treatment plant where it passes over shaking tables, producing a concentrate and a tail. The concentrate pulp is pumped to the grinding plant for secondary grinding. The secondary classifier overflow is pumped back to the treatment plant where it is thickened in two 40xl 5 tray thickeners and agitated in four 20x24 mechanical agitators arranged in series of two (with 2 spares).The table tail pulp is thickened in nine primary 40xl5 tray thickeners (with 2 spares). The thickened table tail pulp and the concentrate agitator discharge pulp are mixed and pumped to four rows of 15x45 pachuca agitators, arranged in series of five tanks.The pachuca discharge pulp is divided. The smaller part goes to two rows of 3-step decantation using 40xl5 tray thickeners (with 2 spare rows). The larger part goes to six 14xl6 drum filters. The repulped primary filter cake and the thickened pulp discharged from the decantation rows are fed to eight 14xl6 drum filters. The secondary filter cake is mixed with water and pumped to the tailing line."

Jan 1, 1968

By Carmen Mihaela Neculita

Sulfate-reducing passive bioreactors efficiently treat acid mine drainage (AMD) by increasing its pH and alkalinity and by removing metals as sulfide precipitates, in addition to hydroxides and carbonates. However, the new mining legislation sets discharge limits based on physicochemical parameters and on the non-toxicity of treated effluent as well. In this study, toxicity was assessed for effluent from 3.5 L column bioreactors filled with mixtures of natural organic carbon sources and operated at two different hydraulic retention times (HRTs) of 7.3 d and 10 d for the treatment of a highly contaminated AMD. Acute toxicity was observed for D. magna, and a toxicity identification evaluation (TIE) procedure was performed to identify potential toxicants. To help the interpretation of the toxicity results, metal partition in the effluent was also determined by filtration and ultrafiltration, while geochemical modeling was performed for metal speciation evaluation. The 10 d HRT effluent was nonacutely lethal for rainbow trout but acutely lethal for D. magna, whereas the 7.3 d HRT effluent was toxic for both organisms. The toxicity to D. magna was eliminated by 2 h of aeration for the 10 d HRT effluent, and the TIE procedure suggested iron as a cause of toxicity. However, the 2 h aeration eliminated the toxicity of the 7.3 d HRT effluent, only if performed after the pH rising to 9.3. Data regarding metal speciation indicated instability of both effluents during aeration and were consistent with the toxicity being caused by iron. The study indicates that design of passive treatment should include sufficient HRT and posttreatment aeration to meet acute toxicity requirements. Keywords: acid mine drainage, sulphate-reducing passive bioreactors, hydraulic retention time, effluent toxicity, metal speciation

Nov 1, 2008

By A Frostick

Ranger uranium mine is surrounded by the World Heritage listed Kakadu National Park and offsite impacts resulting from mining at Ranger would affect the Magela Creek catchment first. Radionuclides, lead isotope ratios and metal concentrations were measured in sediment cores from the Magela Creek catchment around Ranger to investigate the dispersion of mine-derived particulates. The study found a localised area within the Ranger mineral lease, where upper sediments had elevated levels of radionuclides and trace metals and radiogenic Pb isotope ratios. There is no current evidence of sediment dispersion from Ranger further downstream within the Magela Creek catchment.

Aug 1, 2010

By M. Pelser, H. Gous, J. H. W. M. Brits, K. Voogt

The effect of impurities; iron and selenium, on current efficiency and deposit quality in copper electrowinning was investigated using an at-line cell at the Rustenburg Base Metals Refiners’ Copper Tankhouse. The tankhouse is a direct leach, starter sheet electrowinning operation with a unique electrolyte composition. The effect of iron on the current efficiency in the at-line cell was found to be less severe than reported in literature. However, the distinction between the chemical impact of iron on current efficiency and the physical impact of iron on a starter sheet operation needs to be realised. The effect of selenium on current efficiency and deposit quality is also reported. INTRODUCTION Anglo American Platinum’s Rustenburg Base Metals Refiners’ (RBMR) Copper Tankhouse serves a unique function in the refinery flowsheet in that its primary functions are the removal of copper and production of acid for the leach circuit. This function needs to be performed to specification in order to keep the leach circuit operating optimally. Current efficiency is seen as an indicator of tankhouse performance and efficient energy use. The tankhouse plates starter sheets onto titanium blanks which, after stripping and preparation, return into the production cells for further plating. Approximately 78% of tankhouses surveyed (Robinson et al., 2013) use permanent cathode technology whereas older, smaller operations still use starter sheet technology. Drawbacks of the starter sheet technology are that the starters are thin and bend easily leading to increased shorts, they need to be removed after two days of plating to be manually straightened and are then returned to the cells. In addition, the electrolyte levels in the cells need to be adjusted during the plating cycle to improve loop adhesion to the cathode and reduce loop corrosion. The starter sheet technology is therefore very labour intensive and the final quality of the product is dependent on the quality of the starter sheets produced. In comparison, permanent cathode or direct plating technology lends itself to automation, allows for higher current densities and better cathode quality due to more robust geometric stability in the cell. The present barrier to implementation of permanent cathode plating at the RBMR Copper Tankhouse is the capital cost required to convert the tankhouse, i.e. blanks and stripping machine. At present, alternatives to titanium blanks, such as LDX 2101, are being investigated for the unique electrolyte composition of this tankhouse.

Jan 1, 2019

By R. P. Charbonnier

"To understand the operation of jigs and other separating processes based on differences in gravity, it is desirable to analyse the motion of various particles in given media under specific conditions. A summary of the earlier theories shows that they lead to unsatisfactory or excessively complicated formulae in cases such as coal jigging. By adopting a simplified theory of the motion of particles in a fluid, the size and duration of the initial accelerated motion of a particle falling in hindered settling can be calculated. Based on this mathematical analysis, examples referring to actual coal jigging practice show that it is desirable to prolong and/or repeat, as frequently as possible, the optimum conditions corresponding to an effective acceleration of the particles. The separation is then obtained according to gravity rather than to particle dimensions, as would be the case if the velocity of the particles• were allowed to approach its terminal value. To achieve this gravity separation, it is suggested that the suction phase of the jigging cycle be artificially prolonged and the frequency of pulsations increased as much as possible for washing the coals of smaller sizes."

Jan 1, 1959

By Renato Macciotta

The argillaceous formation known as the Opalinus Clay is host to the Mont Terri Underground Research Laboratory. The URL was developed to investigate the hydrogeological, geochemical and geotechnical properties of the formation. Underground excavations at the site are characterized by non-elastic and time-dependent deformations. Excavations trending perpendicular to the strike of the bedding presented relatively minor deformations, however, tunnels parallel to the strike of the bedding showed much larger non-elastic deformations. This paper presents an analysis of the measured response of a tunnel trending parallel to the strike of the bedding in Opalinus Clay. The measured convergence is compared to tunnel experience in weak rocks to assess the tunnel stability.

Jun 1, 2012

By D. A. Brown, C. Lowe

"Abstract -A multi-sensor (electromagnetics, magnetics, gamma-ray spectrometry, and VLF) geophysical survey was recently conducted in the Purcell Basin of southeastern British Columbia in an area encompassing the world-class Sullivan Sedex deposit. The purpose of the survey was to geophysically characterize the known mineral occurrences and to develop strategies to enhance future mineral exploration in the basin. Follow-up investigations of geology and measurements of rock physical properties helped constrain interpretations of the geophysical data sets.Despite extraction of more than 90% of the ore from the Sullivan deposit, it is clear that the mineralization there, as well as at other Sedex occurrences surveyed, is characterized by high conductivity values, strong bedrock conductors, and elevated magnetic anomaly values. The albite-sericite alteration which accompanies this mineralization yields high potassium and low thorium-potassium ratio values. Associated Moyie gabbro is characterized by low electrical conductivity and radioelement concentrations and by high magnetic anomaly values. The study has also aided exploration for other mineral deposit types including stratabound Cu-Ag, placer gold, vein and disseminated basemetal deposits. In addition, the geophysical data have helped constrain geological interpretations of the basin."

Jan 1, 1998

By Ahmad Ghahreman, Denver Cowan, Yeonuk Choi

"Previous testwork by the authors has indicated that both granular coconut shell based activated carbon and a new microporous carbonaceous bead, Lewatit® AF 5, are highly effective in improving the oxidation kinetics and decreasing the elemental sulphur yield when assisting the atmospheric oxidation of pyrite. While the two carbon catalysts show similar results and improvements to the pyrite oxidation, there are differences between them. For example, the activated carbon assisted oxidation showed faster kinetics than the AF 5 assisted oxidation, while elemental sulphur seemed to remain adhered to the surface of the AF 5 after oxidation when much of the elemental sulphur formed during the activated carbon assisted tests was found in the oxidation residue. The objective of this paper is to investigate the differences between the two catalysts through the use of scanning electron microscopy, photography, Raman spectroscopy, surface chemistry analysis, and Brunauer-Emmett-Teller (BET) surface area analysis. With this investigation the authors should be able to achieve a better understanding of the surface chemistry of the two carbon catalysts (i.e. functional groups, charges, acidity etc.) and help explain the strengths and/or the weaknesses of the catalysts when applied to the atmospheric oxidation of pyrite.INTRODUCTION Previous research conducted by the authors, indicated that both a granular coconut shell based activated carbon (activated carbon) and a microporous carbonaceous bead, Lewatit® AF 5 (AF 5), acting as a catalyst enhanced atmospheric pyrite oxidation in ferric sulphate media by approximately 40% when compared to unassisted oxidation and that elemental sulphur yields were increased from 23% to approximately 69 and 63% when assisted by AF 5 and activated carbon respectively (Cowan, Jahromi, & Ghahreman, 2017). The study indicated that although total pyrite oxidation assisted by the two catalysts was similar, differences were noticed in the impact that the two catalysts had on the outcome of the oxidation. Kinetic curves in the Cowan et al. study indicated that atmospheric oxidation of pyrite seemed to reach completion after 48 hours of oxidation when assisted by activated carbon. A similar trend was noticed in the kinetic curves for oxidation when assisted by AF 5, however the oxidation curve did not reach completion until approximately 96 hours (Cowan et al., 2017). This observation indicated that the kinetics of pyrite oxidation are enhanced by both catalysts, oxidation kinetics are much faster when assisted by activated carbon rather than AF 5."

Jan 1, 2017

By G. Robert Guillet

"Ontario enjoys the advantage of its proximity to the large industrial market of northeastern North America. Together with Quebec and fifteen nearby U.S. states, an industrial base of more than half the total for both countries is provided by a population of 100 million.Mineral fillers such as barite, calcium carbonate, kaolin, mica, silica and talc are of major importance in the manufacture of paper, paint, plastics, rubber, vinyl tile, carpet backing, putty, caulks, sealants and adhesives. Ontario has all of these minerals, but domestic production falls short of demand in some cases, and mica and kaolin are not produced at all.This paper reviews the market and outlines the distribution of the major filler-consuming plants of northeastern North America. It also identifies the principal producers of filler minerals in the United States and Canada.IntroductionOf all the uses for industrial minerals, mineral fillers are probably the least understood. Originally used primarily to add low-cost bulk to a variety of commodities without seriously affecting their properties, filler minerals have become essential components that enhance and extend the properties of their hosts. Hence the term ""extender"" which, for all practical purposes, is synonymous with ""filler"". Today, fillers and extenders are not solely low-cost additives, but ingredients that provide essential characteristics to the products they fill.The major markets for mineral fillers are paint, paper, plastics, rubber, adhesives and floor coverings. But there are many others: putty, pesticides, wire and cable coatings, caulking and asphaltic compounds, oilcloth, wall-paper, joint cements, inks, sealants and many others. Barite, calcium carbonate, kaolin, mica, silica and talc are among the more important filler minerals discussed in this paper, but there are at least a dozen others. And the increasing practice of coating the mineral surfaces with certain coupling agents or metal stearates is constantly expanding the range of characteristics and hence the marketing opportunities.Whiteness, particle shape, fineness of grind and cost are the principal considerations in selecting a filler, but other characteristics such as weight, opacity, hardness and chemical inertness may be dominant considerations in specific applications."

Jan 1, 1987

By C. F. Gleeson

"I HAVE BEEN ASSIGNED the pleasant task of adding a last word in way of summary to this session. I think that the trend set in having a geochemical session at the CIM Annual General Meeting is a good one and I hope it will continue so that members of the mining fraternity can get together on a regular basis and discuss geochemical problems peculiar to our Canadian scene. A series of interesting papers have been presented on what has been considered by many in the past as not only a problem area for geochemistry but even an impossible one. The papers should help dispel some of these fears. Lack of development of geochemical techniques in the Shield has been due to a variety of factors, many of which have been mentioned. However, I feel that one of the main reasons for this lack of development has been timidness on behalf of many people to get out and do the necessary work which helps produce answers to many of the questions asked. New developments and applications of geochemical methods in the Shield will not be done by the faint of heart but by those who have the curiosity, interest, imagination, drive, resources and, most of all, the intestinal fortitude to stick to it and not give up because results or conditions do not meet with preconceived notions. Another element that will be needed by those who will successfully apply the techniques and who will be responsible for new breakthroughs is the systematic approach to sampling and analysis and the integration of geochemical, geological and geophysical results. Erratic sampling patterns - a sample here and a sample there - generally lead to more problems than solutions."

Jan 1, 1971

By R. Rosenblum, J. Y. Kim, D. Shink, K. G. Stowe

"A small scale continuous flotation cell was designed and constructed at the Noranda Technology Centre to develop a mobile testing flotation circuit located next to the production circuit. This mobile small scale flotation pilot plant can easily and economically be set up at any circuit in a mill and used to evaluate and predict plant performance against a given flotation condition. The main feature of the cell, beside the provisions of automatic pulp level and feed rate control, is to provide a reduced surface to volume ratio approaching to that of an industrial machine.This paper describes the basic design and the results of various testwork carried out to date including aeration, conditioning, reagents, recirculation, flotation gas, etc. are presented in comparison with the plant production results. INTRODUCTIONStudy of the behaviour of an ore during concentration by flotation process starts, in general, with a series of tests in laboratory cells. Given the discontinuous nature of these tests, the information obtained by such a study is somewhat limited. It has long been recognized that laboratory test results are not always reproducible in pilot plant tests and subsequently, the pilot plant results do not always predict plant performance. The differences are often blamed on batch vs continuous operations, grinding media, size distribution, aging (oxidation), water quality, equipment, operator's skills, etc."

Jan 1, 1992

By Blair Hrabi

Introduction; ?Risks Associated with the Resource Estimation Process; ?Limitations of Available Data; ?Geological and Structural Control of Mineralization Trends; ?Recognizing Multiple Mineralization Controls and Complex Trends; ?Accurate Geodomaining; ?Conclusions.

May 1, 2013

By C. Bourgoin

"The limited development and application of new technologies to narrow-vein mining is currently impeding progress in the industry to a point where productivity, working conditions and operating costs have reached a critical stage in many Canadian mines.In an effort to foster the emergence and application of new technologies, Université Laval, together with other organizations, has carried out a R&D project focussed on the need of identifying new development pos5ibilities and of facilitating the transfer of technology in matters related to underground mining methods and machinery.The first phase of the study consisted of an in-depth assessment of current practices and problems in narrow-vein mining. The second phase was a state-of-the-art survey of mining methods and equipment, together with proposed modifications and innovations. The third phase of the project concerns the experimental development as well as the principal application of various proposed technologies.This paper presents an overview of the organizational, technical and strategic considerations which guided the University group and it serves as well to introduce a newly-published brochure, entitled ""Underground Mining Techniques of Vein Deposits"", which summarizes such work to date. More specifically, a perspective of technological developments applicable to narrow-vein mining is presented."

Jan 1, 1994

By Ross Dimmock

The productivity of a mine has the same goals as the quest to complete a tunnel project both safely and within time. A key aspect to achieving this goal is to excavate and support the working face as quickly as possible, allowing early re-entry to start construction of the next advance in safety. As with tunnels, the use of high performance wet-mix sprayed concrete has increased in mines dramatically over the last few years, and year on year has seen a stepwise increase in shotcrete performance, particularly in the first few hours. Shotcrete early age strength development clearly offers mining teams enormous advantages by allowing early, safe re-entry to the face. This paper illustrates that in the hard rock tunnelling field, modern shotcrete approaches do not cause delays in re-entry, but allow a continuous construction cycle. Furthermore, the paper identifies the most suitable, up-to-date shotcrete admixture and equipment technology available to the mining community to permit early re-entry.

May 1, 2003

By Yuki Tsunazawa, Sho Fukui, Chiharu Tokoro

"A milling process is one of the important unit processes in mineral processing. In the milling process, the control of particle size distribution is necessary for the later physical separation process. Since the milling process is strongly affected by the equipment designs, operation conditions and raw materials, the control for particle size distribution of milling products is mainly based on empirical rules. Recently, the discrete element method (DEM) has been widely used as an effective tool to investigate the behavior of grinding media. Although the behavior of grinding media can be accurately simulated, it is not well-established the way to predict particle size distribution after milling by the DEM simulation.The objective of this study was to develop a new method to predict particle size distribution in ball milling using the DEM simulation. To predict particle size reduction, a correlation between experiments and the simulation was investigated in various operation conditions, that is, wide range of the rotation speed and the filling ratio of grinding media. In experiments, after lime stone was ground, the average and variance of particle size distribution were evaluated using some kinds of distribution equations. In the DEM simulation, direct simulation of particle disintegration in milling process is generally very difficult because of huge computational load. Therefore, collision energy between grinding media balls orgrinding media ball and wall was used as a factor of particle disintegration.As a result, a good correlation between experiments and simulation was obtained. In addition, the particle size distribution in ball milling could be predicted using this correlation. Therefore, these results indicated that the use of the correlation between experiments and the DEM simulation was an effective approach of predicting the particle size distribution in milling process."

Jan 1, 2016

By I. G. L. Sinclair

"ABSTRACTTill geochemistry was the principal exploration tool employed by Selco in the investigation of a zinc prospect in Lanark Township, Ontario. The geochemical studies were undertaken after the discovery at surface of marble boulders containing bands of sphalerite and lesser amounts of galena.The source of these mineralized boulders was sought by means of B-horizon soil geochemistry, followed, where appropriate, by back-hoe trenching and till sampling. The thickness of the tiff in this part of Ontario seldom exceeds 3 m or 4 m and it was generally possible to sample complete profiles of the till and thus to obtain detailed information concerning the dispersion of zinc. Drilling and trenching at the up-ice extremity of the zinc dispersion trains showed the source of the float to be sub-economic bands of sphalerite in the Grenville marble. IntroductionTill geochemistry was used extensively by Selco during the period 1975 through 1980. It was the principal exploration technique employed in a program seeking carbonate-hosted lead-zinc deposits of the Balmat-Edwards type in the Grenville rocks of southeastern Ontario. This paper describes the use of the method to locate the bedrock source of a group of sphalerite boulders found at the surface near the village of Hopetown in Lanark Township (Fig. I) . Geochemical investigations at this locality were the subject of a previous publication by the author (Sinclair, 1979).All of the analytical work related to this project was done by Bondar-Clegg & Co. Ltd. in their Ottawa laboratory. In all cases the -80 mesh fraction of the till was analyzed. The total zinc content was determined by atomic ab sorption spectrometry after hot aqua regia digestion; readily extractable zinc was determined by atomic absorption after samples had been cold-leached for 18 hours in an ammonium citrate hydroxylamine hydrochloride solution, the pH of which was adjusted to 8.5"

Jan 1, 1986

By G. J. Kipouros, P. W. Kingston, W. F. Caley

"Current practice in advanced materials processing has involved the dispersion of a second-phase, commonly SiC, glass, or other non-metallic fibres/whiskers/particulates to enhance physical characteristics of the matrix. These matrices may be either ceramic or metallic but in either case it is important that the physical and chemical properties of each phase be compatible. The dispersed, finely divided second phase is often both expensive and foreign-produced and may exhibit chemical instability at elevated temperatures. Thus, it would appear that certain naturally-occurring industrial mineral fibres/whiskers/particulates might offer an attractive and low-cost alternative. Suitable candidates, such as wollastonite and sillimanite, are available in Canada and may be milled to produce particles whose shape is appropriate for this application. It is suggested that wollastonite, in particular, may be a useful modifying agent in metal matrix composites utilizing AI, Al-Si, Zn-Al or other light metals and alloys. In addition, the possible use of these minerals in ceramic matrix materials is also addressed, with reference to a preliminary thermodynamic analysis of the systems involved.IntroductionBackgroundA composite material may be thought of as consisting of a matrix component, the properties of which are improved by the addition of a second component or dispersed phase which forms a distinct three-dimensional interface with the matrix. This matrix is commonly either metallic or ceramic, which leads to the two most popular types of composites: metal-matrix composites (MMC); and ceramic-matrix composites (CMC).Such materials are by no means new. Early composites include fiberglass, mixtures of plaster and horsehair in wall plastering, as well as use of natural materials such as bamboo in fishing rods. Some of the newer applications include high strength, low-weight bicycle rims, as well as the more familiar role in aircraft and automobile engine components."

Jan 1, 1993

By T. Backers

In fracture mechanics modeling the rock mass is seen as the combination of rock material and fractures, which are evident on different scales. The activation and propagation of fractures dominates the stability and performance issues as well as fluid flow potential. Once fractures have started to grow into each other and create through going separations, stability of rock engineering constructions may be violated. Further, as fracture mechanics inherently incorporated a scaling law in the form of the participating fracture length, scaling issues can be more accurately mirrored. This contribution highlights some examples of application of a fracture mechanics based code on some rock mechanics issues such as numerical laboratory experiments, borehole stability or the creation of fluid pathways due to underground constructions. The results show explicitly the creation of fractures and connection of existing ones, and hence the safety and performance issues. It has to be concluded that the basis for frequent application of fracture mechanics on specialized challenges in Geomechanics is laid, and the needs for future developments are defined.

May 1, 2009

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