Yade

Equations for the contact law "Law2_ScGeom6D_CohFrictPhys_CohesionMoment"

Asked by M.Meidani

Hi,
I am working on direct shear test simulation. I want to use "Law2_ScGeom6D_CohFrictPhys_CohesionMoment" as a contact law and I read a publication about this (Bourrier2013). But it did not contain equations for Kr,Ktw, M max. So I don't know what are parameters like alphaKr or etaRoll.
Does anybody know about the exact formulation of this material type and contact law?

thanks,

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
wangxiaoliang (wangxiaoliang) said : 		#1

you need to calibrate these parameters

Revision history for this message
Bruno Chareyre (bruno-chareyre) said : 		#2

The definition of contact parameters is done by the Ip2 functor, this is where you can find the equations.
See:
https://yade-dem.org/doc/yade.wrapper.html?highlight=ip2_cohfrictmat_cohfrictmat_cohfrictphys#yade.wrapper.Ip2_CohFrictMat_CohFrictMat_CohFrictPhys

Revision history for this message
M.Meidani (masood-meidani) said : 		#3

Thanks for your responses.

I have already seen this link but there is no equation for Maximum moment and also EtaRoll.
Do you know any other documents explaining these equations?

Revision history for this message
Chareyre (bruno-chareyre-9) said : 		#4

You are rigth, this part is missing, I will add it in the doc.
In the meantime, you can read the relationships in the source code
(CohesiveFrictionalContactLaw.cpp:282):

contactPhysics->maxRollPl = min(sdec1->etaRoll*Da,sdec2->etaRoll*Db);
contactPhysics->maxTwistPl = min(sdec1->etaTwist*Da,sdec2->etaTwist*Db);

It is very similar to the equations used for stiffness.

Bruno

Can you help with this problem?

Provide an answer of your own, or ask M.Meidani for more information if necessary.

To post a message you must log in.