yade-users team mailing list archive

[Tina Asia <question630281@xxxxxxxxxxxxxxxxxxxxx>]

New question #630281 on Yade:
https://answers.launchpad.net/yade/+question/630281

Hi,

I have created a spherical packing using randomDensePack without giving spheresInCell, but when I using this packing in my simulation, this packing explodes. 

I add i.phys.unp=i.geom.penetrationDepth, but this method also failed. 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.dt=utils.PWaveTimeStep()
O.usesTimeStepper=True

qt.Controller()
qt.View()

O.run(1,True)
for i in O.interactions:
          i.phys.unp=i.geom.penetrationDepth

BTW: when i using regularHexa packing, there is no error.
Thanks in advances.

Tina


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