Coordination number
Morning Yade community.
I've been working with the eodometric test example. Everything went well implementing the script. However, I would like to estimate the coordination number of the particles at different stages of the test. Right now I'm trying to estimate the coordination number just after the interactions to see how this number change. You're going to find what I've done and if you have any suggestion for me I would appreciate it.
Thanks in advance,
Aurora
$ yade-batch --job-threads=1 03-oedometric-
#
# load parameters from file if run in batch
# default values are used if not run from batch
readParamsFromT
# make rMean, rRelFuzz, maxLoad accessible directly as variables later
from yade.params.table import *
# create box with free top, and ceate loose packing inside the box
from yade import pack, plot
O.bodies.
sp=pack.
sp.makeCloud(
sp.toSimulation()
O.engines=[
ForceResetter(),
# sphere, facet, wall
InsertionSor
InteractionLoop(
# the loading plate is a wall, we need to handle sphere+sphere, sphere+facet, sphere+wall
[
[
[
),
NewtonIntegr
# the label creates an automatic variable referring to this engine
# we use it below to change its attributes from the functions called
PyRunner(
PyRunner(
]
O.dt=.5*
O.step()
def c_number():
for b in O.bodies:
c_num1=
print c_num_1
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply: