yade-users team mailing list archive

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

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

    Status: Open => Answered

Jan Stránský proposed the following answer:
> I am new to Yade and in the available documents there is not much
information about it.

Yade uses Python as an user interface. In python, you have access to almost any data in the simulation (tutorial [1] and user's manual [2] is full of this).
So you have access to positions of all particles, then it is easy to check if they are in a cylinder or not and sum mass of those which are. There are several options e.g.
###
def isBodyInCylinder(body,center,radius):
   d = body.state.pos - center # vector between cylinder center and particle center
   d[2] = 0 # only xy components are relevant
   distance = d.norm() # horizontal distance
   return distance < radius # particle is inside, if distance < raidus
def masses():
   mass1 = mass2 = 0.0 # intially zero
   center1 = Vector3(xp1,yp1,0) # center of 1st cylinder
   center2 = Vector3(-xp2,yp2,0) # center of 2nd cylinder
   for b in O.bodies: # test all particles
      if not isinstance(b.shape,Sphere): continue # skip non-spherical
      if isBodyInCylinder(b,center1,radius): # if a particle is inside 1st cylinder ...
         mass1 += b.state.mass # ... increase mass1 by the particle mass
      if isBodyInCylinder(b,center2,radius): # if a particle is inside 2nd cylinder ...
         mass2 += b.state.mass # ... increase mass2 by the particle mass
   return mass1,mass2 # return computed values
###

cheers
Jan

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