yade-users team mailing list archive

[Luc Scholtès <question432617@xxxxxxxxxxxxxxxxxxxxx>]

Question #432617 on Yade changed:
https://answers.launchpad.net/yade/+question/432617

    Status: Open => Answered

Luc Scholtès proposed the following answer:
Hi Tina,

If you are using cohFrictMat, I guess that you are also using
cohFrictPhys and Law2_ScGeom6D_CohFrictPhys_CohesionMoment (right?).

In this case, you should be able to loop over interactions at any time
of the simulation and count the cohesive/frictional contacts by checking
if the cohesionBroken attribute is false (cohesive contact) or true
(frictional contact) using lines such as:

nbCohBond=0
for i in O.interactions:
  if not i.isReal : continue # not mandatory
  if i.phys.cohesionBroken==False:
    nbCohBond+=1

You could include these into a function that you would run periodically
during your simulation (e.g. using PyrRunner()) and thus get the
evolution of the number of cohesive contacts. By running it before your
simulation starts, you will then be able to compute the ratio of broken
bonds / initial bonds.

I hope it helps

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.