← Back to team overview

yade-users team mailing list archive

Re: [Question #689025]: Difference in unbalancedForce() between yade and yadedaily

 

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

Chareyre posted a new comment:
So, this issue seems to be solved indeed.
Thanks for feedback.
B

On Wed, 4 Mar 2020 at 15:47, Sacha Duverger <
question689025@xxxxxxxxxxxxxxxxxxxxx> wrote:

> 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     : yade-users@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-users
> More help   : https://help.launchpad.net/ListHelp
>


-- 
-- 
_______________
Bruno Chareyre
Associate Professor
ENSE³ - Grenoble INP
Lab. 3SR
BP 53
38041 Grenoble cedex 9
Tél : +33 4 56 52 86 21
________________

Email too brief?
Here's why: email charter
<https://marcuselliott.co.uk/wp-content/uploads/2017/04/emailCharter.jpg>

You received this question notification because your team yade-users is
an answer contact for Yade.