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     

Elastic solutions for stresses in a transversely isotropic half‐space subjected to three‐dimensional buried parabolic rectangular loads

In many areas of engineering practice, applied loads are not uniformly distributed but often concentrated towards the centre of a foundation. Thus, loads are more realistically depicted as distributed as linearly varying or as parabola of revolution. Solutions for stresses in a transversely isotropic half‐space caused by concave and convex parabolic loads that act on a rectangle have not been derived. This work proposes analytical solutions for stresses in a transversely isotropic half‐space, induced by three‐dimensional, buried, linearly varying/uniform/parabolic rectangular loads. Load types include an upwardly and a downwardly linearly varying load, a uniform load, a concave and a convex parabolic load, all distributed over a rectangular area. These solutions are obtained by integrating the point load solutions in a Cartesian co‐ordinate system for a transversely isotropic half‐space. The buried depth, the dimensions of the loaded area, the type and degree of material anisotropy and the loading type for transversely isotropic half‐spaces influence the proposed solutions. An illustrative example is presented to elucidate the effect of the dimensions of the loaded area, the type and degree of rock anisotropy, and the type of loading on the vertical stress in the isotropic/transversely isotropic rocks subjected to a linearly varying/uniform/parabolic rectangular load. Copyright © 2002 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…     

The second order characteristic condition for plastic flow

In this paper the second order characteristic (discontinuous bifurcation) condition is derived for the granular flow (fully plastic) equations. This second order bifurcation equation is shown to be formally identical to the first order localization requirement during steady elastoplastic deformation provided the elastic compliance tensor is substituted for the product of the plastic multiplier with the flow Hessian. For isotropic yield and flow functions the invariant form of the characteristic condition is given in detail, as well as an alternative expression in adapted co‐ordinates. The characteristic condition can be regarded as defining a hardening function which is maximized to identify the critical angles. When the method is applied to 3D Coulomb flow, Mohr's 3D fracture plane conditions are obtained uniquely. Copyright © 2003 John Wiley & Sons, Ltd.     

Magnitude of σ1, σ2, and σ3 at mid-crustal levels in an orogenic belt: Microboudin method applied to an impure metachert from Turkey

We revised an equation for estimating palaeostress magnitude using the microboudin technique by incorporating the influence of time on the fracture strength of minerals. The equation was used to estimate triaxial palaeostresses from a rare sample of metachert from Turkey that contains microboudinaged, columnar tourmaline grains in a wide range of orientations within the foliation plane. The estimated principal palaeostresses are σ₁ = 605 MPa, σ₂ = 598 MPa, and σ₃ = 597 MPa. As the microboudinage is considered to have occurred immediately before the rock encountered the brittle-plastic transition during exhumation, these stress values correspond to conditions at approximately 18 km depth and 300 °C within a Cretaceous orogenic belt.     

江西浒坑钨矿含矿石英脉的地质特征及成矿构造演化

位于赣西武功山地区的江西浒坑钨矿是一个大型石英脉型钨矿床,与赣南典型的“五层楼”模式不同,矿脉主要集中在岩体内接触带。含矿石英脉除了普通的块状构造以外,还普遍存在条带状构造。详细的野外编录和镜下研究表明含矿块状石英脉在形成后遭受不同程度的韧性剪切,从而形成不同成分的分异条带。剪切带以出现新生的面理和线理（拉伸线理）为特征。岩石学研究表明条带状矿石属于糜棱岩,具有丰富的塑性变形显微构造特征。通过对含矿石英脉的构造解析,认为与岩浆岩有关的原生节理系统基本上是同时形成的,沿节理充填的石英脉经历了长期的、多期次的构造运动,从而形成了复杂的石英脉和矿石类型。建立了原生节理形成→石英充填→韧性剪切发生→二期石英充填→脆性碎裂的矿区构造演化序列。韧性剪切作用的动力学背景可能是受后期岩浆侵入的影响,最后的脆性破裂应与区域上存在南北向的压应力有关。