yadeusers team mailing list archive

yadeusers team

Mailing list archive

Message #08240
Re: [Question #235934]: Reference for coeff of restitution in HertzMindlin with nonlinear damping
Question #235934 on Yade changed:
https://answers.launchpad.net/yade/+question/235934
Anton Gladky proposed the following answer:
Hi Vaclav,
probably this one:
http://dx.doi.org/10.1016/00325910(92)88030L
Anton
2013/9/18 Václav Šmilauer <question235934@xxxxxxxxxxxxxxxxxxxxx>:
> New question #235934 on Yade:
> https://answers.launchpad.net/yade/+question/235934
>
> I use the simplesceneplot.py script to test the coefficient of restition for viscous contact laws, only replacing the two lines in the InteractionLoop with
>
> [Ip2_FrictMat_FrictMat_MindlinPhys(en=.5)],
> [Law2_ScGeom_MindlinPhys_Mindlin()],
>
> and setting NewtonIntegrator(damping=0.,...).
>
> The scene is a sphere in free fall onto a box, then rebounding, falling down again etc, until all energy is dissipated. According to the documentation and source code, *en* is the normal coefficient of restitution, so I expect to see the energy quartered (1/4=.5^2) after each rebound. The sphere is originally with v=0 at z=2, entering in contact with the box when z=1.5. Therefore, after one rebound, it should fly up to about 1.625 etc, eventually approaching 1.5.
>
> In the simulation, the sphere peaks at about z=1.55 after the first rebound. Considering that the equilibrium position will be at about z=1.45 due to still weight of the sphere in gravity, I still obtain coefficient of restitution of about sqrt((1.551.45)/(21.45))=.42, not 0.5.
>
> The source (pkg/dem/HertzMindlin.cpp, lines 89 and 320) references some (Tsuji, 1992) paper for computing the viscous coefficient from the coefficient of restitution (cn=alpha*sqrt(mbar), with alpha=sqrt(5/6.)*2*log(en)/sqrt(log(en)^2)+M_PI^2)*sqrt(2*E*sqrt(R) and mbar=(m1*m2)/(m1+m2)). I am not able to identify this paper, can I ask for a full reference (emailing fulltext is highly appreciated)?
>
> How reliable is this code in Yade? Has someone doublechecked that it does what it should?
>
> Cheers, Václav
>
> 
> You received this question notification because you are a member of
> yadeusers, which is an answer contact for Yade.
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yadeusers
> Post to : yadeusers@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yadeusers
> More help : https://help.launchpad.net/ListHelp

You received this question notification because you are a member of
yadeusers, which is an answer contact for Yade.