Hello! As we already said ( http://www.mail-archive.com/[email protected]/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
_______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
