yade-users team mailing list archive

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

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

    Status: Open => Answered

Jan Stránský proposed the following answer:
Something like this?
#####################
r = 1
O.bodies.append(wall(6*r,2))
def inject():
   ss1 = (
      sphere((0,0*2**r,0*r),r),
      sphere((0,1*2**r,2*r),r),
      sphere((0,2*2**r,0*r),r),
      sphere((0,3*2**r,2*r),r),
   )
   #
   ss2 = (
      sphere((0,6*2**r,0*r),r),
      sphere((0,7*2**r,2*r),r),
      sphere((0,8*2**r,0*r),r),
      sphere((0,9*2**r,2*r),r),
   )
   #
   O.bodies.append(ss1)
   cid,bids = O.bodies.appendClumped(ss2)
   v = (1,0,1)
   for s in ss1 + ss2:
      s.state.vel = v
   O.bodies[cid].state.vel = v

O.engines = [
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Wall_Aabb()]),
   InteractionLoop(
      [Ig2_Sphere_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom()],
      [Ip2_FrictMat_FrictMat_FrictPhys()],
      [Law2_ScGeom_FrictPhys_CundallStrack()],
   ),
   NewtonIntegrator(),
   PyRunner(iterPeriod=500000,command='inject()',initRun=True),
]
O.dt = 1e-5
O.step()
#####################
cheers
Jan

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