← Back to team overview

yade-users team mailing list archive

Re: [Question #690066]: Using O.forces along normal and tangential directions at contact

 

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.