Stress anisotropy in shear-jammed packings of frictionless disks

Abstract

We perform computational studies of repulsive, frictionless disks to investigate the development of stress anisotropy in mechanically stable (MS) packings. We focus on two protocols for generating MS packings: 1) isotropic compression and 2) applied simple or pure shear strain $\gamma$ at fixed packing fraction $\varphi$. MS packings of frictionless disks occur as geometric families (i.e. parabolic segments with positive curvature) in the $\varphi$-$\gamma$ plane. MS packings from protocol 1 populate parabolic segments with both signs of the slope, $d\varphi /d\gamma >0$ and $d\varphi /d\gamma <0$. In contrast, MS packings from protocol 2 populate segments with $d\varphi /d\gamma <0$ only. For both simple and pure shear, we derive a relationship between the stress anisotropy and dilatancy $d\varphi /d\gamma$ obeyed by MS packings along geometrical families. We show that for MS packings prepared using isotropic compression, the stress anisotropy distribution is Gaussian centered at zero with a standard deviation that decreases with increasing system size. For shear jammed MS packings, the stress anisotropy distribution is a convolution of Weibull distributions that depend on strain, which has a nonzero average and standard deviation in the large-system limit. We also develop a framework to calculate the stress anisotropy distribution for packings generated via protocol 2 in terms of the stress anisotropy distribution for packings generated via protocol 1. These results emphasize that for repulsive frictionless disks, different packing-generation protocols give rise to different MS packing probabilities, which lead to differences in macroscopic properties of MS packings.

1. change the vertical axis label in Fig. 7 (b), fix the Eq. 13 in the revised manuscript; 2. Add new results of pure shear case in Appendix A; 3. page 3, fix typos, use d$\varphi$ in the revised manuscript.; 4. define the pure shear strain as $\gamma =\ln (L{x}^{\prime }/L{y}^{\prime })$, which makes the relationship between the components of the stress tensor and dilatancy for simple and pure shear the same;