An experimental investigation of the force network ensemble

TL;DR

This study experimentally investigates force networks in a cyclically compressed quasi-2D granular system, visualizing force chains and analyzing their statistical properties to compare with theoretical predictions from the force network ensemble theory.

Contribution

It provides the first detailed experimental visualization and analysis of force network evolution under cyclic loading in a frictionless granular system, validating theoretical models.

Findings

Force chains form and evolve with cyclic compression.

The contact force distribution matches theoretical predictions.

Force network properties are statistically characterized.

Abstract

We present an experiment in which a horizontal quasi-2D granular system with a fixed neighbor network is cyclically compressed and decompressed over 1000 cycles. We remove basal friction by floating the particles on a thin air cushion, so that particles only interact in-plane. As expected for a granular system, the applied load is not distributed uniformly, but is instead concentrated in force chains which form a network throughout the system. To visualize the structure of these networks, we use particles made from photoelastic material. The experimental setup and a new data-processing pipeline allow us to map out the evolution subject to the cyclic compressions. We characterize several statistical properties of the packing, including the probability density function of the contact force, and compare them with theoretical and numerical predictions from the force network ensemble theory.

Equations

$$F=\Big{\langle}\sqrt{F_{N}^{2}+F_{T}^{2}}\leavevmode\nobreak\ \leavevmode\nobreak\ \Big{\rangle}$$

Full text

1 Introduction

Although the positions of particles in a jammed granular system are fixed to a specific geometrical configuration, the particle positions alone are not sufficient to determine the force network that carries the load on that packing. As such, an underdetermined mechanical system, there are many ways in wich force and torque balance on each particle can be statistfied for any given packing geometry and boundary conditions Snoeijer2004 ; Tighe2010 . There are both stable configurations, and unstable configurations and initially stable granular systems can evolve into catastrophic failure. While two packings might have the same occupied volume or internal pressure they might have vastly different bulk material properties Dagois-Bohy2012 . For frictional granular system it remains an open question to determine whether the statistics due to the stress state and the volume of the system can be decoupled from each other.

To make predictions for the physical behavior of granular systems, tools and concepts from statistical physics have become widely used Bi2015 , but what is the correct ensemble to describe jammed granular packings? To get more insight into these jammed packings one needs to look at the structure of force networks that form in loaded granular packings. When the packing is subjected to external load, not all particles share the load equally but the forces are highly localized into force chains. In this work, we present an experiment to look at the distribution of forces in a loaded granular system, while disentangling the effects of configuration from other influences. This is of interest to compare these distributions to predictions from the Force Network Ensemble theory Snoeijer2004 .

The Force Network Ensemble

The Force Network Ensemble (FNE) is a concept introduced by Snoeijer et al. Snoeijer2004 in which they use an ensemble approach for examining the force distribution in static granular packings. Since forces in fixed granular packings are typically underdetermined, the ensemble averages over all microscopic variations of a packing, an approach that goes back to Edwards Edwards1989 . The FNE predicts, among other things Tighe2010 , a fitnite value for $P(F)$ as $F\rightarrow 0$ and a faster than exponential decay of the probability density function (PDF) of the contact forces at large forces, depending on the dimensionality of the packing Snoeijer2004 . For a two dimensional system, $P(F)$ is predicted to have a gaussian tail. For a review see Tighe2010 (and references therein) where it is also discussed that the peak in the PDF should vanish for anisotropic stresses. Saithoh et al. Saitoh2015 , using molecular dynamics simulations to determine transition rates for contact changes, recently also found a master equation that describes the PDF of forces in soft particle packings.

2 An Experimental Approach

While a number of experimental works exist, e.g. Corwin2005 ; Howell1999 , PDFs of contact forces are most often produced from numerical simulations, including a recent pair of papers by Pugnaloni Pugnaloni2016 and Kondic Kondic2016 where they study the structure of force networks in tapped particulate systems of disks superimposed by a gravity force. However, in most experimental studies the external load that probes the force network is not the only force applied to the system, there is additionally a load superimposed by gravity, or basal friction Kovalcinova2016 . Further, there is only few experiments Puckett2013 probing granular ensembles.

In this manuscript, we present an experiment that is designed to enumerate how many force configurations of a single hyperstatic granular arrangement are practically accessible, while at the same time keeping the external load the only force that is beeing applied to the system. To achieve this we prepare a horizontal quasi 2D granular system that is floated on a gentle air cushion, thereby generating an effectively gravity free system without basal friction Puckett-thesis . The particles are confined by a piston that can apply an uniaxial load to the packing. A schematic drawing of the experimental setup is detailed in Fig. 1. By cyclically loading and unloading the packing in a way that will not change the particle configuration (no neighbors changes), the system cycles trough many contact force configurations due to microscopic changes of the exact contact point. For the experiments done here, we compress the packing in steps of constant volume ($\Delta V=0.002869V_{initial}$) and the compression steps are applied quasistatically over 20 substeps of $\Delta x=0.01\leavevmode\nobreak\ \text{mm}$. The initial volume $V_{initial}$ was chosen to be close to the onset of jamming, and the final volume so that the mean contact force rises by more than a factor of 3. We performe the experiment with several random configuration of 29 particles of two different radii ($r_{1}=5.5\text{mm}$ and $r_{2}=7.6\text{mm}$) to prevent crystallization.

In order to extract force information from the experiment, the particles are made of photoelastic material (Vishay PSM-4), which will shift the polarization of light that is shined through it as a function of the applied load. A model of the force modulated light intensity can then be fitted to camera images of the particles Daniels2016 ; Puckett2013 .

3 Results

When we start the load cycling we observe that even after an initial annealing period there are variations in the force network, (see Fig. 3), while the particle configuration is unchanged. This validates our experimental approach and allows us to probe the nature of the FNE.

We run the experiment for $\approx 1000$ cycles and observe strong fluctuations in the contact number, determined by a minimum threshold force ($F_{th}>0.01$ N) and in the number of load bearing particles (the number of particles with one ore more contacts above the threshold force). Figure 2 shows the average contact force

[TABLE]

and its standard deviation over all cycles, as a function of the applied compression. We see that the standard deviations in $F$ grow with the applied compression. For higher packing fractions, the intervals given by the STD($F$) around $F$ begin to overlap for consecutive compression steps.

Figure 4 shows the PDFs of the contact forces for 5 compression steps over all cycles. The PDF exhibits a strong peak, decaying in the limit of both low and high forces. These results are qualitatively similar to the predictions by FNE theory: most notably a stronger than exponential decay in $P(F)$. We also see the peak in $P(F)$, along with the average contact force, move towards higher values as the system gets compressed stronger. In fact, the PDFs collapse onto a single curve when normalized by the average contact force at the corresponding compression step. Although we perform an experiment with anisotropic loading, we nonetheless identify a peak in the distribution, a feature that is suggested to vanish for anisotropic loads Tighe2010 . For small forces we find, $P(F)$ to rise exponentially, as approximately $F^{3/2}$ . Wyart Wyart2012 , showed that the exponent is determined by the pair distribution function g(r). However, in our small system with a fixed particle configuration, g(r) is undersampled, so the question arises what the exponent is set by here.

4 Conculsions and Outlook

We have designed an experiment that can explore the Force Network Ensemble of a two dimensional granular packing, while excluding forces other than the applied load. We find that we can qualitatively reproduce some of the features in contact force distribution predicted by the FNE theory.

Due to the limited number of cycles and initial packing configurations explored here, it is not immediately clear whether some features we observed in the contact force PDFs are due to sample size or the specific initial configuration. Repeating the experiment not only for more compression cycles but also many different packing configurations, would allow us to probe ergodic effects by contrasting the time average statistics to the ensemble average statistics.

Furthermore, future work should try to identify whether the are several subpopulations of networks which differ in their contact force PDFs. These subpopulations can be found by using tools from network theory such as community detection Bassett2015 . More generally, the data generated by this experiment will be useful in relating force network features to macroscopic packing properties Giusti2016 ; Papadopoulos2016 .

Acknowledgements and References

We gratefully acknowledge James Puckett for the design and construction of the air table on which the apparatus is based, and for the inspiration for the new parallelized version of the contact-force code. This research was suppored by the James S. McDonnell Foundation and the NSF through grants DMR-0644743 and DMR-1206808.