yade-users team mailing list archive

[Jan Stránský <question291326@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

Jan Stránský proposed the following answer:
Hi Luis,

I used slightly modified script from #2. The result is 57%.
cheers
Jan

picture result:
https://www.dropbox.com/s/93j2bfhujhcimd7/Screenshot%20from%202016-04-22%2022%3A46%3A19.png?dl=0

script:
#########################################################
#!/usr/bin/python
# -*- coding: utf-8 -*-
#======================================
from yade import pack
import gts
from yade import utils
surf=gts.read(open('tetrahedron.gts'))
pred=pack.inGtsSurface(surf,noPad=True)
vol = surf.volume()
#======================================
# Spheres
O.materials.append(JCFpmMat(type=1,young=1e8,poisson=0.3,frictionAngle=radians(30),density=3000,tensileStrength=1e6,cohesion=1e6,jointNormalStiffness=1e7,jointShearStiffness=1e7,jointCohesion=1e6,jointFrictionAngle=radians(20),jointDilationAngle=0.0,label='sphereso'))
#======================================
O.bodies.append(pack.gtsSurface2Facets(surf,color=(0,0,1)))
sp = pack.randomDensePack(pred,radius=0.08,rRelFuzz=0.01, material =
'sphereo',useOBB=False,returnSpherePack=True,spheresInCell=300)
sp.toSimulation(color=(0.9,0.8,0.6))
#======================================
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
[Ip2_JCFpmMat_JCFpmMat_JCFpmPhys(cohesiveTresholdIteration=1)],
[Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(smoothJoint=False)]
),
GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.2),
NewtonIntegrator(damping=0.9,gravity=[0,0,-9.81]),
]
ve=utils.getSpheresVolume()
print vol,ve,ve/vol
#=======================================
from yade import qt
qt.View()
qt.Controller()
#########################################################

tetrahedron.gts:
4 6 4
-1.63299 -0.942809 -0.666667
0 1.88562 -0.666667
1.63299 -0.942809 -0.666667
0 0 2
2 1
3 2
2 4
1 3
4 1
3 4
1 2 4
5 3 1
3 6 2
6 5 4


2016-04-21 19:13 GMT+02:00 Luis Barbosa <
question291326@xxxxxxxxxxxxxxxxxxxxx>:

> Question #291326 on Yade changed:
> https://answers.launchpad.net/yade/+question/291326
>
>     Status: Answered => Open
>
> Luis Barbosa is still having a problem:
> Hi Jan,
> I didn't know if I got the idea,
> but my volume of spheres is near 39%, and again, is not dense.
> Could you share what you have done?
>
> Here is what I'm using to calc volume of sphere,creating Gts and dense
> pack, respectively:
>
> utils.getSpheresVolume()
>
> pred=pack.inGtsSurface(surf,noPad=True)
>
> pack.randomDensePack(pred,spheresInCell=500,radius=0.08,rRelFuzz=0.01,
> material = 'sphereo',useOBB=False,returnSpherePack=True)
>
> Thanks
> Luis
>
> --
> 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     : yade-users@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-users
> More help   : https://help.launchpad.net/ListHelp
>

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.