yade-users team mailing list archive

[Adel ALBABA <question237037@xxxxxxxxxxxxxxxxxxxxx>]

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

Adel ALBABA posted a new comment:
Christian, jduriez and Bruno.. Thanks for the effective solutions. That
has clarified a lot of YADE capabilities.

Bruno: For free surface problems, I tried the the Tesselation.Wrapper by
looping over the particles and summing the values of Voronoi’s cell of
all spheres (Total Volume) where I can workout the porosity.  The
porosity I am getting is the same as the voxelporosity of a bounding box
(which does not consider the free surface of the problem). I do not know
why  the Tesselation.Wrapper  does not take the irregular surface into
acoount ? :s

Here is the mini-scrip:

# Porosity by TesselationWrapper

TW=TesselationWrapper()
TW.triangulate()
TW.computeVolumes()
volTW=0
for i in O.bodies:
	volTW+=TW.volume(i.id)
TWporo=(1-(solid_volume/volTW))

#------------------------------------
# Porosity by aabbExtrema and voxelporosity

voxelPorosity(resolution=200,start=aabbExtrema()[0],end=aabbExtrema()[1])

Adel

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.