Question #699589 on Yade changed:
https://answers.launchpad.net/yade/+question/699589
Status: Open => Answered
Jan Stránský proposed the following answer:
Hello,
there might be several reasons, difficult to tell without actual data
[1]. For a general case:
A) contact normal vector vs. contact normal force:
> ... contact normal force vector ...
> print("force",i.phys.normalForce + i.phys.shearForce)
you are printing TOTAL force, not normal force.
If i.geom.shearForce is non-zero, then you naturally get non-parallel vector
O.forces.f(i.id1) result is "independent" of constituent interactions
B) branch vector vs. contact normal vector:
i.geom.normal is computed in InteractionLoop.
Body positions are THEN updated by NewtonIntegrator.
So b2.state.pos - b1.state.pos and i.geom.normal (in general) are not parallel.
Cheers
Jan
[1] https://www.yade-dem.org/wiki/Howtoask
--
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
