yade-users team mailing list archive

[Surej Kumar <question672270@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Answered => Open

Surej Kumar is still having a problem:
Thanks Jan for the reply.

I just try to run the packing in 3D but it is taking ages (almost 4
hours sand still not  completed for 400 grains.)

s@s-VirtualBox:~/Desktop/TEST4August2018/Triax_results$ yade prepareRVE.py 
Welcome to Yade 1.20.0 
TCP python prompt on localhost:9000, auth cookie `dsayuc'
XMLRPC info provider on http://localhost:21000
Running script prepareRVE.py
400
3.16023157237e-05

below is my script

sp = pack.SpherePack()
size =.24
sp.makeCloud(minCorner=(0,0,0),maxCorner=(size,size,.05),rMean=.005,rRelFuzz=0.4,num=400,periodic=True,seed=1)
sp.toSimulation()
O.cell.hSize = Matrix3(size,0,0, 0,size,0, 0,0,.1)
print len(O.bodies)
for p in O.bodies:
   p.state.blockedDOFs = 'zXY'
   p.state.mass = 2650 * 0.1 * pi * p.shape.radius**2 # 0.1 = thickness of cylindrical particle
   inertia = 0.5 * p.state.mass * p.shape.radius**2
   p.state.inertia = (.5*inertia,.5*inertia,inertia)

O.dt = utils.PWaveTimeStep()
print O.dt

O.engines = [
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb()]),
   InteractionLoop(
      [Ig2_Sphere_Sphere_ScGeom()],
      [Ip2_FrictMat_FrictMat_FrictPhys()],
      [Law2_ScGeom_FrictPhys_CundallStrack()]
   ),
   PeriTriaxController(
      dynCell=True,
      goal=(-1.e5,-1.e5,0),
      stressMask=3,
      relStressTol=.001,
      maxUnbalanced=.001,
      maxStrainRate=(.5,.5,.0),
      doneHook='term()',
      label='biax'
   ),
   NewtonIntegrator(damping=.1)
]

def term():
   O.engines = O.engines[:3]+O.engines[4:]
   print getStress()
   print O.cell.hSize
   setContactFriction(0.5)
   O.cell.trsf=Matrix3.Identity
   O.cell.velGrad=Matrix3.Zero
   for p in O.bodies:
      p.state.vel = Vector3.Zero
      p.state.angVel = Vector3.Zero
      p.state.refPos = p.state.pos
      p.state.refOri = p.state.ori
   O.save('0.yade.gz')
   O.pause()

O.run();O.wait()

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