Question #690066 on Yade changed:
https://answers.launchpad.net/yade/+question/690066
Status: Open => Answered
Jan Stránský proposed the following answer:
> F_n = f(alfa) * normal_strain_rate + g(alfa)
> F_t = h(alfa) * tangential_strain_rate
thanks for more info, then it is as easy as:
# i is an Interaction instance, e.g. within "for i in O.interactions:"
length = i.geom.refR1 + i.geom.refR2 # or something similar
normal_strain_rate = i.geom.incidentVel / length
tangential_strain_rate = i.geom.shearIncrement / O.dt / length
F_n = f(alfa) * normal_strain_rate + g(alfa)
F_t = h(alfa) * tangential_strain_rate
F = F_n + F_t
# add F (and probably corresponding torque, too) to O.forces for both bodies
(maybe open a new question for this topic if needed)
just a note, all this seems like a task for a Law2 functor (yes, which would
need to be written in C++ and compiled with Yade, but the computation would be
significantly faster....).
So according to your problem, plans, domain size...., consider it as a future
step.
cheers
Jan
--
You received this question notification because your team yade-users is
an answer contact for Yade.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
