yade-users team mailing list archive

[Rioual <question690280@xxxxxxxxxxxxxxxxxxxxx>]

New question #690280 on Yade:
https://answers.launchpad.net/yade/+question/690280

Hello,

This question concerns the way to implement a new tricky contact law that can not be derived from the other classical 
 contact laws from Yade. It follows the previous advices of Jan (https://answers.launchpad.net/yade/+question/690234).
I want to implement in yade this new contact law that can be expressed as:

Force on normal direction of contact:
F_n = f1(alfa,mu, R_0) * (Normal_strain_rate) +  g(alfa,gamma,R_0)

Force on tangential direction:
F_t = f2(alfa,mu,R_0)*(Tangential_strain_rate)

+ 

d(alfa)/dt = f3( mu,R_0, gamma, Normal_strain_rate, alfa) 

with R_0: average radius of the particle ; mu: viscosity of the material of the particle; gamma: surface tension.
f1,f2,f3,g: functions
alfa(t) is a dynamical geometrical parameter of the contact.


As far as I have seen (https://answers.launchpad.net/yade/+question/681016), I must define:
-  a new Law2 functor
- a new IPhys
- a new lp2 functor

Could you give guidance on the way and methodology to implement this new law in Yade ???
(I am a quite recent Yade user and have very little knowledge of c++...) 

Thank you very much for your assistance in this ambitious project, 

All the best!

Vincent,




-- 
You received this question notification because your team yade-users is
an answer contact for Yade.