Initial planar deformation of dilatant granular materials

A theory for the initial planar deformation of dilatant granular materials based on a kinematic proposal of R. Butterfield and R.M. Harkness (1972) is presented. The theory introduces an additional parameter called the angle of dilatancy into the traditional structure of plasticity theories for granular materials and soils. When the angle of dilatancy is zero, the present theory reduces to the theory introduced by A.J.M. Spencer in 1964. When the angle of dilatancy is equal to the angle of internal friction, the present theory reduces to the planar form of the theory introduced by D. C. Drucker and W. Prager in 1952. The properties of the theory presented here include coincidence of the stress and velocity characteristics, realistic energy dissipation predictions, and, in general, non-coincidence of the principal axes of stress and strain-rate. However, the angle of dilatancy is assumed to be a constant in this analysis and it does not decrease to zero with increased monotonic shearing deformation as experiment requires that it should, the theory therefore being limited to the initial deformation of dilatant granular materials.     

颗粒材料内部滑动变形局部化的离散元模拟

颗粒材料在加载过程中表现出非常复杂的变形性质.在微观尺度上,单个颗粒的运动、单个孔隙胞元内的变形都是非均匀的,但也不是完全随机、没有规律可循的,而是呈现一定结构性的分布,如微带、剪切带等.本文用基于孔隙胞元的离散元方法对颗粒体进行双轴加载数值试验,模拟了以滑动变形表征的变形局部化现象.数值结果发现,对应加载过程中的不同...     

Micro-mechanical analysis of deformation characteristics of three-dimensional granular materials

The deformation characteristics of idealized granular materials have been studied from the micro-mechanical viewpoint, using Bagi’s three-dimensional micro-mechanical formulation for the strain tensor [Bagi, K., 1996. Mechanics of Materials 22, 165–177]. This formulation is based on the Delaunay tessellation of space into tetrahedra. The set of edges of the tetrahedra can be divided into physical contacts and virtual contacts between particles. Bagi’s formulation expresses the continuum, macro-scale strain as an average over all edges, of their relative displacements (between two successive states) and the complementary-area vectors. This latter vector is a geometrical quantity determined from the set of edges, i.e. from the structure of the particle packing.Results from Discrete Element Method simulations of isotropic and triaxial loading of a three-dimensional polydisperse packing of spheres have been used to investigate statistics of the branch vectors and complementary-area vectors of edges (subdivided into physical and virtual contacts) and of the relative displacements of edges. The investigated statistics are probability density functions and averages over groups of edges with the same orientation. It is shown that these averages can be represented by second-order Fourier series in edge orientation.Edge orientations are distributed isotropically, contrary to contact orientations. The average lengths of the branch vectors and the normal component of the complementary-area vectors are distributed isotropically (with respect to the edge orientation) and their average values are related to each other and to the volume fraction of the assembly. The other two components of the complementary-area vector are zero on average.The total deformation of the assembly, as given by the average of the relative displacements of the edges of the Delaunay tessellation follows the uniform-strain prediction. However, neither the deformation of the physical contact network nor of the virtual contact network has this property. The average relative displacement of physical edges in the normal direction (determined by the branch vector) is smaller than that according to the uniform-strain assumption, while that of virtual contacts is larger. This is caused by the high interparticle stiffness that hinders compression. The reverse observation holds for the tangential component of the relative displacement vector. The contribution of the deformation of the empty space between physical contacts to the continuum, macro-scale strain tensor is therefore very important for the understanding and the prediction of the macro-scale deformation of granular materials.     

Influence of porosity on the non-linear deformation of granular materials

An algorithm is suggested to determine non-linear deformative properties of granular composite materials with non-linearly deformed porous components. Stochastic differential equations of physically non-linear elasticity theory with the following use of conditional averaging are assumed as a basis. Solution of non-linear equations with respect to conditional moment is constructed by the method of successive approximation. The volume content of the components and pores in the components have been studied for their effect on the character of the strain curves of porous granular composite.     

Micromechanical definition of an entropy for quasi-static deformation of granular materials

A micromechanical theory is formulated for quasi-static deformation of granular materials, which is based on information theory. A reasoning is presented that leads to the definition of an information entropy that is appropriate for quasi-static deformation of granular materials. This definition is based on the hypothesis that relative displacements at contacts with similar orientations are independent realisations of a random variable. This hypothesis is made plausible based on the results of Discrete Element simulations. The developed theory is then used to predict the elastic behaviour of granular materials in terms of micromechanical quantities. The case considered is that of two-dimensional assemblies consisting of non-rotating particles with an elastic contact constitutive relation. Applications of this case are the initial elastic (small-strain) deformation of granular materials. Theoretical results for the elastic moduli, relative displacements, energy distribution and probability density functions are compared with results obtained from the Discrete Element simulations for isotropic assemblies with various average numbers of contacts per particle and various ratios of tangential to normal contact stiffness. This comparison shows that the developed information theory is valid for loose systems, while a theory based on the uniform-strain assumption is appropriate for dense systems.     

Experimental analysis of crack tip fields in rubber materials under large deformation

A three-nested-deformation model is proposed to describe crack-tip fields in rubber-like materials with large deformation.The model is inspired by the distribution of the measured in-plane and out-of-plane deformation.The inplane displacement of crack-tip fields under both Mode I and mixed-mode(Mode I-II) fracture conditions is measured by using the digital Moire’ method.The deformation characteristics and experimental sector division mode are investigated by comparing the measured displacement fields under different fracture modes.The out-of-plane displacement field near the crack tip is measured using the three-dimensional digital speckle correlation method.     

On phase transformations in shape memory alloy materials and large deformation generalized plasticity

A new version of rate-independent generalized plasticity, suitable for the derivation of general thermomechanical constitutive laws for materials undergoing phase transformations, is proposed within a finite deformation framework. More specifically, by assuming an additive decomposition of the finite strain tensor into elastic and inelastic (transformation induced) parts and by considering the fractions of the various material phases as internal variables, a multi-phase formulation of the theory is developed. The concepts presented are applied for the derivation of a three-dimensional thermomechanical model for shape memory alloy materials. The ability of the model in simulating several patterns of the extremely complex behavior of these materials, under both monotonic and cyclic loadings, is assessed by representative numerical examples.     

非线性弹性大变形材料梁的抗爆特性

从理论上探讨了非线性弹性大变形材料应用于抗爆结构的可行性,为此,基于等效结构体系的分析原理,将两端固定铰支梁的横向和纵向位移表示为三角级数形式,应用第二类Lagrange方程建立了非线性大变形材料梁的非线性分析方法,并且用ABAQUS有限元软件中的超弹性材料模型验证了所提出的方法的有效性。对典型的爆炸荷载作用下非线性弹性大变形材料梁的抗爆特性进行了分析,讨论了动力放大系数和材料性质及动荷载之间的关系。结果表明:与线弹性小变形材料相比,非线性弹性大变形材料具有优良的抗爆特性,结构的抗爆能力随结构变形的增大而显著提高。     

Entropy productions in granular materials

Granular materials display more abundant dissipation phenomena than ordinary materials. In this paper, a brief energy flow path with irreversible processes is illustrated, where the concept of granular temperature T_g, initially proposed for dilute systems, is extended to dense systems in order to quantify disordered force chain configurations. Additionally, we develop the concept of conjugate granular entropy s_g and its production equation. Our analyses find out that the granular entropy significantly undermined the elastic contact between particles, seriously affecting the transport coefficients in granular materials and creating new transport processes.     

Change of scale in granular materials

This paper analyzes relationships allowing global variables to be deduced from local ones in granular materials. This change of scale is now clear for static variables but rather confusing for kinematic variables. This paper compares the various established formulations for kinematic variables and also proposes original formulations. By running numerical simulations all the formulations in question are compared and the basic premises to define the change of scale for kinematic variables in granular media put forward.     

Incorporating the effects of fabric in the dilatant double shearing model for planar deformation of granular materials

The objective of this paper is to incorporate the effects of fabric and its evolution into the Dilatant Double Shearing Model [Mehrabadi, M.M., Cowin, S.C., 1978. Initial planar deformation of dilatant granular materials. J. Mech. Phys. Solids 26, 269–284] for granular materials in order to capture the anisotropic behavior and the complex response of granular materials in cyclic shear loading. An important consequence of considering the fabric is that one can have unequal shearing rates along the two slip directions. This property leads to the non-coaxiality of the principal axes of stress and strain rate, which is more appropriate for a material that exhibits initial and induced anisotropy. In addition, we employ a fabric-dependent elasticity tensor with orthotropic symmetry. The model developed in this paper also predicts one of the experimentally observed characteristics of granular materials: the gradual concentration of the contact normals towards the maximum principal stress direction.We implement the constitutive equations into ABAQUS/Explicit by writing a user material subroutine in order to predict the strength anisotropy of granular materials in a plane strain biaxial compression test and investigate the mechanical behavior of granular materials under the cyclic shear loading conditions. The predictions from this model show good quantitative agreement with the experiments of [Park, C.S., 1990. Anisotropy in deformation and strength properties of sands in plane strain compression, Masters Thesis, University of Tokyo; Park, C.S., Tatsuoka, F., 1994. Anisotropic strength and deformation of sands in plane strain compression. In: XIII ICSMFE, New Delhi, India; Okada, N., 1992. Energy dissipation in inelastic flow of cohesionless granular media. Ph.D. Thesis, University of California, San Diego].     

A return-map algorithm for general associative isotropic elasto-plastic materials in large deformation regimes

The present paper addresses a flexible solution algorithm for associative isotropic elasto-plastic materials, i.e. for materials whose elastic and plastic behaviors are described through an isotropic free-energy function, an isotropic yield function and an associative flow rule. The discussion is relative to a large deformation regime, while no hardening mechanisms are included. The algorithm is based on a combination of the operator split method and a return map scheme. Both the algorithm linearization and the requirements for the yield criterion convexity are discussed in detail. Finally, to show the algorithm flexibility and performance, the discussion is specialized to three yield criteria and some test problems are studied.     

Effects of large deformation and material nonlinearity on spherical indentation of hyperelastic soft materials

Up to now, the indentation of hyperelastic soft materials has not been completely understood. In this paper, the spherical indentation on hyperelastic soft solids was systematically investigated through theoretical analysis and finite element method (FEM). The validation and application of the Hertzian load-displacement relation for indentation of hyperelastic soft materials were clarified, the effects of large deformation and material nonlinearity on spherical indentation of hyperelastic soft materials were analyzed and discussed. It was found that the complicated indentation behaviors of hyperelastic soft solids mainly depended on the coupling interactions of large deformation and material nonlinearity. Besides, we proposed two new nonlinear elastic contact models to separate the effects of large deformation and material nonlinearity on spherical indentation responses of hyperelastic soft solids. Our efforts might help to enhance the understanding of hyperelastic indentation problems and provided necessary instructions for the mechanical characterization of hyperelastic soft materials.     

Kinematic hardening of granular materials

Summary Some recently developed constitutive equations (yield function, loading criteria, flow rule) for kinematic hardening of granular materials are discussed and some relevant material dependent parameters are determined on the basis of test results.

---

Kinematische Verfestigung granularer Stoffe

Übersicht Einige in jüngster Zeit entwickelte konstitutive Beziehungen (Fließbedingung, Belastungsbedingung und Fließregel) für kinematisch verfestigende granulare Stoffe werden diskutiert und einige relevante materialabhängige Parameter werden auf der Grundlage experimenteller Ergebnisse bestimmt.

     

Couette flow of granular materials

The flow of granular materials between rotating cylinders is studied using a continuum model proposed by Rajagopal and Massoudi (A method for measuring material moduli for granular materials: flow in an orthogonal rheometer, DOE/PETC/TR90/3, 1990). For a steady, fully developed condition, the governing equations are reduced to a system of coupled non-linear ordinary differential equations. The resulting boundary value problem is non-dimensionalized and is then solved numerically. The effect of material parameters, i.e., dimensionless numbers on the volume fraction and the velocity fields are studied.     

On flows of granular materials

This article reviews the behavior of materials made up of a large assemblage of solid particles under rapid and quasi static deformations. The focus is on flows at relatively high concentrations and for conditions when the interstitial fluid plays an insignificant role. The momentum and energy exchange processes are then primarily governed by interparticle collisions and Coulomb-type frictional contact. We first discuss some physical behavior —dilatancy, internal friction, fluidization and particle segregation — that are typical to the understanding of granular flows. Bagnold's seminal Couette flow experiments and his simple stress analysis are then used to motivate the first constitutive theories that use a microstructural variable — the fluctuation energy or granular temperature — governing the subscale fluctuating motion. The kinetic theories formalize the derivation of the field equations of bulk mass, momentum and energy, and permit derivation of constitutive relations for stress, flux of fluctuation energy and its dissipation rate for simple particle assemblages and when frictional rubbing contact can be ignored. These statistical considerations also show that formulation of boundary conditions needs special attention. The frictional-collisional constitutive behavior in which both Coulomb-type rubbing contact and collisional encounters are significant are discussed. There is as yet no rigorous formulation. We finally present a phenomenological approach that describes rapid flows of granular materials under simultaneous transport of heat and close with a summary of stability analyses of the basic flow down an inclined plane.Dedicated to Professor Dr.-Ing. Franz Gustav Kollmann on the occasion of his sixtieth brithday     

Analysis of local behaviour in granular materials

A local scale, called the meso-scale, has recently been introduced to the multi-scale approach for 2D granular materials. This local scale is defined at the level of meso-domains enclosed by particles in contact. Stress and strain have been defined at this local scale, and their relation with the local structure has been studied. The purpose of this paper is to analyse the behaviour of granular materials at the meso-scale, i.e. the stress–strain–structure relationship at this scale. Analyses are performed on a 2D numerical granular sample subjected to a biaxial compression test and simulated with the Discrete Element Method (DEM). The sample is quite dense and it is loaded at a relatively low strain rate so that the state of the sample can be considered as being quasi-static. The size of sub-domains in the sample varies largely from 3 to 12 particles. It is shown that the evolution of the internal state of the sample corresponds, at the meso-scale, to a clear evolution of the quantity of meso-domains oriented in different directions. In addition, the behaviour of meso-domains is highly governed by their orientation rather than their density, especially for the strongly elongated meso-domains: the meso-domains oriented in the compression (resp. extension) direction behave like a dense (resp. loose) granular material.