yade-dev team mailing list archive

[chiara modenese <c.modenese@xxxxxxxxx>]

Hello!

As we already said (
http://www.mail-archive.com/yade-users@xxxxxxxxxxxxxxxxxxx/msg02133.html), I
think there is something wrong/weird in the computation of
plasticDissipation() counter (inside Law2_ScGeom_FrictPhys_Basic functor).
Now, I was wondering if is there any other way to compute the plastic
dissipation in the code, I mean looping over all the contacts and over time
since it is incremental. I ask this since the following line is formally
correct

plasticDissipation +=
((1/currentContactPhysics->ks)*(trialForce-shearForce)).dot(shearForce); //
line 87 of ElasticContactLaw.cpp

and so maybe the mistake is due to the OpenMPAccumulator or something else I
ignore atm. Just let me know if you have any idea how to code this in a
different way. Or also any idea of debugging it as it is now would be
appreciated.

cheers! Chiara