Integration algorithms of elastoplasticity for ceramic powder compaction

Inelastic deformation of ceramic powders (and of a broad class of rock-like and granular materials), can be described with the yield function proposed by Bigoni and Piccolroaz (2004, Yield criteria for quasibrittle and frictional materials. Int. J. Solids and Structures, 41, 2855-2878). This yield function is not defined outside the yield locus, so that 'gradient-based' integration algorithms of elastoplasticity cannot be directly employed. Therefore, we propose two ad hoc algorithms: (i.) an explicit integration scheme based on a forward Euler technique with a 'centre-of-mass' return correction and (ii.) an implicit integration scheme based on a 'cutoff-substepping' return algorithm. Iso-error maps and comparisons of the results provided by the two algorithms with two exact solutions (the compaction of a ceramic powder against a rigid spherical cup and the expansion of a thick spherical shell made up of a green body), show that both the proposed algorithms perform correctly and accurately.Comment: 21 pages. Journal of the European Ceramic Society, 201

How granular materials deform in quasistatic conditions

Based on numerical simulations of quasistatic deformation of model granular materials, two rheological regimes are distinguished, according to whether macroscopic strains merely reflect microscopic material strains within the grains in their contact regions (type I strains), or result from instabilities and contact network rearrangements at the microscopic level (type II strains). We discuss the occurrence of regimes I and II in simulations of model materials made of disks (2D) or spheres (3D). The transition from regime I to regime II in monotonic tests such as triaxial compression is different from both the elastic limit and from the yield threshold. The distinction between both types of response is shown to be crucial for the sensitivity to contact-level mechanics, the relevant variables and scales to be considered in micromechanical approaches, the energy balance and the possible occurrence of macroscopic instabilitie

An elastoplastic framework for granular materials becoming cohesive through mechanical densification. Part I - small strain formulation

Mechanical densification of granular bodies is a process in which a loose material becomes increasingly cohesive as the applied pressure increases. A constitutive description of this process faces the formidable problem that granular and dense materials have completely different mechanical behaviours (nonlinear elastic properties, yield limit, plastic flow and hardening laws), which must both be, in a sense, included in the formulation. A treatment of this problem is provided here, so that a new phenomenological, elastoplastic constitutive model is formulated, calibrated by experimental data, implemented and tested, that is capable of describing the transition between granular and fully dense states of a given material. The formulation involves a novel use of elastoplastic coupling to describe the dependence of cohesion and elastic properties on the plastic strain. The treatment falls within small strain theory, which is thought to be appropriate in several situations; however, a generalization of the model to large strain is provided in Part II of this paper.Comment: 42 pages, 27 figure

Development of Stresses in Cohesionless Poured Sand

The pressure distribution beneath a conical sandpile, created by pouring sand from a point source onto a rough rigid support, shows a pronounced minimum below the apex (`the dip'). Recent work of the authors has attempted to explain this phenomenon by invoking local rules for stress propagation that depend on the local geometry, and hence on the construction history, of the medium. We discuss the fundamental difference between such approaches, which lead to hyperbolic differential equations, and elastoplastic models, for which the equations are elliptic within any elastic zones present .... This displacement field appears to be either ill-defined, or defined relative to a reference state whose physical existence is in doubt. Insofar as their predictions depend on physical factors unknown and outside experimental control, such elastoplastic models predict that the observations should be intrinsically irreproducible .... Our hyperbolic models are based instead on a physical picture of the material, in which (a) the load is supported by a skeletal network of force chains ("stress paths") whose geometry depends on construction history; (b) this network is `fragile' or marginally stable, in a sense that we define. .... We point out that our hyperbolic models can nonetheless be reconciled with elastoplastic ideas by taking the limit of an extremely anisotropic yield condition.Comment: 25 pages, latex RS.tex with rspublic.sty, 7 figures in Rsfig.ps. Philosophical Transactions A, Royal Society, submitted 02/9

Jamming and Stress Propagation in Particulate Matter

We present simple models of particulate materials whose mechanical integrity arises from a jamming process. We argue that such media are generically "fragile", that is, they are unable to support certain types of incremental loading without plastic rearrangement. In such models, fragility is naturally linked to the marginal stability of force chain networks (granular skeletons) within the material. Fragile matter exhibits novel mechanical responses that may be relevant to both jammed colloids and cohesionless assemblies of poured, rigid grains.Comment: LATEX, 3 Figures, elsart.cls style file, 11 page

Quasistatic rheology and the origins of strain

Features of rheological laws applied to solid-like granular materials are recalled and confronted to microscopic approaches via discrete numerical simulations. We give examples of model systems with very similar equilibrium stress transport properties -- the much-studied force chains and force distribution -- but qualitatively different strain responses to stress increments. Results on the stability of elastoplastic contact networks lead to the definition of two different rheological regimes, according to whether a macroscopic fragility property (propensity to rearrange under arbitrary small stress increments in the thermodynamic limit) applies. Possible consequences are discussed.Comment: Published in special issue of "Comptes-Rendus Physique" on granular material

Stiffness pathologies in discrete granular systems: bifurcation, neutral equilibrium, and instability in the presence of kinematic constraints

The paper develops the stiffness relationship between the movements and forces among a system of discrete interacting grains. The approach is similar to that used in structural analysis, but the stiffness matrix of granular material is inherently non-symmetric because of the geometrics of particle interactions and of the frictional behavior of the contacts. Internal geometric constraints are imposed by the particles' shapes, in particular, by the surface curvatures of the particles at their points of contact. Moreover, the stiffness relationship is incrementally non-linear, and even small assemblies require the analysis of multiple stiffness branches, with each branch region being a pointed convex cone in displacement-space. These aspects of the particle-level stiffness relationship gives rise to three types of micro-scale failure: neutral equilibrium, bifurcation and path instability, and instability of equilibrium. These three pathologies are defined in the context of four types of displacement constraints, which can be readily analyzed with certain generalized inverses. That is, instability and non-uniqueness are investigated in the presence of kinematic constraints. Bifurcation paths can be either stable or unstable, as determined with the Hill-Bazant-Petryk criterion. Examples of simple granular systems of three, sixteen, and sixty four disks are analyzed. With each system, multiple contacts were assumed to be at the friction limit. Even with these small systems, micro-scale failure is expressed in many different forms, with some systems having hundreds of micro-scale failure modes. The examples suggest that micro-scale failure is pervasive within granular materials, with particle arrangements being in a nearly continual state of instability

Numerical simulations of materials with micro-structure : limit analysis and homogenization techniques

Continuum-based numerical methods have played a leading role in the numerical solution of problems in soil and rock mechanics. However, for stratified soils and fractured rocks, a continuum assumption often leads to difficult parameters to define and over-simplified geometry to be realistic. In such cases, approaches that consider the micro-structure of the material can be adopted. In this paper, two of such approaches are detailed, namely limit analysis incorporating fractures and individual blocks, and elastoplastic homogenization of layered soils