Re: [Yade-dev] Contact damping

2010-06-20 Thread Bruno Chareyre
Hi Chiara, I have a better suggestion that would not need a new class. You could add a variant of ScGeom::rotateAndGetShear() that would return incident velocity instead of shear increment. It will not be a big deal to compute normal and shear increments in your constitutive law if you have

Re: [Yade-dev] static_cast in GSTS

2010-06-20 Thread chiara modenese
On 18 May 2010 16:55, Bruno Chareyre bruno.chare...@hmg.inpg.fr wrote: 2) Why do we need to continue if fn==0? IIRC Bruno told me that was to avoid contacts that transmit no force, even if they have non-zero nominal stiffness. Exact zero will be there for interactions that really had no

Re: [Yade-dev] Contact damping

2010-06-20 Thread chiara modenese
Hi guys, just another quick question concerning damping before I commit my changes. ATM VisoelasticPM requires the two constants cn and cs to work out the viscous components. However, in HM I preferred to code the possibility to enter the beta coefficients (betan and betas) through which is