a non-linear contact law type Hertz-Mindlin

Asked by emnaemna

Hii,

I want to use a non linear contact law type Hertz-Mindlin.
I want to ask if the following command are true . And if these command translate a non liner contact law type Hertz-Mindlin
InteractionLoop(
  [Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()],
  [Ip2_FrictMat_FrictMat_MindlinPhys(en=0.1,es=0.1)],
  [Law2_ScGeom_MindlinPhys_Mindlin()]
 ),

Also I want to know the mathematical model that modelise this Hertz-Mindlin contact law . To be more accuracy, , i want to know the mathematical equation thay Yade utilize to cmpute the contact force .

Best Regards

Question information

Language:
English Edit question
Status:
Solved
For:
Yade Edit question
Assignee:
No assignee Edit question
Solved by:
Klaus Thoeni
Solved:
Last query:
Last reply:
Revision history for this message
Best Klaus Thoeni (klaus.thoeni) said :
#1

Hi,

the beauty of open-source is that you can check the equations in the source code :-)

If you don't have the source code on your computer you can find all the equations used for Law2_ScGeom_MindlinPhys_Mindlin here [1]. Some variations of the Hertz-Mindlin contact law are available (e.g. Law2_ScGeom_MindlinPhys_MindlinDeresiewitz, Law2_ScGeom_MindlinPhys_HertzWithLinearShear).

Your IntersectionLoop looks right to me.

HTH
Klaus

[1] https://github.com/yade/trunk/blob/master/pkg/dem/HertzMindlin.cpp

Revision history for this message
Jérôme Duriez (jduriez) said :
#2

Hi,

We may add that your interaction loop, provided that your bodies have "FrictMat" as material (which is the default behavior).

Jérôme

Revision history for this message
Jérôme Duriez (jduriez) said :
#3

"that your interaction loop IS OK", sorry...

Revision history for this message
emnaemna (emna-masmoudi92) said :
#4

Thanks a lot for all of you . Now I'm sure that my code is correct :)

Revision history for this message
emnaemna (emna-masmoudi92) said :
#5

Thanks Klaus Thoeni, that solved my question.