Question #291326 on Yade changed:
https://answers.launchpad.net/yade/+question/291326
Status: Open => Answered
Jan Stránský proposed the following answer:
Hi Luis,
thanks for the files. For some reason, the "smooth" surface, created by a
wall in triaxial compression, is placed wrongly. There are several
solutions, all parameters of randomDensePack, see docs [1]:
spheresInCell=500 (or another number). It does internally periodic triaxial
compression, overcoming the problem with fixed wall position
cropLayers=3, it simply enlarge the compressed volume
dim, according to the documentation, it could also help
my suggestion is 1), since it also reduces time needed for preparation.
cheers
Jan
[1] https://yade-dem.org/doc/yade.pack.html#yade.pack.randomDensePack
2016-04-19 2:17 GMT+02:00 Luis Barbosa <question291...@answers.launchpad.net
>:
> Question #291326 on Yade changed:
> https://answers.launchpad.net/yade/+question/291326
>
> Status: Answered => Open
>
> Luis Barbosa is still having a problem:
> Hi Jan,
>
> Sure I can.
>
> #!/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)
> 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)
> 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
> #=======================================
> from yade import qt
> qt.View()
> qt.Controller()
> #########################################################
>
> And here follows the gts surface...
> https://drive.google.com/open?id=0B9LhrwAf_vasanZ4bVc4b0FHVUE
>
> 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@lists.launchpad.net
> 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.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
