Search Documents

Flotation works very well for particles that are typically in the range 20 to 150 µmin diameter, for base metal ores. In this range, it is possible to get quite high recoveries in conventional flotation machines. Outside this range, the recoveries decline progressively, whether it be with the very fine or the coarse end of the size spectrum. In this paper, the reasons for the drop-off in recovery are briefly reviewed, together with strategies for future work. Recent progress in the author's laboratory will be described. These include the Concorde Cell, which takes advantage of the high dissipation rates in the shock wave created when a gas-liquid mixture reaches supersonic velocities. The environment thus created is very beneficial for the flotation of ultrafines. For coarse particles, the reasons for the drop-off in recovery can be related to the highly turbulent nature of the pulp in a conventional flotation cell. To improve recovery, it is necessary to find a way of bringing particles and bubbles into contact in a quiescent environment. A new process for coarse particle flotation will be described, in which a fluidized bed is created in the flotation cell. The flow conditions are very gentle, and the high solids concentration leads to rapid rates of capture of the particles. Experimental results are presented. It is apparent that the maximum floatable size for coal and minerals can be increased by a factor often over the limits found in current practice. Advances

Jan 1, 2009

By F. R. Nogas

"Because of the elements of uncertainty and instabilitp associated with single-enterprise communities, there is increasing interest in the alternatives to developing additional new communities of this nature. Gulf Minerals Canada Limited has chosen a unique alternative means to accommodate employees at a new uranium mine/mill complex at Rabbit Lake in Northern Saskatchewan. Gulf has adopted a fly-in program which moves their staff from a 450-mile radius, in and out of the on-site facilities, on a 7-day rotation."

Jan 1, 1976

By Kyle Paulsen

The Cameco Key Lake mill began processing U-bearing waste rock from the McArthur River mine in 2000. This ore often contains elevated levels of calcium (>4% by mass), resulting in less than favourable conditions in parts of the Key Lake mill, yet still contains U3O8 at a concentration of about 0.3%. To recover the U3O8 and reject as much calcium as possible, modifications were made to the grinding circuit in late 2007 to include a gravity concentrator. The concentrator allows this mineralized waste to be processed separate from the regular blend material, rejecting calcium and concentrating U3O8, thus improving the operating conditions in areas such as SX. This will be presented in greater detail along with results and further work.

Aug 1, 2010

By J. G. Lang

"This paper gives an outline of the principles used in E.M. wire rope testing, showing the practical application of the device and its value in industry and in mining operations.It presents: case histories involving deterioration, such as corrosion, peening and broken wires; anomalies indicating conditions peculiar to drum hoists and friction hoists; tests of locked coil ropes; low-frequency testing; and experiments and development."

Jan 1, 1969

By M. Rafieazad, P. Fathi, M. Mohammadi, A. M. Nasiri

This study examines the impact of low-temperature heat-treatment cycles on microstructure and corrosion performance of Direct Metal Laser Sintered (DMLS) AlSi10Mg_200C alloy. It is known that the morphology and the size of eutectic silicon are the two most critical factors affecting mechanical properties of the AlSi10Mg alloy. Although high-temperature heat-treatment cycles are common practice for the strengthening of the cast Al-Si-Mg alloys, in the case of DMLS fabricated AlSi10Mg alloy, high-temperature solutionizing and low temperature annealing cycles result in mechanical properties degradation primarily due to a dramatic decrease in the solubility of Si in the super-saturated a-Al matrix of the as-printed alloy. In this study, low-temperature heat-treatment cycles that only promote precipitation of Mg2Si accompanied by the interruption of the silicon network were investigated. SEM and XRD analysis were utilized to observe the evolution of microstructure after different annealing heat-treatments at 200, 250, 300, and 350°C for 3 h. Increasing the heat-treatment temperature contributed to the Si networks interruption and an increased precipitation of Mg2Si phase. Low-temperature annealing from 200 to 350°C was found to promote the homogeneity of the microstructure characterized by a uniform distribution of eutectic Si in ??-Al matrix. Additionally, in order to investigate the impact of different heat-treatment cycles on corrosion resistivity of DMLS-AlSi10Mg_200C, the potentiodynamic polarization testing and electrochemical impedance spectroscopy were performed in a 3.5 wt.% NaCl solution. The results revealed more uniformly distributed pitting attack on the corroded surface by further increase in the annealing temperature, which was attributed to the uniformity of Si dispersion. In comparison, low temperature annealing at 200°C led to a penetrating selective corrosion attack along the melt pool boundaries, resulting in a higher corrosion rate associated with the lower absolute value of impedance of the protective passive film on the surface.

Jan 1, 2018

By Jack Young

U. S. Steel Canada Inc.’s Lake Erie Works is the last green field integrated plant to be built in North America. The plant was designed using the best available environmental technology to develop an environmentally friendly production facility. A summary of the technology and measures used in the construction of the plant are summarized in this paper.

Jan 1, 2015

By D. J. Zuliani

A mass and energy balance for a modern electric arc furnace ("EAF") indicates that about 70% of the process energy losses are contained in the off-gas as chemical energy in the form of uncombusted CO and H2 and as thermal energy. Given the magnitude of the off-gas energy losses, improved EAF energy efficiency needs to focus first on controlling furnace practice to maximize "in-EAF" energy utilization and second, on recovering residual energy in the off-gas existing in the fume system. This paper focuses on the application of off-gas analyzers, sensors and process models to control and optimize "in-EAF" energy utilization efficiency. The paper will also discuss the application of heat recovery technology to EAF steelmaking.

Jan 1, 2015

By Lars A. Bergstrom

High-voltage direct -current transmission is a relatively new technology. The successful development and operation of AC/DC converters, capable of handling large power blocks, have made HVDC transmission generally accepted by utility engineer s and planners around the world. HVDC technology still has significant development potential and, as a result, it is expected that equipment costs will escalate at a substantially slower rate than those of conventional AC equipment. This development will also result in still more compact terminal stations. This paper reviews two modern HVDC installations. In Zaire in Africa. a 1700-km overhead DC line rated at 500 kv, will supply power to the copper mining district in the Shaba region. Another HVDC link connects Norway and Denmark through a submarine cable rated at 250 kv. This is the longest HVDC cable in the world, 127 km, and the water depth , maximum 530 m, is also much greater than for any previous similar cable. The reasons for selecting HVDC are discussed, and the design and layouts are presented.

Jan 1, 1978

By Yongsheng Song, Yong Chen, Wenjuan Li

"Detailed characterization and improved recovery of a micro-fine, low-grade refractory gold ore from a beneficiation plant in Yunnan Province, China, was investigated. Different analytical techniques, such as videlicet size analysis, scanning electron microscopy, and mineral analysis by a Mineral Liberation Analyser, were used for mineralogical characterization of the ore. This paper describes an experimental study on the effect of the grinding fineness and the use of various flotation reagent regimes involving collectors, depressors, and activators on the flotation behaviour of micro-fine auriferous pyrite. Two-stage flotation resulted in high gold grade (32.71 g/t) and gold recovery (80.02%) of the bulk gold concentrate. This process opens the possibilities of maximizing Au recoveries in a micro-fine low-grade refractory gold ore where low recoveries are common.INTRODUCTIONWith the rapid development of the gold industry, refractory gold ores are becoming increasingly important resources for gold production (O’Connor & Dunne, 1994; Agorhom, Skinner, & Zanin, 2012). Among these refractory gold ores, sulfur and micro-grained native gold ores represent a significant challenge, since they are not amenable to gold recovery by direct cyanidation. The gold contained in this type of ore is very finely disseminated and encapsulated in sulfides, which makes it inaccessible to the cyanide solution during direct cyanidation (Bulatovic, 1997; Allan & Woodcock, 2001). Autoclave pressure oxidation (AC POX) of the sulfide ore and subsequent cyanidation is a common technology used to recover the gold in such types of ore. Flotation is often used as a pre-treatment method to refractory gold ore (Monte Lins, & Oliveira, 1997).As the gold-carrier minerals are usually finely disseminated (sometimes under 10 µm) and intimately associated with gangue minerals, a fine grind is necessary to get a high degree of liberation. At the same time, the fine grind usually causes the production of significant quantities of slimes, which comprise easy-to-float gangue minerals such as mica, dolomite, and calcite. These mineral phases are naturally hydrophobic and float readily, competing with the value minerals in the froth phase and resulting in lower concentrate grades."

Jan 1, 2016

By H. Lundberg

Introduction The world has been suffering from an excessive production of all kinds of raw materials, and the metals and petroleum are no exception. As a consequence, at the moment, the urge to find new mineral deposits is almost non-existent and prospecting languishes. But from previous experience we know that sooner or later new supplies will be needed. Mineral deposits become exhausted, and to keep up the supplies of our metals and petroleum, new deposits must be found. Even now, a supply of petroleum of domestic production would not be unwelcome in Canada. It is not, therefore, so entirely malapropos to discuss means of finding ores and oil as it would seem. Almost any physical property connected in definite manner with minerals and rocks may be made the basis of geophysical exploration. This is particularly true of electrical resistivity, and its utilization in electrical prospecting has been the deciding factor in the discovery of such big ore-bodies as Buchans mine, Newfoundland, and the large copper ore-bodies in northern Sweden, etc., as well as in the solution of many foundation and structural problems. The importance of electrical prospecting has been growing steadily, and already there exists hardly any exploration project involving geological problems which may not profit by its application. Electrical prospecting is marching on, and neither temporary failures, nor incompetence, nor over-optimism, nor ill-will and pessimism, nor bad times and business depression, will be able to halt its progress.

Jan 1, 1931

By Claude Hubert, Robert Sansfaçon

"The Malartic gold camp includes the Malartic Gold Fields, Barnat, East-Malartic, Canadian Malartic and Sladen mines. From 1935 to 1984, these mines have yielded 212 403 kg of gold. They occur within a 900 m wide belt of highly sheared, dislocated and faulted volcanic and sedimentary rocks (Piché and part of the Pontiac groups) herein named the Malartic Tectonic Znne (MTZ). The zone trends WNW-ESE and separates two lithotectonic sedimentary domains, the Cadillac to the north and the Pontiac to the south. The sedimentary domains are characterized by their coherent fold generations consisting of open to tight reclined and recumbent F1 folds that have been refolded into series of asymmetric upright F2 folds. A flat-lying schistosity (S1 ) and a vertical pressure-solution cleavage (S2) are ubiquitous throughout. Discrete sets of brittle and brittleductile faults are synchronous and parallel to each of these two fabrics.The MTZ trends N320° in the northwestern sector of the town of Malartie, changing to N2800 at its easternmost extremity. The Malartie and Sladen faults constitute its northern and southern boundaries respectively. These two faults also bound another set of faults oriented N305°. Throughout the MTZ, the faults at N280° and N320° and those at N305° outline a series of anastomosed motifs that constitute the principal structural feature of the Malartic district. Locally, a fourth set of minor faults at N065° combines with the other three prominent sets. All these faults have disrupted the gold-bearing volcanic, sedimentary and intrusive rocks encountered throughout the MTZ.The MTZ comprises numerous lozenge-shaped blocks of komatiites and basalts or greywackes and black shales which are cut by disjointed dikes, sills and small intrusions of felsic and mafic rocks. The latter includes gabbro, diorite, lamprophyre, monzonite and pegmatite. They frequently fill conjugate sets of fractures. In some areas, they are clustered into coalescent intrusions. The present morphology of these intrusive bodies( cube, chocolate-tablet boudin, parallel epiped and cigar) is not representative of their primary forms but the result of late shearing and faulting indicative of a major tectonic event subsequent to their emplacement. The intrusions were injected during a brittle stage of deformation related to the F1 folding but prior to the development of the F2 folds and subsequent tearing and shearing by late brittle and brittle-ductile faults . The ultramafic and mafic rocks are intensely altered to serpentine, tremolite, actinolite, carbonate, biotite, hematite and pyrite according to the nature of their proto/it h. In the vicinity of the MTZ, all lithologies exhibit mineral assemblages indicative of the upper greenschist or lower amphibolite facies of metamorphism.The structural evolution of the MTZ involved superimposed folding, shearing, transposition, mylonitization, brittle and brittle-ductile faulting. Each lenticular motif results from the combined effect of conjugate sets of shear zones that outline lozenge-shaped blocks where strain markers indicate different strain states from cores to margins. Bedding in ultramafic flows and/or sedimentary rocks is either NW-SE or NE-SW within the blocks but, in the vicinity of the Malartie and S laden faults, it has been rotated to become subparallel to the faults.Ore shoots are most commonly associated with altered and deformed bodies of diorite and monzonite, or some mafic and ultramafic schists and, at the Canadian Malartic, Sladen and East Malartic mines, with metasedimentary rocks. The shoots occur in close proximity to the various faults of the MTZ but the gold itself is intimately associated with pyrite, biotite and carbonate in halos bordering fractures crisscrossing the different host rocks. The density of fractures and the width of their altered zones have been significant parameters in the production of economic ore shoots. Grains of gold are visible in some veins of quartz a"

Jan 1, 1990

By J. Lobel, M. Bourassa, J. A. Kozinski, B. Peng

Non-isothermal kinetics studies of the reduction phenomenon of electric arc furnace (EAF) dust generated in stainless steelmaking have been carried out. The dust was mixed with carbon as a reducing agent and dolomite as a binder to form pellets. A Thermo-Gravimetric Analyzer (TGA) was used to determine the weight loss of pellets heated at an average rate of 40 K/min. up to a temperature of 1500°C. The degree of metal reduction was calculated by considering the pellet weight loss, and evaporation of moisture, dust, zinc and lead at high temperature. The reduction process was divided into three different steps. The first covers the 300°C to 760°C temperature range (from 7 to 18 min.), the second step ranges between 760°C and 980°C (18 to 24 min.), and the third covers temperatures above 980°C. Each step is characterized by different chemical reactions. Non-isothermal reduction kinetics equations were set up to describe each step. The kinetics parameters such as apparent activation energies and frequency factors were established. It was found that the first step was chemically controlled, while the second step was diffusion-controlled. Reduction behaviour in the third step was strongly dependent on the initial carbon content in the pellet. Chromium oxide (CrO) could only be reduced in this last step. High temperatures and higher initial carbon content favoured the reduction of CrO.

Jan 1, 2001

By A. J. Williams

"IntroductionAN EXCELLENT paper entitled Potash Discoveries in Western Canada (l) was presented to the 1948 Annual Meeting of the C.l.M. by L. Heber Cole, of the Department of Mines and Resources, Ottawa. In the intervening three and a half years, the search for oil and gas in Saskatchewan has added considerably to the information available on the occurrence of potash. It is the purpose of this paper to summarize the additional data. For greater clarity, it may be useful, first, to recapitulate some of the highlights of Mr. Cole's paper.Underlying the western plains of Canada there is a large deposit of sodium chloride, believed to be of early Devonian age. Associated with the sodium chloride are certain potash salts, primarily 'sylvinite,' a mechanical mixture of sodium chloride (halite) and potassium chloride (sylvite) crystals. The potash salts, when present, are usually near the top of the salt formation, though in at least one well they were found near the central part of the deposit. Cole reported on some seven wells in Saskatchewan which had encountered appreciable amounts of potash; these were the Ogema, Radville, Davidson Nos. 1 and 2, Prairie Salt Nos. 1 and 2, and Verbata wells. Depth of potash from surface ranged from 7,653 feet in the Radville well to a minimum of 3,466 feet in the Verbata well. Such depths do not compare favourably with those at which the potash is found in mines in the United States, France, and Germany, where about 1,000 feet is the average. The greatest concentration of sylvite encountered was in the Verbata well, where a bed 11 feet thick carried the equivalent of 21.64 per cent K20."

Jan 1, 1952

By Bernhard Schwab

A popular rule in the mining industry is that a property should only be developed if average unit production cost is less than a given proportion of the estimated long-run price of the mineral. This paper provides a comparative analysis of the average cost approach (ACA) with the use of net present value (NPV) as an investment criterion. It is concluded that the ACA is an easily understood and easily applied rule of thumb which may serve as a rough first screening device, for example in the early stages of deposit exploration. Its operational convenience is at the expense of accuracy and validity, and the NPV provides a sounder method of making decisions, especially in the case of marginal deposits.

Jan 1, 1978

By J. Ambrose, R. W. Zimmerman

"Shales are highly anisotropic in their mechanical behavior. The strength ( ??1) of anisotropic shales depends not only on the magnitude of the other two principal stresses ( ??2, ??3), but also on the bedding plane orientations (?) relative to the principal stresses. In this study, the role of the intermediate stress ( ??2) on the strength of shale is evaluated using new data from triaxial compression and extension tests. To understand the failure behavior of shales, we have carried out triaxial compression tests ( ??1 > ??2 = ??3) on two organic rich mudstones, at different confining stresses and loading angles. We first fit the triaxial compression data with Jaeger’s plane-of-weakness model, and find that the model is able to represent the strength behavior of anisotropic shales reasonably well. We then carried out five triaxial extension tests ( ??1 = ??2 > ??3) for a high strength-anisotropy shale, and two extension tests for a low strength-anisotropy shale. Results from the compression and extension tests for failure along the weak plane shows that the role of ??2 is not significant for weak plane failure. However, for the low strength-anisotropy shale, the failure strength of the intact rock (i.e., at ??= 0o and 90o) for compression and extension tests were significantly different. This indicates that for intact rock failure, the role of ??2 is important and cannot be ignored.INTRODUCTIONFor many years, various researchers have studied the roles of ??2 for failure of isotropic rocks. (Mogi, 1967 and 2007; Colmenares and Zoback, 2002; Al-Ajmi and Zimmerman, 2005 and 2006; Haimson, 2009; Ma and Haimson, 2012). The findings from these studies show that the effect of ??2 on the strength of isotropic rocks is generally significant, and cannot be ignored.Since the role of ??2 in the failure of isotropic rocks was recognized, other researchers attempted to extend some of these isotropic failure criteria to anisotropic rocks, but with limited success (e.g., Tiwari & Rao, 2007; Singh et al., 1998; Zhang and Zhu, 2007; Pei, 2008). Smith and Cheatham (1980), on the other hand, conducted true-triaxial experiments on the Green River Shale and used a J2-I1 type relationship to predict the strength behavior, but without considering the effect of the three-dimensional orientation of the weak plane on strength. Therefore, unlike for the case of isotropic rocks, additional studies need to be done to determine the role of 2 on the failure of anisotropic rocks."

Jan 1, 2015

By E. Haccuria

Due to the growing share of electrical cars, the use of lithium ion batteries significantly increased. Recycling of these batteries is essential to safely dispose hazardous materials as well as to recover valuables (Cu, Ni, Co, etc.). In order to optimise the recycling process, accurate information is required on the slag phase equilibria for battery smelting. The slag is a multi-component system containing alumina from the battery cases, silica and lime from fluxes, and manganese and lithium from the electrode materials. The present study accurately determined the phase equilibria in the "MnO"-Al2O3-SiO2 sub-ternary system at metal alloy saturation. Particular attention was given to the accuracy and reliability of the final results by highlighting the different reaction pathways, mass transfer- and reaction mechanisms taking place in the system during experiments. This attention enabled improved experimental methodology and measurements. Based on the present study, further measurements will be undertaken in the quaternary "MnO"-Al2O3-CaO-SiO2 system at low "MnO" concentrations, which closely represents the future recycling industrial slags.

Jan 1, 2014

By V. A. G. Preto, W. H. White, J. E. Harakal

Potassium-argon age dating performed recently on intrusive rocks from Copper Mountain indicates that the satellite Smelter Lakes and Voigt stocks have radio-metric ages which are indistinguishable from that of the Copper Mountain stock, as determined by Sinclair and White (1967). The complex of intrusive rocks known as the Lost Horse Intrusions also yields essentially the same radiometric age, as do biotite-sulphide-pegmatite veins on Copper Mountain.

Jan 1, 1971

By X. G. Zhao, M. Cai, L. K. Ma, P. F. Li

"Crack initiation stress is one of the characteristic stress thresholds in the rock deformation process under compression. By analyzing AE responses of brittle rocks in uniaxial compression, a new method which is called cumulative AE hit (CAEH) method is proposed for determining the crack initiation stress. An advantage of the CAEH method is that the method removes the user's subjective judgment. To verify the CAEH method, uniaxial compression tests with AE measurement using 77 crystalline rock specimens were conducted. The test specimens cover eight rock classes with different mineral grain distributions. Crack initiation stresses of these specimens were determined by the proposed CAEH method and the LSR method, and the resulting mean crack initiation stresses determined using the two methods were comparable. In addition, taking one rock class with 20 specimens for example, the mean crack initiation stresses determined using the CAEH method and those from several established strain-based and AE-based methods were similar. Hence, the CAEH method can be used as an alternative method for crack initiation stress determination.INTRODUCTIONExperimental observations indicate that the failure process of brittle rocks is associated with crack initiation, crack propagation, crack damage and final strain localization (Brace et al., 1966; Bieniawski, 1967; Lockner, 1993; Wawersik & Fairhurst, 1970). Compression tests on cylindrical specimen are often used to characterize the crack development process. The crack initiation stress (sci), crack damage stress (scd), and peak stress (sc) represent important stress thresholds during the failure process of brittle rocks. Among the three characteristic stress magnitudes, sci is difficult to be determined objectively, directly from examining the stress–strain curves. Some researchers argued that sci obtained from laboratory compression tests is an important parameter for estimating the spalling failure observed in brittle rocks around underground excavations (Andersson et al., 2009; Diederichs, 2007; Martin, 1997). sci can also be used to estimate rock tensile strength and Hoek-Brown failure criterion parameter mi (Cai, 2010). Over the last 40 years, various strain-based methods have been proposed to determine the crack initiation stress in compression tests (Brace et al., 1966; Diederichs, 2007; Eberhardt et al., 1998; Lajtai, 1974; Martin & Chandler, 1994; Stacey, 1981). Nicksiar and Martin (2012) proposed a Lateral Strain Response (LSR) method for establishing sci for low-porosity rocks. Meanwhile, they also used different strain-based methods to evaluate sci of Äspö diorite, and the results showed that all strain-based methods gave acceptable crack initiation stress values. Subsequently, they utilized the LSR method to further investigate the crack initiation stresses in igneous, metamorphic, and sedimentary rocks in compression tests (Nicksiar & Martin, 2013). Compared with previously proposed strain-based methods, an advantage of the LSR method is that it removes the user judgment when interpreting sci. However, the use of the LSR method relies on accurate determination of scd and on having a polynomial fitting equation to find the maximum LSR value."

Jan 1, 2015

By I-H. Jung, J. Kim, B. Konar

"The usage of FeNdB permanent magnet, containing Rare earth elements (REEs), is steeply increasing with its demand for wind-turbine and electrical vehicle. The improvement and recycling of the permanent magnet require good understanding of chemistry of all constituent elements. In order to understand the chemical interactions between REEs and Fe, all Fe-RE binary systems have been critically assessed and optimized based on reliable literature data. The Modified Quasichemical Model and Compound energy formalism is employed to model the liquid phase and the solid solutions present in the binary systems respectively. The enthalpy of the liquid and the intermetallic phases showed a systematic trend along the lanthanide series. This helped to mitigate the discrepancies in literature data and also to predict the binary interactions between some less well investigated Fe-RE binary systems. INTRODUCTIONRare Earth elements (REEs) are a group of seventeen elements (including Scandium, Yttrium and the lanthanide series) that exhibit same characteristic chemical and physical properties. These properties make them valuable to the electronic, appliance, green technology, weapon and medical device manufacturing industries. Compounds of REE with other metals (like Co, Fe, Ni etc.) have unique properties of luminescence, thermal and electrical conductivity, magnetism and ability to act as catalysts and polishing compounds (Gschneidner ,2011).They are not only indispensable for some electronic gadgets of today, like computers, music systems, mobile devices, televisions and MRI machines, but also holds a key to more greener technology of the future. They contribute to the development of more efficient and environment friendly version of everyday products such as cars and light bulbs. REE compounds revolutionized the magnet industry by making them smaller and more powerful with increased longevity, and high temperature workability, which are used in motors and generators that power electronics, telecommunications, electric automobiles and turbines(Humphries, 2011)."

Jan 1, 2012

By M. L. F. Paullo, S. Cristian Mejia, D. Roehl

"Rock fracture simulation most often presents non-linear behaviour due to loss of material integrity and strain concentration in the fracture zone. These strains generate local instabilities that propagate to the global system of equilibrium equations, generating convergence problems of the numerical solution. Often, in such cases, analyses in the post-critical regime are unfeasible. To overcome the numerical difficulties associated with problems of stiffness loss in the post-critical regime control methods, also called continuation methods, are employed. Among them are the arc-length method, the energy control method and the indirect displacement control method. These methods are employed in conjunction with a Newton- Raphson scheme for the solution of non-linear systems of equations. However, depending on the system instability, convergence in not guaranteed even with very small increments. In the present work these control techniques are investigated in combination with cohesive elements to simulate rock fracture. The constitutive model employed with the cohesive elements represents stiffness degradation through a damage law, which leads to serious convergence difficulties as reported in the literature. Here the continuation methods mentioned are applied to a Mode I rock fracture problem and to hydraulic fracture simulations. The effectiveness of the different continuation methods is compared.INTRODUCTIONThe process of rock fracture induced by the fluid pressure is called hydraulic fracturing. The propagation of a hydraulic fracture is a complex process, which is basically defined by the mechanical deformation of the rock, the fluid flow within the fracture and the fracture propagation itself (Mokryakov, 2011). In the works of Bendezu et al. (2013), Chen (2012) and Carrier et al. (2012), vertical hydraulic fracture propagation is modeled by cohesive elements with traction-separation law. The softening behavior of rock in post-critical regime is demonstrated by experimental testing and is characterized mainly by high strains and low stresses (Crowder and Bawden, 2004). This nonlinear behaviour can arise due to material nonlinearity, such as micro-cracking and damage, leading to softening behaviour (Chandrakant, S.D., 2012). The cohesive fracture model considers in its formulation the softening behaviour that most rocks present in post-critical regime with propagation criteria."

Jan 1, 2015