Yade

How to define a material using “Coefficient of rolling elastic limit”?

Asked by Amauri Aguiar de Freitas

Hello every one,
I am trying to reproduce the article “A spherical discrete element model: calibration procedure and incremental response”, but I can’t find a combination of Law2 and Ip2 that permit to use a material with the parameter “Coefficient of rolling elastic limit”. How can I solve this? What is the best combination of Law2 and Ip2 to reproduce the macro behavior of Labene Sand?

Regards,
Amauri

Question information

Language:
English Edit question
Status:
Answered
For:
Yade Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Jérôme Duriez (jduriez) said : 		#1

Hi,

Ip2_CohFrictMat_CohFrictMat_CohFrictPhys and Law2_ScGeom6D_CohFrictPhys_CohesionMoment should be what you want, I guess.

This law makes a distinction between twisting and bending moments though, contrary to the above-mentioned paper, it seems (after a quick look). As such, you have two parameters: CohFrictMat.etaRoll [1] and CohFrictMat.etaTwist [2] that have a similar meaning than the "Coefficient of rolling elastic limit" you mention.

It may be still possible to establish some equivalence with Plassiard et al.'s paper (using maybe CohFrictMat.etaRoll = CohFrictMat.etaTwist = etaYADE, and "Coefficient of rolling elastic limit" = sqrt ( 2 eta^2) ? ), but I would advice you to dig into the corresponding equations [3] to be sure about that, or to wait for other answers.

Jérôme

[1] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.CohFrictMat.etaRoll
[2] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.CohFrictMat.etaTwist
[3] https://github.com/yade/trunk/blob/e4e757f2e98a620e3177b7a36a1d10f69f6a6a28/pkg/dem/CohesiveFrictionalContactLaw.cpp

Can you help with this problem?

Provide an answer of your own, or ask Amauri Aguiar de Freitas for more information if necessary.

To post a message you must log in.