tracking a sphere in a pack

Hello Mohammad,

for both your questions there is number of solutions. If you want to track
position of one particle, you can do it with plot module [1, 2]

b = O.bodies[idOfParticleYouWant]
def plotAddData():
  pos = b.state.pos
  plot.addData( i=O.iter, x=pos[0], y=pos[1], z=pos[2] )

...

O.engines = [
  ...
  PyRunner(iterPeriod = 10, command="plotAddData()"), # change iterPeriod
according to your requirements
]

after simulation, plot.data will be a dictionary looking like this:
{
  'i': [10,20,30,...],
  'x':[xPositionAtIter10, xPositionAtIter20, xPositionAtIter30, ... ],
  'y':[yPositionAtIter10, yPositionAtIter20, yPositionAtIter30, ... ],
  'z':[zPositionAtIter10, zPositionAtIter20, zPositionAtIter30, ... ],
}

Concerning number of interactions, if you want to get number of
interactions of one particle, you can do it like

nIntrs = len( b.intrs() ) # [3]

or for all particles

nIntrs = numIntrsOfEachBody() # [4]

you can add this information also to the plotAddData function:

plot.addData( i=O.iter, x=pos[0], y=pos[1], z=pos[2] ,nIntrs=len(b.intrs())
)

cheers
Jan

[1] https://yade-dem.org/doc/yade.plot.html
[2] https://yade-dem.org/doc/user.html#tracking-variables
[3] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.Body.intrs
[4] https://yade-dem.org/doc/yade.utils.html#yade._utils.numIntrsOfEachBody

2013/2/22 MOHAMMAD MOHSENI <email address hidden>

> New question #222563 on Yade:
> https://answers.launchpad.net/yade/+question/222563
>
> Hi
>
> I have a sphere packing which falls into a cylinder whose base is
> rotating. I want to find a way in which I can track only one sphere from
> the packing and find number of interactions. Any suggestions?
>
> --
> You received this question notification because you are a member of
> yade-users, which is an answer contact for Yade.
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-users
> Post to : <email address hidden>
> Unsubscribe : https://launchpad.net/~yade-users
> More help : https://help.launchpad.net/ListHelp
>

many Thanks for your comprehensive respond Jan , hope it works .

Dear Jan

I changed the code as you said but still have errors , can you please kindley help me with this

from yade import pack, plot

O.bodies.append(utils.geom.facetCylinder((.5,.5,.5),radius=0.5,height=1,segmentsNumber=30, wallMask=4))
cylIds=O.bodies.append(utils.geom.facetCylinder((.5,.5,.5),radius=0.6,height=0,wallMask=1))
sp=pack.SpherePack()
O.materials.append(FrictMat(density = 4000, frictionAngle=90,young=1e9,poisson=.11))
sp.makeCloud((0.3,0.3,-0.5),(0.75,0.75,-0.75),rMean=.05,rRelFuzz=0, num=20)
b=O.bodies.append(utils.sphere(center=(0.35,0.35,-0.1),radius=.10,fixed=False))

sp.toSimulation()

O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_L3Geom(),Ig2_Facet_Sphere_L3Geom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_L3Geom_FrictPhys_ElPerfPl()]

),
GravityEngine(gravity=(0,0,9.81)),
NewtonIntegrator(damping=0.4),
RotationEngine(rotateAroundZero=True,zeroPoint=(0.5,0.5,0.5),rotationAxis=(0,0,1),angularVelocity=1000*(2*pi/60),ids=cylIds),
PyRunner(iterPeriod = 10, command="plotAddData()")

]
O.dt=.5*utils.PWaveTimeStep()

def plotAddData():
pos = b.state.pos
plot.addData( i=O.iter, x=pos[0], y=pos[1], z=pos[2] )

nIntrs = len( b.intrs() )
nIntrs = numIntrsOfEachBody()
plot.addData( i=O.iter, x=pos[0], y=pos[1], z=pos[2] ,nIntrs=len(b.intrs()

The Error we are getting is

1. For tracking the partcle p

NameError Traceback (most recent call last)
<string> in <module>()

NameError: name 'plotAddData' is not defined

2. For number of interactions

Traceback (most recent call last):
File "/usr/bin/yade", line 172, in runScript
execfile(script,globals())
File "plotandint.py", line 37

^
SyntaxError: invalid syntax

Many thanks

hi,

from your script i see some problems.
i am not sure, if this solves all errors ...

for 1. :

wrong:

def plotAddData():
pos = b.state.pos
plot.addData( i=O.iter, x=pos[0], y=pos[1], z=pos[2] )

correct:

def plotAddData():
  pos = b.state.pos
  plot.addData( i=O.iter, x=pos[0], y=pos[1], z=pos[2] )

-> you forgot indents

for 2. :

wrong:

plot.addData( i=O.iter, x=pos[0], y=pos[1], z=pos[2] ,nIntrs=len(b.intrs()

correct:

plot.addData( i=O.iter, x=pos[0], y=pos[1], z=pos[2] ,nIntrs=len(b.intrs())

-> you forgot ")"

regards,

christian

Hi Mohammad,

if you wrote to the email complete code, the second error is because of
missing final braces:
plot.addData( i=O.iter, x=pos[0], y=pos[1], z=pos[2] ,nIntrs=len(b.intrs()
should be
plot.addData( i=O.iter, x=pos[0], y=pos[1], z=pos[2] ,nIntrs=len(b.intrs()
) )

with the first error I don't know, I will try to find the solution.

cheers
Jan

2013/2/28 MOHAMMAD MOHSENI <email address hidden>

> Question #222563 on Yade changed:
> https://answers.launchpad.net/yade/+question/222563
>
> Status: Answered => Open
>
> MOHAMMAD MOHSENI is still having a problem:
> Dear Jan
>
> I changed the code as you said but still have errors , can you please
> kindley help me with this
>
> from yade import pack, plot
>
> O.bodies.append(utils.geom.facetCylinder((.5,.5,.5),radius=0.5,height=1,segmentsNumber=30,
> wallMask=4))
>
> cylIds=O.bodies.append(utils.geom.facetCylinder((.5,.5,.5),radius=0.6,height=0,wallMask=1))
> sp=pack.SpherePack()
> O.materials.append(FrictMat(density = 4000,
> frictionAngle=90,young=1e9,poisson=.11))
> sp.makeCloud((0.3,0.3,-0.5),(0.75,0.75,-0.75),rMean=.05,rRelFuzz=0, num=20)
>
> b=O.bodies.append(utils.sphere(center=(0.35,0.35,-0.1),radius=.10,fixed=False))
>
> sp.toSimulation()
>
> O.engines=[
> ForceResetter(),
> InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
> InteractionLoop(
> [Ig2_Sphere_Sphere_L3Geom(),Ig2_Facet_Sphere_L3Geom()],
> [Ip2_FrictMat_FrictMat_FrictPhys()],
> [Law2_L3Geom_FrictPhys_ElPerfPl()]
>
> ),
> GravityEngine(gravity=(0,0,9.81)),
> NewtonIntegrator(damping=0.4),
>
> RotationEngine(rotateAroundZero=True,zeroPoint=(0.5,0.5,0.5),rotationAxis=(0,0,1),angularVelocity=1000*(2*pi/60),ids=cylIds),
> PyRunner(iterPeriod = 10, command="plotAddData()")
>
> ]
> O.dt=.5*utils.PWaveTimeStep()
>
> def plotAddData():
> pos = b.state.pos
> plot.addData( i=O.iter, x=pos[0], y=pos[1], z=pos[2] )
>
> nIntrs = len( b.intrs() )
> nIntrs = numIntrsOfEachBody()
> plot.addData( i=O.iter, x=pos[0], y=pos[1], z=pos[2] ,nIntrs=len(b.intrs()
>
> The Error we are getting is
>
> 1. For tracking the partcle p
>
> NameError Traceback (most recent call last)
> <string> in <module>()
>
> NameError: name 'plotAddData' is not defined
>
> 2. For number of interactions
>
> Traceback (most recent call last):
> File "/usr/bin/yade", line 172, in runScript
> execfile(script,globals())
> File "plotandint.py", line 37
>
> ^
> SyntaxError: invalid syntax
>
>
> Many thanks
>
> --
> You received this question notification because you are a member of
> yade-users, which is an answer contact for Yade.
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-users
> Post to : <email address hidden>
> Unsubscribe : https://launchpad.net/~yade-users
> More help : https://help.launchpad.net/ListHelp
>

I supposed that indents were missing from copy-paste process.. the function
plotAddData work and did not result in SyntaxError, but I don't know yet
how to solve the NameError..
Jan

2013/2/28 Christian Jakob <email address hidden>

> Question #222563 on Yade changed:
> https://answers.launchpad.net/yade/+question/222563
>
> Christian Jakob posted a new comment:
> hi,
>
> from your script i see some problems.
> i am not sure, if this solves all errors ...
>
> for 1. :
>
> wrong:
>
> def plotAddData():
> pos = b.state.pos
> plot.addData( i=O.iter, x=pos[0], y=pos[1], z=pos[2] )
>
> correct:
>
> def plotAddData():
> pos = b.state.pos
> plot.addData( i=O.iter, x=pos[0], y=pos[1], z=pos[2] )
>
> -> you forgot indents
>
>
> for 2. :
>
> wrong:
>
> plot.addData( i=O.iter, x=pos[0], y=pos[1], z=pos[2]
> ,nIntrs=len(b.intrs()
>
> correct:
>
> plot.addData( i=O.iter, x=pos[0], y=pos[1], z=pos[2]
> ,nIntrs=len(b.intrs())
>
> -> you forgot ")"
>
>
> regards,
>
> christian
>
> --
> You received this question notification because you are a member of
> yade-users, which is an answer contact for Yade.
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-users
> Post to : <email address hidden>
> Unsubscribe : https://launchpad.net/~yade-users
> More help : https://help.launchpad.net/ListHelp
>

Ok , I will send the whole code to your Email (<email address hidden> )
Thanks for your helps

Thanks Jan Stránský, that solved my question.