No Arabic abstract
We have made experimental observations of the force networks within a two-dimensional granular silo similar to the classical system of Janssen. Models like that of Janssen predict that pressure within a silo saturates with depth as the result of vertical forces being redirected to the walls of the silo where they can then be carried by friction. By averaging ensembles of experimentally-obtained force networks in different ways, we compare the observed behavior with various predictions for granular silos. We identify several differences between the mean behavior in our system and that predicted by Janssen-like models: We find that the redirection parameter describing how the force network transfers vertical forces to the walls varies with depth. We find that changes in the preparation of the material can cause the pressure within the silo to either saturate or to continue building with depth. Most strikingly, we observe a non-linear response to overloads applied to the top of the material in the silo. For larger overloads we observe the previously reported giant overshoot effect where overload pressure decays only after an initial increase [G. Ovarlez et al., Phys. Rev. E 67, 060302(R) (2003)]. For smaller overloads we find that additional pressure propagates to great depth. This effect depends on the particle stiffness, as given for instance by the Youngs modulus, E, of the material from which the particles are made. Important measures include E, the unscreened hydrostatic pressure, and the applied load. These experiments suggest that when the load and the particle weight are comparable, particle elasticity acts to stabilize the force network, allowing non-linear network effects to be seen in the mean behavior.
A remarkable feature of static granular matter is the distribution of force along intricate networks. Even regular inter-particle contact networks produce wildly inhomogeneous force networks where certain chains of particles carry forces far larger than the mean. In this paper, we briefly review past theoretical approaches to understanding the geometry of force networks. We then investigate the structure of experimentally-obtained granular force networks using a simple algorithm to obtain corresponding graphs. We compare our observations with the results of geometric models, including random bond percolation, which show similar spatial distributions without enforcing vector force balance. Our findings suggest that some aspects of the mean geometry of granular force networks may be captured by these simple descriptions.
Disordered biopolymer gels have striking mechanical properties including strong nonlinearities. In the case of athermal gels (such as collagen-I) the nonlinearity has long been associated with a crossover from a bending dominated to a stretching dominated regime of elasticity. The physics of this crossover is related to the existence of a central-force isostatic point and to the fact that for most gels the bending modulus is small. This crossover induces scaling behavior for the elastic moduli. In particular, for linear elasticity such a scaling law has been demonstrated [Broedersz et al. Nature Physics, 2011 7, 983]. In this work we generalize the scaling to the nonlinear regime with a two-parameter scaling law involving three critical exponents. We test the scaling law numerically for two disordered lattice models, and find a good scaling collapse for the shear modulus in both the linear and nonlinear regimes. We compute all the critical exponents for the two lattice models and discuss the applicability of our results to real systems.
Starting from a microscopic model of randomly cross-linked particles with quenched disorder, we calculate the Laudau-Wilson free energy S for arbitrary cross-link densities. Considering pure shear deformations, S takes the form of the elastic energy of an isotropic amorphous solid state, from which the shear modulus can be identified. It is found to be an universal quantity, not depending on any microscopic length-scales of the model.
We carry out a direct comparison of experimental and numerical realizations of the exact same granular system as it undergoes shear jamming. We adjust the numerical methods used to optimally represent the experimental settings and outcomes up to microscopic contact force dynamics. Measures presented here range form microscopic, through mesoscopic to system-wide characteristics of the system. Topological properties of the mesoscopic force networks provide a key link between micro and macro scales. We report two main findings: the number of particles in the packing that have at least two contacts is a good predictor for the mechanical state of the system, regardless of strain history and packing density. All measures explored in both experiments and numerics, including stress tensor derived measures and contact numbers depend in a universal manner on the fraction of non-rattler particles, $f_{NR}$. The force network topology also tends to show this universality, yet the shape of the master curve depends much more on the details of the numerical simulations. In particular we show that adding force noise to the numerical data set can significantly alter the topological features in the data. We conclude that both $f_{NR}$ and topological metrics are useful measures to consider when quantifying the state of a granular system.
The response to a localized force provides a sensitive test for different models of stress transmission in granular solids. The elasto-plastic models traditionally used by engineers have been challenged by theoretical and experimental results which suggest a wave-like (hyperbolic) propagation of the stress, as opposed to the elliptic equations of static elasticity. Numerical simulations of two-dimensional granular systems subject to a localized external force are employed to examine the nature of stress transmission in these systems as a function of the magnitude of the applied force, the frictional parameters and the disorder (polydispersity). The results indicate that in large systems (typically considered by engineers), the response is close to that predicted by isotropic elasticity whereas the response of small systems (or when sufficiently large forces are applied) is strongly anisotropic. In the latter case the applied force induces changes in the contact network accompanied by frictional sliding. The larger the coefficient of static friction, the more extended is the range of forces for which the response is elastic and the smaller the anisotropy. Increasing the degree of polydispersity (for the range studied, up to 25%) decreases the range of elastic response. This article is an extension of a previously published letter [1].