yade-users team mailing list archive

[lingran <question234129@xxxxxxxxxxxxxxxxxxxxx>]

Question #234129 on Yade changed:
https://answers.launchpad.net/yade/+question/234129

    Status: Answered => Open

lingran is still having a problem:
Thanks a lot, Jan, you are saving my life. 
But I am still a little confused. I did a small test, it turns out i.phys.shearForce is equal to Ft rather than Ftt. 
Did you miss something after O.step()?
see:

i=O.interactions[975,2702]
Ft=i.phys.shearForce  # -- Ft=(-263.88, -780.43,0)
O.step() # now you use the same value of Ft as in c++
du=i.geom.shearInc
kt=i.phys.ks
Ftt=Ft-kt*du   # --Ftt=(-274.83, -812.82, 0)
maxFs=Fn.norm()*i.phys.tangensOfFrictionAngle  #--maxFs=223.55
If (Ftt.norm() > maxFs):
      ratio = maxFs / Ftt.norm()
      trialForce=Ftt
      Ftt *= ratio  # Ftt=(-71.61, -211.77,0)
      e=((1/kt)*(trialForce-Ftt)).dot(Ftt)
      i.phys.shearForce     # -- i.phys.shearForce =(-263.88, -780.43,0)


Best regards:
Lingran

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.