A novel coupled FDM-DEM modelling method for flexible membrane boundary in laboratory tests

Difficulties are involved in discrete element method (DEM) modelling of the flexible boundary, that is, the membranes covering the soil sample, which can be commonly found in contemporary laboratory soil tests. In this paper, a novel method is proposed wherein the finite difference method (FDM) and DEM are coupled to simulate the rubber membrane and soil body, respectively. Numerical plane strain and triaxial tests, served by the flexible membrane, are implemented and analysed later. The effect of the membrane modulus on the measurement accuracy is considered, with analytical formulae derived to judge the significance of this effect. Based on an analysis of stress-strain responses and the grain rotation field, the mechanical performances produced by the flexible and rigid lateral boundaries are compared for the plane strain test. The results show that (1) the effect of the membrane on the test result becomes more significant at larger strain level because the membrane applies additional lateral confining pressure to the soil body; (2) the tested models reproduce typical stress and volumetric paths for specimens with shear bands; (3) for the plane strain test, the rigid lateral boundary derives a much higher peak strength and larger bulk dilatation, but a similar residual strength, compared with the flexible boundary. The latter produces a more uniform (or ‘diffuse') rotation field and more mobilised local kinematics than does the former. All simulations show that the proposed FDM-DEM coupling method is able to simulate laboratory tests with a flexible boundary membrane.     

An approximate method for evaluating the shear band thickness in saturated soils

A formula for the thickness of a shear band formed in saturated soils under a simple shear or a combined stress state has been proposed. It is shown that the shear band thickness is dependent on the pore pressure properties of the material and the dilatancy rate, but is independent of the details of the combined stress state. This is in accordance with some separate experimental observations. Copyright © 2004 John Wiley & Sons, Ltd.     

Strain localization in sand: an overview of the experimental results obtained in Grenoble using stereophotogrammetry

Experimental results are presented from the extensive program of drained plane strain compression tests on sand carried out in Grenoble over the last two decades. Systematic analysis of photographs of the deforming specimen allowed for measuring deformations and determining strain fields throughout the test, that is: prior to, at, and after the onset of strain localization. The principles, details and accuracy of the procedure are described, as well as its suitability to properly depict the patterns of deformation. Findings concerning the occurrence and progression of strain localization are discussed. The issues of shear band orientation and thickness are addressed, as well as temporary and persistent complex localization patterns, and the volumetric behaviour inside a band after its formation. The influence of such variables as initial state of the sand (effective stress and relative density), specimen size and slenderness, as well as grain size, is discussed. Copyright © 2004 John Wiley & Sons, Ltd     

Paleoproterozoic Extensional Structure and Its Controlling on Mineralization in the East of Liaoning Province

IntroductionThe area of eastern Liaoning is an importantmetal and nonmetal metallogenetic district in China,and the Liaohe group is one of the most importantstrata that hosts Pb, Zn, Au, B and Mg etcstratabound deposits. Up to now many geo1ogistssuch as Z…     

Partial derivatives of surface wave phase velocities for flat anisotropic models

Summary. Surface wave behaviour in flat anisotropic structures is first illustrated by performing an exact computation on a simple two-layer model. The variational procedure of Smith & Dahlen is then used to compute the partial derivatives of surface wave phase velocities with respect to the elastic parameters in more realistic earth models. Linear relationships between the partial derivatives for a general anisotropic structure and those for a transversely isotropic structure are derived. When considering waves propagating in a fixed direction, there are only four independent derivatives for Rayleigh waves, and two for Love waves. To avoid the lack of resolution in an inverse method, we propose to use physically constrained models. These results are illustrated by using a model with hexagonal symmetry and a symmetry axis oriented either vertically or horizontally. Quasi-Love- and quasi-Rayleigh-wave partial derivatives are computed for both axis orientations. Modes up to the second overtone and periods ranging between 45 and 130 s have been considered. Finally, anomalies of phase velocity are computed in an oceanic model made of 1/6 oriented olivine crystals with horizontal or vertical preferred orientations of the a -axis.     

Inversion of multimode surface wave data: evidence for sub-lithospheric anisotropy

Summary. Multiparameter inversions of multimode dispersion data are performed for two large regions: the Pacific Ocean and North America. Anisotropy is taken into account by considering transversely isotropic structures with a vertical axis of symmetry. Two fundamental questions are studied in detail: (1) how to make the inverted models consistent when using different sets of parameters, (2) what is the significance of transversely isotropic inversion for the actual Earth's structure? It is proved that full consistency of the inverted models can be achieved by properly taking into account some a priori informations on the model and it is shown that the use of transversely isotropic models with vertical axis of symmetry does not cause severe limitations when interpreting the data. The models we have obtained are discussed in the light of these investigations. Considering an olivine-rich upper mantle, we make a tentative interpretation of these models in terms of preferred orientation of the a -axis of the crystals in one fixed horizontal direction.     

Numerical modelling of the impact crater depth-diameter dependence in an acoustically fluidized target

2D numerical modelling of impact cratering has been utilized to quantify an important depth-diameter relationship for different crater morphologies, simple and complex. It is generally accepted that the final crater shape is the result of a gravity-driven collapse of the transient crater, which is formed immediately after the impact. Numerical models allow a quantification of the formation of simple craters, which are bowl-shaped depressions with a lens of rock debris inside, and complex craters, which are characterized by a structural uplift. The computation of the cratering process starts with the first contact of the impactor and the planetary surface and ends with the morphology of the final crater. Using different rheological models for the sub-crater rocks, we quantify the influence on crater mechanics. To explain the formation of complex craters in accordance to the threshold diameter between simple and complex craters, we utilize the Acoustic Fluidization model. We carried out a series of simulations over a broad parameter range with the goal to fit the observed depth/diameter relationships as well as the observed threshold diameters on the Moon, Earth and Venus.