yade-users team mailing list archive

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

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

    Status: Open => Answered

Jan Stránský proposed the following answer:
Hi Nguyen


>
> Nguyen N.G. Hien is still having a problem:
> Hello Jan,
> So acording to your suggestion, if I want to avoid the boudary effect by
> counting only the contact point of a zone inside the specimen, for example
> hereby I consider only a smaller box which each edge has size 85% of the
> box, can I use the contact point position instead of bodies position as the
> criteria? For example:
>

I tkink conact point is also good :-)


>
> [...]
>   for i in O.interactions:
>     if (abs(i.geom.contactPoint[0])<=0.85*(width*0.5) and
> abs(i.geom.contactPoint[1])<=0.85*(height*0.5) and
> abs(i.geom.contactPoint[2])<=0.85*(depth*0.5)): continue
> # width height and depth if the dimension of the specimen
>

the operation should be opposite, i.e.

  if (abs(i.geom.contactPoint[0])>0.85*0.5*width and ...): continue #
continue means skip interactions for those the > operation is true is true
(i.e. skip interactions that are too far away from the center of the
specimen)
  doStressComputation()

or

  if (abs(i.geom.contactPoint[0])<0.85*0.5*width and ...): # consider only
interactions near the center
    doStressComputatuion


>
> And I realize that it should be f = i.phys.normalForce +
> i.phys.shearForce (remove the bracket) sothat the Vector3 is callable
> when we run the function :)
>

sorry, it was quick code from scratch.. :-)

cheers
Jan

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