Analysis of stresses in two-dimensional isostatic granular systems

We analyze a recently proposed continuous model for stress fields that develop in two-dimensional purely isostatic granular systems. We present a reformulation of the field equations, as a linear first-order hyperbolic system, and show that it is very convenient both for analysis and for numerical computations. Our analysis allows us to predict quantitatively the formation and directions of stress paths and, from these, trajectories and magnitudes of force chains, given the structure in terms of a particular fabric tensor. We further predict quantitatively changes of stresses along the paths, as well as leakage and branching of stress from the main paths into the cones that they make in terms of the fabric tensor. Numerical computations in both Cartesian and cylindrical coordinates verify the analytic results and illustrate the rich behavior discovered. All the phenomena predicted by our solutions have been observed experimentally, suggesting that stresses in isostatic systems can form a base model for a more developed stress theory in granular materials.