Boundary element analysis for effective stiffness tensors: effect of fabric tensor determination method

Second-rank fabric tensors have been extensively used to describe structural anisotropy and to predict orthotropic elastic constants. However, there are many different definitions of, and approaches to, determining the fabric tensor. Most commonly used is a fabric tensor based on mean intercept length measurements, but star volume distribution and star length distribution are commonly used, particularly in studies of trabecular bone. Here, we investigate the effect of the fabric tensor definition on elastic constant predictions using both synthetic, idealized microstructures as well as a micrograph of a porous ceramic. We use an efficient implantation of a symmetric Galerkin boundary element method to model the mechanical response of the various microstructures, and also use a boundary element approach to calculate the necessary volume averages of stress and strain to obtain the effective properties of the media.     

Stress-induced anisotropy in sand under cyclic loading

The anisotropy of a granular material’s structure will influence its response to applied loads and deformations. Anisotropy can be either inherent (e.g. due to depositional process) or induced as a consequence of the applied stresses or strains. Discrete element simulations allow the interactions between individual particles to be explicitly simulated and the fabric can be quantified using a fabric tensor. The eigenvalues of this fabric tensor then give a measure of the anisotropy of the fabric. The coordination number is a particle scale scalar measure of the packing density of the material. The current study examines the evolution of the fabric of a granular material subject to cyclic loading, using two-dimensional discrete element method (DEM) simulations. Isotropic consolidation modifies and reduces the inherent anisotropy, but anisotropic consolidation can accentuate anisotropy. The ratio of the normal to shear spring stiffness at the particle contacts in the DEM model affects the evolution of anisotropy. Higher ratios reduce the degree of anisotropy induced by anisotropic consolidation. The anisotropy induced by cyclic loading depends on the amplitude of the loading cycles and the initial anisotropy.     

Parametric hypoplasticity as continuum model for granular media: from Stokesium to Mohr-Coulombium and beyond

Fluid-particle systems, in which internal forces arise only from viscosity or intergranular friction, represent an important special case of strictly dissipative materials defined by a history-dependent 4th-rank viscosity tensor. In a recently proposed simplification, this history dependence is represented by a symmetric 2nd-rank fabric tensor with evolution determined by a given homogeneous deformation. That work suggests an essential physical link between idealized suspensions (“Stokesium”) and granular media (“Mohr-Coulombium”) along with possible models for the visco-plasticity of fluid-saturated and dry granular media. The present paper deals with the elastoplasticity of dilatant non-cohesive granular media composed of nearly rigid, frictional particles. Based on the underlying physics and past modeling by others, a continuum model based on parametric hypoplasticity is proposed, which involves a set of rate-independent ODEs in the state-space of stress, void ratio and fabric. As with the standard theory of hypoplasticity, the present model does not rely on plastic potentials but, in contrast to that theory, it is based explicitly on positive-definite elastic and plastic moduli. The present model allows for elastic loading or unloading within a dissipative yield surface and also provides a systematic treatment of Reynolds dilatancy as a kinematic constraint. Some explicit forms are proposed and comparisons are made to previous hypoplastic models of granular media.     

Macroscopic model with anisotropy based on micro–macro information

Physical experiments can characterize the elastic response of granular materials in terms of macroscopic state variables, namely volume (packing) fraction and stress, while the microstructure is not accessible and thus neglected. Here, by means of numerical simulations, we analyze dense, frictionless granular assemblies with the final goal to relate the elastic moduli to the fabric state, i.e., to microstructural averaged contact network features as contact number density and anisotropy. The particle samples are first isotropically compressed and then quasi-statically sheared under constant volume (undrained conditions). From various static, relaxed configurations at different shear strains, infinitesimal strain steps are applied to “measure” the effective elastic response; we quantify the strain needed so that no contact and structure rearrangements, i.e. plasticity, happen. Because of the anisotropy induced by shear, volumetric and deviatoric stresses and strains are cross-coupled via a single anisotropy modulus, which is proportional to the product of deviatoric fabric and bulk modulus (i.e., the isotropic fabric). Interestingly, the shear modulus of the material depends also on the actual deviatoric stress state, along with the contact configuration anisotropy. Finally, a constitutive model based on incremental evolution equations for stress and fabric is introduced. By using the previously measured dependence of the stiffness tensor (elastic moduli) on the microstructure, the theory is able to predict with good agreement the evolution of pressure, shear stress and deviatoric fabric (anisotropy) for an independent undrained cyclic shear test, including the response to reversal of strain.     

三维连续体结构仿生拓扑优化新方法

Wolff法则是指骨骼通过重建/生长,保证骨小梁方向趋于与主应力方向一致以不断地适应它的力学环境。根据Wolff法则,建立了一种新的拓扑优化的准则法。该方法的基本思想是:(1)将待优化的结构看作是一块遵从Wolff法则生长的骨骼,骨骼的重建过程作为三维连续体结构寻找最优拓扑的过程;(2)用构造张量描述正交各向异性材料的弹性本构;(3)重建规律为结构中材料的更新规律。通过引入参考应变区间,材料更新规律可解释为:设计域内一点处主应变的绝对值不在该区间时,该点处构造张量出现变化;否则,构造张量不变化,该点处于生长平衡状态。(4)当设计域内所有点都处于生长平衡状态时,结构拓扑优化结束。采用各向同性本构模型,即令二阶构造张量与二阶单位张量成比例,分析三维结构拓扑优化。实例进一步验证基于Wolf法则的连续体结构优化方法的正确性和可行性。     

On micromechanical characteristics of the critical state of two-dimensional granular materials

In micromechanics of quasi-static deformation of granular materials, relationships are investigated between the macro-scale, continuum-mechanical characteristics, and the micro-scale characteristics at the particle and interparticle contact level. An important micromechanical quantity is the fabric tensor that reflects the distribution of contact orientations. It also contains information on the coordination number, i.e. the average number of contacts per particle. Here, the focus is on characteristics of the critical state in the two-dimensional case. Critical state soil mechanics is reviewed from the micromechanical viewpoint. Two-dimensional discrete element method (DEM) simulations have been performed with discs from a fairly narrow particle-size distribution. Various values for the interparticle friction coefficient and for the confining pressure have been considered to investigate the effect of these quantities on critical state characteristics (shear strength, packing fraction, coordination number and fabric anisotropy). Results from these DEM simulations show that a limiting fabric state exists at the critical state, which is geometrical in origin. The contact network tessellates the assembly into loops that are formed by contacts. For each loop, a symmetrical loop tensor is defined, based on its contact normals. This loop tensor reflects the shape of the loop. An orientation is associated with each loop, based on its loop tensor. At the critical state, the frequencies with which loops with different number of sides occur depend on the coordination number. At the critical state, these loops have, on average, the following universal characteristics, i.e. independent of the coordination number: (1) loops with the same number of sides and orientation have identical anisotropy of the loop tensor, (2) the anisotropy of the loop tensor depends linearly on the number of sides of the loop, (3) the distribution of loop orientations is identical, (4) Lewis’s law for the loop areas, which is a linear relation between the number of sides of loops and their area, is satisfied (not exclusively at the critical state) and (5) the areas of the loops do not depend on their orientation.     

Influence of particle shape on sheared dense granular media

We study by means of molecular dynamics simulations of periodic shear cells, the influence of particle shape on the global mechanical behavior of dense granular media. At large shear deformation samples with elongated particles, independent of their initial orientation, reach the same stationary value for both shear force and void ratio. At the micro-mechanical level the stress, the fabric and the inertia tensors of the particles are used to study the evolution of the media. In the case of isotropic particles the direction of the principal axis of the fabric tensor is aligned with the one of the principal stress, while for elongated particles the fabric orientation is strongly dependent on the orientation of the particles. The shear band width is shown to depend on the particle shape due to the tendency of elongated particles to preferential orientations and less rotation.     

A 3D fourth order fabric tensor approach of anisotropy in granular media

We propose a micromechanical approach for granular media, with a particular account of the texture-induced anisotropy and of the strain localization rule. The approach is mainly based on the consideration of a fourth order fabric tensor able to capture general anisotropy which can be induced by complex distribution of contacts. Incorporation of this fourth order fabric tensor in a suitable homogenization scheme allows to determine the corresponding macroscopic elastic properties of the granular material. For this purpose, in addition to the classical Voigt upper bound, a new kinematics-based localization rule is proposed. It generalizes the one formulated by Cambou et al. [B. Cambou, Ph. Dubujet, F. Emeriault, F. Sidoroff, Eur. J. Mech. A/Solids 14 (1995) 225–276] in the case of an isotropic contact distribution. The results of the complete model compare well to numerical simulations results when available [C.S. Chang, C.L. Liao, Appl. Mech. Rev. 47 (1 Part 2) (1994) 197–207] (case of isotropic distribution of contacts). Finally, the interest of the fourth order fabric tensor based approach combined with the proposed localization rule is shown for different distributions of contacts by comparing its predictions to those given by a second order fabric tensor approach.     

Macroscopic material properties from quasi-static, microscopic simulations of a two-dimensional shear-cell

One of the essential questions in the area of granular matter is, how to obtain macroscopic tensorial quantities like stress and strain from “microscopic” quantities like the contact forces in a granular assembly. Different averaging strategies are introduced, tested, and used to obtain volume fractions, coordination numbers, and fabric properties. We derive anew the non-trivial relation for the stress tensor that allows a straightforward calculation of the mean stress from discrete element simulations and comment on the applicability. Furthermore, we derive “elastic” (reversible) mean displacement gradient, based on a best-fit hypothesis. Finally, different combinations of the tensorial quantities are used to compute some material properties. The bulk modulus, i.e. the stiffness of the granulate, is a linear function of the trace of the fabric tensor which itself is proportional to the density and the coordination number. The fabric, the stress and strain tensors are not co-linear so that a more refined analysis than a classical elasticity theory is required. Received: 23 July 1999     

Stress-chain based micromechanics of sand with grain shape effect

The mechanical behaviors of granular media are controlled by grain properties and microstructure. The primary property of granular media is denoted by its grain shape, grain size distribution, stiffness, and interparticle friction. The grain shape itself is of particular importance. Microstructures are formed in the connection paths of contact points between grains. In this paper, the deformation of granular materials with different grain shapes was simulated using two-dimensional DEM under different stress-levels and densities. After analyzing the results, the authors investigated fabric changes. The evolution rule of stress-induced anisotropy and its limitation as well as the existence of a critical state of fabric are revealed.     

Micromechanical behavior of granular materials with inherent anisotropy under cyclic loading using 2D DEM

This paper presents the micromechanical behavior of granular materials due to different initial inherent anisotropic conditions during cyclic loading using the discrete element method (DEM). Oval particles were used to model the samples. Three samples, with three different inherent anisotropic conditions based on the particle’s bedding direction, were prepared and subjected to biaxial cyclic loading. The differences in the inherent anisotropic conditions of the samples affect the stress–strain-dilative behavior of granular materials. The width of the stress–strain cyclic loops decreases as the preferred bedding angle changes from vertical to horizontal. Contact fabric evolution is found to be dependent on the inherent anisotropic fabric of the sample during loading and unloading. The fabric anisotropy is dominant for horizontal particle bedding at the end of loading and for vertical particle bedding at the end of unloading. A change in fabric anisotropy is observed only for the first few loading–unloading cycles for the given conditions depicted in the present study.     

A flow rule incorporating the fabric and non-coaxiality in granular materials

Deficiencies of constitutive models in prediction of dilatancy are often attributed to simplifications associated with flow rules such as assumptions of isotropy and coaxiality. It is thus proposed here to develop a comprehensive flow rule for granular materials by including the effect of fabric and without the assumption of coaxiality. A second-order tensor is introduced as a fabric for the distribution of contact normals and contact forces. By using the energy principle in micro-mechanical scale and a suitable dissipation mechanism in granular materials, a stress-dilatancy relation is obtained. Fabric plays a “bridge-like” role in the dilatancy and non-coaxiality. Non-coaxialities between stress-strain-fabric are attributed to the non-coaxiality between stress-fabric and strain-fabric. In this formulation the constants for modeling fabric depend on non-coaxiality of the system rather than the history that determines such a state. Ability of this stress-fabric-dilatancy for modeling the non-coaxiality shows that this relation can predict the behavior of granular materials in the presence of the rotation of principal stress axes.     

3D experimental quantification of fabric and fabric evolution of sheared granular materials using synchrotron micro-computed tomography

Many experimental studies have demonstrated that mechanical response of granular materials is highly influenced by micro-structural fabric and its evolution. In the current literature, quantification of fabric and its evolution has been developed based on micro-structural observations using Discrete Element Method or 2D experiments with simple particle shapes. The emergence of X-ray computed tomography technique has made quantification of such experimental micro-structural properties possible using 3D high-resolution images. In this paper, synchrotron micro-computed tomography was used to acquire 3D images during in-situ conventional triaxial compression experiments on granular materials with different morphologies. 3D images were processed to quantify fabric and its evolution based on experimental measurements of contact normal vectors between particles. Overall, the directional distribution of contact normals exhibited the highest degree of isotropy at initial state (i.e., zero global axial strain). As compression progressed, contact normals evolved in the direction of loading until reaching a constant fabric when experiments approached the critical state condition. Further assessment of the influence of confining pressure, initial density state, and particle-level morphology on fabric and its evolution was formed. Results show that initial density state and applied confining pressure significantly influence the fabric-induced internal anisotropy of tested specimens at initial states. Relatively, a higher applied confining pressure and a looser initial density state resulted in a higher degree of fabric-induced internal anisotropy. Influence of particle-level morphology was also found to be significant particularly on fabric evolution.     

Macro deformation and micro structure of 3D granular assemblies subjected to rotation of principal stress axes

This paper presents a numerical investigation on the behavior of three dimensional granular materials during continuous rotation of principal stress axes using the discrete element method. A dense specimen has been prepared as a representative element using the deposition method and subjected to stress rotation at different deviatoric stress levels. Significant plastic deformation has been observed despite that the principal stresses are kept constant. This contradicts the classical plasticity theory, but is in agreement with previous laboratory observations on sand and glass beads. Typical deformation characteristics, including volume contraction, deformation non-coaxiality, have been successfully reproduced. After a larger number of rotational cycles, the sample approaches the ultimate state with constant void ratio and follows a periodic strain path. The internal structure anisotropy has been quantified in terms of the contact-based fabric tensor. Rotation of principal stress axes densifies the packing, and leads to the increase in coordination numbers. A cyclic rotation in material anisotropy has been observed. The larger the stress ratio, the structure becomes more anisotropic. A larger fabric trajectory suggests more significant structure re-organization when rotating and explains the occurrence of more significant strain rate. The trajectory of the contact-normal based fabric is not centered in the origin, due to the anisotropy in particle orientation generated during sample generation which is persistent throughout the shearing process. The sample sheared at a lower intermediate principal stress ratio \((b=0.0)\) has been observed to approach a smaller strain trajectory as compared to the case \(b=0.5\), consistent with a smaller fabric trajectory and less significant structural re-organisation. It also experiences less volume contraction with the out-of plane strain component being dilative.     

Dependence of elastic constants of an anisotropic porous material upon porosity and fabric

The elastic properties of an anisotropic porous material can be represented as functions of the material's solid volume fraction (or porosity) and the principal diameters of the material's fabric ellipsoid. The fabric ellipsoid is a measure of the anisotropy of the microstructure of a material. The definitions and measurement techniques for fabric ellipsoids in granular materials, foams, cancellous bone, and rocks are discussed. The principal results presented in this work are algebraic expressions for the dependence of the orthotropic elastic constants upon both solid volume fraction and the fabric ellipsoid.     

The influence of rolling resistance on the stress-dilatancy and fabric anisotropy of granular materials

The effects of rolling resistance on the stress-dilatancy behavior and fabric anisotropy of granular materials were investigated through a three-dimensional discrete element method (DEM). A rolling resistance model was incorporated into the DEM code PFC3D and triaxial DEM simulations under simulated drained and undrained conditions were carried out. The results show that there existed a threshold value of the rolling friction. When the rolling friction was smaller than this value, the mechanical behavior of granular materials under both drained and undrained conditions were substantially influenced by the rolling friction, but the influence diminished when it was larger than the threshold value. A linear relationship has been observed between the dilatancy coefficient and the natural logarithm of the rolling-friction coefficient when it was smaller than the threshold value. An increase in the rolling friction led to an increase in the fabric anisotropy of all strong contacts under both testing conditions until the threshold value was attained. The investigation on the effect of rolling friction on the microstructure of granular materials revealed that the rolling friction enhanced the stability of force chains, which resulted in the difference in the stress-dilatancy behavior. Finally, the relationship between the stress ratio q/p\(^{\prime }\) and the fabric measure at strong contacts \(\hbox {H}_{\mathrm{d}}^{\mathrm{s}} /\hbox {H}_{\mathrm{m}}^{\mathrm{s}}\) was found independent of the inter-particle friction, rolling friction and testing conditions.     

Effect of fabric anisotropy on shear localization in sand during plane strain compression

Summary The paper focuses on the effect of fabric anisotropy on shear localization in cohesionless granular materials. For the numerical simulation, a hypoplastic constitutive model was used. In order to take into account a characteristic length of the micro-structure, the constitutive model was extended to include the second gradient of the Euclidian norm of the deformation rate. The hypoplastic model captures the salient features of granular bodies in a wide range of density and pressure with a single set of parameters. Transversal isotropy is described by the dyadic product of the normal vector of the space orientation of the plane of symmetry. FE-simulations of plane strain compression under constant lateral pressure were carried out with a medium dense specimen for both uniform and stochastic distribution of the initial void ratio. The effect of the direction of the bedding plane and the initial void ratio distribution on the load-deformation behavior was investigated. Moreover, the location, thickness and inclination of the shear zone were also analyzed.     

Characterization of microstructural anisotropy in orthotropic materials using a second rank tensor

A second rank symmetric tensor which describes the degree of orientation in orthotropic materials is presented and shown to reflect accurately patterns of experimental data. The use of this tensor to describe microstructural anisotropy is compared to currently accepted methods and is found to be more useful and accurate in experimental studies. A method for determining the anisotropy tensor in a material is given, based on measurements on any three mutually perpendicular planes, and the fundamental restriction of this method to orthotropic materials is discussed. Experimentally determined anisotropy tensors in five specimens of cancellous bone from five different human bones are given.     

Procedures for stereological estimation of structural anisotropy

Actual procedures are given for detecting structural anisotropy by the stereological method, i.e. by observing only cross-sections in the material, on the basis of the theoretical formulation of the previous article [1]. First, computation by the Monte Carlo method is described explicitly in terms of observed data, determining the distribution density completely. Next, alternative procedures are given by the use of restricted cross-sections on the assumption that the anisotropy is “weak,” facilitating actual experiments a great deal. Observations are made only on planes parallel to the three coordinate planes or cylindrical surfaces around the three coordinate axes. “Fabric tensors” describing the anisotropy and the “equivalent strain” are given explicitly in terms of observed data.