questions about Kn and Ks

Asked by Fu zuoguang

Dear Jan Stránsky :
     I am sorry for providing a unclear question, so I have give some more details zbout it:
(1). According to my simulation requirement, all the particls in specimen must have the same
Kn and Ks which are all equal to 4e5.
(2). From these equations:
     Kn=2.*E*R1*E*R2/(E*R1+E*R2) (1)
     Ks/Kn=P (2)
I can obviously understand that all particles can not have the same values of "young" for that the radii of them are not the same. But I failed to get a correct input of "young" because the order may force me to do in this way(one type material can only be assigned only one value of "young"):

# spheremat = O.materials.append(FrictMat(poisson=0.5,density=1500,young=1.48e8,frictionAngle=0.7))

How can I do to solve this problem?
Seeking your help!

Question information

Language:
English Edit question
Status:
Solved
For:
Yade Edit question
Assignee:
No assignee Edit question
Solved by:
Bruno Chareyre
Solved:
Last query:
Last reply:
Revision history for this message
Jan Stránský (honzik) said :
#1

> (1). According to my simulation requirement, all the particls in specimen
> must have the same
> Kn and Ks which are all equal to 4e5.
>

Kn and Ks are attributes of interactions. when new interaction is created,
this Kn and Ks are evaluated according to formulas from your previous email.

> (2). From these equations:
> Kn=2.*E*R1*E*R2/(E*R1+E*R2) (1)
> Ks/Kn=P (2)
> I can obviously understand that all particles can not have the same values
> of "young" for that the radii of them are not the same. But I failed to get
> a correct input of "young" because the order may force me to do in this
> way(one type material can only be assigned only one value of "young"):
>
> # spheremat =
> O.materials.append(FrictMat(poisson=0.5,density=1500,young=1.48e8,frictionAngle=0.7))
>
>
in your script, tou have unifrrom radius, so you can identify E and P from
equations (1) and (2). If you have nonuniform particle radii and you want
Kn and Ks, this is quite non-standard requirement (btw, why do you want is
so?). Just now I don't know how to solve it with existing code, I will
think about it..

cheers
Jan

Revision history for this message
Jan Stránský (honzik) said :
#2

fixed typo from my previous answer:
In your script, you have unifrrom radius, so you can identify E and P
from equations
(1) and (2). If you have nonuniform particle radii and you want smae Kn and
Ks for all interactions, this is quite non-standard requirement (btw, why
do you want it so?).

Revision history for this message
tranviet (tranviet2302) said :
#3

Hi,

I think there is another way, but I really don't recommend it.

If you are using the source code, you can go to the Ip2 functor, for instance: Ip2_FrictMat_FrictMat_FrictPhys.cpp and set:

contactPhysic->kn = your value
contactPhysic->ks = your value

Viet

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

>contactPhysic->kn = your value
>contactPhysic->ks = your value

Since "4e5" is hardcoded in the question, this is a smart answer actually.
I vote for it.

Revision history for this message
Fu zuoguang (zgfu1985) said :
#5

Thanks Bruno Chareyre, that solved my question.