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.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
