Coordination Number

Hello,

roughly:

1) select the relevant spheres, whose center is inside a box
sphs = [b for b in O.bodies if isInsideRelevantBox(b)]

2) compute number of interactions of each sphere using b.intrs() [1]
you can also use information about shape to exclude sphere-nonsphere contact if it is not desirable

3) averaging is trivial

cheers
Jan

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

Thanks Jan for answering.
By the way, I can't understand how I could define the condition of being in the box I want. Moreover, which is the difference between intrs() and avgNumInteractions?
Thanks

> I can't understand how I could define the condition of being in the box I want.

please be more specific what you can't understand (principles, how to code it, how to use values from yade, ...).
It also pretty much depends on what "box" actually is.

> which is the difference between intrs() and avgNumInteractions?

documentation [2,3] is clear IMO. The main difference is returned type:
- b.intrs() returns a **list of interactions**
- avgNumInteractions() returns average **number** of interactions per particle

the other difference is domain considered:
- b.intrs() deals with one body
- avgNumInteractions() deals with all bodies (hence not suitable for your task)

cheers
Jan

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

My box is a facetbox with its top open. Let's suppose that I'd like to determine the averege coordination number of the spheres which are inside:

geom.facetBox((0,0,0),(d,l,s),wallMask=31)

You suggested me to detect the relevant spheres whose center is inside the facetbox

sphs = [b for b in O.bodies if isInsideRelevantBox(b)]

I can't understande how to code the isInsideRelevantBox condition.

I think it could be something related to the x,y,z relative position between a center and the box. Basically a sphere is inside it

if -d/2<b.state.pos[0]<d/2 and -l/2<b.state.pos[1]<l/2 and -h/2<b.state.pos[2]<h/2

Hi, I confirm your self answer.
You could put this in a function to be used in the condition.
Bruno

Le ven. 6 mars. 2020 14:17, Paolo <email address hidden> a
écrit :

> Question #689086 on Yade changed:
> https://answers.launchpad.net/yade/+question/689086
>
> Status: Answered => Open
>
> Paolo is still having a problem:
> My box is a facetbox with its top open. Let's suppose that I'd like to
> determine the averege coordination number of the spheres which are
> inside:
>
> geom.facetBox((0,0,0),(d,l,s),wallMask=31)
>
> You suggested me to detect the relevant spheres whose center is inside
> the facetbox
>
> sphs = [b for b in O.bodies if isInsideRelevantBox(b)]
>
> I can't understande how to code the isInsideRelevantBox condition.
>
> I think it could be something related to the x,y,z relative position
> between a center and the box. Basically a sphere is inside it
>
> if -d/2<b.state.pos[0]<d/2 and -l/2<b.state.pos[1]<l/2 and
> -h/2<b.state.pos[2]<h/2
>
> --
> 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 : <email address hidden>
> Unsubscribe : https://launchpad.net/~yade-users
> More help : https://help.launchpad.net/ListHelp
>
>
>