Shear zones in granular mixtures of hard and soft particles with high and low friction

Granular materials show inhomogeneous flows characterized by strain localization. When strain is localized in a sheared granular material, rigid regions of a nearly undeformed state are separated by shear bands, where the material yields and flows. The characteristics of the shear bands are determined by the geometry of the system, the micromechanical material properties, and the kinematics at the particle level. For a split-bottom shear cell, recent experimental work has shown that mixtures of hard, frictional and soft, nearly frictionless particles exhibit wider shear zones than samples with only one of the two components. To explain this finding, we investigate the shear zone properties and the stress response of granular mixtures using discrete element simulations. We show that both interparticle friction and elastic modulus determine the shear-band properties of granular mixtures of various mixing ratios, but their stress response depends strongly on the interparticle friction. Our study provides a fundamental understanding of the micromechanics of shear band formation in granular mixtures.


Introduction
Strain localization is a characteristic effect of sheared granular matter, which is not observed in the flow of other materials, such as liquids.Once strain localizes within a narrow shear zone (also referred to as the shear band), the stress, strain, and strain rate of the material assume high values within that zone.On the contrary, outside of it, the material remains nearly undeformed.The knowledge of the properties of shear zones and the factors that influence them is crucial for understanding the intrinsic properties of granular materials under shear.
In many systems where particles come into contact with boundaries, narrow shear zones appear at the interface because of particle-boundary interactions.These shear zones are triggered by boundary effects, which hinder the accurate characterization of bulk material properties based on the shear zone characteristics.In our work, we present a unique setup called a split-bottom shear cell, which is a modified version of a Couette cell.This setup involves splitting the cell's bottom at a specific radius (R s ), resulting in two parts consisting of the inner disk and outer cylinder components (see Figure 1).After loading the cell with grains to a certain height (H), the inner disk is set to rotational motion.This design allows granular materials to be sheared to develop wide shear zones inside the bulk of the materials, far away from the rigid boundaries, which allows for a comprehensive characteriza-tion of their yielding properties and flow behavior at a macromechanical level.This is quite different from the usual narrow shear bands that we observe close to the boundary, which are formed due to particle-boundary interactions.This shear cell setup has been used in experiments 1 and simulations 2 to explore the bulk properties of granular systems.
The universal geometrical properties of shear bands in the splitbottom shear cell have been studied extensively in the literature [3][4][5] .The absolute value of the macroscopic normalized angular velocity of the granular material was measured as a function of the radial position, and it has been found to follow closely the shape of an error function, which is characterized by two parameters: the width W and the center position R c of the shear zone 4,5 .The width of the shear band increases with increasing height within the bulk of the material, while the distance of its central position R c from the center of rotation decreases (that is, the shear band moves inward with increasing height).The width of the shear band is affected by the size, shape, and interaction properties of the particles, such as friction between the particles, cohesion, and the shear rate [6][7][8][9] .The position of the shear band, R c , at the surface of the granular material in the shear cell, depends on two geometrical parameters, the filling height H and split radius R s 10 .The bulk profile of the shear band as a function of height, R c (z), can be linked to a variational problem.The solution of the function R c (z) that minimizes the integral representing the dissipation rate, is related to the minimal mechanical torque and thus depends on the material properties.Therefore, the shape of the shear band within the bulk is not only influenced by geometric factors such as the split position (R s ) and the filling height (H) but is also closely connected to the material properties, e.g., the shear stress, and the packing fraction of the material, as highlighted in previous studies 5,10 .Interparticle friction and particle softness play a key role in influencing material bulk flows, which have been widely studied.Wang et al. 11 experimentally characterized the shear zones of binary mixtures made of hard, frictional and soft, slippery particles.Their study identified distinct shear band characteristics specific to the mixture, distinguishing them from those in pure-species systems.Luding 8 explored the impact of interparticle friction on shear zones in split-bottom shear cells through numerical simulations.Increasing friction resulted in narrower shear bands that shifted towards the center.Ashour et al. 12 found that particle softness adds qualitatively new features to the dynamics and the packing characteristics of silo flows.Soft, deformable particles, such as hydrogel particles, can adapt to local stress conditions, leading to alterations in pressure profiles and flow fields.Hong et al. 13 highlighted the importance of particle softness in hopper flows, where gravitational forces must be less than repulsive forces from particle stiffness to induce clogging.Harth et al. 14 observed longer time scales for internal rearrangement in soft hydrogel particles during silo orifice flow compared to hard grain systems, where shorter time scales have been reported by Unac et al. 15 .This was attributed to mass flow rate differences, roundness, and increased viscous damping of hydrogel grains.The rate-dependent frictional regimes observed in hydrogel particle suspensions have been attributed to their high deformability and low friction 16 .Wang et al. 17 observed experimentally that, unlike granular materials with hard particles, soft slippery particles in a silo flow do not follow Beverloo's 18 law, according to which the flow rate is independent of the overburden pressure.Adding just 5% hard particles to soft silo flows restored this behavior and significantly altered system flow properties.Based on the micromechanical observations, Saitoh et al. 19 highlighted that the evolution of the system behavior in soft particle packings extends beyond just particle contacts.Inter-particle gaps between nearest neighbors are also crucial and affect the stochastic evolution of the probability distribution function of forces.Recently, Luding et al. 20 studied the contact and force networks of frictionless soft granular materials, revealing that the evolution of the pressure in the system was correlated with the existence of loops formed in the force network during particle rearrangements.Götz et al. 21observed that soft particles with larger deformations and increased contact areas result in larger frictional forces and improved holding capacity of a granular gripper.In summary, both particle friction and softness are crucial factors influencing bulk flow behavior, and the intricate interplay between these properties and their combined effects on bulk behavior remains a subject of ongoing investigation.
The findings of Wang et al. 11 established that mixtures of hard and soft particles exhibit wider shear zones than assemblies made of only one of the individual materials (only hard or only soft particles).However, the underlying reason for this behavior remains unclear.Our objective here is to gain insight into the combined effect of microscopic parameters, namely particle stiffness and inter-particle friction, on the macroscopic shear zone character- The bottom plate (dark-red color) rotates around the vertical axis at constant rotation frequency ω o , while the outer wall (marked in grey) stays stationary.The differential motion triggers the formation of a shear band (marked in orange color) starting from the split radius R s , which widens with increasing height.In contrast, the material outside the shear band (marked in yellow) stays undeformed.
istics for ensembles of soft, slippery particles and hard, frictional particles sheared in a split-bottom shear cell.

Numerical setup and materials
We use an open-source Discrete Element Method (DEM) code MercuryDPM 22 to simulate sheared systems of hard and soft particles.The numerical setup used for this study is a splitbottom shear cell, as shown in Figure 1, where the inner bottom plate rotates around a central axis at a constant angular velocity, ω o = 0.52 rad/s, while the outer boundary of the cell remains stationary.The granular material located in the outer parts remains static, while the material near the center of the cell rotates with the bottom disk.The shear in this system is localized in an intermediate region, defining a wide shear band.The position and the width of the shear zone depend on the ratio H/R s , where H is the filling height and R s is the split radius of the rotating bottom plate 4,23 .For small and moderate filling height (H < 0.6R s ), the shear zone starts from the bottom split position, evolves through the entire height of the granular layer, and ends up at the surface 3 , as shown in Figure 1, while the inner part of the granular material co-rotates synchronously with the bottom plate.In our numerical setup, the filling height is H = 20 mm and H/R s = 0.42, which is considered a moderate filling height.This numerical setup is similar to that employed in the experimental work by Wang et al. 11 , with the difference that the experimental configuration incorporates an additional small rod in the center of the cell, where the shear load was imposed.

Material parameters
We consider a mixture of soft low-frictional particles and hard frictional particles.In this Section, we describe the generation of binary granular mixtures of such particles.We take the proper- ties of hydrogel spheres as the soft, low-frictional material component, and mustard seeds represent hard, frictional particles.
A reference sample of pure soft hydrogel particles consists of N 0 A = 26795 polydisperse particles of species A with mean diameter, d A = 2.5 mm.Starting from this baseline sample volume, we create mixtures by substituting N 0 A − N T A particles of species A with particles of species B of diameter d B = 3 mm, such that the same volume of material A is replaced by B, i.e.
For both species A and B, we use particles of homogeneous size distribution of mean particle diameter d A and d B , respectively, with polydispersity of ±8%.
To characterize the mixtures with different ratios of species B, we define the mixing ratio We analyze systems with X ∈ {0, 0.05, 0.25, 0.50, 0.75, 0.85, 1.0}.Table 1 presents the number of particles of the two species in each sample.Figure 2 shows some examples of different mixtures.
The properties of soft particles of species A correspond to that of hydrogel spheres with friction coefficient µ A = 0.005 and elastic modulus E A = 50 kPa.The frictional properties of hard particles correspond to mustard seeds with friction coefficient µ B = 0.20 and the elastic modulus E B = 5000 kPa.All relevant material and interaction properties of the two species are given in Table 2.

Contact model
In the soft-particle DEM, the evolution of particle positions and velocities is calculated by integrating Newton's equation of motion.When two particles i and j interact, they can slightly overlap.Only viscoelastic interactions are considered, i.e., after separating the contacting particles, they recover their initial shape.The compression of contacting particles ξ is written as: where R i , R j are the particle radii and ⃗ r i , ⃗ r j are the particle positions.The normal force is calculated as in 24 , where an elastic Hertzian normal force is considered in conjunction with a viscous damper of velocity-dependent coefficient of restitution. with where E * is the effective elastic modulus, R * = R i R j / R i + R j is the effective radius, ⃗ e n = (⃗ r i −⃗ r j )/|⃗ r i −⃗ r j | is the normal unit vector, and A n = 4.1 × 10 −6 s is the normal dissipative parameter, calculated as in 25 , considering a coefficient of restitution of 0.70 for the tangential translational velocity of the rotating bottom plate at the split radius.The effective elastic modulus E * is given by where E i , E j and ν i , ν j are the elastic moduli and Poisson ratios of particles i, j, respectively.
We model the tangential viscoelastic forces considering a path-dependent calculation, following the no-slip solution of Mindlin 26 for the elastic part and Parteli and Pöschel 27 for the tangential dissipative parameter A t ≈ 2A n E * , which are capped by the static friction force between two particles.The tangential force is thereby given by where µ is the friction coefficient, ⃗ e t is the unit vector in the tangential contact direction and G * the effective shear modulus, given by The shear modulus G i = E i 2(1+ν i ) is calculated assuming isotropic material.For the interaction between different species, the effective friction coefficient is considered as . It is typical for many DEM formulations to consider the minimum coefficient of friction when two particles of different materials come into contact.Since the harmonic average of two values favors the smaller one, considering this calculation of the friction coefficient for contacts between high and low friction coefficients gives values close to considering the minimum value.This assumption is chosen as a milder approach to simply assuming the minimum friction coefficient of the two, and more experimental evidence is required to understand the friction between hard and soft particles.Considering the harmonic average here, we have three scenarios of possible contacts between frictional-frictional, slippery-slippery and frictional-slippery (with friction very close to slippery-slippery) particles.

Segregation effects
When granular mixtures with particles differing in size and/or density are subjected to external forces such as shaking, stirring, or shearing, they often segregate to form complex patterns.These patterns of segregation can lead to different and inhomogeneous macroscopic behavior.In many applications where a homogeneous packing is sought, this is undesirable, and significant effort is often invested in attempts to avoid or control segregation during processing and handling of granular materials.Shearinduced percolation segregation that occurs due to kinetic sieving and squeeze expulsion is the dominant mechanism in dense granular flows.In this study, we study a bi-disperse mixture of two species that differ slightly in terms of particle size.Therefore, in this section, we examine the extent to which particles segregate in this system of hard frictional and soft slippery particles, sheared in the split-bottom shear cell.First, we analyze the probability distribution of the particles of each species in the radial r and vertical z directions for a given snapshot after 200 s of simulation time, as shown in Figure 3(a) and (b), respectively.The distribution increases in the radial direction for both species due to the increased measuring volume (the probability is calculated using co-centric cylinders of increasing radius, aligned with the center of the cell).Note that the distribution of the two species overlaps with each other in the radial r direction, indicating that there is no segregation effect of the two species horizontally.The distribution in the vertical z direction for species B is nearly constant.However, the distribution of species A decreases in vertical z direction from 0.065 to 0.035, indicating a higher probability of species A near the bottom of the shear cell than near the free surface.This is a signature of the segregation of species A in the vertical direction due to kinetic sieving, and is attributed to the smaller particle size compared to species B. To further analyze this effect, Figure 4(a) and (b) illustrate the mean radial positions r and perpendicular positions z of the two species as a function of shear displacement λ .Here, it is to be noted that the mean position of Species A is about 0.4 mm less in the r direction and 1.2 mm less in the z direction than that of Species B. However, the differences are negligibly small and of less than one particle diameter.

Micro-macro transition
The current work aims to derive appropriate macroscopic fields commensurate with the shear band features and macromechanical stress analysis, e.g., strain rate, velocity profile, and stress, based on the given micromechanical properties.The flow that rapidly relaxed to a steady state is purely azimuthal and is proportional to the driving rate of the shear cell, ω o .The averaging is thus performed in toroidal volumes over many snapshots of time, leading to generic continuum fields Q(r, z) as a function of the radial and vertical positions.We choose the representative elementary volume according to the 28 as equal to the minimum particle diameter d p = 1.37 mm.
Shearing of a granular system leads to volume dilation and build-up of shear stress and anisotropy in fabric.Since we are interested in the steady-state flow profile, we probe for the minimum strain required to reach the steady-state flow regime.We define the shear displacement λ = ωt rot through which the system rotates (where t rot is the simulation time) as the shear control parameter.As reported by Singh et al. 29 , both the global quantities of kinetic energy and average number of contacts reach a steady state very fast (λ ∼ 5).We analyze the relaxation of local quantities, specifically the local velocity profiles, and conclude that λ ∼ 20 − 30 is the amount of shear displacement required to form a stable shear band, which is in agreement with Ries et al. 30 .
The data depicting the initial, transient evolution of local shear band properties is not presented here, as our focus is solely on the steady-state characteristics of the shear band.Consequently, we perform the local averaging over almost 600 snapshots distributed over λ ∼ 80 to 100.
The consistency of the local averaged quantities also depends on the accumulated local shear strain during the averaging time.We concentrate our interest in the region where the system can be considered to reach a critical state.The critical state is a unique steady state reached after extensive shear, where the material deforms with applied strain without any further change in normal stress, shear stress, and volume fraction.This is a state where the system forgets its loading, packing, or other sample-preparation history 31 .In accordance with previous experimental and numerical results in the same setup 32,33 , the criterion of large strain rate identifies the shear band region, e.g., higher than a critical strain rate of γc = 0.08 s −1 .Here, the shear band center region is defined by strain rates higher than 80% of the maximum strain rate for different heights in the shear cell.Therefore, the critical strain rate here is a function of height in the shear cell and is denoted as γc (z) = 0.8 γmax (z).

Flow profile of granular mixtures
Once the flow of the mixtures relaxes to steady state shearing, the azimuthal velocity is proportional to the driving rate ω o .Figure 5(a) shows the scaled angular velocity profile ω(r, z)/ω o .Within the limits of assumptions in our simulations, we do not see a clear difference in the flow profiles of the three cases of pure particles and mixtures for X = 0, 0.50, and 1.However, we observe that the shear zone of the mixture (X = 0.5) is slightly more diffused compared to the pure case simulations involving only soft (X = 0) or only hard particles (X = 1).This indicates a wider shear band for the mixture, although it is not clearly concluded from this figure.A closer and more precise introspection on the shear band properties is therefore essential to conclude the shear zone properties of the mixture.On this ground, we analyze the dimensionless ratio of the measured angular velocity ω/ω o as a function of the radial coordinate r.The dimensionless velocity profile in a split bottom shear cell is well approximated by an error function 5,10 where R c and W are the position and width of the shear band developed above the split at R s , respectively.Both the position and the width of the shear band depend on the height z in the system.Therefore, we extract the position and the width of the shear band at layers of varying height along the filling height in the shear cell by fitting the scaled ω(r) against the radial position r data.With the goal to extract quantitative data for the shear band features, Figure 5(b) shows values of non-dimensional angular velocity ω(r)/ω o at different heights in the bulk of granular materials against the radial distance.
A quantitative study of the velocity gradient leads to the strain rate γ, which is calculated from the shear-free sheet strain tensor components d 1 and d 2 34 .The shear intensity in the shear plane is  thus expressed in terms of velocity gradients ∂ v θ /∂ r and ∂ v θ /∂ z as: Figure 6(a) shows the strain rate γ(r, z) evaluated from Equation 9and plotted as a function of radial and vertical position for the sample with X = 0.50.The dashed lines overlaid on top of the contour graph indicate the center R c and the edges R c + W and R c − W of the shear band, indicated in the figure, as obtained from the fit function in Equation 8.
An important objective of this present study is to characterize the shear zone profiles of granular materials for ensembles of soft slippery and hard frictional materials.Figure 6(b) shows the shear band profiles for samples with X = 0, 0.50 and 1.As observed from the figure, the width of shear band corresponding to a mixture sample X = 0.50 is wider than that of the pure mixtures, corresponding to X = 0 and 1.Also, the position of the shear band corresponding to a mixture sample X = 0.50 is shifted closer to the center than the pure mixtures.Note that the profile for the scaled width W /W top as a function of z/H shown in Figure 6(c) does not follow a quadrant shape, as described by the analysis of Ries et al. 30 , for any of the samples.The reason behind this deviation needs further exploration.However, a possible reason for this deviation can be due to the low W /d and z/d ratios, where d is the mean particle diameter.The normalized position of the shear band for the different mixtures as a function of bulk height is shown in Figure 6(d), and the trends agree with that shown by Unger et al. 10 .From both Figure 6(c) and (d), we observe a marked deviation in the trend for X = 0.50 compared to the pure samples with X = 0 and 1. Next, we analyze the width and position of the shear band for a horizontal slice through the axis of the shear cell at z = z top = 16 mm, as represented by the dashed line in Figure 6(b), for all samples of different fraction of hard particles X and show them in Figure 6(e) and (f), respectively.Here, we find that the shear band gets wider as well as moves inwards in a mixture of soft and hard spheres up to a fraction X = 0.50, which corresponds to an intermediate mixture of 50% soft -50% hard particles.With a further increase of X, the shear band narrows and moves outwards in the shear cell.The non-monotonic trend of the width of the shear band with increasing X was also observed by Wang et al. 11 .
From the above observations, we conclude that the ratio of hard spheres in a mixture of hard frictional and soft slippery particles significantly affects the flow profiles at the bulk scale.For up to intermediate ratios of hard particles, the shear band gets wider by up to 22%, and also the center shifts inwards by up to 2%, compared to the values for pure soft particles.Both of these features are influenced by the effect of inter-particle friction and stiffness, which are dominant for frictional hard particles 8 .Wang et al. 11 observed experimentally a notable difference in the shear band properties between samples made of pure materials and mixtures.It was observed that the shear zones of the mixtures are wider than those of the pure materials.Unlike the experimental observations, which show a significant widening of the shear zone for mixtures, our simulations show a rather small effect of widening the shear zone.However, the distinct nonmonotonic trend of the shear band width as a function of the mixing ratio is the feature observed in experiments, which is also captured in our simulations.In addition, we observe the shift in shear band position for the mixtures, which were not reported in the experimental work.It is important to note that the capillary bridges resulting from small amounts of liquid between hydrogel particles potentially have an impact on experiments, despite not being taken into consideration in this numerical study.For larger quantities of liquid in the mixtures, beyond the capillary regime, the material would exhibit a more dispersed shear behavior.
6 Simplified cases to explore the role of interparticle friction and elastic modulus Following the findings of this study and the qualitative agreement of our results with the observations noted by Wang et al. 11 , it becomes evident that mixtures of soft, slippery and hard frictional grains demonstrate shear bands of variable width for increasing mixing ratio, which evolves non-monotonically with increasing values of the ratio.Since the materials of interest have two independent material parameters differing distinctly, namely friction and elastic modulus, it is interesting to look into their individual effects on the shear band characteristics.Hence, to explore the role of inter-particle friction and elastic modulus on the flow profile of the granular mixture systematically and independently of each other, we vary one parameter at a time in additional discrete element simulations of a parametric nature, while keeping the other parameter constant.Therefore, we conducted four case studies, as outlined in Table 3.In these studies, our primary focus is on analyzing the shear band characteristics at a specific height, denoted as z = z top in Figure 6(b) as a function of mixing ratio X.
We then proceed to compare the obtained results with those from the real simulations.

Shear band characteristics of granular mixtures
Luding 8 examined the effects of friction on shear bands via simulations of a split-bottom shear cell.He showed that the presence of inter-particle friction does not affect the qualitative behavior of the shear band.However, the shear band moves inwards and gets narrower in the presence of friction.This is justified by the fact that with increasing friction, there is less relative displacement between the particles, particles tend to move more as a rigid block, and hence, the shear band width decreases.Here, we vary the inter-particle friction of the mixture by varying the fraction of hard particles X, keeping the elastic modulus constant in cases E soft µ real and E hard µ real as implied by the red and blue curves, respectively, in Figure 7.For both the cases, we observe that the shear band width decreases with increasing X in Figure 7(a).Also, the shear band position moves inwards with increasing X as shown in Figure 7(b).Our results agree with the findings of Luding on the effects of friction on shear band characteristics.
On the contrary, varying the elastic modulus and considering same friction coefficient for the two materials leads to a wider shear band formation with increasing fraction of hard particles X,  as demonstrated by cases E real µ soft and E real µ hard .This is implied by the green and orange curves in Figure 7(c).However, for the same friction between the two species, varying only the elastic modulus of the species has only a small effect on the position of the shear band as a function of X as shown by the green and orange trends in Figure 7(d).
Overall, there is a dual effect of friction and elastic modulus that influences the shear band properties.While increasing friction tends to reduce the width of the shear band with increasing X, increasing elastic modulus plays the opposite role in increasing the width for increasing X.For smaller values of X, the friction of soft particles is smaller and hence is less dominant, i.e., it is easier for inter-particle contacts to slide.With increasing X, there are more contacts between hard frictional particles, and hence friction dominates over the effect of elastic modulus.

Macroscopic friction coefficient of granular mixtures
The local shear stress τ(r, z) plotted against the local normal stress P(r, z) is shown in Figure 8(a).Since the system is inhomogeneous in nature, for a given pressure, we observe a wide range of local strain rate values γ, and we find that shear stress τ is higher for increasing strain rate γ.However, for values above the critical strain rate, i.e., γ > γc , τ becomes almost independent of the local strain.This means that τ/p is nearly constant for all data points with strain rate γ > γc .A linear trend is observed for the shear stress as a function of the normal stress, which is fitted well by: where µ is the macroscopic friction coefficient obtained from linear fitting of the local shear stress and normal stress values.
We further analyze the macroscopic friction coefficient µ for different mixtures as a function of X for the real simulations and the simplified cases and show the results in Figure 8(b).The macroscopic friction coefficient µ increases as a quadratic function with X for the simplified cases of E soft µ real and E hard µ real , i.e., same elastic modulus for both materials and varying interparticle friction coefficient.On the contrary, µ increases linearly with X for cases E real µ soft and E real µ hard , i.e. varying elastic modulus and same inter-particle friction coefficient for both materials.Furthermore, the macroscopic friction coefficient µ for the real simulations also increases as a quadratic function with X, which suggests that µ for the real simulations is strongly dominated by the varying inter-particle friction coefficient between the two species.The elastic modulus does not play a significant role in influencing the µ of the bulk behavior.Focusing on the functional form of the quadratic function µ vs X for the real simulations, the black solid line represents the quadratic fitting equation µ = µ soft + aX 2 , where µ soft ≈ 0.14 is the macroscopic friction coefficient corresponding to the pure soft particle mixtures and a = 0.18 here depends on the inter-particle friction coefficient values of pure soft and hard species.

Conclusions
This work investigated the shear zone characteristics in binary mixtures of soft low frictional and hard frictional particles, via discrete element simulations of mixtures inside a split-bottom shear cell.The ratio of the volume of hard particles to the soft particles was varied parametrically to observe its effects on measures of interest, including the shear band width W , shear band position R c and the macroscopic friction coefficient µ.
By comparing the width and position of the shear zone for samples of pure soft low frictional and for samples of pure hard frictional particles, it becomes evident that both systems show somewhat similar shear band features.However, mixtures of the two materials exhibit broader shear bands and a shifting of the shear band position towards the center of the shear cell.Both these effects demonstrate non-monotonicity for increasing fraction of hard particles.Although the changes in shear band characteristics for different mixtures are small, the non-monotonic trend of the shear band features as a function of mixing ratio is an interesting effect of the shear behavior of binary mixtures of this nature.We attribute the effects of shear band widening in these mixtures to an interplay of the role of inter-particle friction and particle stiffness.
Contrary to the shear band characteristics, the macroscopic friction coefficient shows a monotonically increasing trend as a quadratic function with an increasing fraction of hard particles in the mixtures.The macroscopic friction coefficient is strongly influenced by the inter-particle friction coefficient of the species present in the mixture, and the elastic modulus seems to play a less significant role.
This work captures critical features of shear band formation in soft-hard particle mixtures of low-high inter-particle friction, such as the non-monotonic evolution of the shear band width and position for increasing mixing ratios.However, the sensitivity of these features in our simulations is minor as compared to that observed in the experiments by Wang et al. 11 .Some of the additional factors that might play a role in the experiments but are not captured in our simulations include the following: (i) our simulations only consider spherical particles and do not account for non-spherical shapes or the actual particle size distribution, (ii) the real stiffness of the mustard seeds is not taken into account; a sufficiently stiffer material is chosen for the mustard seeds instead, and (iii) capillary bridges due to small amounts of liquid between hydrogel particles might play a role in experiments which are not considered here.The material would exhibit a more dispersed shear behavior for larger quantities of water in the mixtures beyond the

Fig. 1
Fig. 1 Schematic of a split bottom shear cell of radius R o with granular material at filling height H, illustrating the main system components.The bottom plate (dark-red color) rotates around the vertical axis at constant rotation frequency ω o , while the outer wall (marked in grey) stays stationary.The differential motion triggers the formation of a shear band (marked in orange color) starting from the split radius R s , which widens with increasing height.In contrast, the material outside the shear band (marked in yellow) stays undeformed.

Fig. 2 Fig. 3
Fig. 2 Top view of the shear cell for granular mixtures of different percentages of hard particles X, where the blue and red particles represent species A (soft) and B (hard), respectively.

Fig. 4 Fig. 5 (
Fig. 4 Average (a) radial distance and (b) distance from the bottom surface for the species A and B as a function of shear displacement λ = ωt rot .

Fig. 6 (
Fig. 6 (a) Contour plot of strain rate γ(r, z) as a function of the radial and vertical position for the sample with X = 0.50.The lines overlaid on the figure indicate the center R c and the edges R c + W and R c − W of the shear band as obtained from the fit function in Equation 8. (b) Shear band center R c and edges R c +W and R C −W for different fractions of hard particles X, (c) scaled shear band width W /W top as a function of scaled height z/H and (d) relative position of shear band R c /R s as a function of scaled height z/R s for samples with X = 0, 0.50 and 1. (e) Shear band width W as a function of the fraction of hard particles X and (f) shear band position R c as a function of the fraction of hard particles X at z = z top .The error bars in the graph illustrate the fitting error, which was determined by fitting the numerical data of ω/ω o vs. r to the equation presented in Equation 8.The blue ⃝, black ♢ and red △ markers highlight the shear band width and position corresponding to X = 0, 0.50 and 1, respectively.

Fig. 7
Fig. 7 Shear band (a) width W and (b) position R c for fixed elastic modulus and real inter-particle friction, represented by cases E soft µ real and E hard µ real , as compared with the real simulations, and shear band (c) width W and (d) position R c for real elastic modulus and fixed inter-particle friction, represented by cases E real µ soft and E real µ hard , as compared with the real simulations, where all the features are plotted as a function of fraction of hard particles X.The width and position of the shear band are obtained at z top = 16 mm.The error bars represent the fitting error, determined by fitting the numerical data for angular velocity to Equation 8.

Fig. 8 (
Fig. 8 (a) Shear stress τ plotted against normal stress P. The lines represent linear regression fits given by τ = µP.(b) Macroscopic friction coefficient µ as a function of the fraction of hard particles X for the real simulations and the simplified cases.The solid lines indicate linear and quadratic regression fits for the different cases.The error bars in the graph illustrate the fitting error, which was determined by fitting the numerical data of µ Vs P to the equation presented in Equation 10.

Table 1
Number of particles of the species in each sample

Table 2
Microscopic material parameters for the model

Table 3
Simplified cases to explore the effects of friction and elastic modulus on shear band characteristics separately E soft µ real E = 50 kPa for both species µ soft = 0.005 | µ hard = 0.20 E hard µ real E = 5000 kPa for both species µ soft = 0.005 | µ hard = 0.20 E real µ soft E soft = 50 kPa | E hard = 5000 kPa µ = 0.005 for both species E real µ hard E soft = 50 kPa | E hard = 5000 kPa µ = 0.20 for both species