Micromechanical study of elastic moduli of loose granular materials

In micromechanics of the elastic behaviour of granular materials, the macro-scale continuum elastic moduli are expressed in terms of micro-scale parameters, such as coordination number (the average number of contacts per particle) and interparticle contact stiffnesses in normal and tangential directions. It is well-known that mean-field theory gives inaccurate micromechanical predictions of the elastic moduli, especially for loose systems with low coordination number. Improved predictions of the moduli are obtained here for loose two-dimensional, isotropic assemblies. This is achieved by determining approximate displacement and rotation fields from the force and moment equilibrium conditions for small sub-assemblies of various sizes. It is assumed that the outer particles of these sub-assemblies move according to the mean field. From the particle displacement and rotation fields thus obtained, approximate elastic moduli are determined. The resulting predictions are compared with the true moduli, as determined from the discrete element method simulations for low coordination numbers and for various values of the tangential stiffness (at fixed value of the normal stiffness). Using this approach, accurate predictions of the moduli are obtained, especially when larger sub-assemblies are considered. As a step towards an analytical formulation of the present approach, it is investigated whether it is possible to replace the local contact stiffness matrices by a suitable average stiffness matrix. It is found that this generally leads to a deterioration of the accuracy of the predictions. Many micromechanical studies predict that the macroscopic bulk modulus is hardly influenced by the value of the tangential stiffness. It is shown here from the discrete element method simulations of hydrostatic compression that for loose systems, the bulk modulus strongly depends on the stiffness ratio for small stiffness ratios.     

Effective elastic moduli of two dimensional solids with distributed cohesive microcracks

Effective elastic properties of a defected solid with distributed cohesive micro-cracks are estimated based on homogenization of the Dugdale–Bilby–Cottrell–Swinden (Dugdale–BCS) type micro-cracks in a two dimensional elastic representative volume element (RVE).Since the cohesive micro-crack model mimics various realistic bond forces at micro-scale, a statistical average of cohesive defects can effectively represent the overall properties of the material due to bond breaking or crack surface separation in small scale. The newly proposed model is distinctive in the fact that the resulting effective moduli are found to be pressure sensitive.     

Nonlinear granular micromechanics model for multi-axial rate-dependent behavior

Constitutive equations for class of materials that possess granular microstructure can be effectively derived using granular micromechanics approach. The stress–strain behavior of such materials depends upon the underlying grain scale mechanisms that are modeled by using appropriate rate-dependent inter-granular force–displacement relationships. These force–displacement functions are nonlinear and implicit evolutions equations. The numerical solution of such equation under applied overall stress or strain loading can entail significant computational expense. To address the computations issue, an efficient explicit time-integration scheme has been derived. The developed model is then utilized to predict primary, secondary and tertiary creep as well as rate-dependent response under tensile and compressive loads for hot mix asphalt. Further, the capability of the derived model to describe multi-axial behavior is demonstrated through generations of biaxial time-to-creep failure envelopes and rate-dependent failure envelopes under monotonic biaxial and triaxial loading. The advantage of the approach presented here is that we can predict the multi-axial effects without resorting to complex phenomenological modeling.     

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.     

On the quasistatic effective elastic moduli for elastic waves in three-dimensional phononic crystals

Effective elastic moduli for 3D solid–solid phononic crystals of arbitrary anisotropy and oblique lattice structure are formulated analytically using the plane-wave expansion (PWE) method and the recently proposed monodromy-matrix (MM) method. The latter approach employs Fourier series in two dimensions with direct numerical integration along the third direction. As a result, the MM method converges much quicker to the exact moduli in comparison with the PWE as the number of Fourier coefficients increases. The MM method yields a more explicit formula than previous results, enabling a closed-form upper bound on the effective Christoffel tensor. The MM approach significantly improves the efficiency and accuracy of evaluating effective wave speeds for high-contrast composites and for configurations of closely spaced inclusions, as demonstrated by three-dimensional examples.     

Breakage mechanics—Part II: Modelling granular materials

The compression of granular materials has been traditionally modelled with the limitations of classical elasto-plasticity. The energy was implicitly assumed to dissipate from the frictional interaction of particles. However, the fact that brittle granular materials crush suggests that energy must also be dissipated from the fracturing of the grains, as in fracture mechanics. The concept of breakage as a thermomechanical internal variable was introduced in Part I [Einav, I., 2006. Breakage mechanics—Part I: theory. J. Mech. Phys. Solids 00,000-000] to describe the fracturing mechanisms. The theory allows to treat ideal theoretical materials that undergo dissipation purely from breakage with no other mechanism allowed for the energy consumption. However, as accounted for in elasto-plasticity, dissipation must also occur from the frictional rearrangement of grains. The combination of the two dissipative mechanisms of breakage and plasticity must therefore be investigated, as we do in this paper. Those two mechanisms are generally coupled, in the sense that one inevitably appears when the other develops. Plastic dissipation emerges as a by-product of breakage dissipation because after grains crush, local rearrangement must occur. This scenario may be termed an ‘active breakage mechanism’, and typifies compression deformations. In shear the plastic dissipation is dominant but breakage appears inevitably from grains abrasion. This scenario may be called a ‘passive breakage mechanism’. Based on the coupling assumption, models are developed for granular materials. In particular, we show that in compression isotropic hardening of sands may appear without involving plastic strains, i.e., independent of frictional dissipation. This interpretation of hardening is different from the one used in classical critical state soil mechanics. However, frictional dissipation leads to plastic straining that are necessary for the models to be predictive in unloading.     

Estimate of effective elastic moduli with microcrack interaction effects

A method is developed for estimating the effects of microdefect interaction on the effective elastic properties of heterogeneous solids. An effective medium is defined to calculate the global effective elastic moduli of brittle materials weakened by distributed microcracks. Each microcrack is assumed to be embedded in an effective medium, the compliance of which is obtained from the dilute concentration method without accounting for interaction. The present scheme requires no iteration; it can account for microcrack interaction with sufficient accuracy. Analytical solutions are given for several two- and three-dimension problems with and without anisotropy.     

Effective elastic properties of porous materials: Homogenization schemes vs experimental data

In this paper, we focus on the prediction of elastic moduli of isotropic porous materials made of a solid matrix having a Poisson's ratio vm of 0.2. We derive simple analytical formulae for these effective moduli based on well-known Mean-Field Eshelby-based Homogenization schemes. For each scheme, we find that the normalized bulk, shear and Young's moduli are given by the same form depending only on the porosity p. The various predictions are then confronted with experimental results for the Young's modulus of expanded polystyrene (EPS) concrete. The latter can be seen as an idealized porous material since it is made of a bulk cement matrix, with Poisson's ratio 0.2, containing spherical mono dispersed EPS beads. The Differential method predictions are found to give a very good agreement with experimental results. Thus, we conclude that when vm=0.2, the normalized effective bulk, shear and Young's modulus of isotropic porous materials can be well predicted by the simple form (1 − p)² for a large range of porosity p ranging between 0 and 0.56.     

Micromechanical analysis of friction anisotropy in rough elastic contacts

Computational contact homogenization approach is applied to study friction anisotropy resulting from asperity interaction in elastic contacts. Contact of rough surfaces with anisotropic roughness is considered with asperity contact at the micro scale being governed by the isotropic Coulomb friction model. Application of a micro-to-macro scale transition scheme yields a macroscopic friction model with orientation- and pressure-dependent macroscopic friction coefficient. The macroscopic slip rule is found to exhibit a weak non-associativity in the tangential plane, although the slip rule at the microscale is associated in the tangential plane. Counterintuitive effects are observed for compressible materials, in particular, for auxetic materials.     

Effective elastic moduli of granular media

Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 91–96, January–February, 1990.     

On the definition of the stress tensor in granular media

This paper investigates the definition of the stress tensor within a granular assembly, when inertial effects are likely to occur. It is shown that the stress tensor can be expressed as a sum of two terms. A first term corresponds to the standard definition of the stress, according to the Love–Weber formula; this term is related to the contact forces existing within adjoining particles. A second term accounts for dynamic effects related to rotation velocities and accelerations of the particles. These results are checked from discrete numerical simulations in order to examine in which context the contribution of inertial effects should not be omitted. With this aim, the simulation of a granular specimen collapse and then a silo discharge is considered.     

Crushing of particles in idealised granular assemblies

Four idealised assemblies of equally sized spherical particles are subjected to a range of macroscopic compressive principal stresses and the contact forces on individual particles are determined. For each set of contact forces the stress fields within individual particles are studied. A failure criterion for brittle materials is imposed and indicates that crushing (or rupture) occurs when the maximum contact force reaches a threshold particle strength value, irrespective of the presence and magnitude of other lesser contact forces acting on the particle and the material properties of the particle. Combining the crushing mechanism with an assembly instability mechanism enables failure surfaces to be drawn in the three-dimensional stress space. A simple spatial averaging technique has been applied to the failure surfaces to remove the effects of assembly anisotropies. Sections of the failure surfaces on π planes have similarities to those commonly used in sand modelling.     

A micromechanics-based elastoplastic damage model for granular materials at low confining pressure

This paper is devoted to the formulation of a micromechanics-based constitutive model for granular materials under relatively low confining pressure. The constitutive formulation is performed within the general framework of homogenization for granular materials. However, new rigorous stress localization laws are proposed. Some local constitutive relations are established under the consideration of irreversible thermodynamics. Macroscopic plastic deformation is obtained by considering local plastic sliding in a limit number of families of contact planes. The plastic sliding at each contact plane is described by a non-associated plastic flow rule, taking into account pressure sensitivity and normal dilatancy. Nonlinear elastic deformation related to progressive compaction of contacts is also taken into account. Material softening is described by involving damage process related to degradation of microstructure fabric. The proposed model is applied to some typical granular materials (sands). The numerical predictions are compared with experimental data.     

Shock tube study of the dynamical behavior of granular materials

The purpose of the present study is to clarify both the compression phenomenon and the gas filtration effect that take place inside a granular medium when it is dynamically loaded by a shock wave. In order to measure the pore pressure and the total stress at different locations along the granular medium, pressure transducers were placed along the side-wall and at the end-wall of the shock tube test section, which was filled with the granular material. In order to elucidate the gas filtration effect, the results of two experiments with identical granular media but with and without filtration were compared. The gas filtration was eliminated by means of a thin plastic film, which was placed at the front edge of the granular medium. Based on this comparison quantitative information on the gas filtration and its role in the stress formation inside granular media of different material and length was obtained. Furthermore, curves of the dynamic compression and the Young moduli of the granular medium for the range of the operating conditions were reconstructed.     

Effective elastic moduli of an inhomogeneous medium with cracks

In the present paper, the effective elastic moduli of an inhomogeneous medium with cracks are derived and obtained by taking into account its microstructural properties which involve the shape, size and distribution of cracks and the interaction between cracks. Numerical results for the periodic microstructure of different dimensions are presented. From the results obtained, it can be found that the distribution of cracks has a significant effect on the effective elastic moduli of the material. The project supported by the National Education Committee for Doctor     

Effective elastic moduli of inhomogeneous solids by embedded cell model

An embedded cell model is presented to obtain the effective elastic moduli for three-dimensional two-phase composites which is an exact analytic formula without any simplified approximation and can be expressed in an explicit form. For the different cells such as spherical inclusions and cracks surrounded by sphere and oblate ellipsoidal matrix, the effective elastic moduli are evaluated and the results are compared with those from various micromechanics models. These results show that the present model is direct, simple and efficient to deal with three-dimensional two-phase composites. The project supported by the National Natural Science Foundation of China (No. 19704100) and the National Natural Science Foundation of Chinese Academy of Sciences (No. KJ951-1-201)     

Assessing continuum postulates in simulations of granular flow

Continuum mechanics relies on the fundamental notion of a mesoscopic volume “element” in which properties averaged over discrete particles obey deterministic relationships. Recent work on granular materials suggests that a continuum law may be inapplicable, revealing inhomogeneities at the particle level, such as force chains and slow cage breaking. Here, we analyze large-scale three-dimensional discrete-element method (DEM) simulations of different granular flows and show that an approximate “granular element” defined at the scale of observed dynamical correlations (roughly three to five particle diameters) has a reasonable continuum interpretation. By viewing all the simulations as an ensemble of granular elements which deform and move with the flow, we can track material evolution at a local level. Our results confirm some of the hypotheses of classical plasticity theory while contradicting others and suggest a subtle physical picture of granular failure, combining liquid-like dependence on deformation rate and solid-like dependence on strain. Our computational methods and results can be used to guide the development of more realistic continuum models, based on observed local relationships between average variables.     

Micromechanical investigation of granular ratcheting using a discrete model of polygonal particles

We use a two-dimensional model of polygonal particles to investigate granular ratcheting. Ratcheting is a long-term response of granular materials under cyclic loading, where the same amount of permanent deformation is accumulated after each cycle. We report on ratcheting for low frequencies and extremely small loading amplitudes. The evolution of the sub-network of sliding contacts allows us to understand the micromechanics of ratcheting. We show that the contact network evolves almost periodically under cyclic loading as the sub-network of the sliding contacts reaches different stages of anisotropy in each cycle. Sliding contacts lead to a monotonic accumulation of permanent deformation per cycle in each particle. The distribution of these deformations appears to be correlated in form of vortices inside the granular assembly.     

Micromechanical investigation of granular ratcheting using a discrete model of polygonal particles

We use a two-dimensional model of polygonal particles to investigate granular ratcheting. Ratcheting is a long-term response of granular materials under cyclic loading, where the same amount of permanent deformation is accumulated after each cycle. We report on ratcheting for low frequencies and extremely small loading amplitudes. The evolution of the sub-network of sliding contacts allows us to understand the micromechanics of ratcheting. We show that the contact network evolves almost periodically under cyclic loading as the sub-network of the sliding contacts reaches different stages of anisotropy in each cycle. Sliding contacts lead to a monotonic accumulation of permanent deformation per cycle in each particle. The distribution of these deformations appears to be correlated in form of vortices inside the granular assembly.