yade-users team mailing list archive

[Jan Stránský <question293669@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

Jan Stránský proposed the following answer:
Hi Behzad,

      const Vector3r relativeVelocity = (de1.vel+de1.angVel.cross(c1x)) -
> (de2.vel+de2.angVel.cross(c2x)) + shiftVel;
>       const Vector3r normalVelocity = geom.normal.dot(relativeVelocity)*
> geom.normal ;
>       const Vector3r shearVelocity = relativeVelocity-normalVelocity;


these values should be available in IGeom class, you don't need to compute
it "manually". Code duplication is always possible source of problems.

anyway, I think now the best way how to find the problem is to compare your
law to some existing law, which works ok (CohFrictMat?).

Another option is some debugging, checking actual values and comparing them
to expected/correct/analytical values)..

cheers
Jan


2016-05-14 2:17 GMT+02:00 behzad <question293669@xxxxxxxxxxxxxxxxxxxxx>:

> Question #293669 on Yade changed:
> https://answers.launchpad.net/yade/+question/293669
>
> behzad gave more information on the question:
>
> second update:
>
> the complete normal situation is when I use both:
>
> geom->rotate(shearForce)
>
> and
>
> geom->rotate(normalForce)
>
> Is this normal?
>
> --
> 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     : yade-users@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-users
> More help   : https://help.launchpad.net/ListHelp
>

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.