Start-Up Shear Of Spherocylinder Packings: Effect Of Friction

We study the response to shear deformations of packings of long spherocylindrical particles that interact via frictional forces with friction coefficient mu. The packings are produced and deformed with the help of molecular dynamics simulations combined with minimization techniques performed on a GPU. We calculate the linear shear modulus g(infinity), which is orders of magnitude larger than the modulus g(0) in the corresponding frictionless system. The motion of the particles responsible for these large frictional forces is governed by and increases with the length l of the spherocylinders. One consequence of this motion is that the shear modulus g(infinity) approaches a finite value in the limit l -> infinity, even though the density of the packings vanishes, rho proportional to l(-2). By way of contrast, the frictionless modulus decreases to zero, g(0) similar to l(-2), in accordance with the behavior of density. Increasing the strain beyond a value gamma(c) similar to mu, the packing strain weakens from the large frictional to the smaller frictionless modulus when contacts saturate at the Coulomb inequality and start to slide. In this regime, sliding friction contributes a "yield stress" sigma(y) = g(infinity)gamma(c) and the stress behaves as sigma = sigma(y) + g(0)gamma. The interplay between static and sliding friction gives rise to hysteresis in oscillatory shear simulations.