    Next: Vibration induced compaction and Up: Metastable granular materials Previous: Metastable granular materials   Contents

## Janssen effect and the distribution of internal stresses

In 1895 H. A. Janssen  discovered that in a vertical cylinder the pressure measured at the bottom does not depend upon the height of the filling, i.e. it does not follow the Stevin law which is valid for Newtonian fluids at rest : (1.1)

where is the vertical pressure, the density of the fluid, the gravity acceleration and the height of the column of fluid above the level of measurement. The pressure in the granular material follows instead a different law, which accounts for saturation: (1.2)

where is of the order of the radius of the cylinder. This guarantees the flow rate in a hourglass to be constant. Moreover, this law is very important in the framework of silos building, as the difference between ordinary hydrostatic and granular hydrostatic is mainly due to the presence of anomalous side pressure, i.e. force exerted against the walls of the cylinder. It happens that the use of a fluid-like estimate of the horizontal and vertical pressure leads to an underestimating of the side pressure and, consequently, to unexpected (and dramatic) explosions of silos .

The first interpretation of the law has been given by Janssen in his paper, in terms of a simplified model with the following assumptions:

1. The vertical pressure is constant in the horizontal plane.

2. The horizontal pressure is proportional to the vertical pressure where is constant in space.

3. The wall friction (where is the static friction coefficient) sustains the vertical load at contact with the wall.

4. The density of the material is constant over all depths.

In particular the first assumption is not true (the pressure depends also upon the distance from the central axis of the cylinder) but is not essential in this model (as it is formulated as a one-dimensional problem), while the second assumption should be obtained by means of constitutive relations, i.e. it requires a microscopic justification.

Imposing the mechanical equilibrium of a disk of granular material of height and radius (the radius of the container) the following equation is obtained: (1.3)

which becomes: (1.4)

where . This equation is exactly solved by the function (1.2). The particular behavior of the vertical pressure in granular materials is mainly due to the anomalies in the stress propagation. The configuration of the grains in the container is random and the weight can be sustained in many different ways: every grain discharges its load to other grains underlying it or at its sides, creating big arches and therefore propagating the stress in unexpected directions. Moreover, arching is not only a source of randomness, but also of strong fluctuations, i.e. disorder: in a granular assembly some force chains can be very long and span the size of the entire system, posing doubts on the validity of (local) mean field modeling.

Further interesting phenomena have been experimentally observed in the statics of granular materials:

• the fluctuations of the pressure at the bottom of a silo are large, they can change of more than in repeated pouring of grains in the same container , and in a single pouring the distribution of stresses, measured deep inside or at the bottom of the silo, show an exponential tail [141,168,46].

• the vertical pressure below conical sandpiles does not follow the height of the material, but rather has a minimum underneath the apex of the pile 

Different models have been proposed and debated in the last years, in order to understand the problem of the distribution of forces in a silo or in a granular heap.

The q-model has been introduced in 1995 (remarkably a century after the work of Janssen) by Liu et al. [141,66] in order to reproduce the stress probability distribution observed in experiments. The model consists of a regular lattice of sites each with a particle of mass unity. Each site in layer is connected to exactly sites in layer . Only the vertical components of the forces are considered explicitly: a random fraction of the total weight supported by particle in layer is transmitted to particle in layer . Thus the weight supported by the particle in layer D at the -th site, , satisfies a stochastic equation: (1.5)

The random variables are taken independent except for the constraint (1.6)

which enforces the condition of force balance on each particle. Given a distribution of 's, it can be obtained the probability distribution of finding a site that bears a weight on layer . By means of mean field calculations, exact calculations and numerical solutions, the authors conclude that (apart of some limiting cases) a generic continuous distributions of 's lead to a distribution of weights that, normalized to the mean, is independent of depth at large and which decays exponentially at large weights. They find also a good agreement with molecular dynamics simulations of the packing of hard spheres. The q-model has many limits:

• it is a scalar model, i.e. it takes into account only one component of the internal stress (this was solved by P. G. de Gennes  who introduced a vectorial version of the q-model, obtaining a more realistic propagation of forces)

• it does not reproduce the minimum (the dip'') of the pressure measured under central axis at the bottom of a sand heap

• it does not reproduce the Janssen law ; this problem appears as a consequence of the diffusive nature of the q-model solution: the saturation depth where the stress distribution becomes independent of depth scales with the silo width as , at odds with the Janssen observation that predicts .

A more refined version of the Janssen model has been introduced by Bouchaud et al. : the authors have considered a local version of the Janssen assumption on the proportionality between horizontal and vertical stresses: which lead to the linear equation: (1.8)

This equation for the vertical stress is hyperbolic and therefore differs from the equivalent equation for an elastic medium, which is elliptic , and from the q-model equation that is parabolic (as a diffusion equation): it is equivalent to the equation for the wave propagation with as the time'' variable and as the inverse of the propagation velocity. This model well reproduces the dip in the measure of the pressure under the bottom of the conical heap . A cellular automaton was introduced by Hemmingsson  which was capable of reproduce the dip under the heap as well as the correct Janssen law (with the linear scaling). In the framework of the study of force networks in the bulk of a static arrangement of grains a key role was played by the experiments on the propagation of sound. The inhomogeneities present inside a granular medium can drastically change the propagation of mechanical perturbations. Liu and Nagel [139,140,138] have addressed this issue in several experiments. They have discovered  that the in the bulk of a granular medium perturbed by a harmonic force ( Hz) the fluctuations could be very large, measuring power-law spectra of the kind with . Then they have seen  that the sound group velocity can reach times the phase velocity and that a change in the amplitude of vibration can result in a hysteretic behavior (due to a rearrangement of force chains). They have also measured  a variation of the sound transmission as a consequence of a very small (compared to the size of the grains) thermal expansion of a little carbon resistor substituted to a grain of the granular medium. This sensitivity to perturbation is another signature of the strong disorder (arching and chain forces) in the bulk of the medium.    Next: Vibration induced compaction and Up: Metastable granular materials Previous: Metastable granular materials   Contents
Andrea Puglisi 2001-11-14