# Highly strained mixtures of bidimensional soft and rigid grains: an experimental approach from the local scale

Jonathan Barés, Manuel Cárdenas-Barrantes, David Cantor, Emilien Azéma,
Mathieu Renouf

<figure><img src="image_1.png" /></figure>

## To cite this version:

Jonathan Barés, Manuel Cárdenas-Barrantes, David Cantor, Emilien Azéma, Mathieu Renouf. Highly strained mixtures of bidimensional soft and rigid grains: an experimental approach from the local scale. Powders & Grains, Jul 2021, Buenos Aires (virtual), Argentina. pp.05004, 10.1051/epj-conf/202124905004. hal-03259819

HAL Id: hal-03259819

https://hal.science/hal-03259819

Submitted on 14 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
# Highly strained mixtures of bidimensional soft and rigid grains: an experimental approach from the local scale

Jonathan Barés<fnref n="1" />,*, Manuel Cárdenas-Barrantes<fnref n="1" />, David Cantor<fnref n="2" />, Émilien Azéma<fnref n="1" />,<fnref n="3" />, and Mathieu Renouf<fnref n="1" />

**Abstract.** Granular systems are not always homogeneous and can be composed of grains with very different mechanical properties. To improve our understanding of the behavior of real granular systems, in this experimental study, we compress 2D bidisperse systems made of both soft and rigid grains. By means of a recently developed experimental set-up, from the measurement of the displacement field we can follow all the mechanical observables of this granular medium from the inside of each particle up-to the whole system scale. We are able to detect the jamming transition from these observables and study their evolution deep in the jammed state for packing fractions as high as 0.915. We show the uniqueness of the behavior of such a system, in which way it is similar to purely soft or rigid systems and how it is different from them. This study constitutes the first step toward a better understanding of mechanical behavior of granular materials that are polydisperse in terms of grain rheology.

## 1 Introduction

Loaded granular materials made of highly deformable particles are ubiquitous in nature and industry from cell monolayer growth [1, 2] to foam and emulsion compression [3–5] and even metal or plastic powder sintering [6, 7]. In some cases, because of impurities or for technical and recycling purposes, rigid, undeformable particles are mixed among soft ones [8, 9]. This is the case for example in sand-rubber granular layers used as low-cost seismic isolations [10].

These kinds of mixed soft-rigid granular systems display singular behaviors different from purely rigid or purely soft systems. Among other purposes, modifying the mixing ratio and the particle nature permits to tune the mechanical properties of a granular medium. In order to optimize their industrial applications and to control their behavior in natural phenomena, it is needed to know the connection between their macroscopic behaviors and their microscale composition and organisation. Because granular matter is naturally multiscale, a clear understanding of its global behavior necessitates to figure out and to measure what happens at the local/grain scale and at all the intermediate scales up-to the global one. Due to the already complex process of measuring local mechanical properties of each grain for fewly deformed particles [11] the challenge is even more important of measuring them for soft highly deformed particles.

Very recently this challenge has been overcome [12, 13] using an experimental method coupling Digital Image Correlation (DIC) and very acurate imaging. In this paper,

we extend this method to the so called case of mixtures of soft and rigid particles. We show how this experimental set-up is modified to compress these mixtures and how the image post-processing is modified accordingly. A selection of resulting experimental data from the inside of the particle up-to the system size is presented. This allows us to understand better the mechanical behavior of these materials and will permit to improve their compression models [14].

## 2 Experimental method

In Fig. 1 we show the experimental set-up already presented in [13]. It is composed of a bidimensional piston of initial dimension 270 × 202 mm<fnref n="2" /> which permits to compress uniaxially the system. A stepper motor rotates a screw which translate the moving edge of the piston inward in an accurate manner (25 µm per motor step). Two force sensors are attached to this latter, they record the global stress evolution (σ) when compressing a mixed granular system. The global stress is measured continuously with a frequency of 100 Hz while the system is compressed stepwisely. At each step the piston moves of 0.5 mm at a speed of 2 mm/min to make it sure to stay in the quasistatic regime. Then, the system rests during 1 min to relax and is imaged from below with a flatbed [15] with a resolution of 2400 dpi (10.6 µm/px).

The granular systems are composed of soft hyperelastic cylinders made of silicone [13, 16] and of rigid cylinders made of PVC (see Fig. 1). Their Young moduli are 0.45 MPa [18] and 1.2 GPa, respectively. Both soft

<footnote n="1">¹LMGC, Université de Montpellier, CNRS, Montpellier, France</footnote>

<footnote n="2">²Department of Civil, Geological and Mining Engineering, Polytechnique Montréal, Québec, Canada</footnote>

<footnote n="3">³Institut Universitaire de France (IUF), Paris, France</footnote>

<footnote>*e-mail: jb@jonathan-bares.eu</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Experimental set-up. Top: 3D drawing of the compression device and particle packing laying on a photographic flatbed scanner [15]. Bottom-left: squeezed soft silicone particle [13, 16]. The top shiny aspect is due to a layer of thin metallic glitter [17]. Bottom-right: rigid PVC particle. The top shiny aspect is due to silver paint. During experiments particles are used with their shiny part in contact with the scanner.</figcaption></figure>

and rigid granular packings are bidisperse cylinders of diameters 20 mm and 30 mm and of height 15 mm. For the experiments presented in this paper, 96 particles are used, 1/2 of them being soft and 2/3 of them being small. This last ratio is equally spread among soft and rigid particles, this avoids phase crystallization. The ratio of soft to rigid particles, $\kappa = n_{soft}/n_{total}$ is varied from 0.2 to 0.8. To avoid basal friction the scanner glass is coated with oil which also lubricate the contacts making particles almost frictionless. Namely the measured friction coefficient is $0.05 \pm 0.02$.

To measure displacement fields on the particle bottom faces, the DIC method requires a thin random pattern with a high optical contrast. In the case of soft particles, this is obtained by coating the silicone with a very thin layer of shiny glitter [17]. For rigid particles, PVC is simply coated with silver metallic paint. In both cases, the glitter characteristic size is around 25 µm. This creates a random pattern with a correlation length of about 50 µm on the bottom face of each particle.

For each experiment a set of about 100 pictures is obtained with this experimental set-up by compressing the system over 15%. Pictures are then post-processed with an algorithm modified from [13, 19]. First, particles are detected individually by thresholding the initial/undeformed picture. Each particule is then tracked in the full set of pictures by using image correlation. Then, from the measurement of the particle rotation and translation, subsets of im-

ages are extracted following each particle and correcting its solid rigid motion. For soft particles a DIC algorithm aimed for large deformations already presented in [18] is used to obtain the displacement field (u) inside each particle. For the rigid particles a more classical DIC algorithm is used. This latter correlates all images with the initial one as classically done for systems in the small deformation assumption. In both cases the correlation cell size is self-adjusted to take random pattern inhomogeneities into account as already presented in [13].

Then, for soft particles the deformation gradient tensor is computed from u [20]:

$$F = \nabla u + I \quad (1)$$

with I being the second order identity tensor. The right Cauchy-Green strain tensor is finally derived [20]:

$$C = F^T F \quad (2)$$

For rigid particles, the infinitesimal strain tensor is just directly computed ($\epsilon = 1/2((\nabla u)^T + \nabla u)$). From u the evolution of the exact particle boundaries is also deduced which permits to get all geometric observables for both particles and voids, but also contacts, at each compression step. Contacts are detected by considering first the proximity of each point of the particles boundaries with the boundaries of neighbor particles. Then, for each point closer than a certain threshold the von Mises strain around this point is measured and the contact is assumed to be real if this value is above a certain threshold as detailed in [13]. For the sake of conciseness only a meaningful selection of the different mechanical and geometrical observables are presented here.

## 3 Results and discussions

As presented in Fig. 2, a granular system with $\kappa = 0.5$ is compressed for a packing fraction ($\phi$) going from 0.795 to 0.915. While loading, the global stress first stays close to 0 MPa and then increases progressively from the jamming point [21]. This occurs for $\phi$ close to $\phi_j = 0.84$, as already observed in many other experimental studies [22]. It also corresponds with a coordination value slightly above 3 as predicted theoretically [23] since particles are slightly frictional. In Fig. 2 we also plot the evolution of the specific contact length, which is the average of the ratio of the particle perimeter in contact with another particle or with the system boundaries ($l_{contact}$) over the particle perimeter: $l = \langle l_{contact}/\text{perimeter} \rangle_{particles}$. This quantity increases linearly with $\phi$, which is not in agreement with the Hertz's contact law in 2D [24]. This last point is not surprising since soft particles are far from the small deformation assumption in our study.

Figure 2 also shows that the average value of the von Mises strain ($|C|$) in soft particles decreases parabolically all along the compression while the average soft particle asphericity first gently increases and then sharply rises after $\phi = 0.88$. This quantity is computed as the ratio between the actual particle area and the area deduced from
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Global observables. Variation of the global stress ($\sigma$), coordination ($Z$), specific contact length ($l$), average value of the von Mises of the right Cauchy-Green strain tensor ($|C|$) of soft particles, particle asphericity ($a$) and void area ($A$) as a function of the packing fraction ($\phi$) for a compression experiment made with 50% of rigid particles ($\kappa = 0.5$).</figcaption></figure>

the perimeter assuming the particle would be circular. The packing fraction $\phi = 0.88$ also corresponds with the maximum value of the total void area which first increases, since void clusters form, and then decreases, since particles are squeezed and reduce these cluster areas, as already observed for fully soft systems [13].

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Local observables. Field of the first invariant of the right Cauchy-Green strain tensor ($C$) for soft particles and of the infinitesimal strain tensor ($\varepsilon$) for rigid particles. These fields are computed for a compression experiment made with 50% of rigid particles at a packing fraction $\phi = 0.89$ far above the jamming point.</figcaption></figure>

In Fig. 3, the first invariant of the strain fields measured in each particle is shown for $\phi = 0.89$. For soft particles the right Cauchy-Green strain tensor ($C$) is used while it is the infinitesimal strain tensor ($\varepsilon$) which is used for rigid particles. The field for soft particles distinctly shows high strain chains with clear areas near the contact points where the deformation is maximum. On the contrary the patterns observed on the rigid particles suggest that the imaging method is not accurate enough to highlight grain deformations. This is not surprising since the level of compression is rather low and the ratio of soft to rigid particle stiffness covers several orders of magnitude: deformation is almost exclusively carried by soft particles. The lines that appear on the rigid particles actually correspond with the small speed inhomogeneities of the imaging method.

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Statistical observation. Main panel: Evolution of the PDF of the value of the von Mises of the right Cauchy-Green strain tensor ($C$) measured in soft particles as a function of the packing fraction for a compression experiment made with 50% of rigid particles. Scale is semi-log. Inset: Evolution of standard deviation of $C$ for different particle mixes ($\kappa$) as function of the distance to the jamming point.</figcaption></figure>

In Fig. 4, we present the evolution of the strain field statistics. More exactly, we plot the probability density functions (PDF) of the von Mises values of the right Cauchy-Green strain tensor ($|C|$) in soft particles for different packing fractions. For gently squeezed states, the PDF is a Gaussian centered around 1 which corresponds with the fact that the particles are mainly undeformed. While the system is loaded, the PDF gets broader and broader and the average value slightly shifts to the left which corresponds with the global compression of the system. For densely packed states, the PDFs display an exponential tail with a coefficient decreasing with the packing fraction. This exponential tail is reminiscent of what is observed for force values in packing of rigid granular systems [25] and remains unchanged from observations made on purely soft systems [13]. In the inset of Fig. 4, the evolution of the standard deviation of $C$ is shown for different ratios between soft and rigid particles. When the system is mainly
composed of rigid particles ($\kappa = 0.2$), this value rapidly increases right after the jamming point. In this case lot of particles undergo large deformation and display more extreme strain values. On the contrary, when the system is composed of more soft particles ($\kappa = 0.8$) the deformation is shared over more particles and less extreme strain values are observed so the standard deviation increases less rapidly.

## 4 Conclusion

In summary, we have developed an experimental tool that can permit to measure all mechanical and geometrical observables from the displacement field at the particle scale for bidimensional granular systems under uniaxial compression, whatever the particle mechanical properties. We believe this measurement method at the *sub-micro scale* will permits to model granular matter in a more efficient way to predict its behavior. It is applied here to the case of a mixed granular media made of soft squeezable particles and rigid ones. The system is loaded step by step and imaged very accurately, while the global force applied to it is recorded continuously. By mean of DIC, the evolution of mechanical fields inside each particle are measured which permits to deduce observables as diverse as packing fraction, coordination, contact length, particle asphericity, void area and strain fields to name a few.

More specifically we have shown that our system crosses the jamming transition for a packing fraction and a coordination value close from what has been observed in the case of purely rigid systems. This point validate our experimental approach. We have also shown that the average von Mises strain in soft particles decreases with the exact same tendency as the global stress increases in the system. Since rigid particle deformation is negligible, this suggests that rigid particles do not carry more load than the soft ones. The loading seems to be balanced between soft and rigid particles. The stress vs. packing fraction curve is also similar to what has been observed in a recent numerical study [14].

Deep in the jammed state the PDF of the von Mises strain tensor displays an exponential tail as already observed for purely soft systems and similarly to what is observed for the intergranular force PDF of purely rigid systems. Also, the standard deviation of this quantity increases more rapidly for systems composed of more rigid particles. This is explained by the fact that deformation is mainly carried by soft particles and that they are individually more deformed when they are less for a given deformation level.

Beyond the study of mixed soft-rigid granular systems, we believe that this experimental tool offers an opportunity to catch experimentally the behavior of breakable and polydisperse (in size, shape and material) particles under different loading geometries.

Gille Camp, Stéphan Devic and Rémy Mozul are greatly thanked for their technical support. We warmly thank Vincent

Huon for fruitful discussions.

## References

[1] T.E. Angelini, E. Hannezo, X. Trepat, M. Marquez, J.J. Fredberg, D.A. Weitz, PNAS **108**, 4714 (2011)

[2] P. Martin, S.M. Parkhurst, Dev. **131**, 3021 (2004)

[3] F. Bolton, D. Weaire, Phys. Rev. Lett. **65**, 3449 (1990)

[4] G. Katgert, M. van Hecke, Europhys. Lett. **92**, 34002 (2010)

[5] J. Brujić, C. Song, P. Wang, C. Briscoe, G. Marty, H.A. Makse, Phys. Rev. Lett. **98**, 248001 (2007)

[6] K.T. Kim, M.M. Carroll, Int. J. of Plast. **3**, 63 (1987)

[7] J. Barés, M.L. Hattali, D. Dalmas, D. Bonamy, Phys. Rev. Lett. **113**, 264301 (2014)

[8] W.Y. Tsoi, K.M. Lee, Géotech. **61**, 133 (2011)

[9] H. Khatami, A. Deng, M. Jaksa, Geosyn. Intern. **27**, 432 (2020)

[10] A. Tsiavos, N.A. Alexander, A. Diambra, E. Ibraim, P.J. Vardanega, A. Gonzalez-Buelga, A. Sextos, Soil Dyn. and Earth. Eng. **125**, 105731 (2019)

[11] A.A. Zadeh, J. Barés, T. Brzinski, K.E. Daniels, J. Dijksman, N. Docqiuer, H. Everitt, J. Kollmer, O. Lantsoght, D. Wang et al., GM **21**, 83 (2019)

[12] T.L. Vu, J. Barés, S. Mora, S. Nezamabadi, Phys. Rev. E **99**, 062903 (2019)

[13] T.L. Vu, J. Barés, Phys. Rev. E **100**, 042907 (2019)

[14] M. Cárdenas-Barrantes, D. Cantor, J. Barés, M. Renouf, E. Azéma, Phys. Rev. E **102**, 032904 (2020)

[15] Canon™ canoscan 9000f mark ii

[16] Moldstar™ 15 slow, castable silicone from smooth-on, https://www.smooth-on.com/products/mold-star-15-slow/

[17] Cast magic™, silver thin metallic glitter from smooth-on, https://www.smooth-on.com/products/cast-magic/

[18] T.L. Vu, J. Barés, S. Mora, S. Nezamabadi, Expe. Mech. **59**, 453 (2019)

[19] Python codes to post-process experimental images (particle tracking, dic, fields computation, etc.), https://git-xen.lmgc.univ-montp2.fr/Image/granular dic

[20] L.A. Taber, *Nonlinear theory of elasticity: applications in biomechanics* (World Scientific, 2004)

[21] A.J. Liu, S.R. Nagel, Nat. **396**, 21 (1998)

[22] R.P. Behringer, B. Chakraborty, Rep. on Prog. in Phys. **82**, 012601 (2018)

[23] M. van Hecke, J. of Phys. CM **22**, 033101 (2009)

[24] K.L. Johnson, *Contact Mechanics* (Cambridge University Press, Cambridge, 1985)

[25] F. F. Radjai, M. Jean, J.J. Moreau, S. Roux, Phys. Rev. Lett. **77**, 274 (1996)