yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #22937
[Question #690280]: Implementing a new contact law very different from the usual ones
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.