how to use clumps

Asked by Cyurs on 2011-08-27

When I run the following script

from yade import pack
spheremat=O.materials.append(FrictMat(density=2600,young=150e6,poisson=0.5,frictionAngle=30))
a=[]
for i in range(3):
 a.append(pack.SpherePack())
 r=0.03*i+0.02
 a[i].add((0,0,0),r)
 a[i].add((0,0,1.5*r),r)

sp=pack.SpherePack()
sp.makeClumpCloud((0,0,0),(2,2,10),a)
sp.toSimulation(material=spheremat)
wallmat=O.materials.append(FrictMat(density=1000,young=150e6,poisson=0,frictionAngle=0))
O.bodies.append(utils.facetBox((1,1,5),(1,1,5),wallMask=31,material=wallmat))

O.engines=[
 ForceResetter(),
 InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
 InteractionLoop(
  [Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
  [Ip2_FrictMat_FrictMat_FrictPhys()],
  [Law2_ScGeom_FrictPhys_CundallStrack()]
 ),
 GravityEngine(gravity=(0,0,-9.81)),
 NewtonIntegrator(damping=0.4),
]
O.dt=.5*utils.PWaveTimeStep()

I get
Body #2570 has velocity==NaN!
What's wrong

Thanks

Question information

Language:
English Edit question
Status:
Solved
For:
Yade Edit question
Assignee:
No assignee Edit question
Solved by:
Anton Gladky
Solved:
2011-09-10
Last query:
2011-09-10
Last reply:
2011-09-07
Best Anton Gladky (gladky-anton) said : #1

I have ran your script for about 50000 iterations and no errors appeared.
Try to use the latest trunk version.

Cyurs (hexuzhen-1988) said : #2

Thanks Anton Gladky, that solved my question.