yade-users team mailing list archive

[Jan Stránský <question238289@xxxxxxxxxxxxxxxxxxxxx>]

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

Jan Stránský proposed the following answer:
Hello Hicham,

by playing with your script (thanks :-), I have found two problems:


negIds,posIds,axis,crossSectionArea=UX['negIds'],UX['posIds'],UX['axis'],UX['area']
>

returned crossSectionArea = 0.0, which then is not very good for
stress=force/area computation and might be source of NaNs. Set
crossSectionArea = pi*r*r # r = O.bodies[i].shape.radius
or something like this

Another problem is, that there are no interactions in your simulation,
therefore no interaction forces. Either use the intRadius > 1.0 or create
the interaction "by hand" [1]

let us know if it works or not :-)
cheers
Jan

[1] https://yade-
dem.org/doc/yade.utils.html#yade._utils.createInteraction

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.