Question #689025 on Yade changed: https://answers.launchpad.net/yade/+question/689025
Sacha Duverger posted a new comment: I took off the box and printed the clump's final position in the MWE of the first post (cf new MWE). Yade 2018.02b and yadedaily 20200301-3527~67223fc~bionic1 give the exact same final position, so it looks like the trajectories are the same. As for O.interactions.has(0,1), I understood with Jérôme's message that it should be "False" since no contact is created among clump members of a same clump ([*]). If it is "True" there is an interaction in O.interactions that is not supposed to be there, the computation of unbalancedForce() could then be wrong. [*] : https://yade-dem.org/doc/formulation.html#clumps-rigid-aggregates New MWE: ######################## c_id = O.bodies.appendClumped([sphere((0,0,0), 1), sphere((1.5,0,0), 1)])[0] O.engines=[ForceResetter() ,InsertionSortCollider([Bo1_Sphere_Aabb(), Bo1_Box_Aabb()]) ,InteractionLoop( [Ig2_Sphere_Sphere_ScGeom(), Ig2_Box_Sphere_ScGeom6D()], [Ip2_FrictMat_FrictMat_FrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack()]) ,NewtonIntegrator(damping=0, gravity=(0,-10,0)) ] O.dt = 1e-6 O.run(1000, wait=True) print(O.bodies[c_id].state.pos) -- 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 : [email protected] Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
