Why the sum of TW.volume is not equal to the volume of pack
Dear all,
I use TesselationWrapper to compute the volume of spheres, but I find that the sum of TW.volume is not equal to the volume of pack. Does the Voronoi cell of each sphere overlap each other?
I run the MWE and I got
volume=
TW_Volume=
Thanks for any suggestion.
####MWE####
import numpy as np
O.periodic=True
O.cell.
num_spheres=1000 # number of spheres
den_ball=2600 # density of particles from the experimental test
young=1e8
iso_pressure=-1e4
compFricDegree = 30 # initial contact friction during the confining phase
O.materials.
sp=pack.
sp.makeCloud(
sp.toSimulation
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
NewtonIntegrat
PeriTriaxContr
]
O.dt=PWaveTimeS
O.run();O.wait()
TW=TesselationW
TW.triangulate()
TW.setState()
TW.computeVolumes()
TW_Volume=0
for i in range(1000):
TW_Volume += TW.volume(i)
print('
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Chareyre
- Solved:
- Last query:
- Last reply: