Constitutive relations for the isotropic deformation of frictionless packings of polydisperse spheres

The isotropic compression of polydisperse packings of frictionless spheres is modeled with the Discrete Element Method (DEM). The evolution of coordination number, fraction of rattlers, isotropic fabric, and pressure (isotropic stress) is reported as function of volume fraction for different system parameters. The power law relationship, with power $\approx 1/2$, between coordination number and volume fraction is confirmed in the jammed state for a broad range of volume fractions and for different (moderate) polydispersities. The polydispersity in the packing causes a shift of the critical volume fraction, i.e., more heterogeneous packings jam at higher volume fractions. Close to jamming, the coordination number and the jamming volume fraction itself depend on both history and rate. At larger densities, neither the deformation history nor the loading rate have a significant effect on the evolution of the coordination number.Concerning the fabric tensor, comparing our DEM results to theoretical predictions, good agreement for different polydispersities is observed. An analytical expression for the pressure as function of isotropic (volumetric) strain is proposed for polydisperse packings, based on the assumption of uniform deformation. We note that, besides the implicit proportionality to contact number density (or fabric), no single power-law is evidenced in the relation between pressure and isotropic strain. However, starting from zero pressure at the jamming point, a linear term with a quadratic correction describes the stress evolution rather well for a broad range of densities and for various polydispersities. Finally, an incremental evolution equation is proposed for both fabric and stress, as function of isotropic strain, and involving the coordination number and the fraction of rattlers, as starting point for further studies involving anisotropic deformations.     

基于球形颗粒几何排列的离散元试样高效生成方法

在球体离散元数值模拟中,颗粒的初始排列状态是影响计算效率和计算结果的重要环节。本文采用前进面几何构造算法,提出了一种基于网格搜索的球形颗粒随机排列高效算法。通过求解空间三边方程,满足了粒径设置的任意大小的颗粒依次置入前进面的外侧,并与构成前进面的三个颗粒相互接触。为获得高体积分数的颗粒簇,该算法允许颗粒改变其粒径大小。采用颗粒网格化方法可以简化前进面的搜索,并由此提高排列效率。通过计算平均配位数、体积分数和二阶结构张量的特征值,对不同粒径比下得到的立方体试样进行了分析,得到试样配位数及体积分数均随着粒径比的增大而增大,且得到的试样为各向同性。此外,空间网格的大小和初始颗粒的生成点对随机排列的效率均会产生显著的影响。最后,对非规则铁路道砟进行了精细构造及压碎模拟,发现DEM模拟得到的应力-应变曲线与试验结果基本吻合,验证了该算法得到的颗粒试样在模拟道砟裂纹起裂、扩展等过程的有效性。     

Packing properties of binary mixtures in disordered sphere systems

Mixtures of binary spheres are numerically simulated using a relaxation algorithm to investigate the effects of volume fraction and size ratio. A complete profile of the packing properties of binary spheres is given. The density curve with respect to the volume fraction has a triangular shape with a peak at 70% large spheres. The density of the mixture increases with the size ratio, but the growth becomes slow in the case of a large size disparity. The volume fraction and size ratio effects are reflected in the height and movement, respectively, of specific peaks in the radial distribution functions. The structure of the mixture is further analyzed in terms of contact types, and the mean coordination number is demonstrated to be primarily affected by “large–small” contacts. A novel method for estimating the average relative excluded volume for binary spheres by weighting the percentages of contact types is proposed and extended to polydisperse packings of certain size distributions. The method can be applied to explain the density trends of polydisperse mixtures in disordered sphere systems.     

Packing properties of binary mixtures in disordered sphere systems

Mixtures of binary spheres are numerically simulated using a relaxation algorithm to investigate the effects of volume fraction and size ratio, A complete profile of the packing properties of binary spheres is given. The density curve with respect to the volume fraction has a triangular shape with a peak at 70% large spheres. The density of the mixture increases with the size ratio, but the growth becomes slow in the case of a large size disparity, The volume fraction and size ratio effects are reflected in the height and movement, respectively, of specific peaks in the radial distribution functions. The structure of the mixture is further analyzed in terms of contact types, and the mean coordination number is demonstrated to be primarily affected by ＂large-small＂ contacts. A novel method for estimating the average relative excluded volume for binary spheres by weighting the percentages of contact types is proposed and extended to polydisperse packings of certain size distributions. The method can be applied to explain the density trends of polydisperse mixtures in disordered sphere systems,     

Effect of particle size distribution on micro- and macromechanical response of granular packings under compression

The role of particle size heterogeneity on micro- and macromechanical properties of assemblies of spherical particles was studied using DEM simulations. The response to an imposed load of a granular material composed of non-uniformly sized spheres subjected to uniaxial confined compression was investigated. A range of geometrical and micro-mechanical properties of granular packings (e.g., void fraction, contact force distribution, average coordination number and degree of mobilisation of friction at contacts between particles) were examined, and provided a more accurate interpretation of the macroscopic behaviour of mixtures than has previously been available. The macromechanical study included stress transmission, stiffness and angle of internal friction of the granular assemblies.The degree of polydispersity showed slight effect on both, the void fraction and the elastic properties of the system. The tendency for increase in the lateral-to-vertical pressure ratios was observed with an increasing degree of particle size heterogeneity; however, the different pressure ratios calculated for samples with various degrees of polydispersity lay within the range of data scatter.     

Discrete element modelling of pebble beds: With application to uniaxial compression tests of ceramic breeder pebble beds

In this paper, a discrete element simulation scheme for pebble beds in fusion blankets is presented. Each individual pebble is considered as one element obeying equilibrium conditions under contact forces. We study not only the rearrangement of particles but also the overall behaviour of an assembly under the action of macroscopic compressive stresses. Using random close packing as initial configurations, the discrete element simulation of the uniaxial compression test has been quantitatively compared to experiments. This method yields the distribution of the inter-particle contact forces. Moreover, the micro-macro relations have been investigated to relate the microscopic information, such as the maximum contact force and the coordination number inside the assembly, to the macroscopic stress variables.     

Linear viscoelastic behavior of copper phthalocyanine dispersions used in printing inks

In this paper we describe the linear viscoelastic properties of copper phthalocyanine (CuPCN) dispersions that are used in the manufacturing of offset lithographic printing inks. Transmission electron microscopy shows that the primary pigment particles are rod-like and have sizes in the range of 10 to 300 nm. Steady shear measurements show that the dispersions are Newtonian at a pigment volume fraction of 0.073 and become increasingly shear thinning as the pigment volume fraction is increased. The strong shear-thinning nature of these dispersions can be attributed to the highly flocculated nature of the dispersions, which is due to interparticle attractions. The structural complexity of the dispersions also results in an unexpected linear viscoelastic response. While at low frequencies (0.1 and 1.0 Hz) the ex tent of the linear region decreases with increasing pigment concentration, at a higher frequency (10 Hz) the extent of the linear region increases with increasing pigment concentration. This increase in the linear region with increasing pigment concentration suggests that at higher frequencies the dispersion is less brittle, and that the rheological behavior is dominated by intra-aggregate associations. In addition, frequency sweeps show that the dispersions behave like a viscoelastic liquid at low pigment concentrations. However, at higher pigment concentrations (yet significantly lower than the maximum packing fraction) the dispersions behave like a cross-linking polymer at its gel point.     

Normal stress effects in the gravity driven flow of granular materials

In this paper, we study the fully developed gravity-driven flow of granular materials between two inclined plates. We assume that the granular materials can be represented by a modified form of the second grade fluid where the viscosity depends on the shear rate and volume fraction and the normal stress coefficients depend on the volume fraction. We also propose a new isotropic (spherical) part of the stress tensor which can be related to the compactness of the (rigid) particles. This new term ensures that the rigid solid particles cannot be compacted beyond a point, namely when the volume fraction has reached the critical/maximum packing value. The numerical results indicate that the newly proposed stress tensor has obvious and physically meaningful effects on both the velocity and the volume fraction fields.     

Statistics of contact force network in dense granular matter

We investigate the distribution of contact forces in a static granular system and in annular shear granular flow, using the discrete element method, and considering the influences of both packing fraction and friction coefficient. We find the existence of a critical packing fraction. If the packing fraction is lower than this critical value, all contact forces in granular system vanish. For shear granular flow, the critical packing fraction is significantly smaller than that for static granular system. The ...     

Lattice gas simulations of two-dimensional liquid foams

Liquid foam is a dense random packing of gas bubbles in a small amount of immiscible liquid containing surfactants. The liquid within the Plateau borders, although small in volume, causes considerable difficulties to investigations of the physical properties of foams, and the situation becomes even more complicated if the flow of the liquid through the foam is considered too. Here we propose a fresh approach to tackling these issues by introducing a discrete two-dimensional hybrid lattice gas model of liquid foams. While lattice gas models have been used to model two-phase liquids in the past, their application to the study of liquid foams is novel and proves promising. We represent bubble surfaces by a finite number of nodes, and model the surrounding liquid as a lattice gas (with a finite number of liquid particles). The gas in the bubbles is treated as an ideal gas at constant temperature. The model is tested by choosing an arbitrarily shaped bubble that evolves into a circular shape in agreement with Laplaces law. The model is then employed to simulate periodic ordered and disordered dry and wet foams. Since our model is specifically designed to handle wet foams up to a critical liquid fraction of 0.16 (void fraction of random packing of disks), we are able to compute the variation in coordination number (average number of neighbours of a bubble) over the whole range of liquid fractions, and we find it to be a linear function of the shear modulus.This paper was presented at the first Annual European Rheology Conference (AERC) held in Guimarães, Portugal, 11–13 September 2003.     

Physical characterisation and yield stress of a concentrated Miscanthus suspension

We present yield stress measurements of a concentrated suspension of Miscanthus x giganteus, a biomass substrate for the production of biofuels and chemicals. The particle porosity and the apparent solid packing fraction are quantified, and a simple relationship between volume and mass fractions proposed. The relationship between yield stress and solid concentration was measured with the maximum torque method of Dzuy and Boger (J Rheol 27: 321–349, 1983) for concentrations between 11 and 15 %, just below maximum packing. In the range of concentrations considered, the measured dependence of the yield stress on concentration appears to be significantly stronger than that reported for corn stover, a plant for which reliable data is available: the apparent power law exponent is found to be larger than 10, while the corresponding exponent reported for corn is in the range 4–6. We argue that differences in particle rigidity between corn and Miscanthus may play a role in explaining this result.     

转动驱动圆角立方体颗粒有序堆积的离散元模拟

宏观颗粒体系的有序堆积研究可以为热系统中微观粒子的自组装提供研究模型,也有助于工业生产中提高颗粒材料的填充率.实验发现圆筒中的圆角立方体颗粒受容器往复旋转剪切作用会实现有序堆积.为探究旋转圆筒中立方体颗粒有序堆积过程的内部结构演化过程和动力学机理,采用超二次曲面方程构造了圆角立方体颗粒,基于离散元方法对旋转剪切作用下圆...     

Fabric and connectivity as field descriptors for deformations in granular media

Granular materials involve microphysics across the various scales giving rise to distinct behaviours of geomaterials, such as steady states, plastic limit states, non-associativity of plastic and yield flow, as well as instability of homogeneous deformations through strain localization. Incorporating such micro-scale characteristics is one of the biggest challenges in the constitutive modelling of granular materials, especially when micro-variables may be interdependent. With this motivation, we use two micro-variables such as coordination number and fabric anisotropy computed from tessellation of the granular material to describe its state at the macroscopic level. In order to capture functional dependencies between micro-variables, the correlation between coordination number and fabric anisotropy limits is herein formulated at the particle level rather than on an average sense. This is the essence of the proposed work which investigates the evolutions of coordination number distribution (connectivity) and anisotropy (contact normal) distribution curves with deformation history and their inter-dependencies through discrete element modelling in two dimensions. These results enter as probability distribution functions into homogenization expressions during upscaling to a continuum constitutive model using tessellation as an abstract representation of the granular system. The end product is a micro-mechanically inspired continuum model with both coordination number and fabric anisotropy as underlying micro-variables incorporated into a plasticity flow rule. The derived plastic potential bears striking resemblance to cam–clay or stress–dilatancy-type yield surfaces used in soil mechanics.     

Shear viscosity of slightly-elastic concentrated suspensions at low and high shear rates

A quasi-static asymptotic analysis is employed to investigate the elastic effects of fluids on the shear viscosity of highly concentrated suspensions at low and high shear rates. First a brief discussion is presented on the difference between a quasi-static analysis and the periodic-dynamic approach. The critical point is based on the different order-of-contact time between particles. By considering the motions between a particle withN near contact point particles in a two-dimensional “cell” structure and incorporating the concept of shear-dependent maximum packing fraction reveals the structural evolution of the suspension under shear and a newly asymptotic framework is devised. In order to separate the influence of different elastic mechanisms, the second-order Rivlin-Ericksen fluid assumption for describing normal-stress coefficients at low shear rates and Harnoy's constitutive equation for accounting for the stress relaxation mechanism at high shear rates are employed. The derived formulation shows that the relative shear viscosity is characterized by a recoverable shear strain,S _R at low shear rates if the second normal-stress difference can be neglected, and Deborah number,De, at high shear rates. The predicted values of the viscosities increase withS _R , but decrease withDe. The role ofS _R in the matrix is more pronounced than that ofDe. These tendencies are significant when the maximum packing fraction is considered to be shear-dependent. The results are consistent with that of Frankel and Acrivos in the case of a Newtonian suspension, except for when the different divergent threshhold is given as [1 ? (Φ/Φ _m )^1/2] ? 1.     

Determining a representative element volume for DEM simulations of samples with non-circular particles

Numerical studies on the number of particles or system size required to attain a representative element volume (REV) for discrete element method (DEM) simulations of granular materials have almost always considered samples with spherical or circular particles. This study considers how many particles are needed to attain a REV for 2D samples of 2-disc cluster particles where the particle aspect ratio (AR) was systematically varied. Dense and loose assemblies of particles were simulated. The minimum REV was assessed both by considering the repeatability of static packing characteristics and the shearing behaviour in biaxial compression tests, and by investigating the effect of sample size on the measured characteristics and observed shearing behaviour. The repeatability of the data considered generally improved with increasing sample size. The packing characteristics of the dense samples were more repeatable suggesting that the minimum REV reduces with increasing packing density. The minimum REV was observed to be sensitive to the characteristic measured. Although the overall responses of the samples during shear deformation were similar irrespective of the sample sizes, the smaller the sample size, the higher the fluctuations observed in the responses. Analysis of the coefficient of variation of the fluctuations around the critical state stress ratio can provide insight as to whether a REV is attained. The particle AR influences the effect of sample size on shearing characteristics and thus the minimum number of particles required to attain a REV; this can be explained by the influence of AR on the number of contacts within the samples.     

Three-dimensional numerical simulation of drops suspended in simple shear flow at finite Reynolds numbers

A three-dimensional study of suspension of drops in simple shear flow has been performed at finite Reynolds numbers. Results are obtained using a finite difference/front tracking method in a periodic domain. The effects of the Reynolds number and the Capillary number are addressed at two volume fractions: 0.195 and 0.34. It is observed that suspensions of deformable drops exhibit a shear-thinning behavior. Similar to the motion of a single drop, drops migrate away from the walls. The effective viscosity, the first and the second normal stress differences oscillate around a mean value in all cases. The first normal stress difference increases with the Capillary number, the Reynolds number and the volume fraction. Results show that drops deform more and orient more in the flow direction as the Capillary number or the volume fraction is increased. Also, the average size of clusters is smaller than for suspension of rigid particles. The radial dependence of the pair distribution function across the channel has been studied. This dependency shows that the tendency to form clusters is reduced as the Capillary number increases or the volume fraction decreases.     

Random packing of tetrahedral particles using the polyhedral discrete element method

The random packing of tetrahedral particles is studied by applying the discrete element method (DEM), which simulates the effects of friction, height ratio, and eccentricity. The model predictions are analyzed in terms of packing density and coordination number (CN). It is demonstrated that friction has the maximal effect on packing density and mean CN among the three parameters. The packing density of the regular tetrahedron is 0.71 when extrapolated to a zero friction effect. The shape effects of height ratio and eccentricity show that the regular tetrahedron has the highest packing density in the family of tetrahedra, which is consistent with what has been reported in the literature. Compared with geometry-based packing algorithms, the DEM packing density is much lower. This demonstrates that the inter-particle mechanical forces have a considerable effect on packing. The DEM results agree with the published experimental results, indicating that the polyhedral DEM model is suitable for simulating the random packing of tetrahedral particles.     

Stress redistribution in notched specimens during fatigue cycling

Most materials exhibit a change in stress-strain relationship when subjected to fatigue stresses. In this work, the effect of this change on the stress distribution across the throat of notched-plate specimens of mild steel is examined. Using a set of strain gages, the strain distribution across the specimens was determined under dynamic conditions for various numbers of cycles. Tests of unnotched specimens were used to obtain the cyclic stress-strain properties for corresponding numbers of cycles, and from these data the stress distribution in the notched specimens was determined. Tests in which the strain amplitude at the notch root was held constant revealed a decreasing maximum stress with fatigue cycles. In another series of tests, in which the load amplitude was constant, the maximum stress amplitude was observed to decrease with number of fatigue cycles, despite an increasing strain amplitude. In both types of tests, the stress-concentration factor was observed to decrease with increasing number of fatigue cycles.     

浓悬浮体的屈服应力和最大填充率

对剪切稀化的浓悬浮体，假设在剪切状态下结构参数变化的速率服从一级动力学关系，在此基础上导出浓悬浮体最大填充率与剪应力、屈服应力与填充率之间关系的数学模型．用硅粉－甘油水溶液及可可粉－可可脂两种悬浮体的实验确认了模型对工业微米级颗粒悬浮体的适用性     

Gelation of a radiation crosslinked model polyethylene

A radiation crosslinked model linear low-density polyethylene (LLDPE) exhibits power-law relaxation,G(t) =St ^–n at its gel point (GP). The relaxation exponent has a value of about 0.46. The relaxation behavior is dominated by power laws, not only directly at GP, but in a very broad vicinity of GP and in a frequency window, which narrows with distance from the gel point. The power law exponent decreases with increasing radiation dose (increasing extent of crosslinking). Independent measurements of the gel fraction and the molecular-weight distribution of the radiated samples' soluble fraction support the rheological observations.Delivered as a Keynote Lecture at the Golden Jubilee Conference of the British Society of Rheology and Third European Rheology Conference, Edinburgh, 3–7 September, 1990.