A numerical study of flexural buckling of foliated rock slopes

The occurrence of foliated rock masses is common in mining environment. Methods employing continuum approximation in describing the deformation of such rock masses possess a clear advantage over methods where each rock layer and each inter-layer interface (joint) is explicitly modelled. In devising such a continuum model it is imperative that moment (couple) stresses and internal rotations associated with the bending of the rock layers be properly incorporated in the model formulation. Such an approach will lead to a Cosserat-type theory. In the present model, the behaviour of the intact rock layer is assumed to be linearly elastic and the joints are assumed to be elastic-perfectly plastic. Condition of slip at the interfaces are determined by a Mohr-Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformations. The model is incorporated into the finite element program AFENA and validated against an analytical solution of elementary buckling problems of a layered medium under gravity loading. A design chart suitable for assessing the stability of slopes in foliated rock masses against flexural buckling failure has been developed. The design chart is easy to use and provides a quick estimate of critical loading factors for slopes in foliated rock masses. It is shown that the model based on Euler's buckling theory as proposed by Cavers (Rock Mechanics and Rock Engineering 1981; 14:87-104) substantially overestimates the critical heights for a vertical slope and underestimates the same for sub-vertical slopes. Copyright (C) 2001 John Wiley & Sons, Ltd.

Finite element modelling of rock alteration and metamorphic process in hydrothermal systems

We use the finite element method to simulate the rock alteration and metamorphic process in hydrothermal systems. In particular, we consider the fluid-rock interaction problems in pore-fluid saturated porous rocks. Since the fluid rock interaction takes place at the contact interface between the pore-fluid and solid minerals, it is governed by the chemical reaction which usually takes place very slowly at this contact interface, from the geochemical point of view. Due to the relative slowness of the rate of the chemical reaction to the velocity of the pore-fluid flow in the hydrothermal system to be considered, there exists a retardation zone, in which the conventional static theory in geochemistry does not hold true. Since this issue is often overlooked by some purely numerical modellers, it is emphasized in this paper. The related results from a typical rock alteration and metamorphic problem in a hydrothermal system have shown not only the detailed rock alteration and metamorphic process, but also the size of the retardation zone in the hydrothermal system. Copyright (C) 2001 John Wiley & Sons, Ltd.

A director theory for visco-elastic folding instabilities in multilayered rock

A model for finely layered visco-elastic rock proposed by us in previous papers is revisited and generalized to include couple stresses. We begin with an outline of the governing equations for the standard continuum case and apply a computational simulation scheme suitable for problems involving very large deformations. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered beam under compression. We analyse folding up to 40% shortening. The standard continuum solution becomes unstable for extreme values of the shear/normal viscosity ratio. The instability is a consequence of the neglect of the bending stiffness/viscosity in the standard continuum model. We suggest considering these effects within the framework of a couple stress theory. Couple stress theories involve second order spatial derivatives of the velocities/displacements in the virtual work principle. To avoid C-1 continuity in the finite element formulation we introduce the spin of the cross sections of the individual layers as an independent variable and enforce equality to the spin of the unit normal vector to the layers (-the director of the layer system-) by means of a penalty method. We illustrate the convergence of the penalty method by means of numerical solutions of simple shears of an infinite layer for increasing values of the penalty parameter. For the shear problem we present solutions assuming that the internal layering is oriented orthogonal to the surfaces of the shear layer initially. For high values of the ratio of the normal-to the shear viscosity the deformation concentrates in thin bands around to the layer surfaces. The effect of couple stresses on the evolution of folds in layered structures is also investigated. (C) 2002 Elsevier Science Ltd. All rights reserved.

Large amplitude folding in finely layered viscoelastic rock structures

We analyze folding phenomena in finely layered viscoelastic rock. Fine is meant in the sense that the thickness of each layer is considerably smaller than characteristic structural dimensions. For this purpose we derive constitutive relations and apply a computational simulation scheme (a finite-element based particle advection scheme; see MORESI et al., 2001) suitable for problems involving very large deformations of layered viscous and viscoelastic rocks. An algorithm for the time integration of the governing equations as well as details of the finite-element implementation is also given. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered plate. The domain is compressed along the layer axis by prescribing velocities along the sides. First, for the viscous limit we consider the response to a series of harmonic perturbations of the director orientation. The Fourier spectra of the initial folding velocity are compared for different viscosity ratios. Turning to the nonlinear regime we analyze viscoelastic folding histories up to 40% shortening. The effect of layering manifests itself in that appreciable buckling instabilities are obtained at much lower viscosity ratios (1:10) as is required for the buckling of isotropic plates (1:500). The wavelength induced by the initial harmonic perturbation of the director orientation seems to be persistent. In the section of the parameter space considered here elasticity seems to delay or inhibit the occurrence of a second, larger wavelength. Finally, in a linear instability analysis we undertake a brief excursion into the potential role of couple stresses on the folding process. The linear instability analysis also provides insight into the expected modes of deformation at the onset of instability, and the different regimes of behavior one might expect to observe.

Environmental influence on the stress corrosion cracking of rock bolts

In order to understand rock bolt Stress Corrosion Cracking (SCC), a series of experiments have been performed in Linearly Increasing Stress Test (LIST) apparatus. One series of experiments determined the threshold stress of various bolt metallurgies (900 MPa for Steel A, and 800 MPa for Steel B and C). The high values of threshold stress suggest that SCC begins in rock bolts when they are sheared by moving rock strata. Typical crack velocity values have been measured to be 2.5 x 10(-8) m s(-1), indicating that there is not much benefit for rock bolt steel of higher fracture toughness. Another series of experiments were performed to understand the environmental conditions causing SCC of steel A and galvanised Steel A rock bolt steel. SCC only occurred for environmental conditions for which produce hydrogen on the sample surface, leading to hydrogen embrittlement and SCC. Fracture surfaces of LIST samples failed by SCC were found to display the same fracture regions as fracture surfaces of rock bolts failed in service by SCC: Tearing Topography Surface (TTS), Corrugated Irregular Surface (CIS), quasi Micro Void Coalescence (qMVC) and Fast Fracture Surface (FFS). Water chemistry analysis were carried out on samples collected from various Australian mines in order to compare laboratory electrolyte conditions to those found in underground mines.

Material influence on the stress corrosion cracking of rock bolts

Rock bolt stress corrosion cracking (SCC) has been investigated using the linearly increasing stress test (LIST). One series of experiments determined the threshold stress of various bolt metallurgies (900 MPa for 1355AXRC, and 800 MPa for MAC and MA840B steels). The high values of threshold stress suggest that SCC begins in rock bolts when they are sheared by moving rock strata. SCC only occurred for environmental conditions which produce hydrogen on the sample surface, leading to hydrogen embrittlement and SCC. Different threshold potentials were determined for a range of metallurgies. Cold work was shown to increase the resistance of the steel to SCC. Rock bolt rib geometry does not have a direct impact on the SCC resistance properties of the bolt, although the process by which the ribs are produced can introduce tensile stresses into the bolt which lower its resistance to SCC. (C) 2004 Elsevier Ltd. All rights reserved.

Numerical simulation of double-diffusion driven convective flow and rock alteration in three-dimensional fluid-saturated geological fault zones

Numerical methods are used to simulate the double-diffusion driven convective pore-fluid flow and rock alteration in three-dimensional fluid-saturated geological fault zones. The double diffusion is caused by a combination of both the positive upward temperature gradient and the positive downward salinity concentration gradient within a three-dimensional fluid-saturated geological fault zone, which is assumed to be more permeable than its surrounding rocks. In order to ensure the physical meaningfulness of the obtained numerical solutions, the numerical method used in this study is validated by a benchmark problem, for which the analytical solution to the critical Rayleigh number of the system is available. The theoretical value of the critical Rayleigh number of a three-dimensional fluid-saturated geological fault zone system can be used to judge whether or not the double-diffusion driven convective pore-fluid flow can take place within the system. After the possibility of triggering the double-diffusion driven convective pore-fluid flow is theoretically validated for the numerical model of a three-dimensional fluid-saturated geological fault zone system, the corresponding numerical solutions for the convective flow and temperature are directly coupled with a geochemical system. Through the numerical simulation of the coupled system between the convective fluid flow, heat transfer, mass transport and chemical reactions, we have investigated the effect of the double-diffusion driven convective pore-fluid flow on the rock alteration, which is the direct consequence of mineral redistribution due to its dissolution, transportation and precipitation, within the three-dimensional fluid-saturated geological fault zone system. (c) 2005 Elsevier B.V. All rights reserved.

Prevalence and correlates of the injection of methadone syrup in Sydney, Australia

A sample of 312 heroin users was interviewed on their injection of methadone syrup. Methadone injecting was widespread, with 52% of subjects having injected methadone syrup, 29% in the preceding six months. Males and females were equally likely to report methadone injecting. Forty per cent Of current methadone injectors reported weekly or more frequent methadone injecting over the preceding six months. A history of methadone injecting was;associated with abscesses and infections in injection sites, having been diagnosed with a venous thrombosis and a history of heroin overdose. Current methadone injectors were in poorer general health, had more injection-related symptoms, higher levels of psychological distress, were more likely to have recently passed on used injecting equipment and to have recently committed criminal acts. Implications for the reduction in the prevalence of methadone injecting and associated harm are discussed.

Transitions between routes of benzodiazepine administration among heroin users in Sydney

A sample of 312 heroin users were interviewed regarding their benzodiazepine use. The majority (94%) had used benzodiazepines, 72% in the 6 months prior to interview. Benzodiazepine injecting was common, with 28% of the sample having injected these drugs, 13% in the 6 months preceding interview. Current benzodiazepine injectors showed greater polydrug use, injection-related HIV risk-taking behaviour, criminal involvement, psychological distress and injection-related health problems, as well as poorer general health, and an increased risk of having overdosed than other users of benzodiazepines. Of those subjects who had injected benzodiazepines, 55% were no longer current benzodiazepine injectors. Concern for general health emerged as the most common reason for having made a transition away from injecting, and for being likely to make such a transition.