yade-users team mailing list archive

[lida <question237753@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Answered => Open

lida is still having a problem:
Hi Alex @Alexander Eulitz,

Thank you for your help. That is  great! However, since I don't have any
background in yade. I don't know how to define the velocity and
timestep. What does "timestep=O.dt # in seconds per iteratino" mean?

I can show you my poor command:
from yade import pack
    
pred=pack.inCylinder(centerBottom=(0,0,0),centerTop=(0,0,0.012),radius=0.002))
  
O.materials.append(FrictMat(young=200e9,poisson=.3,frictionAngle=.5,label="steel")


spheres=pack.randomDensePack(pred,radius=0.001,color=None,material="steel",sphereInCell=200,memoDbg=True,memoizeDb=None)

O.bodies.append(spheres)

spheres.toSimulation()

O.engines=[
   
   ForceResetter(),

InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb(),Bo1_Wall_Aabb()]),

   InteractionLoop(

[Ig2_Sphere_Sphere_L3Geom(),Ig2_Facet_Sphere_L3Geom(),Ig2_Wall_Sphere_L3Geom()],

      [Ip2_FrictMat_FrictMat_FrictPhys()],
      [Law2_L3Geom_FrictPhys_ElPerfPl()]
   ),

   GravityEngine(gravity=(0,0,-9.81)),

   NewtonIntegrator(damping=0.4),

   PyRunner(command='velocity()'),

]

O.dt=.5*util.PWaveTimeStep()

def velocity():

   velocity=0,0001

   timeNeed=0.006/velocity

   iterNeed=timeNeed/
O.run()
O.stopAtlter(iterNeeded)
plot.saveGnuplot(baseName=plotpath)

Could you give me any specific informations to finish this model?

Regards,

Lida

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