yade-dev team mailing list archive

Thread
Date

significat timestep improvement with non-incremental shear computation

To: Yade Development Group <yade-dev@xxxxxxxxxxxxxxxxxxx>
From: Václav Šmilauer <eudoxos@xxxxxxxx>
Date: Fri, 19 Sep 2008 12:48:17 +0200
User-agent: Mozilla/5.0 (X11; U; Linux x86_64; en-US; rv:1.8.1.16) Gecko/20080724 Lightning/0.8 Thunderbird/2.0.0.16 Mnenhy/0.7.5.666

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