yade-users team mailing list archive

[Jingchi Yu <question699459@xxxxxxxxxxxxxxxxxxxxx>]

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

Dear all,

I tried to use I tried to use mat to simulate the connection behavior of 400 particles, but when running I found that the calculation could not reach the goal and  finished (I believe there was an error), and the particles bounced off each other at the beginning of the run instead of coming together quickly (I believe the bond was not successfully connected). Attached is my code. What are the possible reasons for this, please? 

code:
O.materials.append(CohFrictMat(young=6e8,poisson=.8,density=2650,isCohesive=True,frictionAngle=26.565,normalCohesion=1e5,shearCohesion=1e5))#density=2e3

sp = pack.SpherePack()
size = .24
sp.makeCloud(minCorner=(0,0,.05),maxCorner=(size,size,.05),rMean=.005,rRelFuzz=.4,num=400,periodic=True,seed=1)
sp.toSimulation()
O.cell.hSize = Matrix3(size,0,0, 0,size,0, 0,0,.1)
print(len(O.bodies))
for p in O.bodies:
   p.state.blockedDOFs = 'zXY'
   p.state.mass = 2650 * 0.1 * pi * p.shape.radius**2 # 0.1 = thickness of cylindrical particle
   inertia = 0.5 * p.state.mass * p.shape.radius**2
   p.state.inertia = (.5*inertia,.5*inertia,inertia)

O.dt = utils.PWaveTimeStep()
print(O.dt)

O.engines = [
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb()]),
   InteractionLoop(
      [Ig2_Sphere_Sphere_ScGeom6D()],
      [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys()],
      [Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]
   ),
   PeriTriaxController(
      dynCell=True,
      goal=(-1.e5,-1.e5,0),
      stressMask=3,
      relStressTol=.001,
      maxUnbalanced=.001,
      maxStrainRate=(.5,.5,.0),
      doneHook='term()',
      label='biax'
   ),
   NewtonIntegrator(damping=.1)

]

Thank you all,
YU

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