yade-users team mailing list archive

[feda <question404455@xxxxxxxxxxxxxxxxxxxxx>]

New question #404455 on Yade:
https://answers.launchpad.net/yade/+question/404455

Hello,
I am simulating a periodic isotropic compression of well graded material (Cu>3) to reach the stable loosest state (target porosity).
I use ViscElMat to define the material (Kn,Ks,cn,cs,frictionangle,density):
-The friction angle is set to 0.7 (the porosity will increase by increasing the friction)
-cn and cs are in functions of a predefined critical damping ratio, kn and the mean radius of particles. 
-The damping in Newton Integrator is set to 0. The energy dissipation is then introduced in the system by means of a viscous damping defined at contacts.
This set of parameters allows obtaining porosity smaller than the target porosity (a more dense state) and if I increase the values of cn and cs by increasing the critical damping ratio, I obtain the same value of porosity.

Now, if I set cn and cs to 0 and I use a global numerical damping instead, I can reach the target porosity (obtain a loosest state).
How we can interpret this difference, and which method is better to use with viscElMat?  
Thank you
feda


-- 
You received this question notification because your team yade-users is
an answer contact for Yade.