yade-users team mailing list archive

[Sacha Duverger <question689025@xxxxxxxxxxxxxxxxxxxxx>]

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.