The cohesion of granular materials reflects their ability to resist external tensile stresses. This property enhances also the shear strength that is used to assess the cohesion of granular soils from shear tests performed at different confining pressures. The microscopic origin of cohesion is the presence of attraction forces between particles. These interactions, to which we refer as cohesive interactions, have different physical and chemical origins. They manifest themselves between two particles by resistance to separation, shear or rolling. The phenomenology of simple interactions can be recovered by integrating the electromagnetic forces acting on very smaller scales. This upscaling procedure is, however, rather complex. Moreover, most properties of granular media emerge from their many-body feature and the dissipative nature of their interactions. For this reason, the relevant scale for modeling granular materials is often that of the particles and the cohesive interactions need to be considered at this scale.

DEM algorithms can easily be extended to account for cohesive interactions which are often simply supplemented to the repulsive elastic and frictional interactions of cohesionless materials. We present here a general framework for the implementation of cohesion in a discrete-element framework. Generally, simple models of interaction laws are privileged in DEM so as to reduce time-consuming operations arising from too many elementary operations. In any case, one has to check the impact of the simplifications on the resulting global behavior of the material.

In the case of a non-frictional contact, the force acting between two particles is normal to the contact plane and can be additively decomposed into three forces:

| (1) |

where f_{n}^{e} is the repulsive force at contact, f_{n}^{d} a viscous force, and f_{n}^{c} the cohesion
force.

(a) (b)

For spherical particles, the contact surface takes the form of a disk of radius a at which the pressure is not uniformly distributed; figure 1. We make the following assumptions about the contact:

- particles are perfectly smooth;
- the behavior is elastic and isotropic;
- the tangential component of the force does not affect the normal component ;
- the contact deflexion is small (δ
_{n}<< a).

Under these assumptions, the contact force is given by the Hertz contact law (Maugis [2000]):

| (2) |

avec = + where E_{1} et E_{2} are the Young modulus of the two particles
and ν_{1} and ν_{2} the poisson ratio and R^{⋆} is the harmonic mean of the particles
radii.

For cohesive granular materials, the effect of cohesion often dominates that of confining stresses. The macroscopic behavior is much more dictated by the cohesive interactions than by the non-linear behavior at the contact. Therefore, a linear approximation of the contact force is relevant:

| (3) |

where k_{n} is the contact normal stiffness.

The mechanical approach presented here in a 2D configuration (Delenne et al. [2002]) is a Lagrangian approach where the kinematics and deformations are described from a reference state. The presence of cohesion limits the motion so that the displacements are more relevant than velocity as kinematic variables.

The kinematics is described in a global reference frame . We introduce
two material points I_{1} and I_{2} belonging respectively to particle 1 (center C_{1} and
radius R_{1}) and particle 2 (center C_{2} and radius R_{2}); figure 2a and 2b. The reference
state is characterized by the coordinates of the centers of particles and their radii.
The point of cohesion between two particles 1 and 2 is located by the angle α_{12}.

At the reference state, the points I_{1} and I_{2} coincide; figure 2a. Subsequently, these
two points are tracked in their movements through the movements and the rotations (θ_{1}
and θ_{2}) of the two particles. The normal δ_{n} and tangential δ_{t} displacements and the
rotation θ around the cohesion point is calculated in the local frame (where
= ∕
and is a unit vector orthogonal to ):

| (4) |

| (5) |

| (6) |

The cohesion is taken into account through local relationships between the
displacements (δ_{n},δ_{t},θ) and the contact force (f_{n},f_{t},M) where M is the torque at
contact. In general, these relationships can be expressed as a cohesion law
ψ:

| (7) |

A local failure criterion governs the transition between the cohesion state and the sliding-contact state. This criterion is expressed through a yield function:

| (8) |

Different failure criteria may be used depending on the nature of the interface. Figure 4a shows the case of a failure criterion only based on in tensile strength. In this case, we consider that no rupture occurs in pure shear or torsion. This extreme case is similar to the case of normal adhesion for which the action of cohesion is only normal to the contact.

In the case of cemented particles, the failure may occur in tension, but also under shearing and bending. Moreover, for particles that are not in contact because of the presence of an interposed thick binder, fracture may occur in compression. In this case it may be necessary to use different failure thresholds for the different degrees of freedom; figure 4b. It may also be necessary to take into account the cross effects due to combined loading such as compression plus shearing or tension plus torque. This can be done by using more complex criteria of elliptical shape, for example; figure 4c.

(b)

(c)

To test the relevance of this approach, the cohesion law and failure criteria were determined experimentally for a 2D model material consisting of aluminum cylinders glued together with an epoxy resin (Delenne et al. [2004]). Tensile, compression and shear forces as well as torques were applied on a pair of two cylinders with a cohesive bond. For this model medium, these experiments show that the compressive failure threshold is very high compared to the threshold in tension so that the threshold in compression may be set to infinity (Figure 5). A paraboloidal yield surface was found to fit the data and used in the numerical simulations:

| (9) |

For κ < 0, there is cohesion bond, otherwise the bond fails and the contact becomes frictional without cohesion. In the case of an irreversible failure (no rejoining), one can use a parameter η which describes the damaged state of the bond (Delenne et al. [2008]). Another possibility is to modify the yield surface, corresponding to transition from elastic behavior to plastic behavior with a hardening parameter.

Luding et al. proposed a model for plastic contact as shown in figure 6. This model takes into account the normal direction but can easily be extended to more complex loading paths:

| (10) |

where k_{1} < k_{2}. Note that for an initial load, the force increases linearly with the
deflexion δ to a maximum value δ_{max}. The maximum value is stored as a memory
variable. The line of slope k_{1} is the limit value of force for a given value of δ. During
unloading, the force declines along the line of slope k_{2} down to 0 for values of δ
which vary from δ_{max} to δ_{0} = δ_{max}. The value δ_{0} represents a
plastic deformation at contact. Upon reloading, the force increases following
the line of slope k_{2} to the maximum value already reached δ_{max} and the
follows the line of slope k_{1}. Unloading for values below f_{0} (defined for δ = δ_{0})
leads to an attraction force. This force may decrease until f_{min} for a value
δ_{min} = δ_{max}∕ along the slope k_{2}. For δ < δ_{min}, the force follows
the slope -k_{c}.

J. Y. Delenne, M. S. El Youssoufi, and J. C. Bénet. Comportement mécanique et rupture de milieux granulaires cohésifs. CRAS mécanique, 330: 475–482, 2002.

Jean-Yves Delenne, Moulay Saïd El Youssoufi, Fabien Cherblanc, and Jean-Claude Benet. Mechanical behaviour and failure of cohesive granular materials. International Journal for Numerical and Analytical Methods in Geomechanics, 28(15):1577–1594, 2004. URL http://dx.doi.org/10.1002/nag.401.

J.Y. Delenne, Y. Haddad, J.C. Bénet, and J. Abecassis. Use of mechanics of cohesive granular media for analysis of hardness and vitreousness of wheat endosperm. Journal of Cereal Science, 47(3):438–444, 2008.

D. Maugis. Contact, Adhesion and Rupture of Elastic Solids. Springer, 2000.