← Back to team overview

yade-users team mailing list archive

Modification in Cundall damping


Hello all,

There is few modifications to be made in the equations of damping (force 
and momentum).
The current equations are of the type : for i=1,...,3 { f[i] *= 1- 
damping * sign(f[i]*vel[i]); }

in fact it should be : f *= 1- damping * sign(f.vel);  (where *= denotes 
product between vector and real and "." is the scalar product of two 

It can make a difference in some cases.

If everybody agree, I'll try to correct this and commit the change to 
the SVN repository. In the meantime I suggest that everybody doing 
calculation with Yade change the corresponding lines in force and 
momentum damping engines. (Luc : does this little sphere roll to the 
bottom of the big sphere with these new equations?)

Yade-users mailing list

Follow ups