How to assign the values of kn ks of particles in Yade directly

Hello,

sorry, your question is too ambiguous, please be more specific (especially what is wrong)

> the contacts can have the same kn and ks, what should I do?

simple answer, you should do nothing special :-)

yes, contacts can have the same kn and ks (depending on meaning):
- kn=ks would be for poisson=1.0
- kn_interaction1=kn_interaction2 would be if their material properties and radii related to interactions have mutually special properties.

Is "same kn and ks" (whatever it means) a problem? why? why would you want it different?

thanks
Jan

Dear Jan,

Thanks for your suggestions. My question is actually the second one you mentioned that all the interactions have the same kn values and also the same ks values.

 Actually, you know I used PFC before, and in PFC we can directly assign the values of kn, ks in a simple linear contact model. But when I use Yade to do the same simulation to compare both results, I found that in Yade the particle materials properties like young, and poisson need to be assigned when using the p2_FrictMat_FrictMat_FrictPhys() functor. And the normal and shear stiffness of particles are determined based on the particle radius and also the materials properties. So I just wonder how to directly assign their values like kn=1.0e8, ks=1.0e8*0.3? I think it needs to modify some source code maybe? I have little experience about it.

Regards,
jiayan

Dear Jan,
Thanks. I have tried c), and I found that is was not useful. But in option a), after some minor improvements of source code, that solved my problem.

Thanks Jan Stránský, that solved my question.

Hello,
for future reference. could you please describe why c) was not useful?
thanks
Jan

Hi,

> after some minor improvements

Should read "after some minor regressions" IMO. :)
I think c) was a smart answer and I'm also wondering why it didn't help.
Bruno

For future reference, a MWE how to achieve this without modification of source code:

###
kn = 1e6
radii = [0.7, 0.8, 0.9, 1.0]
centers = [
   (0,0,0),
   (1,0,0),
   (0,1,0),
   (0,0,1),
]
youngs = [kn/r for r in radii]
materials = [FrictMat(young=young,poisson=0.3) for young in youngs]
O.bodies.append([sphere(c,r,material=mat,fixed=True) for c,r,mat in zip(centers,radii,materials)])
O.step()
for i in O.interactions:
   print i.phys.kn,i.phys.ks # prints "1000000.0 300000.0" for 6 interactions
###

cheers
Jan

Hi, Jan,

Thanks for your MWE. It functions well. but indeed when using this method, for a fixed kn, the sphere with a smaller radius will have a larger young, which results in a small critical time step. That's really time-consuming when conduct simulations with a gap-graded specimen.

Maybe for a simple linear contact model, like in PFC, we can directly assign the same values of normal and shear stiffness kn ,ks for all the contacts, and the critical time step is depends on the minimum of math.sqrt((mass/kn)) of the particle. I do not know if Yade has the same contact model since I am a novice. So I modified some source code to achieve this simple contact model.

Hi, Bruno, I just modified some source code to achieve the above mentioned contact model. Thanks.

Hi, Bruno, I just modified some source code to achieve the above mentioned contact model. Thanks.

"the sphere with a smaller radius will have a larger young, which results in a small critical time step"

How do you define your time step ? I would imagine that behavior may happen if you use PWaveTimeStep (?) but I'm expecting it would not with GlobalStiffnessTimeStepper (which is the way to go), that directly uses the kn.

> for a fixed kn, the sphere with a smaller radius will have a larger young, which results in a small critical time step.

Absolutely.
And that's exactly what happens in PFC when you assign a unique kn.
Got it?

Bruno

p.s. @Jérôme, the effect is exactly the same with both PWaveTS or GSTS