二维颗粒体系单轴压缩形成的力链结构

从接触力、能量分布和接触网络结构特点出发,提出了强力链的力大小(Fc)判据(Fc大于平均接触力〈F〉)和角度θc判据(θc=180/〈Z〉,其中〈Z〉是平均配位数),指出强力链和弱力链是本质不同的两类结构存在于颗粒体系中,其中强力链网络与体系的宏观性质直接相关.以二维颗粒体系的单轴压缩为例,计算发现了强力链长度的幂率分布规律,分析了侧向压力系数与相应强力链网络结构的关系:当内部强力链网络充分发育而不再变化时,侧向压力系数趋于稳定数值.     

Stresses in isostatic granular systems and emergence of force chains

Progress is reported on several questions that bedevil understanding of granular systems: (i) Are the stress equations elliptic, parabolic, or hyperbolic? (ii) How can the often-observed force chains be predicted from a first-principles continuous theory? (iii) How do we relate insight from isostatic systems to general packings? Explicit equations are derived for the stress components in two dimensions including the dependence on the local structure. The equations are shown to be hyperbolic and their general solutions, as well as the Green function, are found. It is shown that the solutions give rise to force chains, and the explicit dependence of the force chains trajectories and magnitudes on the local geometry is predicted. Direct experimental tests of the predictions are proposed. Finally, a framework is proposed to relate the analysis to nonisostatic and more realistic granular assemblies.     

Elasticity from the force network ensemble in granular media

Transmission of forces in static granular materials is studied within the framework of the force network ensemble, by numerically evaluating the mechanical response of hexagonal packings of frictionless grains and rectangular packings of frictional grains. In both cases, close to the point of application of the overload, the response is nonlinear and displays two peaks, while at larger length scales, it is linear and elasticlike. The crossover between these two behaviors occurs at a depth that increases with the magnitude of the overload and decreases with increasing friction.     

重力作用下颗粒介质应力链的离散元模拟

用离散元的方法模拟了仅有重力作用的二维颗粒系统内部力的分布情况,并根据力的大小得到颗粒之间的应力链.模拟结果与颗粒介质研究中的两个著名模型q模型和α模型作了对比,并与光弹实验的结果作了比较.对比结果表明,模拟结果与实验相似,而与两个概率模型有一定的差异.另外计算结果还表明,颗粒介质中力大小的概率分布极为不均匀,较大的力概率呈指数衰减,应力链的分布具有分形特征. 关键词： 颗粒介质 离散元 应力链 光弹实验     

颗粒链在振动条件下的行为研究

采用实验方法研究打结颗粒链在外部激振下解结过程与激振源以及颗粒链长度的关系.实验结果表明颗粒链结解开所需的时间与振动频率有关,颗粒链解结存在截止频率和最佳振动频率;颗粒链解结时间与链长之间呈非线性关系.同时采用Monte Carlo时步模拟的方法,对珠链解结的动力学模型进行仿真.     

Fluctuations in granular media

Dense slowly evolving or static granular materials exhibit strong force fluctuations even though the spatial disorder of the grains is relatively weak. Typically, forces are carried preferentially along a network of "force chains." These consist of linearly aligned grains with larger-than-average force. A growing body of work has explored the nature of these fluctuations. We first briefly review recent work concerning stress fluctuations. We then focus on a series of experiments in both two- and three-dimension [(2D) and (3D)] to characterize force fluctuations in slowly sheared systems. Both sets of experiments show strong temporal fluctuations in the local stress/force; the length scales of these fluctuations extend up to 10(2) grains. In 2D, we use photoelastic disks that permit visualization of the internal force structure. From this we can make comparisons to recent models and calculations that predict the distributions of forces. Typically, these models indicate that the distributions should fall off exponentially at large force. We find in the experiments that the force distributions change systematically as we change the mean packing fraction, gamma. For gamma's typical of dense packings of nondeformable grains, we see distributions that are consistent with an exponential decrease at large forces. For both lower and higher gamma, the observed force distributions appear to differ from this prediction, with a more Gaussian distribution at larger gamma and perhaps a power law at lower gamma. For high gamma, the distributions differ from this prediction because the grains begin to deform, allowing more grains to carry the applied force, and causing the distributions to have a local maximum at nonzero force. It is less clear why the distributions differ from the models at lower gamma. An exploration in gamma has led to the discovery of an interesting continuous or "critical" transition (the strengthening/softening transition) in which the mean stress is the order parameter, and the mean packing fraction, gamma, must be adjusted to a value gamma(c) to reach the "critical point." We also follow the motion of individual disks and obtain detailed statistical information on the kinematics, including velocities and particle rotations or spin. Distributions for the azimuthal velocity, V(theta), and spin, S, of the particles are nearly rate invariant, which is consistent with conventional wisdom. Near gamma(c), the grain motion becomes intermittent causing the mean velocity of grains to slow down. Also, the length of stress chains grows as gamma-->gamma(c). The 3D experiments show statistical rate invariance for the stress in the sense that when the power spectra and spectral frequencies of the stress time series are appropriately scaled by the shear rate, Omega, all spectra collapse onto a single curve for given particle and sample sizes. The frequency dependence of the spectra can be characterized by two different power laws, P proportional, variant omega(-alpha), in the high and low frequency regimes: alpha approximately 2 at high omega; alpha<2 at low omega. The force distributions computed from the 3D stress time series are at least qualitatively consistent with exponential fall-off at large stresses. (c) 1999 American Institute of Physics.     

Numerical simulation of two-dimensional granular shearing flows and the friction force of a moving slab on the granular media

This paper studies some interesting features of two-dimensional granular shearing flow by using molecular dynamic approach for a specific granular system.The obtained results show that the probability distribution function of velocities of particles is Gaussian at the central part,but diverts from Gaussian distribution nearby the wall.The macroscopic stress along the vertical direction has large fluctuation around a constant value,the non-zero average velocity occurs mainly near the moving wall,which forms a shearing zone.In the shearing movement,the volume of the granular material behaves in a random manner.The equivalent friction coefficient between moving slab and granular material correlates with the moving speed at low velocity,and approaches constant as the velocity is large enough.     

Strain localisation in granular media

This paper discusses strain localisation in granular media by presenting experimental, full-field analysis of mechanical tests on sand, both at a continuum level, as well as at the grain scale. At the continuum level, the development of structures of localised strain can be studied. Even at this scale, the characteristic size of the phenomena observed is in the order of a few grains. In the second part of this paper, therefore, the development of shear bands within specimen of different sands is studied at the level of the individual grains, measuring grains kinematics with x-ray tomography. The link between grain angularity and grain rotation within shear bands is shown, allowing a grain-scale explanation of the difference in macroscopic residual stresses for materials with different grain shapes. Finally, rarely described precursors of localisation, emerging well before the stress peak are observed and commented.     

Modeling force transmission in granular materials

The probability density function of contact forces in granular materials has been extensively studied and modeled as an outstanding signature of granular microstructure. Arguing that particle environments play a fundamental role in force transmission, we analyze the effects of steric constraints with respect to force balance condition and show that each force may be considered as resulting from a balance between lower and larger forces in proportions that mainly depend on steric effects. This idea leads to a general model that predicts an analytical expression of force density with a single free parameter. This expression fits well our simulation data and generically predicts the exponential fall-off of strong forces, a small peak below the mean force and the non-zero probability of vanishingly small forces.     

Physics of humid granular media

A very small amount of liquid added or condensed from a vapor in a granular heap can induce dramatic changes of its static properties. In this paper we review recent advances in humid granular media. We discuss the first approaches for describing the cohesion forces acting between spherical rough beads, and their effect on the maximum avalanche angle of a granular heap. We also discuss the time dependency of these cohesive forces leading to ageing effects in the properties of the medium. To cite this article: L. Bocquet et al., C. R. Physique 3 (2002) 207–215.     

Interrupted shear of granular media

The response of a granular material during a stop-and-go shear experiment is investigated using an annular shear cell and silicagel powders of different particle sizes. The experimental results are examined on the basis of the Dieterich-Rice-Ruina model for solid friction. In addition to making this analogy with solid friction, we describe a new instability that is observed when restarting shear, where the powder bed is found to slip and compact for short hold times but only dilates for long hold times. The minimum hold time to restore a non-slip behaviour has been investigated for different size particles and normal loadings. The observed dependencies show analogies between this behaviour and the sliding rearrangements seen above the stick-slip threshold.     

Knots and random walks in vibrated granular chains

We study experimentally statistical properties of the opening times of knots in vertically vibrated granular chains. Our measurements are in good qualitative and quantitative agreement with a theoretical model involving three random walks interacting via hard-core exclusion in one spatial dimension. In particular, the knot survival probability follows a universal scaling function which is independent of the chain length, with a corresponding diffusive characteristic time scale. Both the large-exit-time and the small-exit-time tails of the distribution are suppressed exponentially, and the corresponding decay coefficients are in excellent agreement with theoretical values.     

Velocity statistics in excited granular media

We present an experimental study of velocity statistics for a partial layer of inelastic colliding beads driven by a vertically oscillating boundary. Over a wide range of parameters (accelerations 3-8 times the gravitational acceleration), the probability distribution P(v) deviates measurably from a Gaussian for the two horizontal velocity components. It can be described by P(v) approximately exp(-mid R:v/v(c)mid R:(1.5)), in agreement with a recent theory. The characteristic velocity v(c) is proportional to the peak velocity of the boundary. The granular temperature, defined as the mean square particle velocity, varies with particle density and exhibits a maximum at intermediate densities. On the other hand, for free cooling in the absence of excitation, we find an exponential velocity distribution. Finally, we examine the sharing of energy between particles of different mass. The more massive particles are found to have greater kinetic energy. (c) 1999 American Institute of Physics.     

Drag induced lift in granular media

Laboratory experiments and numerical simulation reveal that a submerged intruder dragged horizontally at a constant velocity within a granular medium experiences a lift force whose sign and magnitude depend on the intruder shape. Comparing the stress on a flat plate at varied inclination angle with the local surface stress on the intruders at regions with the same orientation demonstrates that intruder lift forces are well approximated as the sum of contributions from flat-plate elements. The plate stress is deduced from the force balance on the flowing media near the plate.     

Impact craters in loose granular media

The craters formed by the impact of steel balls with a loose sand bed are experimentally studied. Both crater size and morphology depend strongly on the impact angle. Two scaling laws, corresponding to length and width, respectively, are proposed.     

Strain stiffening in random packings of entangled granular chains

Random packings of granular chains are presented as a model system to investigate the contribution of entanglements to strain stiffening. The chain packings are sheared in uniaxial compression experiments. For short chain lengths, these packings yield when the shear stress exceeds the scale of the confining pressure, similar to granular packings of unconnected particles. In contrast, packings of chains which are long enough to form loops exhibit strain stiffening, in which the effective stiffness of the material increases with strain, similar to many polymer materials. The latter packings can sustain stresses orders-of-magnitude greater than the confining pressure, and do not yield until the chain links break. X-ray tomography measurements reveal that the strain-stiffening packings contain system-spanning clusters of entangled chains.     

Compaction and force propagation in granular packings

Motivated by experimental results in granular media, and some recently proposed scalar force models for these materials, some aspects of packing properties and their force distributions are studied in a frustrated Ising lattice gas.     

Acoustic virtual mass of granular media

Mechanical coupling between grains in a randomly packed unconsolidated granular medium is shown to cause an increase in the effective inertia, hence, a reduction in sound and shear wave speeds, relative to predictions by the standard expressions for a uniform elastic solid. The effect may be represented as a virtual mass term, and directly related to the scintillation index of the grain-to-grain contact stiffness.