Search Documents

By P. van der Heever

The genesis of Simrac?s technology transfer practices and future initiatives are discussed in this report. Simrac produces a relatively large volume of research-based health and safety outputs that must be communicated to its tripartite stakeholders and to other interested and affected parties. The wide spectrum and complexity of the transferables, together with the diversity of end-user needs and skills of the recipient population, necessitate customizing of technology transfer practices. In an attempt to optimize this process, a variety of practices are employed. Information is transferred by means of publication and dissemination of a variety of reports (paper, website and CD).Knowledge transfer is directed to health and safety practitioners by means of workshops, seminars, best practice booklets, handbooks and textbooks. Promotion of recently completed projects is directed at targeted recipients by means of project launches. Only isolated cases of piloting and implementation of technologies have been undertaken due to lack of resources and past organizational focus. It is anticipated that Simrac?s new organizational structure and management system, together with end-user involvement by industry champions, from project inception to implementation, will enhance the effectiveness of the technology transfer process.

Jan 1, 2002

By A. M. Hearn, J. G. Roos

Energy consumption, being a large proportion of the input of the ferroalloy production process, can be measured in terms of the efficiency with which the energy is utilized and the effectiveness with which it is used. The efficiency of the usage of energy revolves around the specific consumption of the energy for the production of the alloy, which encompasses the thermodynamic aspects of the process and is more of a metallurgical function that will not be addressed further in this paper. The effectiveness with which the energy is utilized is determined by the optimisation of the electricity supply from the utility (in this case Eskom). The South African electricity supply industry is on the brink of considerable restructuring, implying changes in ownership, structure and regulation. In other countries, however, restructuring has often seen higher prices. Owing to the diversity of growth prospects in South African markets, Eskom is turning to new market management and generation strategies. In future, Eskom will develop pricing strategies to increase demand from attractive customer segments such as paper & pulp, iron & steel, aluminium, and other metallurgical industries. There will be a broader spectrum of product offerings, including dynamic pricing options and demand-side management incentives. In a future electricity market that will become more dynamic and demanding, it will become increasingly important for electricity consumers to optimise their electricity supply. There are many facets in optimising the electricity supply, which include optimal utilisation of the tariff, integrating energy cost management into the production planning process, implementation of demand-side management actions, modelling, scenario planning (.what-ifs?.), optimisation, maintenance scheduling coordinated with the tariff periods, etc. This paper gives consideration to the status of the development of technology at the manganese alloy works of Samancor at Meyerton, South Africa in order to increase the effectiveness with which energy is utilized.

Jan 1, 2004

By H. Willenberg, J. S. Coggan

ABSTRACT: Although individually both continuum and discontinuum numerical analyses provide useful means to analyze rock slope stability problems, complex failures typically involve mechanisms related to both existing discontinuities (i.e. discontinuum) and deformation and brittle fracture of intact rock (i.e. continuum). One such example is the 1991 Randa rockslide in the southwestern Swiss Alps, where failure involved several complex mechanisms, which require the consideration of progressive failure and rock mass strength degradation. This paper presents results from a series of hybrid finite-/discrete-element models directed towards the explicit modelling of brittle fracturing in natural rock slopes, enabling the modelling of progressive failure, from initiation, through to the development of catastrophic failure, transport and deposition. RÉSUMÉ: Si les modélisations numériques de type continu ou discontinu fournissent des outils utiles pour l’analyse de problèmes de stabilité de pentes rocheuses, elles montrent de sérieuses limitations lorsque des mécanismes de rupture complexes impliquent des processus liés autant aux discontinuités préexistantes qu’au développement progressif de nouvelles fractures. Ceci se retrouve typiquement dans l’exemple de l’éboulement de 1991 à Randa dans les Alpes suisses du sud-ouest, où la rupture a impliqué plusieurs mécanismes complexes nécessitant la prise en considération de la fracturation progressive ainsi que de la dégradation de la résistance des masses rocheuses. Cet article présente des résultats provisoires d'une série de modèles hybrides orientés vers la modélisation explicite de la fracturation de pentes naturelles, permettant la modélisation de fractures progressive de leur initiation jusqu’à leur développement, leur rupture catastrophique, leur transport et leur déposition. ZUSAMMENFASSUNG: Obwohl numerische Kontinuum- und Diskontinuum-Berechnungen einzeln nützliche Ergebnisse zur Stabilitätsanalyse von Felshängen liefern, sind komplexe Felsrutschungen typischerweise mit Mechanismen verbunden, die sowohl durch bestehende Trennflächen (s.a. Diskontinuum-Ansatz) als auch durch die Deformation und das Versagen intakter Felsbereiche (s.a. Kontinuum-Ansatz) bedingt werden. Ein Beispiel für einen solchen Sachverhalt ist der Bergsturz von Randa 1991 im südwestlichen Teil der Schweizer Alpen, bei dem die Rutschung mehrere komplexe Mechanismen umfasste, so dass “progressive failure” und Verminderung der Gebirgsfestigkeit in Betracht gezogen werden müssen. Dieser Artikel beschreibt vorläufige Ergebnisse einer Reihe von hybriden Modellrechnungen (Finite-/Diskrete-Elemente). Diese haben zum Ziel, explizit das bruchhafte Versagen natürlicher Felshänge zu modellieren, wobei das progressives Versagen im ganzen modelliert wird, von der Anfangsphase über die Entwicklung zum katastrophalen Versagen und den Transport und die Ablagerung der Rutschmassen.

Jan 1, 2003

By M. Castelli, C. Scavia, S. Marello, A. Allodi

ABSTRACT – This paper investigates the possibility of interpreting the shear band evolution in hard soils and soft rocks as the result of shear propagation from pre-existing natural defects. Thesis done through the application of the principles of Fracture Mechanics, a Slip-Weakening Model(SWM) being used to simulate the non-linear zone at the tips of the discontinuity. Experimental observations have been carried out on prismatic samples of Beaucaire marl through biaxial compression tests in plane strain conditions. A numerical simulation confirms the validity of the proposed approach. RÉSUMÉ – Cet article présente la possibilité d’interpréter l’évolution des bandes de cisaillement dans les roches tendres et les sols raides comme le résultat de la propagation pour cisaillement à partir d’un défaut préexistant. Le Slip-Weakening Model (SWM) a été employé pour modéliser la zone non linéaire a l’extrémité d’une discontinuité. Nous présentons des observations expérimentaux obtenues à travers des essais de compression en déformation plane sur des échantillons prismatiques de marne. La simulation numérique d’un essai confirme la validité de la méthode proposée. KURZFASSUNG – Die vorliegende Arbeit versucht die Ausbreitung von Scherzonen in(bindigen) Böden und im Lockergestein anhand der Entwicklung von Scherrissen ausgehend vonexistierenden Schwachstellen zu erklären. Die Umsetzung erfolgte mit Hilfe der Bruchmechanik,für die Modellierung der nichtlinearen Zone an den Rissenden wurde ein Slip-Weakening-Modell(SWM) eingesetzt. Es wurden zweiaxiale Druckversuche an prismatischen Probekörpern vonBeaucaire Mergel unter der Bedingung eines ebenen Verzerrungszustandes durchgeführt. Dienumerische Untersuchungen bestätigten die Annahmen der untersuchten Theorie.

Jan 1, 2003

By S. M. Spottiswoode

Much as one hopes to predict whether damaging seismicity might occur somewhere in the next shift, this is not possible at present. Recently, an analysis of seismicity at two mines has shown that a small increase in seismicity occasionally occurs before both large and small events. This analysis is extended here to five mines, with the same results. Earthquake seismologists forecast hazard in terms of the average incidence of earthquakes over the last 50 years, adjusted on an hourly basis according to the possibility of aftershocks of events that have happened. The rate of seismicity in any area remains approximately constant when expressed as seismicity per area mined, or better still, per unit of elastic strain energy released. It is suggested that seismic hazard in mines be estimated using historical seismicity and forward modelling, adjusted according to the possibility of aftershocks. In summary, medium-term forecasting of seismicity is viable; short-term prediction of large events is not.

Jan 1, 2010

By O. Rondon

Localized uniform conditioning is a technique developed to spatially locate selective mining unit grades that have been derived using uniform conditioning for the assessment of recoverable resources. The technique has the advantage of producing selective mining unit estimates conforming to the uniform conditioning panel-specific grade-tonnage curve while introducing spatial information at the scale of the selective mining units. This paper describes an alternative technique to localized uniform conditioning which does not explicitly require the uniform conditioning panel-specific grade-tonnage curve to localize the selective mining unit estimates. The technique can therefore be implemented in mining software where uniform conditioning is not available.

Jan 1, 2014

By J. N. Van der Merwe

Due to the approaching depletion of reserves in the Witbank Coalfield, and the fact that the Waterberg Coalfield is as yet under developed, there is increasing pressure to extend the lives of the operating coal mines in the Witbank Coalfield. One potential source of coal is the old pillars, often small and at shallow depth, that were not left with secondary extraction in mind. At face value, some of these pillars appear to be suitable for secondary mining due to their high safety factors. However, it is known that over time, these pillars have scaled and the current sizes are smaller than the as-mined dimensions. It is also known that at shallow depth, the overburden is often less likely to fail during secondary mining, resulting in high abutment loads on the unmined pillars. The paper proposes a systematic method to pre-evaluate those old pillars for the possibility of stooping. The method consists of elements of fundamental methods and newly developed technology. In essence, it revolves around using empirical methods to estimate current pillar dimensions, followed by numerical modelling to investigate the probability of progressive pillar failure and then fundamental methods to determine the limits of applicability of the numerical model.

Jan 1, 2012

By H. Möhle

Premier Mine is a large underground diamond mine, situated at Cullinan in the Gauteng Province of South Africa. The mine hoists over 3 millions tons of diamond bearing kimberlite per year, from a modern, mechanized panel retreat block cave mining operation. Troughs sited at the base of the cave are developed into the overlying caved ore, allowing loading access to the ore through draw points. Trough opening is therefore a vital aspect of the mining operation. Most of the production is dependent on fast, reliable and safe trough opening. Successful trough openings are achieved under difficult mining conditions due to restricted void, long up-holes and the inherent properties of kimberlite?soft, elastic rock with a high swell factor. Early in 1998, Premier Mine undertook a joint development project (with an electronic initiation system supplier) to determine whether the use of a fully programmable electronic detonator system would enhance the mining cycle in the trough opening process. Blasting commenced in February 1998, initially with blasts of a similar size and duration to those that had been previously blasted with shock-tube systems, to confirm the operational capabilities of the blasting system. The accuracy, reliability and flexibility of the electronic delay detonators (EDDs) minimized associated risks for re-engineering of the trough opening process. Further, as blast timing was optimized, larger blasts were successfully initiated. To date over 20 troughs have been opened successfully with the use of EDDs. Not only did it take less time to open troughs, but also drilling and blasting operations have been optimized, as well as an overall improvement in safety. This paper details the benefits that EDDs have delivered to Premier Mine in this area. The logic followed is: ?Review of the original trough opening practices ?Re-engineering of the trough opening process ?Implementation of EDDs ?Evaluation of improvements and conclusions.

Jan 1, 2000

By M. B. Shongwe

Pt-based alloys are being developed for high-temperature applications with the aim of replacing some of the currently used Ni-based superalloys (NBSAs) in the highest temperature applications. The Pt-based alloys have a similar structure to the NBSAs, and since Pt is more chemically inert than nickel and has a higher melting point, they can potentially be used at higher temperatures, up to 1300°C, and in more aggressive environments. Several experimental Pt-based alloys were studied at Mintek, and an optimum composition was found to be Pt84:Ahi:Ru2:Cr3 (at.%). On the basis of hardness and microstructure, a later study identified a new optimum: Ptr8:A111:Ru5:Cr6 (at.%). There are at least two Pt3A1 allotropes, and the high-temperature cubic structure has better properties than the lower temperature tetragonal form, and so needs to be stabilized. This work describes the latest results obtained in transmission electron microscopy (TEM) studies of the quaternary Pt-based superalloys. These results are both interesting and important, because the samples have a higher precipitate density compared to those from earlier work. The precipitate morphology is mainly cubic, with minor areas having irregular-shaped precipitates. The high volume fraction of the precipitates is a major breakthrough, since the objective of this work is to improve the alloys. A prior disadvantage was that the proportion of the precipitates was too low in samples before this work, especially compared with the work from Germany on Pt-Al-Cr-Ni¬-based alloys as well as the NBSAs. TEM -Pt3A1 diffraction patterns were studied, and for each diffraction pattern, many possible lattice point combinations were tried, with the angle between the lattice spots as well as interplanar spacings being calculated for each phase (cubic or tetragonal) to match the measured results. An overall analysis of the diffraction results indicates that the cubic phase fitted the experimental lattice points with much lower errors compared to the tetragonal phase. Thus, with the close match achieved with the cubic structure, the structure of -Pt3A1 precipitates is likely to be cubic. X-ray diffraction has been carried out on selected samples, and the results confirmed the presence of cubic -Pt3A1 and (Pt). Different X¬ray diffractometers were used to further verify the results, and the results were also compared with those from TEM.

Jan 1, 2012

By C. Scheffer

This paper describes a model that can simulate the behaviour of a pneumatic intensifier that is generally referred to as an ?air pump? (AP). The AP is used in many South African underground mines, and utilizes compressed air with pressures between 0.3?0.6 MPa to pump water up to pressures of 12?18 MPa. The high pressure water is used in numerous underground mines to pressurize props and roof bolts that act as safety and roof support devices. The simulation model reproduces the intensifier's operation under various conditions and therefore can be used to investigate the AP's sensitivity towards environmental conditions, design parameters (such as spring stiffness, plunger diameter, nozzle diameter etc.) and operating conditions. Model validation by comparison with laboratory experiments showed good correlation. The model can subsequently be used to investigate design changes aimed at improving the AP's performance and reliability over a wider range of operating conditions.

Jan 1, 2005

The gold concentration profiles which develop inside the activated carbon particles during the adsorption process were measured, so as to obtain additional direct proof of intra-particle mass transfer rate-control and to determine how individual carbon particles in a population vary in their gold uptake behaviour. The observed gold concentration profiles inside the activated carbon particles were typical of situations where an intra-particle mass transfer mechanism is rate-controlling. Gold di-cyanide penetrated the carbon particles with a wave-like front until it was uniformly distributed throughout the particle. This corresponds with the profiles predicted by the film and surface diffusion model. However, the rate with which gold di-cyanide was transported inside the individual carbon particles varied within the population of carbon particles. This variation was characterized by a lognormal distribution in surface diffusivity. Because the estimated variances of these surface diffusivity distributions were small, the dynamic adsorption behaviour of the carbons used in the present study did not deviate significantly from that predicted by a model based on film mass transfer and a single surface diffusivity. Van Deventer (1984, 1985 and 1986a) reported data that could not be adequately described by a model based on film mass transfer and a single surface diffusivity. In this study a model based on film mass transfer and a wide lognormal distribution of surface diffusivity values was applied to his data. It was shown that a population of carbon particles with a wide distribution in surface diffusivity values could offer an alternative explanation to the mechanism proposed by Van Deventer which was based on film mass transfer, with partitioning of the carbon into micro- and macropores with surface diffusion in the macropores and slow mass transfer into the micropores.

Jan 1, 1997

By M. W. Hildyard

A blast was engineered close to a tunnel in a deep-level mine with the purpose of studying the wave interaction with the tunnel. Numerical modelling of seismic wave propagation was used in both the forward planning and the back-analysis of this experiment. This paper presents the forward analysis. A likely model for a propagating blast was investigated and developed. This model predicts the influence of the velocity of detonation of the blast and the rise-time of the blast on the radiation pattern, and hence on the expected positions of maximum velocity at the tunnel surface. An attempt was made to quantify the parameters of the blast model, using waveform data from two other blast experiments. The blast model caused particle velocities which were an order of magnitude lower than these measured velocities. Large values for blast pressure were required to match the amplitudes of the recorded velocities. An alternative ?effective? source with a similar mechanism but increased diameter was proposed to account for the higher velocities. It was established that distant velocities scale with the square of the borehole diameter. A calibration blast at the experimental site was also modelled. Agreement in the waveforms was not sufficient for conclusive statements to be made about the validity of the source model. However, the modelled waveforms show an encouraging correspondence with some measured waveforms. This gave confidence that the model could be used to examine wave propagation around the tunnel, and allowed a scaled source to be developed for forward analysis of the main experiment.

Jan 1, 2001

By O. Celep, T. Uslu, Savas

This study investigates the enrichment of a low-grade colemanite concentrate (-3 mm) using a Knelson centrifugal gravity concentrator. Due to its low boron content, the concentrate is unsaleable and has to be stored under appropriate conditions to avoid potential environmental problems. The low-grade colemanite concentrate was comminuted to size fractions of-1 mm, -0.5 mm, and -0.15 mm before treatment in the Knelson Concentrator. The effects of particle size, fluidizing water velocity, and bowl speed on the enrichment process were examined. The B2O3content of the concentrate was increased from 33.96% to a maximum of 45.52%. B2O3recovery increased with increasing bowl speed and particle size, and decreased with increasing fluidizing water velocity. The enrichment process also rejected arsenic and iron to some extent, with a maximum reduction of arsenic from 1360 g/t to 765 g/t and iron from 0.88% to0.33%.

Jan 1, 2015

By Frik Badenhorst, René Roux, Peter Terbrugge

This paper presents a case study where economics dictated an aggressive slope angle and where an appropriate slope management programme allowed successful completion of mining to 190m below surface. The Navachab Gold Mine is located in Namibia and is one of the few gold mines in that country. Mining by open pit methods commenced in 1990 with economics dictating an aggressive slope design to minimise waste stripping and make the mine viable. With steep slope angles in the design, and a steeply dipping footwall geology, it was necessary to introduce a comprehensive slope management programme. This required good mining practice in limit blast techniques with pre-splits and trims, and good housekeeping to effect a safe and economical operation. The slope management programme in addition required in pit mapping and logging procedures, regular survey monitoring evolving to a robotic Geomos system, as well as slope instrumentation to monitor stability. Two relatively large scale failures occurred at the mine, the second of which impacted on the east wall access ramp, while the first was confined to an operating stack. The paper describes the design and slope management techniques developed for the east wall.

Jan 1, 2006

Rock engineering has now been practised as a discipline for about50 years, and many sophisticated methods of analysis for use in rock engineering design have been developed. However, analysis is only a tool for use in engineering design, and rock masses are often so complex that realistic modelling of them, even with sophisticated methods of analysis, is impossible. At the other end of the scale, rock mass classification methods have been developed to quantify the rock mass. These ?numbers? are now widely used for assessment of stability of underground and surface excavations in rock, for the determination of rock support requirements, and for the estimation of rock mass strength and deformation properties. In considering these areas, the two concepts?design, and prediction of behaviour?will be dealt with. The discussion will be introduced by a consideration of the characteristics of the rock mass, which will highlight the variability in its properties and hence the uncertainty in the behavioural parameters. Examples of real projects will be used to illustrate the discussion.

Jan 1, 2003

By Merete Tangstad, Kai Tang

The primary aim of the present investigation is to develop a model which permits to calculate/predict the viscosities of ferromanganese slags, typically focus on the MnO-SiO2-Al2O3-CaO-MgO-FeO multicomponent system. The present viscosity model emphasizes the ability of prediction, i.e. using the binary parameters to extrapolate the viscosities in ternary and higher order systems. To this end, a physical sound model is essential. Temperature dependence of viscosity is described by the Eyring?s formalism. The well-known cell model is used to define the SiO44- tetrahedra in the silicate melts. Composition dependence is expressed in term of viscous activation energy, which is characterized by the binary structural units. The viscosities in ternary and higher order systems are hence predicted using the binary activation parameters. The model has been used to calculate the viscosities in the whole liquid composition range of ferromanganese slag.

Jan 1, 2007

By R. Trammell, O. Eterigho-Ikelegbe, S. O. Bada

In response to the enormous amounts of coal discard generated during coal mining and preparation, the development of an efficient and sustainable strategic use of this resource is essential. Furthermore, the rising urban population over the next decades is confronted with the depletion of quality raw materials for building components. To this end, this study reports new information on the morphology, water absorption, and flexural strength properties of ceramic composites produced from three different coal discards and polysiloxane pre-ceramic polymer (PCP) resin. In addition, test results relating to the continuous operating temperature, chemical resistance, and efflorescence potential of the composites are presented. The results show that the water absorption and flexural strength of the coal composites, up to 1.94% and 36.46 MPa respectively, exceed the requirements for ceramic and clay roof tiles. The continuous operating temperature of composites is found to be more thermally stable than conventional roofing tiles (concrete and ceramic) between ambient temperature and 600°C. In addition, the excellent chemical resistance of the composites (94.43%−99.98%) compared to conventional roofing tiles (67.82%−99.97%) eliminates the need for additional external coatings. The interesting results documented so far suggest that this technique could be used to produce low-temperature application building products such as bricks, panels, roofing tiles, etc. This new recycling technique offers an excellent opportunity to eliminate enormous volumes of coal discard and to advance the circular economy in the coal industry.

Aug 2, 2022

By K. Lombard

"‘Roof support awaiting time’ (RSAT) is a term used at Goedehoop Colliery’s Simunye Shaft to describe the potential production time lost due to the continuous miner (CM) standing idle waiting for roof support to catch up. Investigations revealed that in 2013, Simunye Shaft had approximately1400 hours of RSAT, which suggests that the mine could have potentially produced an additional 280 000 t of coal. This project consisted of two parts. Firstly, the causes of the high RSAT and means to improve the situation were investigated. Secondly, as insisted by mine management, the CM cutting sequence was investigated as a possible cause of high RSAT. Machine-related challenges due to the roofbolter installing support too slowly, geological conditions (mostly hard roof conditions and slips),logistical challenges pertaining to the CM cutting sequence, man-related challenges related to operator fatigue, re-support, operator inexperience, and the absence of support targets were identified as main contributors to RSAT. Furthermore, results showed that the roofbolters in the sections at Simunye Shaft are slower than the CMs. A target of 28% reduction in RSAT was set. Experts from Kennametal and Fletcher were consulted to find solutions for the identified causes. In total, eight solutions for RSAT were identified, but the solution that contributed most significantly to reducing RSAT was to use hard roof drill bits as a standard product at Simunye Shaft. Calculations showed that by using hard roof drill bits, RSAT can be reduced by 43%, which is more than the specified 28% target. The cutting sequences at Kriel, Greenside, and Simunye Shaft, togetherwith three newly developed cutting sequences, were simulated using the UCMS (Underground Coal Mining Simulation) program. A re-aligningprinciple was incorporated into the newly developed cutting sequences toalign the cutting sequences to general support work and to reduce RSAT. A decision matrix revealed that a cutting sequence in which boxing takesplace in R3 (third road to the right of the belt road) and in which the re-aligning principle has been incorporated will be the best option for Simunye Shaft. The recommended cutting sequence will lead to a 5% increase in production.The work presented in this paper was carried out as partial fullfilment for the degree BEng (Mining Engineering)"

Jan 1, 2015

The influence of varying amounts of magnesia, silica, and alumina in iron ore sinter on its mineralogy, reducibility index (RI), reduction disintegration index (RDI), physical breakdown (AI and TI), and production properties (coke breeze rate) were examined. It was found that the mineralogy of iron sinter can more easily be predicted from its chemical composition than from the RI, RDI, AI or TI. Anticipating the consequence that varying amounts of MgO and SiO2would have on sinter properties is complex, and not necessarily predictable. High concentrations of Al2O3in the sinter result in high concentrations of the SFCA phase, but with drastically deteriorated properties. This study also confirmed that the form in which fluxes are added to the raw material sinter mixture affects the mineralogy and properties of the produced sinter.

Jan 1, 2010

By T. Nishimura, H. Kiyama, K. Wakimoto

Discussion on the load bearing capacity of ground and the stability of tunnel is presented. We derive a simplified solution of a circular lined-tunnel in an elastic ground. The solution results in a better understanding of characteristic curves. Use of DEM allows us to consider effects of possible yield of the overstressed zone around the tunnel and inherent discontinuity of ground on the curves. Diskussion auf der Last Tragfähigkeit Bodens und die Stabilität des Tunnels wird dargestellt. Wir leiten eine vereinfachte Lösung eines kreisförmigen Ausbautunnels in einem elastischen Boden ab. Die Lösung ergibt ein besseres Verständnis der Kennlinien. Um Gebrauch DEM erlaubt uns, zu betrachten Effekte des möglichen Ergebnisses der überbetonten Zone den Tunnel und zugehörige Unstimmigkeit des Bodens auf den Kurven. La discussion sur la portance de charge de la terre et de la stabilité du tunnel est présentée. Nous dérivons une solution simplifiée d'un revêtement-tunnel circulaire dans une terre élastique. La solution a comme conséquence une meilleure compréhension des courbes caractéristiques. L'utilisation de DEM nous permet de considérer des effets du rendement possible de la zone surchargée autour du tunnel et de la discontinuité inhérente de la terre sur les courbes.

Jan 1, 2003