Bond rolling resistance and its effect on yielding of bonded granulates by DEM analyses

A discrete element modelling of bonded granulates and investigation on the bond effect on their behaviour are very important to geomechanics. This paper presents a two‐dimensional (2‐D) discrete element theory for bonded granulates with bond rolling resistance and provides a numerical investigation into the effect of bond rolling resistance on the yielding of bonded granulates. The model consists of mechanical contact models and equations governing the motion of bonded particles. The key point of the theory is that the assumption in the original bond contact model previously proposed by the authors (55th CSCE‐ASCE Conference, Hamilton, Ont., Canada, 2002; 313–320; J. Eng. Mech. (ASCE) 2005; 131 (11):1209–1213) that bonded particles are in contact at discrete points, is here replaced by a more reliable assumption that bonded particles are in contact over a width. By making the idealization that the bond contact width is continuously distributed with the normal/tangential basic elements (BE) (each BE is composed of spring, dashpot, bond, slider or divider), we establish a bond rolling contact model together with bond normal/tangential contact models, and also relate the governing equations to local equilibrium. Only one physical parameter β needs to be introduced in the theory in comparison to the original bond discrete element model. The model has been implemented into a 2‐D distinct element method code, NS2D. Using the NS2D, a total of 86 1‐D, constant stress ratio, and biaxial compressions tests have been carried out on the bonded granular samples of different densities, bonding strengths and rolling resistances. The numerical results show that: (i) the new theory predicts a larger internal friction angle, a larger yielding stress, more brittle behaviour and larger final broken contact ratio than the original bond model; (ii) the yielding stress increases nonlinearly with the increasing value of β, and (iii) the first‐yield curve (initiation of bond breakage), which define a zone of none bond breakage and which shape and size are affected by the material density, is amplified by the bond rolling resistance in analogous to that predicted by the original bond model. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical manifold method for dynamic nonlinear analysis of saturated porous media

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.     

A major sulphur isotope event at c. 510 Ma: a possible anoxia–extinction–volcanism connection during the Early–Middle Cambrian transition?

A new approach to constraining seawater δ³⁴S and sulphate concentration using francolite‐bound sulphate reveals an abrupt increase in δ³⁴S to +50‰ around the Early–Middle Cambrian boundary. Such high δ³⁴S values are best explained by increased rates of pyrite burial due to ocean anoxia coupled with an increased sensitivity of the ocean sulphate reservoir to perturbations due to low sulphate concentrations of 500–700 μgL^?1. We argue that the spread of anoxic waters at this time was partly the result of greenhouse warming related to the eruption of the Kalkarindji Large Igneous Province of northern Australia and that it triggered the collapse of early metazoan reef ecosystems during the latest Early Cambrian. Mass extinctions of the last 260 Myr have all coincided with enhanced volcanic activity, while several are also associated with positive shifts in seawater δ³⁴S. Extending this correlation back in time further implicates volcanically induced climate change as a major determining factor in biosphere evolution. Terra Nova, 18, 257–263, 2006     

Mezen Bay—a key area for understanding Weichselian glaciations in northern Russia

Sediment successions in coastal cliffs around Mezen Bay, southeastern White Sea, record an unusually detailed history of former glaciations, interstadial marine and fluvial events from the Weichselian. A regional glaciation model for the Weichselian is based on new data from the Mezen Bay area and previously published data from adjacent areas. Following the Mikulinian (Eemian) interglacial a shelf‐centred glaciation in the Kara Sea is reflected in proglacial conditions at 100–90 ka. A local ice‐cap over the Timan ridge existed between 75 and 65 ka. Renewed glaciation in the Kara Sea spread southwestwards around 60 ka only, interrupted by a marine inundation, before it advanced to its maximum position at about 55–50 ka. After a prolonged ice‐free period, the Scandinavian ice‐sheet invaded the area from the west and terminated east of Mezen Bay about 17 ka. The previously published evidence of a large ice‐dammed lake in the central Arkhangelsk region, Lake Komi, finds no support in this study. Copyright © 2003 John Wiley & Sons, Ltd.     

A multi‐level parallelized substructuring‐frontal solution for coupled thermo/hydro/mechanical problems in unsaturated soil

This paper proposes a multi‐level parallelized substructuring–frontal combined algorithm for the analysis of the problem of thermo/hydraulic/mechanical behaviour of unsaturated soil. Temperature, displacement, pore water pressure and pore air pressure are treated as the primary variables in a non‐linear analysis. Details are given firstly of the substructuring–frontal combined approach. The incorporation of the algorithm in a multi‐level parallel strategy is then discussed. The parallel processing can thus be carried out at different substructural levels. The method thus developed impacts, in a positive way, on both computer storage requirement and execution time. Copyright © 2003 John Wiley & Sons, Ltd.     

Multistep matched-field inversion for broad-band data from ASIAEX2001

This paper discusses the results of geoacoustic inversion carried out using explosive charge data from the Asian Seas International Acoustic Experiment (ASIAEX) East China Sea (ECS) Experiment. A multifrequency incoherent matched-field inversion processor and a genetic algorithm (GA) are used for the inversion. A multistep matched field inversion approach is presented, which makes use of the varying sensitivities of wave fields at various frequencies to reduce the inversion problem into a sequence of smaller inversions with fewer unknowns to estimate at each stage. Different parameters are estimated using data at different frequencies according to their sensitivities. Inversion results for different areas in the ECS region are summarized and compared with core data.     

A one-dimensional model for the interaction between continental-scale ice sheets and atmospheric stationary waves

 The great continental ice sheets of the Pleistocene represented a significant topographic obstacle to the westerly winds in northern midlatitudes. This work explores how consequent changes in the atmospheric stationary wave pattern might have affected the shape and growth of the ice sheets themselves. A one dimensional (1-D) model is developed which permits an examination of the types and magnitudes of the feedbacks that might be expected. When plausible temperature perturbations are introduced at the ice-sheet margin which are proportional to the stationary wave amplitude, the equilibrium shape of the ice sheet is significantly altered, and depends on the sign of the perturbation. The proposed feedback also affects the response of the ice sheet to time-varying climate forcing. The results suggest that the evolution of a continental-scale ice sheet with a land-based margin may be significantly determined by the changes it induces in the atmospheric circulation. Received: 1 October 1999 / Accepted: 17 July 2000     

Small angle neutron scattering by opals

Micro- and non-crystalline opals, chalcedony and flint show diffuse small angle neutron scattering (SANS). Precious opals give rise to two additional intensity maxima at very small scattering angles which are due to Bragg reflections from the closest packed non-crystalline silica spheres. A small angle texture diagram reveals that the closest packing is faulty. Synthetic non-crystalline opals yield much less intense small angle scattering due to lower contrast between silica spheres and interstitial cement or particles; in this case intensity maxima were not observed. The outer part of the scattering curves of opal-CT and microcrystalline quartz deviates from Porod's law. The specific surface of natural non-crystalline opals ranges from 0.006 to 0.018 nm^–1. In microcrystalline opals, the specific surface is about 10 times larger than in non-crystalline opals.