← Back to team overview

yade-users team mailing list archive

Re: [Question #291326]: Packing and predicate dimension

 

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 <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,
>
> 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@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.