yade-users team mailing list archive

[Sacha Duverger <question689025@xxxxxxxxxxxxxxxxxxxxx>]

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

Hello,

While I was performing triaxial tests I noticed that the same script gave different values for unbalancedForce() depending on which version of yade I was using: yade 2018.02b or yadedaily 20191105-2658~35899c7~bionic1. It seems that the problem only occurs with clumps.

Here is a MWE: 

##################

from yade import plot
import numpy as np

O.bodies.appendClumped([sphere((0,0,0), 1), sphere((1.5,0,0), 1)])

wallIds = O.bodies.append(aabbWalls(thickness=1e-10))

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()])
		   ,GlobalStiffnessTimeStepper()
		   ,NewtonIntegrator(damping=0, gravity=(0,-10,0))
		   ,PyRunner(command='yade.plot.addData(uF=unbalancedForce())', iterPeriod=1)
		  ]

O.run(1000, wait=True)

print(np.mean(plot.data['uF']))

##################

yade 2018.02b outputs 0.11062286226025234 .
yadedaily 20191105-2658~35899c7~bionic1 outputs 75.5814592275 .

Which version of Yade should I trust?

Thanks in advance

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