Search Documents

By A. Galibois, H. Monette, J. L. Collins

This paper deals with the problem of caving, which occurs in varying degrees in blast holes before loading. It is shown that the distribution of the amount of caving in materials may be described by a gamma density function. A decision model to determine the optimal depth of over-drilling 'in order to minimize the expected additional cost per hole is described. Two different policies are discussed. In addition, charts which permit the mine operator to realize significant savings have been drawn for different actual situations.

Jan 1, 1971

By Narasimha Mangadoddy, Balraju Vadlakonda

"Bubble size distribution, bubble rise velocity and gas and solids hold-up are the key major hydrodynamic parameters that are extensively used for column flotation design and performance evaluation. The distribution characteristics of mean gas hold-up is studied using electrical resistance tomography (ERT) coupled with pressure transducers (PT) in column flotation. The ERT technique facilitates non-invasive and nonintrusive visualization of different flow regimes even in the presence of particles in a column flotation. The effect of superficial gas velocity, feed flow rate and slurry height in the column on mean gas hold-up and its radial distribution has been analysed for both two-phase and three phase flows. Experimental results show ERT coupled with PT is suitable as an online monitoring tool to provide useful information on the three-phase flow hydrodynamic parameters of the column flotation. The tomography images, which were generated using a modified sensitivity back projection algorithm, are employed to explore the influence of parameters in two phase flow in column flotation. Mean gas and solids hold-up values extracted from ERT have been critically analysed for the column operating in various flow regimes. Gas hold-ups and solid hold-ups results using ERT are in very good agreement with conventional estimation and correlations obtained using pressure transmitter methods. Meanwhile, the results show that the gas hold-ups in the centre region increase constantly with an increase in the superficial gas velocity, whereas solid hold-ups distribution is very heterogeneous for high gas velocity. Furthermore, the results also indicate that ERT is a very powerful tool for diagnosing the 'inside' flow behaviour of gas –liquid–solid three phase flow in columns operating for a wide range of flow regimes.INTRODUCTIONColumn flotation is perhaps the most significant technological innovation in flotation machinery and its successful introduction in mineral processing plants. There has been a lot of attention in understanding the various mechanisms that take place during operation. Due to many reasons column flotation has been found to yield better performance than conventional flotation cells, particularly with fine particles (Honaker & Mohanty, 1996). Due to complex flow behavior in nature, the performance of column flotation influenced by on flow characteristics of the process, (Jena, Biswal, Das, & Reddy, 2008). Previously basic research is carried out to evaluate the effect of the hydrodynamic parameters on performance of column flotation (Finch, J.A.;Dobby, 1991). Many techniques have been developed that could be used to experimentally analyse the hydrodynamics of column flotation. Tomography techniques are well sound to investigate the hydrodynamics for process equipment. Among that, Electrical resistance tomography (ERT) has been used in various industrial investigations for visualization of the concentration profiles and characteristics of the fluid dynamics in gas–liquid two-phase systems (Wang, Dickin, & Mann, 1999; Williams, R.A; Wang, 2000). ERT is an imaging technique to analyse the internal flow behavior non-intrusively. ERT provide very useful tools for measuring and monitoring of column flotation process on line. Due to its high speed capability, low cost measurement, robust sensors, ERT is considered to be the most powerful tool among other tomography techniques such as positron emission tomography (PET), magnetic resonance imaging (MRI), optical and infrared tomography. Application of ERT has extended in process systems due to its improved speed of data acquisition, sensitivity, flexibility and noise immunity. In multiphase flow of bubble columns, gas hold-up and solid hold-up are explored by using the coupling of ERT and difference pressure method (Jin, Han, Yang, & He, 2010)."

Jan 1, 2016

By K. Himbeault, D. Campbell, N. C. Garisto, D. A. Orava, A. Rosaasen

"The level of scientific knowledge about the effects of contaminants is increasing and changes in environmental monitoring are likely to allow for the assessment of environment effects particularly to the aquatic environment. This paper provides an overview of recent developments in this area which include the environmental effects monitoring program for pulp and paper operations, the assessment of the Aquatic Effects of Mining (AQUAMIN) program, the Aquatic Effects Technology Evaluation (AETE) program, and the review of the federal Metal Mining Liquid Effluent Regulations (MMLER). The paper also provides a look forward at the direction of environmental monitoring of mine projects in Canada. IntroductionEnvironmental monitoring programs provide a wide range of data to environmental management systems, regulatory authorities, and stakeholders. The data obtained are used to demonstrate that operations are in compliance with regulatory environmental protection requirements, and provide the scientific basis for assessing the effectiveness of environmental protection measures. Environmental monitoring requirements for mine projects are established based on the type of operation, site specific issues and concerns, and the requirements of regulators and mine proponents.The level of scientific knowledge regarding the potential effects of contaminants on the receiving environment is increasing, and in particular with regards to the receiving aquatic environment. Regulatory controls regarding the protection of fish, fish habitat, and fishery resources are being reviewed with regards to their adequacy, and the environmental effects monitoring program on the pulp and paper industry is expected to provide useful data relevant to those types of operations. These developments are of interest to mine proponents as it is presently expected that an environment effects monitoring program will be developed specifically for the metal mining industry under a regulated environmental effects monitoring program as part of the Fisheries Act. This paper reviews recent developments in this area, and provides a look forward on the future directions of environmental monitoring programs."

Jan 1, 2000

By A. Karzhau, B. Aldamzharov, Y. Tussupbekov

"This case study produced for the metal mine operating in Kazakhstan, which has ventilation controlled manually and working continuously at maximum speed to keep air flowing through each part of the mine 24 hours a day. The mine auxiliary ventilation has been designed with the maximum length of ducting systems with extra air volumes to compensate for potential leakage. The mine provided auxiliary ventilation to all production stopes even to those, which for various operating reasons become temporarily inactive. However, mine safety has to consider risks associated with the deficit of fresh air for mine gases dilution and removal. The mine needs additional investments for re-designing ventilation network and therefore needs to estimate financial risks associated with it. The aim of this study is to conduct risk analysis for ventilation network optimisation processes at the operational metal mine in Kazakhstan. Mine ventilation surveys were undertaken for validation of the computerised ventilation model. The mine ventilation network modelling has been undertaken to be compliant with Kazakhstan mining regulations. Results of the mine ventilation optimisation processes were analysed in terms of financial risks for investments into reconstruction of the ventilation shafts, extension of the ventilation network, variable speed fans and automated regulators. Savings were resulted from improving ventilation system efficiency by reduce airflow in areas without mining activity and redistributing existing ventilation capacity. Based on savings derived from the optimised ventilation network an additional NPV has been received. Stochastic modelling use Monte Carlo simulation producing representation of possible outcomes for the additional NPV considering electricity pricing as a major variability. As result probability distribution for the optimised ventilation network demonstrates the additional NPV in a range of 1.6 - 2.2 million US$ in the corresponding defined value of 90%. Moreover, there is only 5% chance that the additional NPV to be less than 1.5 million US$. Conducted risk analysis for ventilation network optimisation processes gives considerations on volume of investments into the mine ventilation system.INTRODUCTIONCosts for mine ventilation are typically increasing when mines expand to deep levels. Some underground mines in Kazakhstan still use continuously operating fans at its maximum capacity. This is not helping to save energy and costs that are in many cases very significant. This case study produced for the metal mine operating in Kazakhstan, which has ventilation controlled manually and working continuously at maximum speed to keep air flowing through each part of a mine 24 hours a day. The mine auxiliary ventilation has been designed with maximum length of ducting systems and extra air volumes to compensate for potential leakage. The mine provided auxiliary ventilation to all production stopes even to those, which for various operating reasons become temporarily inactive. However, mine safety has to consider risks associated with the deficit of fresh air for mine gases dilution and removal from the mine workings (Sabanov, 2019). The mine needs additional investments for re-designing ventilation networks and therefore need to estimate financial risks associated with it. The aim of this study is to conduct risk analysis for ventilation networks optimisation processes at the operational mine in Kazakhstan. Mine ventilation surveys for validation of the computerised ventilation model have been undertaken"

Jan 1, 2019

By Jan Butra

Longstanding mining of copper deposits in the Legnica-Glogow Mining Basin, Poland has accumulated broad experience concerning the impact of geological and mining conditions on seismic and rockburst hazard. Methods of rockburst prevention for deep copper mines have been developed based on gathered data. Methods of rockburst hazard evaluation and rockburst control constitute a set of rockburst prevention measures. Due to their nature and location within the production process, prevention methods can be passive or active. Active methods of rockburst prevention consist of the liquidation of stresses in the areas where they are concentrated through controlled generation of rock-mass tremors and rockburst by blasting works (group extraction blasting, release blasting in the rock mass and pillars able to accumulate the stress, as well as in the roof and floor). Group blasting works as a basic method are very effective in copper mine conditions. Provocation of dynamic events makes it possible to control the timing of their occurrence. In Polish copper mines, all blasting works are made during the absence of miners in the mining fields. Daily distribution of both the number of recorded strong seismic events and volume of energy released, confirm the advisability of using this kind of rockburst hazard prevention. The rockburst hazard prevention methods developed over 40 years of copper deposit extraction are economically reasonable and offer sufficient effectiveness under high rockburst hazard conditions.

Aug 1, 2013

By Kazutoshi Haga, Batnasan Altansukh, Yasushi Takasaki, Mohd Yazid Bin Abdul Razak, Atsushi Shibayama

"In this study, the volatilization behaviour of arsenic (As) and antimony (Sb) from arsenicantimony bearing minerals such as tennantite (Cu12As4S13)/tetrahedrite (Cu12Sb4S13) in copper ore was investigated using a roasting process. Specifically, the thermal decomposition of tennantite and the arsenic volatilization were investigated under nitrogen and air atmosphere. The roasting experiments were performed at temperatures ranging from 500 to 1,200oC and different retention times using a laboratory electric furnace. The results indicated that over 95% volatilization was achieved at 700oC. On the other hand, when the temperature increased to 1,200oC, volatilization rate of Sb reached to 90%. Therefore, volatilization of Sb requires high temperature compared to As. It suggests that As and Sb can be separated using temperature difference. Selective removal of As is possible at 700oC, subsequently Sb can be recovered as a volatilized product by increasing the temperature to 1,200oC. Meanwhile, volatilization of Cu and Fe was not observed on these conditions. The content of As and Sb in the residues obtained from roasting at 700 and 1,200oC was less than 0.3 and 0.75 mass%, respectively. The Cu content in solid residue at 1,200oC increased to 34 mass%. The clean copper concentrate with lower amounts of arsenic and antimony was obtained from the process. INTRODUCTIONRecently, arsenic (As) and antimony (Sb) contents in copper (Cu) concentrates are increasing year by year. Arsenic and antimony are present with Cu minerals such as enargite (Cu3AsS4), tennantite (Cu12As4S13) and tetrahedrite (Cu12Sb4S13) that associate with copper sulfides [chalcopyrite (CuFeS2) and bornite (Cu5FeS4)]. They are normally concentrated by flotation and treated in copper smelters. However, treatment of copper concentrates containing high amounts of arsenic and antimony by a smelting process is difficult due to environmental implications of toxic arsenic and antimony. As demand for copper continues to grow, the need to source copper from several types of resources including those that contain high levels of As/Sb becomes necessary. Therefore, a novel processes for the treatment of As/Sb bearing Cu concentrates are required. Reviews on copper resources such as enargite, tennantite and tetrahedrite containing high amounts of arsenic and antimony are covered widely in literatures by many researchers (Padilla et al., 1998; Welham, 2001; Vinals et al., 2003; Guo & Yen, 2005; Balaž & Achimovicova, 2006a; Filippou et al., 2007; Mihajlovic et al., 2007; Lattanzi et al., 2008)."

Jan 1, 2016

By T. Plikas, L. Gunnewiek

Computational fluid dynamics (CFD) modeling was used to evaluate the performance and optimize the design of a potroom ventilation system in an aluminium smelter facility. The ventilation system must ensure the inflow of fresh air to the workplace by natural draft and the removal of heat and hazardous emissions generated from the pot cell process. A full-scale, three dimensional CFD model of the potroom was developed to predict the air temperature distribution and contaminant concentrations throughout the potroom for the purpose of verifying that the industrial workplace hygiene standards will be met by the proposed ventilation system design. The modeling work was used to optimize the size of the floor gratings around the pot cell and roof ventilator, as well as determine damper settings in the roof ventilator. The effects of outside air temperature and external wind on the performance of the ventilation system was evaluated to ensure the design performs adequately in response to seasonal variations in climatological conditions.

Jan 1, 2006

By M. E. Holter

"Limestone-bearing strata are exposed in seven areas of Alberta, each situated adjacent to rail and road facilities. Six of these centers are located along the eastern edge of the Foothills and Front Ranges of the Rocky Mountains in western Alberta. The other is in northeastern Alberta, along the eastern margin of the Alberta Basin. High-grade limestones are present in Middle Cambrian, Upper Devonian and Mississippian strata. Middle Cambrian (Eldon Formation) limestones have been exploited only in the Exshaw area west of Calgary, where they formerly were used for lime production. Limestone beds of the Upper Devonian Palliser Formation are exploited at Exshaw and at Cadomin in west-central Alberta for the manufacture of cement. Reserves of high-grade limestone also are available in the Palliser Formation near Crowsnest Lake in southwestern Alberta and near BrUle, west of Edmonton. Upper Devonian limestones of the Waterways Formation are widespread in the Fort McMurray area of northeastern Alberta, but available data indicate they are of inferior quality to the limestone deposits of the Rocky Mountains regions. The Mississippian Livingstone Formation contains large reserves of high-grade limestone at several localities in the Rocky Mountains and Foothills. Commercial development of the unit is now underway at Crowsnest Lake and at Exshaw for lime production, and future exploitation of these beds is feasible at Blairmore and Nordegg."

Jan 1, 1973

By Iain G. Bruce

The need to improve supply reliability while improving our environmental record during the past half decade of low metal prices has required the mining industry to prioritise expenditures. Quantitative Risk Assessments using third party reviewers and specialist risk management teams have assisted in identifying hazards that can both disrupt supply and impact the receiving environment. Graphical displays using risk ellipses can simplify the presentation of the results in a Probability-Consequence (PC) graph. Development of risk acceptability envelopes using readily available financial data superimposed onto the same PC space assists management to prioritise their risks and focus on risk reduction measures. Development of this system and examples of the successful implementation of the system from several recent projects are described.

May 1, 2003

By V. Aprahamian, G. B. Harris, C. W. White, M. Tangen

"In 2005, the Royal Canadian Mint (RCM) commissioned a new silver electrorefining circuit. This circuit was intended to operate with a minimal electrolyte bleed for the control of copper. However, by the time of start-up, the levels of copper-rich silver doré to electrorefining had increased, and despite partial pyro-conversion in a small Top Blown Rotary Converter (TBRC), copper concentration in the silver anode still averaged 7.5%. This necessitated frequent replenishment of the electrolyte with fresh silver nitrate since copper displaced dissolved silver during electrolysis, and therefore requiring a robust copper bleed treatment step in order to maintain cathode silver flake purity at 99.99%. The originally-installed bleed treatment system could not meet the RCM’s needs. A thorough review of known technologies, however, did not indicate one capable of handling the copper load, and resulted, therefore, in the development of a novel copper hydrolysis approach for removal of copper as an insoluble basic nitrate salt, and recovery of the associated nitric acid for re-use in the circuit. Laboratory-scale trials showed reproducible results, and existing equipment was subsequently modified.INTRODUCTIONIn 2005, the introduction of a new silver electrorefining circuit at the Royal Canadian Mint (RCM) initiated changes in the way customer silver, collected from the gold doré circuit, was processed. These circuit changes were discussed in a recent paper (Aprahamian, 2015). A wastewater treatment plant had not been budgeted at the time, and this represented a significant constraint when considering the introduction of a wet processing circuit. After a review of various options, a turnkey silver electrorefining circuit was purchased from a European firm. Initially, enough cells were ordered to produce 2 million oz of fine silver (999), with the option to increase to 5 million oz. The selected equipment was based on a Moebius cell design combined with a silver grain pressure dissolver to produce make-up silver nitrate solution. Additionally, an evaporator/hydrolysis unit included in the package appealed to the Mint Refinery, on the basis that it offered a simple way to bleed dissolved copper in electrolyte without generating effluent."

Jan 1, 2016

The problem of thin and very thin coal seam extraction in a long perspective of Ukraine?s fuel and energy industry development was determined. Technical and economic benefits of coal seams selective mining technique including beneficial environmental impact were substantiated. The scientific and practical aspects of development of borehole underground coal gasification (BUCG) technology are considered. The technological scheme of underground gas generators is offered on the basis of the authors? criteria of fitness and adaptability to the mine and geological conditions. These charts foresee the preparation of gas generators from earth surface and directly from the present mine working by inclined horizontal boreholes. Basic research is directed at abandoned supplies extraction in coal mines. Concepts of technological decisions and technical implementation of the underground part of the gasification system are worked out with analytical models, with subsequent test verification on special options and a mine experimental gas generator. Forming the technological scheme of gasification is founded on the basis of tying up spatial and sentinel factosr of coal seam gasification taking into account sufficient impermeability of the system a «gas generator ? rock massif» and to dirigibility by thermochemical processes that flow in a rock massif. The offered solutions contain the possibility of creation of the power and chemical enterprises of gasification, taking into account the complex use of BUCG products in the ecologically reserved technological cycle. Gas generator efficiency is evaluated from an economic standpoint.

Aug 1, 2013

By H. S. Fowler

The Consolidated Mining and Smelting Company has been in the business of mining and exploration since its incorporation in 1906. The Company has reports written as early as 1890. Consequently, its files contain a very large number of reports, and also of maps, of ail descriptions. REPORTS Before describing the filing of maps, a brie{ note on our general system of filing reports will perhaps clarify the discussion. Reports of a non-periodic nature, primarily those based upon examinations by field engineers, are separated from all other data and filed in what we know as the Mine Series System. Each property on which reports are written is given a serial number, known as the Mine Series Number. This serial number is assigned at the time the first report is filed and it remains permanently identified with the property, regardless of its state of activity. Our mine series numbers now total over 5,000. The serial number for a property may be found by referring to an alphabetical index of properties. MAPS We found it necessary to divide maps into two groups: (1) Maps of properties (2) Maps of areas

Jan 1, 1944

By S. Broda

Waste rock piles (WRP) are made of heterogeneous, coarse grained waste rock, with inclined layers due to segregation along the slopes and haul traffic compacted layers, creating highly contrasted hydrogeological conditions. The internal structure of the pile has a major impact on water and oxygen movement and pore water pressure distribution. In the presence of reactive minerals, this influences the production of acid mine drainage (AMD) or contaminated neutral drainage (CND) with metal leaching from the WRP. Water flow and distribution may also affect slope stability, particularly following major infiltration events (rain/snow melts).The prediction and interpretation of water distribution in WRP is a challenging problem. Numerous numerical investigations of short and long term forecasts have been conducted, and many of these were found of limited reliability. The shortcomings of existing approaches are typically due to difficulties in representing preferential flow processes in macropores. Such phenomena are in fact often neglected in the assessment of water flow in WRP because these are unattainable with the classical porous media approach typically used in this field. This paper focuses on the use of modeling techniques originally designed for massive fractured rock media to simulate flow in WRP. Waste rock properties used in this study are retrieved from large scale laboratory permeability and water retention tests and from in situ infiltration tests. The effect of WRP configurations on water flow into potentially reactive zones is examined with the numerical 3D fully-integrated surface/subsurface, variably saturated flow model HydroGeoSphere. Material properties retrieved from the laboratory tests are represented by means of stochastic facies distribution generated within a probability framework. Modelled soil water contents and hydraulic head distributions, as well as fluxes within a simple waste rock pile, with compacted layers, are compared with results retrieved from simulations based on the classical equivalent porous medium approach. Results indicate significantly increased flow towards the looser waste rock zone by considering heterogeneous material distribution, causing weak or absent capillary barriers and preferential flow paths. Practical implications for the design of waste rock piles are also discussed.

Aug 1, 2013

By H. Hao

This paper outlines the development of a 2-D axisymmetric mathematical model to predict the evolution of temperature during the Direct Chill (DC) casting of AZ31 magnesium billets. Model predictions of temperature have been compared to measurements made on an industrial DC casting machine at Timminco Ltd, in Haley, Ontario. These measurements have been used to fine-tune the thermal boundary conditions employed in the mathematical model, The model, which is based on the commercial finite-element code ABAQUS, has been developed as part of a strategic project aimed at improving the DC casting process for light metals through a reduction in casting related defects and increased productivity.

Jan 1, 2004

By G. M. Brownell

A DEPOSIT of hematite occurs on Black island, lake Winnipeg, at a point known locally as Red Rock, on the south shore of the island directly across from Clements point. The mail boat from the Canadian Pacific Railway terminus at Riverton, crossing the lake to Manigotogan twice weekly, passes near Red Rock and provides a convenient means of transportation to this part of Black island (Figure 1). HISTORICAL The occurrence of hematite on Black island has been known for a long time. Dowling's report of 1898 states(1): "An iron-ore deposit on the south side has been known' for many years, and on Jeffrey's map of 1762, showing Canada and northern part of Louisiana, Iron island is shown a short distance south of the narrows, evidently the same as the Black island of today".

Jan 1, 1945

By Edgar Dernedde

Cryolitic melts at a temperature of about 960 °C are used in reduction cells to dissolve alumina and to produce primary aluminium. At these temperatures cryolitic melts with an acidic composition have a significant vapour pressure, which contributes substantially to the evolution of fluorides from these cells, The Antoine equation is commonly used for relating vapour pressure data of liquids at different temperatures, It is also applicable to cryolitic melts, Published equations have at least 20 numerical constants, in order to allow for the effects of melt composition. In the present work a simplified equation with only 12 numerical constants was developed and is applicable for multiple additions to molten cryolite up to 30 wt %, The simplified equation was compared with the published equations of Haupin & Kvande and Gustaysen & Ostvold using 472 vapour pressure data of cryolitic melts. Results showed that the simplified equation has the same accuracy as the published equations in the composition range of interest to industry.

Jan 1, 2004

By Edward (Ted) R. Lord

The application of thickened/paste tailings for surface disposal has been practiced for many years but has recently started to gain wider international acceptance. With tailings ponds coming under increasing environmental scrutiny, there is growing momentum towards mine/tailings plans that greatly reduce the size or even eliminate tailings ponds completely. This paper provides examples from the author's own experiences of the application of thickened/paste tailings for surface disposal in different countries, under different climatic conditions and in different industries. The paper will also comment on current research initiatives and future trends.

May 1, 2001

By Claude Bazin

The Quebec province in Canada is a major producer of aluminium through electrolysis of alumina. The alumina that feeds the electrolysis cells is mainly imported from foreign countries, where it is extracted from bauxite. Other alumina sources, such as clays, have sometimes been considered as potential substitutes for bauxite. There is a large clay deposit in the area of Murdochville in Quebec, with sulfuric acid and electricity available nearby the orebody, as well as sea port facilities leading to the major aluminium plants in Quebec. These are favorable factors to assess the economics of an industrial plant to produce alumina from clays. This paper reviews some of the processes considered for the extraction of alumina from clays and presents results of the initial work carried out on the clay available from Murdochville.

May 1, 2001

By Nirmal Weerasekara

Mill liners serve the dual purpose of protecting the shell from wear and transferring the rotary motion of the mill to an effective grinding action. As the liners wear, the efficiency of energy transfer changes, so there is balance between liner profile for wear and for milling efficiency. A tool that can simultaneously predict the evolution of the liner profile and the associated change in milling efficiency should prove most useful in optimal liner design. This paper presents such an approach, based on measured 3D liner profile data informing 3D DEM simulations of charge motion and liner wear modelling.

Jan 1, 2012

By Y. W. Pan, Y. R. Liu

The disturbed natural rock mass tends to equilibrium state accompanying with non-equilibrium evolution process which has significant influence on the long-term stability of engineering structure. The deformation development in the non-equilibrium evolution is determined by two intrinsic processes, trend to equilibrium and damage development, which can be illustrated by evolution of excavated tunnel. Thus the creep and damage must be taken into account in study of non-equilibrium evolution process. Within Rice’s irreversible thermodynamics with internal state variable (ISV), the conjugate force function which equals to zero is a convex surface which divide the stress space into two parts, called equilibrium and nonequilibrium stress space respectively. The concept of effective stress of non-equilibrium has been proposed, in which the stress difference (d - d ) are the effective driving forces of non-equilibrium evolution process with hardening and damage effective. The formulations of inelastic strain rate and energy dissipation rate have been recast in the form of effective stresses. The inelastic complementary energy (ICE) is proposed to measure intensity of non-equilibrium evolution of material system. The specific equations of conjugated forces and kinetic equations of internal state variables are presented. The ICE and rate of energy dissipation are all calculated to be used to evaluate long-term stability of geo-engineering. The results show that rate of energy dissipation may be more convenience and flexible than ICE for internal state variable constitutive equation of rock mass.

Jan 1, 2015