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.
___
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
