Question #657061 on Yade changed: https://answers.launchpad.net/yade/+question/657061
Jan Stránský proposed the following answer: Concerning load, just do ### def plotAddData(): f = O.forces.f(3)[1] d = 2*rayon_g load = f/pow(d,2) plot.addData( dspl = O.bodies[3].state.displ()[1], load = load, ) ### concerning strange results, it is because there are no cohesive interactions in your model. You can easily try ## for i in O.interactions: print "Hi, I am interaction",i i.phys.unp = i.phys.penetrationDepth ## with no print You should put that code after O.step(): ## ... O.bodies.append(sp) ... O.engines = [...] O.step() for i in O.interactions: i.phys.unp = i.geom.penetrationDepth print len([i for i in O.interactions]) ## However, in this case, there are only 18 interactions, which is relatively low number for the packing.. I am nit very familiar with CohFrictMat and how to use it, I suggest to open a new question with more explicit name. cheers Jan -- 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