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