yade-users team mailing list archive

[Luc Scholtès <question689486@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

Luc Scholtès proposed the following answer:
A good and pretty easy way: make sure that the average number of bonds
per particle (coordination number) is the same for all samples.

For that, you can adjust the interaction range:

intR=X # (value to define, >=1)
InsertionSortCollider(Bo1_Sphere_Aabb(aabbEnlargeFactor=intR,label='Saabb')
Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=intR,label='SSgeom')

You can then compute the coordination number with something like:

nbCohIntrs=0
for i in O.interactions:
  if isinstance(O.bodies[i.id1].shape,Sphere) and isinstance(O.bodies[i.id2].shape,Sphere):
      if i.phys.isCohesive:
          nbCohIntrs+=1
print "nbCohIntrs=",nbCohIntrs 

then, at the first iteration, you can reset it (if you want):

O.step();
SSgeom.interactionDetectionFactor=-1.
Saabb.aabbEnlargeFactor=-1.

This will ensure that all samples have similar mechanical properties (if
the samples are all constructed using the same procedure).

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