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.     

Modelling of deformation response and chemo‐mechanical coupling in chalk

A constitutive relation is derived for describing the mechanical response of chalk. The approach is based on a phenomenological framework which employs chemo‐plasticity. The properties of the material are assumed to be affected by the physico‐chemical processes that occur through the interaction between the skeleton and the pore fluid. The underlying mechanism is discussed by invoking a micromechanical analysis. The performance of the framework is illustrated by examining the evolution of mechanical characteristics in the presence of different pore fluids. Copyright © 2006 John Wiley & Sons, Ltd.     

Chemical alteration of tephra in the depositional environment: theoretical stability modelling

The study of the chemical stability of vitreous material in aqueous media is well‐established. There has to date been little consideration of the implications of variations in the chemical durability of tephra in Quaternary tephrochronology. Chemical alteration can take the form of cationic leaching from the matrix, or complete destruction of the silica network, either of which could constrain the ability to chemically identify distal tephra. Here we apply established models of vitreous durability to the published chemical analyses of a large number of Icelandic tephras in order to predict their relative durabilities under equivalent conditions. This suggests that some important tephras have relatively poor chemical stability, and that rhyolitic tephras are, in general, more stable than basaltic. We conclude that tephras should be expected to show predictable differential chemical stability in the post‐depositional environment. 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.     

Evidence for a mid‐Pleistocene shift of ice‐drift pattern in the Nordic seas

Sediment proxy records from a continuous, 1.5 million year long deep‐sea sediment core from a site in the western Norwegian Sea were used to obtain new insights into the nature of palaeoceanographic change in the northern North Atlantic (Nordic seas) during the climatic shift of the Mid‐Pleistocene Revolution (MPR). Red‐green sediment colour and magnetic susceptibility records both reveal significant differences in their mean values when comparing the intervals older than 700 000 yr (700 ka) with those from the past 500 kyr. The timing and duration of these changes indicates that the MPR in the Nordic seas is characterised by a gradual transition lasting about 200 kyr. Together with further sedimentological evidence this suggests that the mid‐Pleistocene climate shift was accompanied by a general change in ice‐drift pattern. It is further proposed that prior to the onset of the major late Pleistocene glaciations in the Northern Hemisphere a significant proportion of the ice in the eastern Nordic seas originated from a southern provenance, whereas later it dominantly came from the surrounding landmasses. Copyright © 2003 John Wiley & Sons, Ltd.     

Descriptors of impulsive signal levels commonly used in underwateracoustics

Impulsive acoustic waveforms are characterized by a different set of derived quantities than are continuous waveforms. This note presents commonly accepted definitions, units, and symbols used to describe the magnitude of impulsive underwater signals     

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.     

Ocean-floor hydrothermal metamorphism in the Limousin ophiolites (western French Massif Central): evidence of a rare preserved Variscan oceanic marker

The Limousin ophiolite is located at the suture zone between two major thrust sheets in the western French Massif Central. This ophiolitic section comprises mantle‐harzburgite, mantle‐dunite, wehrlites, troctolites and layered gabbros. It has recorded a static metamorphic event transforming the gabbros into undeformed amphibolites and the magmatic ultramafites into serpentinites and/or pargasite‐bearing chloritites. With various thermobarometric methods, it is possible to show that the different varieties of amphibole have registered low‐P (c. 0.2 GPa) conditions with temperature ranging from high‐T, late‐magmatic conditions to greenschist–zeolite metamorphic facies. The abundance of undeformed metamorphic rocks (which is typical of the lower oceanic crust), the occurrence of Ca–Al (–Mg) metasomatism illustrated by the growth of Ca–Al silicates in veins or replacing the primary magmatic minerals, the P–T conditions of the metamorphism and the numerous similarities with oceanic crustal rocks from Ocean Drilling Program and worldwide ophiolites are the main arguments for an ocean‐floor hydrothermal metamorphism in the vicinity of a palaeo‐ridge. Among the West‐European Variscan ophiolites, the Limousin ophiolites constitute an extremely rare occurrence that has not been involved in any HP (subduction‐related) or MP (orogenic) metamorphism as observed in other ophiolite occurrences (i.e. France, Spain and Germany).     

Estimation of dispersion in estuaries by a self-consistent method

Longitudinal and transverse dispersion in estuaries may be determined in principle by measuring the change in time of the concentration field of a water soluble dye. In practice, measurements may take the form of several sampling runs crossing the dye patch. The coefficients of dispersion are usually estimated from the second moments of the concentration field; the use of second moments, in particular, introduces errors because measurements at long times or distance are unduly emphasized.A self-consistent method, free of the disadvantages inherent in the use of higher moments, is described in which the distribution of concentration in space is related to the zeroth spatial moment of concentration. In this formulation, based upon sampling runs across the dye patch, the concentration field in the estuary may be directly and simply related to the dispersion coefficients in both the intercept and gradient of a linear graph.The methods are illustrated by the characterization of a series of experiments in the Inner Estuary of Milford Haven.