yade-users team mailing list archive

[Tina Asia <question628632@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Answered => Open

Tina Asia is still having a problem:
Hi Tina,

Thanks for your reply, I tried randomDensePack withous gives
spheresInCell, but I got a very unstable simulation, once it run, all
spheres fly around. The simulation just gave some warnings without
errors. Here is my code:

from yade import pack
from yade import qt,utils

stone=CohFrictMat(young=5.0e10,poisson=0.25,density=2640,frictionAngle=radians(18),isCohesive=True,normalCohesion=4.5e6,shearCohesion=4.5e7,momentRotationLaw=True)
O.materials.append(stone)
steel=CohFrictMat(young=3.06e11,poisson=0.29,density=7861,frictionAngle=0.545,normalCohesion=0,shearCohesion=0)
O.materials.append(steel)

pred=pack.inSphere((0,0,0.5),0.25)
#sphs=pack.regularHexa(pred,radius=0.01,gap=0,material=stone)
sphs=pack.randomDensePack(pred,radius=0.025,rRelFuzz=0,material=stone)
O.bodies.append(sphs)

for i in sphs:
	#global velocity
	velocity=i.state.vel=(0,0,-25)

O.bodies.append(geom.facetBox((0,0,0.2),(0.3,0.3,0.003),material=steel))

O.engines=[
    ForceResetter(),
    InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=1.5),Bo1_Facet_Aabb()]),
    InteractionLoop(
        [Ig2_Sphere_Sphere_ScGeom6D(interactionDetectionFactor=1.5),Ig2_Facet_Sphere_ScGeom6D()],
        [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True)],
        [Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]
        ),
    VTKRecorder(fileName='post/0deg-',recorders=['all'],iterPeriod=50),
    GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=5,timestepSafetyCoefficient=0.8,defaultDt=PWaveTimeStep()),
    NewtonIntegrator(damping=0.3,gravity=(0,0,-9.81)),
    
]

O.trackEnergy=True
O.step()

O.dt=utils.PWaveTimeStep()
O.usesTimeStepper=True

qt.Controller()
qt.View()

But when I tried regularHexa packing, I got a very  stable simulation.


Thanks in advance.

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