Number of neighbouring sphere particles

Hello,

If your distance threshold matches the one used to define interactions in your YADE simulation you may simply use https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.InteractionContainer.withBody (and ask for the number of elements in the list returned by this function)

Jerome

hello,
thank you for the answer. Unfortunately i would require a much bigger radius
Cheers
John

Hi John,
it pretty much depends on the definition of " immediate neighbours (or neighbours within a distance)". Below there is a simple script counting with distance of centers (regardless particle size)
cheers
Jan

##############################
dist = 3.5 # distance to find "neighbors"

# creates bodies
bodies = (
   sphere(( 0,0,0),1),
   sphere((+5,0,0),1),
   sphere((-3,0,0),1),
)
O.bodies.append(bodies)

# save bodies to SpherePack
sp = SpherePack()
sp.fromSimulation()

# use the data in a separate simulation
O.switchScene()
sp.toSimulation()
for b in O.bodies:
   b.shape.radius = .5*dist
O.dt = 0.
O.engines = [
   InsertionSortCollider([Bo1_Sphere_Aabb()],verletDist=0.),
   InteractionLoop(
      [Ig2_Sphere_Sphere_ScGeom()],
      [Ip2_FrictMat_FrictMat_FrictPhys()],
      [Law2_ScGeom_FrictPhys_CundallStrack()],
   )
]
O.step() # create interactions
print len(O.bodies[0].intrs()) # number of interactions (neighbors) [1]

O.switchScene() # return back to original simulation
##############################

[1] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.Body.intrs