Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.     

Dynamic modelling of retrogressive landslides with emphasis on the role of clay sensitivity

This paper presents a detailed numerical study of the retrogressive failure of landslides in sensitive clays. The dynamic modelling of the landslides is carried out using a novel continuum approach, the particle finite element method, complemented with an elastoviscoplastic constitutive model. The multiwedge failure mode in the collapse is captured successfully, and the multiple retrogressive failures that have been widely observed in landslides in sensitive clays are reproduced with the failure mechanism, the kinematics, and the deposition being discussed in detail. Special attention has been paid to the role of the clay sensitivity on each retrogressive failure as well as on the final retrogression distance and the final run‐out distance via parametric studies. Moreover, the effects of the viscosity of sensitive clays on the failure are also investigated for different clay sensitivities.     

Coronal change at the south-west limb observed by Yohkoh on 9 November 1991, and the subsequent interplanetary shock at Pioneer Venus Orbiter

A spectacular change in the lower corona on the south-west limb has been found in solar images taken by the Yohkoh soft X-ray telescope. The event is characterized by a large topological change in magnetic field and a large intensity decrease observed after the X1. 1/1B flare on 9 November, 1991. A coronal mass ejection (CME) was observed by the Mark III K-coronameter (MK3) at the HAO/Mauna Loa Observatory. Both the MK3 (white-light) and soft X-ray observations showed that one leg of this CME was located above the flare site. An interplanetary shock associated with this event was observed by Pioneer Venus Orbiter, and, possibly, by IMP-8.Also Cooperative Institute for Research in the Environmental Sciences (CIRES), University of Colorado, Boulder, CO 80309, U.S.A.     

Spectroscopy of comet P/Schaumasse

Column density profiles for CN, C₃, C₂ and NH have been determined from a long-slit CCD spectrum of periodic comet P/Schaumasse (1992x). Comparisons of these profiles with Haser models indicate that the ratios of the CN, C₃ and C₂ production rates are typical for a short-period comet. Although the scale lengths for NH and its parent species are uncertain, the results indicate that the production rate for NH is much greater than for either C₂ or CN.     

Going underground: in search of Carboniferous coal forests

The development of coal forests during the Carboniferous is one of the best-known episodes in the history of life. Although often reconstructed as steamy tropical rainforests, these ancient ecosystems were a far cry from anything we might encounter in the Amazon today. Bizarre giant club-mosses, horsetails and tree ferns were the dominant plants, not flowering trees as in modern rainforests. At their height, coal forests stretched all the way from Kansas to Kazakhstan, spanning the entire breadth of tropical Pangaea. Most of what we know of their biodiversity and ecology has been quite literally mined out of the ground through two centuries of hard labour. Without coal mining, our knowledge would be greatly impoverished. Over the past few years, we've been exploring underground coal mines in the United States, where entire forested landscapes have been preserved intact over huge areas. Never before have geologists had the opportunity to walk out through mile upon mile of fossilized forest. In this feature article, we describe some of our recent explorations and attempt to shed new light on these old fossils.     

The geochemistry of lead in rivers,estuaries and the continental shelf of the Southeastern United States

Transport by southeastern rivers has insignificant influence on the lead concentration of southeastern shelf waters. If the rate of mobilization and transport of lead by these rivers is representative of uncontaminated fluvial transport during the Pleistocene, only about 5% of the prehistoric output of lead to the North Atlantic can be accounted for by river input.Lead concentrations in southeastern shelf waters are similar to those of North Atlantic Water from the upper 1000 m, which is probably similar to the intrusion source water for the shelf. Atmospheric inputs to the shelf of the same magnitude as observed for the Western North Atlantic are difficult to reconcile given the residence time of shelf waters and their lead concentration unless the rate of loss of lead to shelf sediments is about the same as the atmospheric flux.     

Calibration of the kinematic method of studying solar wind disturbances on the basis of a one-dimensional MHD solution and a simulation study of the heliosphere disturbances between 22 november and 6 december 1977

An empirical kinematic method developed by Hakamada and Akasofu (1982) is calibrated on the basis of a one-dimensional MHD solution. The calibrated results are used to simulate the stream-stream interaction and the background corotating structure in a simple situation and also during 22 November–6 December 1977. The solar wind disturbances caused by solar activities during this period are then introduced into the above background stream in simulating the heliospheric disturbance event which was observed by an aligned set of spacecraft at distances between 0.6 and 1.6 a.u. The observations and the simulated results are satisfactory, and a little more refinement in the simulation could reconstruct reasonably well the data by filling the data gaps in the solar wind speed, the density and the IMF magnitude.     

Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.