Abstract
The physics of granular materials, including rheology and jamming, is strongly influenced by cohesive forces between the constituent grains. Despite significant progress in understanding the mechanical properties of granular materials, it remains unresolved how the range and strength of cohesive interactions influence mechanical failure or avalanches. In this study, we use molecular dynamics simulations to investigate simple shear flows of soft cohesive particles. The particles are coated with thin sticky layers, and both the range and strength of cohesive interactions are determined by the layer thickness. We examine shear strength, force chains, particle displacements, and avalanches, and find that these quantities change drastically even when the thickness of the sticky layers is only 1% of the particle diameter. We also analyze avalanche statistics and find that the avalanche size, maximum stress drop rate, and dimensionless avalanche duration are related by scaling laws. Remarkably, the scaling exponents of the scaling laws are independent of the layer thickness but differ from the predictions of mean-field theory. Furthermore, the power-law exponents for the avalanche size distribution and the distribution of the dimensionless avalanche duration are universal but do not agree with mean-field predictions. We confirm that the exponents estimated from numerical data are mutually consistent. In addition, we show that particle displacements at mechanical failure tend to be localized when the cohesive forces are sufficiently strong.
1 Introduction
Mechanics of granular materials is of great importance in technologies for sands, foods, and pharmaceutical products [1, 2]. Except for well controlled laboratory experiments, granular materials in nature are usually “wet” with water [3]. Wet granular materials consist of sticky particles, where interactions between them are cohesive due to liquid bridges formed at their contact points [4]. It is known that cohesive interactions strongly influence mechanical properties of granular materials [5]; the critical angle and the angle of repose for landslides significantly increase with the increase of amount of water (or layer thickness), the shear strength (stress) increases with the increase of suction1, the critical acceleration for fluidization of vibrated granular beds increases with the increase of liquid content, and segregation is suppressed and hysteresis is enhanced by the cohesive interactions. Furthermore, it has been suggested that, if granular materials are wet, jamming occurs at low packing fractions and inhomogeneity, i.e., localization of particle motions, is more pronounced [3].
One of the fundamental problems of granular matter is mechanical failure or avalanche which can be related to sediment disasters and earthquakes [6–8]. In seismology, the frequency of earthquake magnitude is explained by the celebrated Gutenberg-Richter (GR) law [7, 8]. As the GR law, statistical properties of mechanical failure are of central interest to physicists, where statistics of avalanches have been studied in the context of self-organized criticality (SOC) [9] or non-equilibrium phase transitions [10]. The SOC indicated by power-law distributions of avalanche (cluster) sizes is realized by a cascade of local failure. Associated the cascade of local mechanical failure with the depinning transition [11], power-law scaling of avalanche size distribution was suggested by a mean-field (MF) theory [12–15]. The MF theory also predicts several scaling laws for slip avalanches and its predictions (including the power-law scaling of avalanche size distribution) have been validated by many experiments of, e.g., granular materials under shear [16–19], compressed nano-crystals [20, 21], bulk metallic glasses [22–25], and light flux from a star [26]. Therefore, the statistics of avalanches have been said to be universal, in the sense that scaling exponents for the avalanche size distribution and other quantities do not depend on any details of materials on a microscopic scale.
In addition to experiments, the statistics of avalanches in granular materials have widely been studied by numerical simulations. Nevertheless, the avalanche size distributions extracted from numerical data quantitatively differ from the MF prediction. For example, the power-law exponents for avalanche size distribution found in molecular dynamics (MD) simulations of foams [27] and athermal quasi-static (AQS) simulations of amorphous solids [28–32] are much smaller than the MF prediction. Moreover, different from the MF theory, a mesoscopic elasto-plastic (EP) model was developed on the basis of yielding transition [33]. The EP model includes a “quadrupolar” elastic propagator in its governing equation and predicts a smaller power-law exponent for avalanche size distribution [34]. Thus, there still exist discrepancies in the theories, experiments, and simulations, and researchers have carefully examined the roles of system size [35], strain rate [36–39], particle inertia [40–43], friction [44, 45], and particle shapes [46, 47] in the statistics of avalanches. However, much less attention has been paid to the influence of cohesive interactions, which are crucial to real granular materials.
In this paper, we carry out numerical simulations of soft cohesive particles under shear. The main aim of our simulations is to clarify effects of cohesive forces (between the particles) on the statistics of avalanches. We employ a cohesive contact model which has been used for the studies of rheology [48–52] and jamming [53–55] of two-dimensional cohesive particles. We assume that our system is in a pendular state, i.e., liquid bridges are formed at contact points so that cohesive interactions are pairwise [4]. We show that not only the statistics of avalanches but also mechanical responses, force-chains, and particle rearrangements are affected by the cohesive interactions even if their range is only 1% of particle diameter. In the following, we explain our numerical methods (Section 2), show our results (Section 3), and discuss our findings (Section 4). All the details of our simulations and supporting data are summarized in Supplementary Material (SM).
2 Numerical methods
In this section, we introduce our numerical methods. We study simple shear deformations of soft cohesive particles in two dimensions by MD simulations. Employing MD simulations, we can easily control the range and strength of cohesive forces, and directly calculate stress from numerical data. Thus, in contrast to the EP model [33] and other continuum models [56], the advantage of our method is that the effect of cohesive interactions on avalanches can be unambiguously examined. In the following, we explain our numerical model of soft cohesive particles (Section 2.1) and show how the system is prepared and applied simple shear deformations (Section 2.2).
2.1 Cohesive contact model
In ordinary MD simulations of soft frictionless particles [57], a contact force between the particles, and , is modeled as . Here, with the inter-particle distance, , is the unit vector parallel to the normal direction, where is the position of the particle . In the contact force, represents viscous interaction, where is the viscosity coefficient and is the relative velocity between the particles in contact. The magnitude of elastic force is given by a linear spring as , where is the stiffness and with the sum of particle radii, , represents an overlap between the particles.
In contrast, soft cohesive particles are modeled by coating the soft frictionless particles with sticky layers [48–55]. It is assumed that every particle is covered by a thin sticky layer with the thickness , where is introduced as a small dimensionless parameter. The interaction between the particles is attractive if two layers are overlapped. In our MD simulations, we model the contact force as , where the magnitude of cohesive force scaled by is given by
Here, is the scaled overlap [58]. As shown in Figure 1, the interaction is attractive, i.e., , if the scaled overlap is in the range, , whereas it is repulsive, i.e., , if the scaled overlap is positive, . Note that a potential energy for the cohesive interaction, Eq. (1), is a continuous and smooth function of the inter-particle distance (see SM). In addition, the viscous force acts on the particle if , where either two layers or particles are overlapped.
FIGURE 1
2.2 Simple shear deformations
We prepare our system as a binary mixture of soft cohesive particles, where two kinds of particles have the same mass and different diameters, and [59]. We randomly distribute the particles in a square periodic box such that packing fraction of the particles is given by . Note that the area of sticky layers is not included in .
To apply simple shear deformations to the system, we employ the Lees-Edwards boundary condition. In each time step, we replace every particle position, , with and numerically integrate equations of motion, . We use for the strain increment and for the time increment, where is a unit of time [60]. The shear rate is defined as which we fix throughout the simulation.
In the following, we analyze the system in a steady state, where the amount of shear strain exceeds unity, i.e., . In a steady state, the energy injected by shear is dissipated by the viscous forces between the particles, . Therefore, the injection of energy and the energy dissipation are balanced such that any observables, e.g., shear stress, are steady. We vary the number of particles (the system size) and the packing fraction in the ranges, and , respectively. Furthermore, we scale every time and length by the units, and , respectively.
3 Results
In this section, we show our numerical results of soft cohesive particles under shear. First, we examine how the cohesive forces alter force-chain networks (Section 3.1) and affect macroscopic mechanical responses (Section 3.2). Second, we analyze the effect of cohesive interactions on time-averaged stress (Section 3.3). Third, we introduce slip avalanches and examine their dependence on the cohesive interactions (Section 3.4). Then, we study how scaling laws (Section 3.5) and statistics of avalanches (Section 3.6) are changed by the cohesive forces. In addition, we show that localized non-affine displacements are characteristic of avalanches in soft cohesive particles (Section 3.7).
3.1 Force-chain networks
The structure of force-chain networks of soft cohesive particles under shear is strongly influenced by the range of cohesive interactions. Figure 2A displays snapshots of force-chain networks, where the systems are sheared (as indicated by the horizontal arrows in the top panel) and have reached steady states. In this figure, a small system size, , is used for visualization, where the dimensionless parameter for the cohesive forces is given by (top) and (bottom). The packing fraction of the particles is which is much smaller than the jamming transition density, , for soft frictionless particles in two dimensions [59]. If the thickness of sticky layers is sufficiently small (Figure 2A top), the particles are homogeneously distributed and the force-chain network (solid lines) is weak and homogeneous. However, increasing , we observe that the particles aggregate with each other—as the “gelation” [61, 62] — and strong force-chains are transmitted through the system (Figure 2A bottom). In other words, if the range of cohesive interactions is sufficient, the system is partially jammed and jammed regions percolate through the system. Figures 2B, C show non-affine displacements of the particles and changes of force-chain networks, respectively, which will be discussed in later (Section 3.7).
FIGURE 2
3.2 Macroscopic mechanical responses
We quantify mechanical responses of soft cohesive particles to simple shear deformations by shear stress. We calculate stress tensor of the system according to the Born-Huang expression [63],
Here, is the magnitude of contact force between the particles, and , and each component of the normal unit vector is written as . The shear stress is defined as the average of off-diagonal elements, . Figure 3A displays stress-strain curves, i.e., vs , where the shear stress is scaled by the stiffness . Increasing the strain from zero, we observe that increases from zero and becomes steady when exceeds unity. The steady state shear stress tends to be large if the thickness of sticky layer or increases (as indicated by the vertical arrow in Figure 3A). In SM, we show our results of granular temperature, pressure, and macroscopic friction coefficient [64–68] as functions of . We find that, though the granular temperature and pressure exhibit similar dependence on the layer thickness , the macroscopic friction coefficient defined as is less sensitive to .
FIGURE 3
The influence of cohesive force on the shear stress is intuitively understood as follows. As shown in Figure 3B, we consider typical configurations of the particles under shear. In general, the system under simple shear deformations is compressed along degrees diagonal and decompressed along 45° diagonal [69]. If the particles, and , are aligned in the compressive direction (Figure 3B left), the particles tend to be over-compressed such that , where the force is repulsive, i.e., . In the compressive direction, the product of - and -components of the normal unit vector is negative, , so that such a configuration increases (see Equation 2). On the other hand, if the particles are aligned in the decompressive direction (Figure 3B right), the particles tend to be connected by the sticky layer, where the force is attractive, i.e., . Since in the decompressive direction, such a configuration also increases . Therefore, both cases contribute to the increase of so that the stronger the cohesive interactions are, the larger the shear stress is.
3.3 Mean shear stress and stress fluctuations
To quantify the influence of cohesive forces on the shear stress , we analyze the mean value of in a steady state. As shown in the stress-strain curve (Figure 3A), fluctuates around its mean value if the shear strain exceeds unity, . Thus, we calculate the mean shear stress as the average of in the strain interval, . Though previous works of the rheology of soft frictionless particles focused on the divergence of viscosity near the jamming transition [70–74], we show the dependence of on the range of cohesive interactions and packing fraction of the particles . Figure 4A plots as a function of , where increases as listed in the legend. It is obvious that increases with because the system becomes rigid with the increase of density [2]. To describe the dependence of on , we note that the magnitude of cohesive force is divided into two parts as , where and do not include (see Equation 1). According to the Born-Huang expression, Equation 2, we find that the shear stress is decomposed aswhere and are independent of (see SM). Our numerical results of are well described by (dashed lines in Figure 4A) if we use the coefficients, and , for fitting parameters.
FIGURE 4
In addition to the mean shear stress , we analyze fluctuations of the shear stress in a steady state. We quantify stress fluctuations by the variance, , where the ensemble averages are taken over the data in the strain interval, . Figure 4B shows the variance as a function of , where the packing fraction increases as in Figure 4A. As can be seen, the variance is also a monotonically increasing function of and . To explain the dependence of on , we substitute into . Neglecting the correlation between and as , we find that the variance is quadratic in as (see SM). Our numerical results are well described by (dashed lines in Figure 4B), where and are used for fitting parameters.
3.4 Slip avalanches
In contrast to the mean shear stress and stress fluctuations (Section 3.3), slip avalanches characterize plastic responses of the system to simple shear deformations. Closely looking at the stress-strain curve in a steady state (Figure 3A), one observes that the shear stress increasing with the shear strain suddenly drops to a lower value. Such a stress drop event, or slip avalanche, makes the mean shear stress (in a steady state) constant such that the shear stress fluctuates around it. A stress drop amplitude for a single slip avalanche is introduced aswhere the shear stress starts decreasing at the strain and stops decreasing at . The duration of a slip avalanche is given by dimensionless avalanche duration and the so-called avalanche size is defined as [33].
We calculate an average of avalanche sizes as when the system is in a steady state . While previous studies paid much attention to finite size effects on the mean avalanche size [29, 35, 75, 76, 77, 78, 79, 80], we focus on the dependence of on the range of cohesive interactions and packing fraction of the particles . Figure 4C displays double logarithmic plots of and , where increases as listed in the legend of Figure 4A. In this figure, we average over stress drop events in a steady state. Similar to the mean shear stress and stress fluctuations (Figures 4A, B), is also a monotonically increasing function of and . Substituting the decomposition, Equation 3, into Equation 4, we find that the avalanche size is decomposed as , where and are independent of (see SM). Our numerical results are well fitted to (dashed lines in Figure 4C) if we use and for fitting parameters.
The slip avalanche defined as Equation 4 can be rephrased as the stress drop rate is negative, i.e., , in the strain interval between and . Because the dimensionless avalanche duration is finite, , we can calculate the maximum stress drop rate as for each stress drop event. Figure 4D shows double logarithmic plots of the mean maximum stress drop rate and , where we took ensemble averages of in a steady state and increases as listed in the legend of Figure 4A. Similar to the mean avalanche size (Figure 4C), also monotonously increases with and . Because of Equation 3, the maximum stress drop rate is decomposed as with the -independent coefficients, and (see SM). We find that our numerical results are well fitted to (dashed lines in Figure 4D) if we use and for fitting parameters.
Note that we cannot see a clear trend in the mean value of dimensionless avalanche duration , where slightly increases with (see SM).
3.5 Scaling laws of slip avalanches
It was predicted by the MF theory of slip avalanches that both and scale as the square root of , i.e., and , in the scaling regime [19]. To examine whether the MF predictions are applicable to our system, we make scatter plots of , , and from the data of slip avalanches in a steady state. Figure 5A shows the scatter plots of and (left), and and (right), where each dot corresponds to each slip avalanche. The symbols (circles) are the averages of and in each bin of . As can be seen, both and tend to increase with if the avalanche size is large enough, e.g., or . In our MD simulations, the lower bound of dimensionless avalanche duration is given by the strain increment as . Therefore, if the dimensionless avalanche duration is the smallest, , the maximum stress drop rate is given by , where is proportional to . The solid line in Figure 5A (right) indicates the proportionality, , for sufficiently small avalanches.
FIGURE 5
We show that the scaling laws of , , and can be confirmed even if the cohesive forces exist. Figure 5B displays the averages of (left) and (right) in each bin of , where we increase the dimensionless parameter as listed in the legend (indicated by the arrows). One can see that both and monotonically increase with if the avalanche size is large enough as . To examine the scaling laws of , , and , we plot and as functions of in Figure 5C. We find that all the data in are nicely collapsed if the dimensionless parameter is not infinitesimal, . The solid lines represent scaling laws.where the exponents, and , estimated from the data of are different from the MF prediction, . Therefore, the scaling laws, Eqs. (5) and (6), can be confirmed for , where the exponents, and , are quite insensitive to the range of cohesive interactions. Note that the influence of cohesive forces is included in the averages, , , and (see Figure 4).
In SM, we show that Equations 5, 6 hold for large avalanches, , regardless of the packing fraction , where we estimate the exponents, and , for each value of . We also examine finite size effects on the scaling, Equations 5, 6, and confirm that the scaling laws are well established if the system size is large enough (see SM).
3.6 Statistics of slip avalanches
In contrast to the mean values, and , the occurrence of slip avalanches is quantified by the probability distribution functions (PDFs) of avalanche sizes and dimensionless avalanche duration, i.e., and . The MF theory predicts the power-law decay of the PDFs as and in the scaling regimes, where the power-law exponents are given by and [12–15]. The power-law decay was validated by many experiments of granular materials [16–19], metallic glasses [22–25], etc. However, some simulations disagree with the MF predictions and there has been much debate as mentioned in the Introduction (Section 1). Figure 6 shows double logarithmic plots of the PDFs, (a) and (b) , where the dimensionless parameter increases as listed in the legend of (a) (indicated by the arrows). Increasing , we observe qualitative changes in the shapes of the PDFs, e.g., their tails extend to higher values of and . In addition, our numerical results are clearly different from the MF predictions, and (dashed lines).
FIGURE 6
To analyze the shapes of and , we scale them by the mean values, and , respectively. Figure 7A displays the scaled PDFs of avalanche sizes. If the range of cohesive interactions is not sufficiently small , the scaled PDFs are well collapsed on top of each other. We describe the tails of the scaled PDFs as(solid line), where the cut-off value [32], , is extracted from the simulation with . The power-law exponent in Equation 7 is estimated from the data of as which is smaller than the MF prediction, . Figure 7B shows the scaled PDFs of dimensionless avalanche duration, where all the data are nicely collapsed if . The solid line indicateswith , where we find an extremely small exponent, , from the data of . Thus, the power-law scaling of the PDFs, Equations 7, 8, can be confirmed for , where the exponents, and , are almost independent of the range of cohesive interactions. The influence of cohesive forces is included in the averages, and , and cutoff values, and (data are not shown).
FIGURE 7
Though the power-law exponents, and , in Equations 7, 8 are smaller than the MF predictions, they are consistent with the exponent for the scaling law, Equation 5. Let us derive a relation between the power-law exponents as follows. The probability that the avalanche size is in the range between and is given by . On the other hand, the probability that the dimensionless avalanche duration is in the range between and is . Assuming the one-to-one correspondence between and , we equate these probabilities as . This means that the PDFs are related to each other asBecause the dimensionless avalanche duration obeys the scaling law, Equation 5, in the scaling regime, the derivative in Equation 9 scales as . Thus, substituting the power-law decay, and , into Equation 9, we find . Therefore, the exponent is related to the others asThe MF predictions, i.e., , , and , satisfy Equation 10, while our results, , , and , are also consistent with (Equation 10)2.
In SM, we analyze the effect of packing fraction on the scaled PDFs, where we estimate the exponents, and , and examine the relation, Equation 10, for each value of .
In addition, we examine finite size effects on the tails of the scaled PDFs, where Equations 7, 8 well explain our numerical results unless the system size is extremely small (see SM).
3.7 Particle rearrangements during a slip avalanche
On a microscopic scale, a slip avalanche is triggered by rearrangements of the particles under shear. In our MD simulations, particle displacements (in each strain step) can be decomposed as , where is a unit vector parallel to the -axis. The first term represents an affine displacement, whereas the second term is a non-affine displacement. The non-affine displacements represent the particle rearrangements under shear [81–83] and suppress the increase of shear stress [84]. It is thus expected that the non-affine displacements are relevant to the avalanche size . To analyze the particle rearrangements during a single slip avalanche, we integrate the non-affine displacements over the strain interval between and as . Figure 2B shows spatial distributions of (arrows), where the gray scale indicates the magnitude of . If the range of cohesive interactions is small , the spatial distributions of are homogeneous (Figure 2B top), where we observe “collective rearrangements” of the particles everywhere in the system [81–83]. In contrast, if the range of cohesive interactions is large , the spatial distributions are typically localized (Figure 2B bottom), where characteristic “quadrupole structures” of can be seen in the green circles [28, 29, 33].
Note that particle rearrangements are directly linked to restructuring of force-chain networks [85, 86]. Because we calculate the stress tensor by the Born-Huang expression (Eq. (2)), a stress drop event is a consequence of restructuring of force-chains. Figure 2C visualizes the changes of force-chain networks during a slip avalanche. The red (blue) solid lines represent the increase (decrease) of the repulsive forces, , where the particles, and , are in contact before the slip avalanche. If cohesive interactions are weak (top), the restructuring of force-chains is not localized. However, if cohesive interactions are strong (bottom), the restructuring is localized at which the cumulative non-affine displacements are localized (green circles).
To quantitatively compare the non-affine displacements with the avalanche size , we introduce the mean squared displacement (MSD) during a slip avalanche [87] as
Figure 8A displays double logarithmic plots of the average of MSD and the dimensionless parameter , where we averaged Equation 11 over slip avalanches in a steady state. Similar to the mean values in Figure 4, monotonously increases with and the packing fraction (symbols). We also quantify the localization of particle rearrangements by the participation ratio,
FIGURE 8
Figure 8B shows the mean participation ratio as a function of , where we took ensemble averages of Equation 12 in a steady state. Different from the average of MSD (and other averages in Figure 4), decreases with the increases of and (symbols). If (see Section 3.1 for the jamming transition density ), the mean participation ratio is extremely small, , such that the rearrangements during a slip avalanche are localized regardless of cohesive forces. However, if and the range of cohesive interactions is small, the rearrangements are not localized, e.g., for and . Interestingly, starts decreasing if exceeds and drops to if . Therefore, if the thickness of sticky layer increases to 1% of particle diameter, the rearrangements during a slip avalanche are mostly localized regardless of .
Finally, we show that the MSD is relevant to the avalanche size , while the localization of non-affine displacements is unrelated to . Figure 9A displays double logarithmic plots of the scaled MSD and scaled avalanche size , where each symbol represents an average of in each bin of . As in the cases of dimensionless avalanche duration and maximum stress drop rate (Figure 5), all the data of in are well collapsed if the range of cohesive interactions is not infinitesimal, . From our numerical results, we find the scaling law of MSD (solid line) aswhere the exponent estimated from the data of is . The positive exponent, , means that a large avalanche, , is accompanied by a large amount of particle rearrangements. In addition, Equation 13 may be explained by the scaling law of (Equation 5) as follows. If the magnitude of non-affine displacement in each strain step is given by (on average), the magnitude of non-affine displacement during a slip avalanche is . The MSD is then estimated as . Substituting Equation 5, we find , where the exponent, , is close to . In contrast, the participation ratio is flat over the whole range of if (Figure 9B)3. If , is constant unless the avalanche size is extremely large and slightly increases with extremely large avalanches (Figure 9B). However, is much more sensitive to and its correlation with is not significant.
FIGURE 9
In SM, we confirm that the scaling law, Equation 13, holds in regardless of the packing fraction . We also show that, if , the cohesive interactions are not important and the non-affine displacements during a slip avalanche are mostly localized except for extremely large avalanches, . Furthermore, we examine finite size effects on the scaling law and participation ratio and find that Equation 13 is well established unless the system size is too small (see SM).
4 Discussion and conclusion
In this study, we have examined mechanical responses of soft cohesive particles to simple shear deformations by MD simulations. In our cohesive contact model [48–50, 53, 54], the range of cohesive interactions is controlled by the dimensionless parameter (Equation 1). We found that, if is large enough, the particles are locally jammed or aggregate each other as if the system exhibits “gelation” (Figure 2A). The shear stress in a steady state, i.e., the shear strength, and the averages of avalanche size and maximum stress drop rate are increasing functions of (and the packing fraction of the particles ) (Figure 4). Since the shear stress is given by a linear function of (Equation 3), the avalanche size and maximum stress drop rate are also linear functions of . We showed that the scaling laws of dimensionless avalanche duration and maximum stress drop rate (Equations 5, 6) are well established even if the cohesive forces exist (Figure 5). In addition, another scaling law of MSD (Equation 13) was confirmed as the consequence of the scaling law of dimensionless avalanche duration. We also found that the PDFs of avalanche sizes and dimensionless avalanche duration are well described by the power-law scaling (Equations 7, 8). However, the power-law exponents extracted from the data of PDFs do not agree with the MF predictions but rather are consistent with the scaling law of dimensionless avalanche duration (Equation 10). Note that all the exponents for the scaling laws and the PDFs are independent of the range of cohesive interactions (if ), while the influence of cohesive forces is included in the averages and cutoff values.
One of the characteristic features of soft cohesive particles under shear is the localization of particle rearrangements. We quantified the localization by the participation ratio of non-affine displacements and found that, if the system is less dense as , the particle rearrangements during a slip avalanche are strongly dependent on the range of cohesive interactions (Figure 8B). Interestingly, the participation ratio does not correlate significantly with the avalanche size (Figure 9B) and the non-affine displacements, if localized, exhibit characteristic quadrupole structures (Figure 2B). A further analysis is necessary to unveil the mechanism of the localization, which is left for future work.
In our MD simulations, we assumed that the system is in a pendular state, where the cohesive force between the particles is pairwise. However, if the liquid content increases, the system transitions to a funicular or capillary state, where more than two particles interact through the liquid [3,4]. We did not implement such many body interactions into our model though their effects on avalanches are interesting to know. In addition, in real granular materials, cohesive forces are intrinsically history-dependent [57]. Therefore, the influence of hysteresis in cohesive contacts has to be examined in future. Moreover, the effect of particle shapes [47] and simulations in three dimensions are important for practical applications of this work.
In conclusion, the shear strength, force-chains, and particle rearrangements are strongly affected by cohesive forces if the system is less dense. The statistics of avalanches, such as the scaling laws and power-law distributions, are well established even if the system is cohesive though the scaling exponents are distinct from the MF predictions.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
KS: Conceptualization, Funding acquisition, Investigation, Visualization, Writing–original draft, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was financially supported by the KAKENHI Grant No. 22K03459 from JSPS.
Conflict of interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2025.1548966/full#supplementary-material
Footnotes
1.^The suction is defined as the pressure difference between a liquid bridge and air.
2.^Substituting our estimates, and , into Equation 10, we find which is approximately equal to our result, .
3.^We confirmed that over the whole range of even if and (data are not shown).
References
1.
JeagerHNagelSBehringerR. Granular solids, liquids, and gases. Rev Mod Phys (1996) 68:1259–73. 10.1103/revmodphys.68.1259
2.
BonnDDennMMBerthierLDivouxTMannevilleS. Yield stress materials in soft condensed matter. Rev Mod Phys (2017) 89:035005. 10.1103/revmodphys.89.035005
3.
MitaraiNNoriF. Wet granular materials. Adv Phys (2006) 55:1–45. 10.1080/00018730600626065
4.
MitaraiNNakanishiH. Simple model for wet granular materials with liquid clusters. Eur Phys Lett (2009) 88:64001. 10.1209/0295-5075/88/64001
5.
CastellanosA. The relationship between attractive interparticle forces and bulk behaviour in dry and uncharged fine powders. Adv Phys (2005) 54:263–376. 10.1080/17461390500402657
6.
UhlJTPathakSSchorlemmerDLiuXSwindemanRBrinkmanBAWet alUniversal quake statistics: from compressed nanocrystals to earthquakes. Sci Rep (2015) 5:16493. 10.1038/srep16493
7.
HatanoTNarteauCShebalinP. Common dependence on stress for the statistics of granular avalanches and earthquakes. Sci Rep (2015) 5:12280. 10.1038/srep12280
8.
LherminierSPlanetRVehelVLSimonGVanelLMåløyKJet alContinuously sheared granular matter reproduces in detail seismicity laws. Phys Rev Lett (2019) 122:218501. 10.1103/physrevlett.122.218501
9.
BakPTangCWiesenfeldK. Self-organized criticality: an explanation of the 1/foise. Phys Rev Lett (1987) 59:381–4. 10.1103/physrevlett.59.381
10.
SethnaJPDahmenKAMyersCR. Crackling noise. Nature (2001) 410:242–50. 10.1038/35065675
11.
FisherDS. Collective transport in random media: from superconductors to earthquakes. Phys Rep (1998) 301:113–50. 10.1016/s0370-1573(98)00008-8
12.
DahmenKABen-ZionYUhlJT. A simple analytic theory for the statistics of avalanches in sheared granular materials. Nat Phys (2011) 7:554–7. 10.1038/nphys1957
13.
FisherDSDahmenKRamanathanSBen-ZionY. Statistics of earthquakes in simple models of heterogeneous faults. Phys Rev Lett (1997) 78:4885–8. 10.1103/physrevlett.78.4885
14.
DahmenKABen-ZionYUhlJT. Micromechanical model for deformation in solids with universal predictions for stress-strain curves and slip avalanches. Phys Rev Lett (2009) 102:175501. 10.1103/physrevlett.102.175501
15.
DahmenKAErtaşDBen-ZionY. Gutenberg-Richter and characteristic earthquake behavior in simple mean-field models of heterogeneous faults. Phys Rev E (1998) 58:1494–501. 10.1103/physreve.58.1494
16.
GellerDAEckeREDahmenKABackhausS. Stick-slip behavior in a continuum-granular experiment. Phys Rev E (2015) 92:060201. 10.1103/physreve.92.060201
17.
DenisovDVLörinczKAUhlJTDahmenKASchallP. Universality of slip avalanches in flowing granular matter. Nat Commun (2016) 7:10641. 10.1038/ncomms10641
18.
DenisovDVLörinczKAWrightWJHufnagelTCNawanoAGuXet alUniversal slip dynamics in metallic glasses and granular matter – linking frictional weakening with inertial effects. Sci Rep (2017) 7:43376. 10.1038/srep43376
19.
LongAADenisovDVSchallPHufnagelTCGuXWrightWJet alFrom critical behavior to catastrophic runaways: comparing sheared granular materials with bulk metallic glasses. Granular Matter (2019) 21:99. 10.1007/s10035-019-0946-y
20.
FriedmanNJenningsATTsekenisGKimJYTaoMUhlJTet alStatistics of dislocation slip avalanches in nanosized single crystals show tuned critical behavior predicted by a simple mean field model. Phys Rev Lett (2012) 109:095507. 10.1103/physrevlett.109.095507
21.
CsikorFFMotzCWeygandDZaiserMZapperiS. Dislocation avalanches, strain bursts, and the problem of plastic forming at the micrometer scale. Science (2007) 318:251–4. 10.1126/science.1143719
22.
AntonagliaJWrightWJGuXByerRRHufnagelTCLeBlancMet alBulk metallic glasses deform via slip avalanches. Phys Rev Lett (2014) 112:155501. 10.1103/physrevlett.112.155501
23.
AntonagliaJXieXSchwarzGWraithMQiaoJZhangYet alTuned critical avalanche scaling in bulk metallic glasses. Sci Rep (2014) 4:4382. 10.1038/srep04382
24.
ColemanJPMengFTsuchiyaKBeadsworthJLeBlancMLiawPKet alEffect of annealing on nanoindentation slips in a bulk metallic glass. Phys Rev B (2017) 96:134117. 10.1103/physrevb.96.134117
25.
SunBAYuHBJiaoWBaiHYZhaoDQWangWH. Plasticity of ductile metallic glasses: a self-organized critical state. Phys Rev Lett (2010) 105:035501. 10.1103/physrevlett.105.035501
26.
SheikhMAWeaverRLDahmenKA. Avalanche statistics identify intrinsic stellar processes near criticality in KIC 8462852. Phys Rev Lett (2016) 117:261101. 10.1103/physrevlett.117.261101
27.
TewariSSchiemannDDurianDJKnoblerCMLangerSALiuAJ. Statistics of shear-induced rearrangements in a two-dimensional model foam. Phys Rev E (1999) 60:4385–96. 10.1103/physreve.60.4385
28.
MaloneyCLemaitreA. Subextensive scaling in the athermal, quasistatic limit of amorphous matter in plastic shear flow. Phys Rev Lett (2004) 93:016001. 10.1103/physrevlett.93.016001
29.
MaloneyCELemaitreA. Amorphous systems in athermal, quasistatic shear. Phys Rev E (2006) 74:016118. 10.1103/physreve.74.016118
30.
BaileyNPSchiøtzJLemaitreAJacobsenKW. Avalanche size scaling in sheared three-dimensional amorphous solid. Phys Rev Lett (2007) 98:095501. 10.1103/physrevlett.98.095501
31.
HeussingerCChaudhuriPBarratJL. Fluctuations and correlations during the shear flow of elastic particles near the jamming transition. Soft Matter (2010) 6:3050. 10.1039/b927228c
32.
ShangBGuanPBarratJL. Elastic avalanches reveal marginal behavior in amorphous solids. Proc Natl Acad Sci USA (2020) 117:86–92. 10.1073/pnas.1915070117
33.
NicolasAFerreroEEMartensKBarratJL. Deformation and flow of amorphous solids: insights from elastoplastic models. Rev Mod Phys (2018) 90:045006. 10.1103/revmodphys.90.045006
34.
LinJLernerERossoAWyartM. Scaling description of the yielding transition in soft amorphous solids at zero temperature. Proc Natl Acad Sci USA (2014) 111:14382–7. 10.1073/pnas.1406391111
35.
DubeyAKHentschelHGEProcacciaISinghM. Statistics of plastic events in post-yield strain-controlled amorphous solids. Phys Rev B (2016) 93:224204. 10.1103/physrevb.93.224204
36.
WhiteRADahmenKA. Driving rate effects on crackling noise. Phys Rev Lett (2003) 91:085702. 10.1103/physrevlett.91.085702
37.
BehringerRPBiDChakrabortyBHenkesSHartleyRR. Phys Rev Lett (2008) 101:268301. 10.1103/physrevlett.101.268301
38.
LootensDVanDammeHHébraudP. Giant stress fluctuations at the jamming transition. Phys Rev Lett (2003) 90:178301. 10.1103/physrevlett.90.178301
39.
LiuCFerreroEEPuosiFBarratJLMartensK. Driving rate dependence of avalanche statistics and shapes at the yielding transition. Phys Rev Lett (2016) 116:065501. 10.1103/physrevlett.116.065501
40.
KhfifiMLoulidiM. Scaling properties of a rice-pile model: inertia and friction effects. Phys Rev E (2008) 78:051117. 10.1103/physreve.78.051117
41.
SalernoKMMaloneyCERobbinsMO. Phys Rev Lett (2012) 109:105703. 10.1103/physrevlett.109.105703
42.
SalernoKMRobbinsMO. Effect of inertia on sheared disordered solids: critical scaling of avalanches in two and three dimensions. Phys Rev E (2013) 88:062206. 10.1103/physreve.88.062206
43.
KarimiKFerreroEEBarratJL. Inertia and universality of avalanche statistics: the case of slowly deformed amorphous solids. Phys Rev E (2017) 95:013003. 10.1103/physreve.95.013003
44.
BarésJWangDWangDBertrandTO’HernCSBehringerRP. Local and global avalanches in a two-dimensional sheared granular medium. Phys Rev E (2017) 96:052902. 10.1103/physreve.96.052902
45.
SaitohK. The role of friction in statistics and scaling laws of avalanches. Eur Phys J E (2021) 44:85. 10.1140/epje/s10189-021-00089-8
46.
PeñaAAMcNamaraSLindPGHerrmannHJ. Avalanches in anisotropic sheared granular media. Granular Matter (2009) 11:243–52. 10.1007/s10035-009-0136-4
47.
MurphyKADahmenKAJaegerHM. Transforming mesoscale granular plasticity through particle shape. Phys Rev X (2019) 9:011014. 10.1103/physrevx.9.011014
48.
IraniEChaudhuriPHeussingerC. Impact of attractive interactions on the rheology of dense athermal particles. Phys Rev Lett (2014) 112:188303. 10.1103/physrevlett.112.188303
49.
IraniEChaudhuriPHeussingerC. Athermal rheology of weakly attractive soft particles. Phys Rev E (2016) 94:052608. 10.1103/physreve.94.052608
50.
IraniEChaudhuriPHeussingerC. Discontinuous shear-thinning in adhesive dispersions. Phys Rev Fluids (2019) 4:074307. 10.1103/physrevfluids.4.074307
51.
ChaudhuriPBerthierLBocquetL. Inhomogeneous shear flows in soft jammed materials with tunable attractive forces. Phys Rev E (2012) 85:021503. 10.1103/physreve.85.021503
52.
ZhengWLiuHXuN. Shear-induced solidification of athermal systems with weak attraction. Phys Rev E (2016) 94:062608. 10.1103/physreve.94.062608
53.
KoezeDJTigheBP. Sticky matters: jamming and rigid cluster statistics with attractive particle interactions. Phys Rev Lett (2018) 121:188002. 10.1103/physrevlett.121.188002
54.
KoezeDJHongLKumarATigheBP. Elasticity of jammed packings of sticky disks. Phys Rev Res (2020) 2:032047. 10.1103/physrevresearch.2.032047
55.
LoisGBlawzdziewiczJO’HernCS. Jamming transition and new percolation universality classes in particulate systems with attraction. Phys Rev Lett (2008) 100:028001. 10.1103/physrevlett.100.028001
56.
Herrera-ValenciaEESánchez-VillavicencioMLSoriano-CorreaCBautistaORamírez-TorresLAHernández-AbadVJet alStudy of the electroosmotic flow of a structured fluid with a new generalized rheological model. Rheol Acta (2024) 63:3–32. 10.1007/s00397-023-01418-8
57.
LudingS. Anisotropy in cohesive, frictional granular media. J Phys Condens Matter (2005) 17:S2623–40. 10.1088/0953-8984/17/24/017
58.
van HeckeM. Jamming of soft particles: geometry, mechanics, scaling and isostaticity. J Phys Condens Matter (2010) 22:033101. 10.1088/0953-8984/22/3/033101
59.
O’HernCSSilbertLELiuAJNagelSR. Jamming at zero temperature and zero applied stress: the epitome of disorder. Phys Rev E (2003) 68:011306. 10.1103/physreve.68.011306
60.
RadjaiFRouxS. Turbulentlike fluctuations in quasistatic flow of granular media. Phys Rev Lett (2002) 89:064302. 10.1103/physrevlett.89.064302
61.
MizunoHHachiyaMIkedaA. Structural, mechanical, and vibrational properties of particulate physical gels. J Chem Phys (2021) 155:234502. 10.1063/5.0072863
62.
MizunoHHachiyaMIkedaA. Phonon transport properties of particulate physical gels. J Chem Phys (2022) 156:204505. 10.1063/5.0090233
63.
BornMHuangK. Dynamical theory of crystal lattices. Oxford: Clarendon Press (1988).
64.
MiDiGDR. On dense granular flows. Eur Phys J E (2004) 14:341–65. 10.1140/epje/i2003-10153-0
65.
da CruzFEmamSProchnowMRouxJNChevoirF. Rheophysics of dense granular materials: discrete simulation of plane shear flows. Phys Rev E (2005) 72:021309. 10.1103/physreve.72.021309
66.
HatanoT. Power-law friction in closely packed granular materials. Phys Rev E (2007) 75:060301. 10.1103/physreve.75.060301
67.
AzémaERadjaïF. Internal structure of inertial granular flows. Phys Rev Lett (2014) 112:078001. 10.1103/PhysRevLett.112.078001
68.
RognonPGRouxJNNaaïmMChevoirF. Dense flows of cohesive granular materials. J Fluid Mech (2008) 596:21–47. 10.1017/s0022112007009329
69.
SinghAMagnanimoVSaitohKLudingS. Effect of cohesion on shear banding in quasistatic granular materials. Phys Rev E (2014) 90:022202. 10.1103/physreve.90.022202
70.
OlssonPTeitelS. Herschel-bulkley shearing rheology near the athermal jamming transition. Phys Rev Lett (2012) 109:108001. 10.1103/physrevlett.109.108001
71.
AndreottiBBarratJLHeussingerC. Shear flow of non-brownian suspensions close to jamming. Phys Rev Lett (2012) 109:105901. 10.1103/physrevlett.109.105901
72.
VågbergDOlssonPTeitelS. Universality of jamming criticality in overdamped shear-driven frictionless disks. Phys Rev Lett (2014) 113:148002. 10.1103/physrevlett.113.148002
73.
KawasakiTCoslovichDIkedaABerthierL. Diverging viscosity and soft granular rheology in non-Brownian suspensions. Phys Rev E (2015) 91:012203. 10.1103/physreve.91.012203
74.
OlssonP. Dimensionality and viscosity exponent in shear-driven jamming. Phys Rev Lett (2019) 122:108003. 10.1103/physrevlett.122.108003
75.
LernerEProcacciaI. Locality and nonlocality in elastoplastic responses of amorphous solids. Phys Rev E (2009) 79:066109. 10.1103/physreve.79.066109
76.
KarmakarSLernerEProcacciaI. Statistical physics of the yielding transition in amorphous solids. Phys Rev E (2010) 82:055103. 10.1103/physreve.82.055103
77.
HentschelHGEKarmakarSLernerEProcacciaI. Size of plastic events in strained amorphous solids at finite temperatures. Phys Rev Lett (2010) 104:025501. 10.1103/physrevlett.104.025501
78.
KarmakarSLernerEProcacciaIZylbergJ. Statistical physics of elastoplastic steady states in amorphous solids: finite temperatures and strain rates. Phys Rev E (2010) 82:031301. 10.1103/physreve.82.031301
79.
DasguptaRKarmakarSProcacciaI. Universality of the plastic instability in strained amorphous solids. Phys Rev Lett (2012) 108:075701. 10.1103/physrevlett.108.075701
80.
TalamaliMPetäjäVVandembroucqDRouxS. Avalanches, precursors, and finite-size fluctuations in a mesoscopic model of amorphous plasticity. Phys Rev E (2011) 84:016115. 10.1103/physreve.84.016115
81.
SaitohKMizunoH. Anomalous energy cascades in dense granular materials yielding under simple shear deformations. Soft Matter (2016) 12:1360–7. 10.1039/c5sm02760h
82.
SaitohKMizunoH. Exstrophy cascades in two-dimensional dense granular flows. Phys Rev E (2016) 94:022908. 10.1103/physreve.94.022908
83.
SaitohKMizunoH. Anisotropic decay of the energy spectrum in two-dimensional dense granular flows. Phys Rev E (2017) 96:012903. 10.1103/physreve.96.012903
84.
MizunoHSaitohKSilbertLE. Elastic moduli and vibrational modes in jammed particulate packings. Phys Rev E (2016) 93:062905. 10.1103/physreve.93.062905
85.
SaitohKMagnanimoVLudingS. A Master equation for the probability distribution functions of forces in soft particle packings. Soft Matter (2015) 11:1253–8. 10.1039/c4sm02452d
86.
SaitohKOyamaNOgushiFLudingS. Transition rates for slip-avalanches in soft athermal disks under quasi-static simple shear deformations. Soft Matter (2019) 15:3487–92. 10.1039/c8sm01966e
87.
CaoPDahmenKAKushimaAWrightWJParkHSShortMPet alNanomechanics of slip avalanches in amorphous plasticity. J Mech Phys Sol (2018) 114:158–71. 10.1016/j.jmps.2018.02.012
Summary
Keywords
granular materials, avalanche, plasticity, cohesive interaction, molecular dynamics
Citation
Saitoh K (2025) Shear strength, avalanches, and structures of soft cohesive particles under shear. Front. Phys. 13:1548966. doi: 10.3389/fphy.2025.1548966
Received
20 December 2024
Accepted
20 February 2025
Published
10 March 2025
Volume
13 - 2025
Edited by
Francisco Vega Reyes, University of Extremadura, Spain
Reviewed by
Edtson Emilio Herrera Valencia, National Autonomous University of Mexico, Mexico
Alvaro Rodriguez Rivas, University of Seville, Spain
Updates
Copyright
© 2025 Saitoh.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kuniyasu Saitoh, k.saitoh@cc.kyoto-su.ac.jp
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.