A discrete element analysis of elastic properties of granular materials

The elastic properties of a regular packing of spheres with different tolerances were evaluated using the discrete element method to elucidate the mechanisms behind the discrepancies between laboratory experiments and theoretical predictions of the classic Hertz-Mindlin contact law. The simulations indicate that the elastic modulus of the packing is highly dependent on the coordination number and the magnitude and distribution of contact normal forces, and this dependence is macroscopically reflected as the influence of confining pressure and void ratio. The increase of coordination number and the uniformity of contact normal forces distribution with increasing confining pressure results in the stress exponent $n$ for elastic modulus being higher than 1/3 as predicted by the Hertz-Mindlin law. Furthermore, the simulations show that Poisson’s ratio of a granular packing is not a constant as commonly assumed, but rather it decreases as confining pressure increases. The variation of Poisson’s ratio appears to be a consequence of the increase of the coordination number rather than the increase of contact normal forces with confining pressure.     

On the thermodynamic requirement of elastic stiffness anisotropy in isotropic materials

In general, thermodynamic admissibility requires isotropic materials develop reversible deformation induced anisotropy (RDIA) in their elastic stiffnesses. Taking the elastic potential for an isotropic material to be a function of the strain invariants, isotropy of the elastic stiffness is possible under distortional loading if and only if the bulk modulus is independent of the strain deviator and the shear modulus is constant. Previous investigations of RDIA have been limited to applications in geomechanics where material non-linearity and large deformations are commonly observed. In the current paper, the degree of RDIA in other materials is investigated. It is found that the resultant anisotropy in materials whose strength does not vary appreciably with pressure, such as metals, is negligible, but in materials whose strength does vary with pressure, the degree of RDIA can be significant. Algorithms for incorporating RDIA in a classical elastic-plastic model are provided.     

Moisture influence on the resilient deformation behaviour of unbound granular materials

The influence of moisture on the resilient deformation properties of unbound granular materials was investigated based on repeated load triaxial tests. Results showed that the resilient modulus (M_R) decreased with increasing moisture for a relatively low number of load cycles (N) where the deformation behaviour was mostly resilient with a negligible amount of associated accumulated permanent deformation (PD). Modelling attempts on this behaviour were quite satisfactory. Furthermore, the M_R showed an increasing trend with increasing moisture, up to the optimum, when the N was relatively large with a significant amount of accumulated PD. Above the optimum, the M_R generally decreased. Further investigation suggested that moisture aided the post-compaction (PC) and possible particle rearrangement that resulted in the increased PD and increased M_R. The existing model did not work in this case indicating that the effect of PC on M_R should be considered in modelling.     

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.     

A thickness mode acoustic wave sensor for measuring interface stiffness between two elastic materials

We studied thickness vibration of 2 elastic layers with an elastic interface mounted on a plate piezoelectric resonator. The effect of the interface elasticity on resonant frequencies was examined. The result obtained suggests an acoustic wave sensor for measuring the elastic property of an interface between 2 materials.     

Analysis of the influence of crushing on the behavior of granular materials under shear

The behavior of granular materials mainly depends on the mechanical and engineering properties of particles in its structural matrix. Crushing or breakage of granular materials under compression or shear occurs when the energy available is sufficient to overcome the resistance of the material. Relatively little systematic research has been conducted regarding how to evaluate or quantify particle crushing and how it effects the engineering properties of the granular materials. The aim of this study is to investigate the effect of crushing on the bulk behavior of granular materials by using manufactured granular materials (MGM) rather than using a naturally occurring cohesionless granular material. MGM allow changing only one particle parameter, namely the “crushing strength”. Four different categories of MGM (with different crushing strength) are used to study the effect on the bulk shear strength, stiffness modulus, friction and dilatancy angle “engineering properties”. A substantial influence on the stress–strain behavior and engineering properties of granular materials is observed. Higher confining stress causes some non-uniformity (strong variations/jumps) in volumetric strain and a constant volumetric strain is not always observed under large shear deformations due to crushing, i.e. there is no critical state with flow regime (with constant volumetric strain).     

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.     

Laser-induced surface acoustic waves for evaluation of elastic stiffness of plasma sprayed materials

The elastic properties of plasma sprayed deposits have been evaluated using a laser-excited surface acoustic wave (SAW) technique and an inversion processing analysis. The SAWs including Lamb and Rayleigh waves were generated in plasma sprayed NiCoCrAlY and ZrO₂, respectively, and their group velocity dispersions were used to determine the elastic properties (i.e.Young's modulus, Poison's ratio and density) of the deposits. Estimated elastic moduli from the velocity dispersions of A₀-mode Lamb waves are in the range of 40–140 GPa for the deposits, which are much lower than the values 220–240 GPa of the comparable dense materials. The dramatic reductions in modulus and density of ZrO₂ deposit have been attributed to the presence of high porosity and particularly microcracks. Moreover, this study has emphasized on exploiting the applicability of each kind of the SAWs for the elastic property evaluation of different sprayed materials. Both Lamb and Rayleigh wave dispersions are useful for the estimation of APS and VPS-deposited NiCrAlY, but S₀-Lamb and Rayleigh waves are exceptional for that of sprayed ZrO₂, because of its characterization of high acoustic attenuation and inconsequent displacement across the weak bonded interface of ZrO₂ and substrate.     

DEM analysis of the influence of the intermediate stress ratio on the critical-state behaviour of granular materials

The critical-state response of granular assemblies composed of elastic spheres under generalised three-dimensional loading conditions was investigated using the discrete element method (DEM). Simulations were performed with a simplified Hertz–Mindlin contact model using a modified version of the LAMMPS code. Initially isotropic samples were subjected to three-dimensional stress paths controlled by the intermediate stress ratio, \(b=[(\sigma '_{2}-\sigma '_{3})/\) \((\sigma '_{1}-\sigma '_{3})]\) . Three types of simulation were performed: drained (with \(b\) -value specified), constant volume and constant mean effective stress. In contrast to previous DEM observations, the position of the critical state line is shown to depend on \(b\) . The data also show that, upon shearing, the dilatancy post-peak increases with increasing \(b\) , so that at a given mean effective stress, the void ratio at the critical state increases systematically with \(b\) . Four commonly-used three-dimensional failure criteria are shown to give a better match to the simulation data at the critical state than at the peak state. While the void ratio at critical state is shown to vary with \(b\) , the coordination number showed no dependency on \(b\) . The variation in critical state void ratios at the same \(p'\) value is apparently related to the directional fabric anisotropy which is clearly sensitive to \(b\) .     

Stresses in granular materials

When circularly polarised light is passed through a granular material under boundary stresses patterns—‘light stripes’—are seen in the resulting images which have been traditionally associated with the directions of major principal stresses in the equivalent continuum. In this paper the passage of polarised light through a single spherical particle under stress is studied experimentally and analytically. The effect of placing the particle within a layer of particles, a layer of thickness 2–3 particles, and within a mass of particles is investigated experimentally. The appearance of light stripes is a visual reinforcement of effects seen at the particle level provided the level of stresses in individual particles is low. The implications for quantitative photoelastic interpretation of granular media are discussed.     

On the influence of inter-particle friction and dilatancy in granular materials: a numerical analysis

Mechanical behavior of granular soils is a classic research realm but still yet not completely understood as it can be influenced by a large number of factors, including confining pressure, soil density, loading conditions, and anisotropy of soil etc. Traditionally granular materials are macroscopically regarded as continua and their particulate and discrete nature has not been thoroughly considered although many researches indicate the macro mechanical behavior closely depends on the micro-scale characteristics of particles. This paper presents a DEM (discrete element method)-based micromechanical investigation of inter-particle friction effects on the behavior of granular materials. In this study, biaxial DEM simulations are carried out under both ‘drained’ and ‘undrained’ (constant volume) conditions. The numerical experiments employ samples having similar initial isotropic fabric and density, and the same confining pressure, but with different inter-particle friction coefficient. Test results show that the inter-particle friction has a substantial effect on the stress-strain curve, peak strength and dilatancy characteristics of the granular assembly. Clearly, it is noted that apart from the inter-particle friction, the shear resistance is also contributed to the dilation and the particle packing and arrangements. The corresponding microstructure evolutions and variations in contact properties in the particulate level are also elaborated, to interpret the origin of the different macro-scale response due to variations in the inter-particle friction.     

The influence of particle geometry and the intermediate stress ratio on the shear behavior of granular materials

The behavior of granular materials is very complex in nature and depends on particle shape, stress path, fabric, density, particle size distribution, amongst others. This paper presents a study of the effect of particle geometry (aspect ratio) on the mechanical behaviour of granular materials using the discrete element method (DEM). This study discusses 3D DEM simulations of conventional triaxial and true triaxial tests. The numerical experiments employ samples with different particle aspect ratios and a unique particle size distribution (PSD). Test results show that both particle aspect ratio (AR) and intermediate stress ratio \((b=({\upsigma }_{2}'-{\upsigma }_{3}')/({\upsigma }_{1}'-{\upsigma }_{3}'))\) affect the macro- and micro-scale responses. At the macro-scale, the shear strength decreases with an increase in both aspect ratio and intermediate stress ratio b values. At the micro-scale level, the fabric evolution is also affected by both AR and b. The results from DEM analyses qualitatively agree with available experimental data. The critical state behaviour and failure states are also discussed. It is observed that the position of the critical state loci in the compression \((e-p')\) space is only slightly affected by aspect ratio (AR) while the critical stress ratio is dependent on both AR and b. It is also demonstrated that the influence of the aspect ratio and the intermediate stress can be captured by micro-scale fabric evolutions that can be well understood within the framework of existing critical state theories. It is also found that for a given stress path, a unique critical state fabric norm is dependent on the particle shape but is independent of critical state void ratio.     

Equivalent normal stiffness of the ball in granular dynamics simulations

In the paper a concept of an equivalent normal ball stiffness based on averaging the work done by the nonlinear Hertz force during the impact is introduced and a numerical assessment is made on its efficiency in simulating collinear collisions. The systems of balls investigated, ranging from 2 to 30, are considered to be coupled and conservative. The energy and momentum conservation principles are used to assess the accuracy of simulation results. The effect of the time step for both linear and nonlinear models is investigated and it is shown that the linear model allows an increased time step compared to a nonlinear one while meeting energy and momentum conservation requirements with the same accuracy. In the paper also the pattern of break up for both models and different number of balls is investigated. It is found that for both models the pattern is the same: the balls are disconnected one at the time with constant rate and this rate does not depend on the number of balls. The financial assistance provided by the National Science and Engineering Research Council of Canada is gratefully acknowledged.     

Effect of elastic stiffness on metal forming plasticity

Summary A parametric study is carried out to investigate the influence of the elastic stiffness on the deformation history of a rectangular specimen under tension. Plane-strain condition is assumed, and two solutions differing by their elastic moduli are compared in the finite-strain range. The material behaviour is taken as elastic-viscoplastic and diffuse type necking deformations are triggered by an initial geometric defect. The resulting nonhomogeneous stress states are shown to be sensitive to the value of the elastic stiffness.     

Mixtures of granular materials

Balance laws are given for a mixture of granular materials of a type described by Goodman and Cowin. Constitutive equations are given for the case of two dry granular constituents, and consequences of the entropy principle are found.     

Creep of granular materials

This paper examines the creep of brittle granular materials subjected to one-dimensional compression. One-dimensional creep tests were performed on aggregates of brittle pasta and compared with the behaviour of sand at much higher stress levels. It was found that for both materials, creep strain is proportional to the logarithm of time. One possible mechanism for creep is particle crushing. However, it is usually difficult to measure changes in the particle size distribution during creep because the fines produced are so small, and the mass of fines is too small to measure accurately unless creep is permitted for a very long time. However, for pasta, the particle fragments produced are large, and it is found that particle crushing does occur during creep for 24 hours. This is consistent with the proposition that the behaviour of all brittle granular materials is essentially the same. A micro mechanical argument is then summarised which predicts that creep strain should be proportional to log time.     

A hyperbolic well-posed model for the flow of granular materials

A plasticity model for the flow of granular materials is presented which is derived from a physically based kinematic rule and which is closely related to the double-shearing model, the double-sliding free-rotating model and also to the plastic-potential model. All of these models incorporate various notions of the concept of rotation-rate and the crucial idea behind the model presented here is that it identifies this rotation-rate with a property associated with a Cosserat continuum, namely, the intrinsic spin. As a consequence of this identification, the stress tensor may become asymmetric. For simplicity, in the analysis presented here, the material parameters are assumed to be constant. The central results of the paper are that (a) the model is hyperbolic for two-dimen-Specifically, sional steady-state flows in the inertial regime and (b) the model possesses a domain of linear well-posedness. it is proved that incompressible flows are well-posed.     

A plasticity theory for the kinematics of ideal granular materials

A plasticity theory is formulated for ideal granular materials. The material is assumed to be isotropic and incompressible. All the equations are expressed in the form of 3-dimensional tensor equations. This theory, is based on a new interpretation of the associated flow rule with associated constraints of deformation. The characteristic surface is analyzed to see the possible discontinuity of the deformation. Then the elastic Strain is incorporated and the singular wave propagation is analyzed to determine the possible velocities and directions of propagation.