Bruno, I compared today results of old brefcom that used incremental computation of epsN (roughly the same as in ElasticContactLaw) with new version that does that based on original contact point position (see in SpheresContactGeometry docs). Before, I was using timestep .2*PWaveTimeStep, now I got to 1*PWaveTimeStep with only small change in the strain-stress diagram (the simulation didn't explode). Time per iteration is a bit larger (by about 5-10%) due to more complicated calculations involved. Let me know if you would like to try, adapting ElasticContactLaw would be quite easy.
Vaclav _______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
