yade-users team mailing list archive

[Jan Stránský <honzik.stransky@xxxxxxxxx>]

Hi Behzad,
what version of Yade dou you use? I tried your example without any problem,
after O.run() the particles does not move at all..
cheers
Jan



2014-12-15 13:46 GMT+01:00 behzad <question259096@xxxxxxxxxxxxxxxxxxxxx>:
>
> New question #259096 on Yade:
> https://answers.launchpad.net/yade/+question/259096
>
> Hi,
>
> I'm running the following script. It makes an orthogonal packing in a
> cylindrical shape. I'm wondering why the system has such overlaps which
> gives the particles an explosive velocity. In this case I need to run the
> calm() function to avoid sample collapsing.
> Even with calm() function running at every iteration, if we choose
> O.dt=1.0e-4, particles move.
> Why this is happening?
>
>
> O.reset()
> from yade import utils, plot
> from yade import pack, qt
>
> O.engines=[
> ForceResetter(),
>
> InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=1.0),Bo1_Box_Aabb()]),
> InteractionLoop(
> [Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()],
> [Ip2_FrictMat_FrictMat_FrictPhys()],
> [Law2_ScGeom_FrictPhys_CundallStrack()]),
> NewtonIntegrator(damping=0.7,gravity=[0,0,0])
> ]
>
>
> Co=O.materials.append(FrictMat(young=1e7,poisson=0.3,density=1377,frictionAngle=1))
>
> pred=pack.inCylinder((0.02,0,0),(0.02,0,0.06),0.02)
>
> O.bodies.append(pack.regularOrtho(pred,radius=0.001,gap=0.0, material=Co))
>
>
>
> --
> You received this question notification because you are a member of
> yade-users, which 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
>