yade-users team mailing list archive

[hafsa <question678583@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Needs information => Open

hafsa gave more information on the question:
hi Jan ,
the fist code is :
#material of spheres
Mat2 = O.materials.append(CohFrictMat(poisson=4.,young=2.4e6,frictionAngle=radians(30),density=25.0))

#material of cylinder
Mat1 = O.materials.append(ViscElMat(kn=10.0e4,ks=10.0e4,frictionAngle=0.455,density=2650.0))

# create the big sphere
big_sphr=O.bodies.append(utils.sphere((0,0,16),4.,fixed=False,color=(1,0,0),highlight=True,material=Mat2))

# create sphere packing
sp=pack.SpherePack()
sp.makeCloud((-7.2,-7.2,18.),(7.2,7.2,82.9),rMean=.5,rRelFuzz=0,num=6000,seed=1)
sp.toSimulation(color=(0,1,0),material=Mat2)

# create Cylinder 
Cylinder=O.bodies.append(geom.facetCylinder((0,0,30.),10.0,height=60.0,wallMask=6,segmentsNumber=64,material=Mat1))

# interactions 
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()]
    ),
    NewtonIntegrator(gravity=(0,0,-9.81),damping=0.3),
 ]
O.dt=utils.PWaveTimeStep() 
O.run(11200)
O.save('init_config.yade.gz')

the second code is :

O.load('init_config.yade.gz')

# interactions and oscillation
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()]
    ),
    NewtonIntegrator(gravity=(0,0,-9.81),damping=0.3),
    HarmonicMotionEngine(A=[0,0,0.5], f=[0,0,1.0], fi = [0.0,0.0,0.0],ids =Cylinder)
 ]
O.dt=utils.PWaveTimeStep()
O.run(17200)

the problem is that the simulation in the second code stop at 11200 iteration and i can't add harmonic motion 
thanks in advance.

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