yade-users team mailing list archive

[Jiannan Wang <question690820@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Needs information => Open

Jiannan Wang gave more information on the question:
Hello Robert,

Oops... Sorry about the missing variable, again. The following MWE runs
as expected. Thank you.

Best
Jiannan Wang

######################3
from builtins import range
from yade import pack,utils
import math
from pylab import rand

O.materials.append(JCFpmMat(type=1,density=2500,young=1800,poisson=0.1,frictionAngle=radians(7),tensileStrength=100,cohesion=1000,label='mat'))

mn,mx=Vector3(0,0,0),Vector3(1,2,1)
walls=aabbWalls([mn,mx],thickness=0,material='mat')
wallIds=O.bodies.append(walls)
sps=pack.SpherePack()
sp=pack.randomDensePack(pack.inAlignedBox(mn,mx),spheresInCell=2000,radius=1/20,returnSpherePack=True)
sp.toSimulation(material='mat')

O.engines=[
    ForceResetter(),
    InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=1)]),
    InteractionLoop(
      [Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=1)],
     [Ip2_JCFpmMat_JCFpmMat_JCFpmPhys(cohesiveTresholdIteration=1)],
     [Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(smoothJoint=True,label='interactionLaw', recordCracks=True, recordMoments=True,Key='identifier',momentRadiusFactor=5,neverErase=True),Law2_ScGeom_FrictPhys_CundallStrack()]
    ),
    GlobalStiffnessTimeStepper(active=1,timestepSafetyCoefficient=0.5),
    NewtonIntegrator(damping=0.1),
]

from yade import qt
v=qt.View()

O.engines=O.engines+[VTKRecorder(Key='identifier',label='vtk',iterPeriod=100,initRun=True,fileName='output-',recorders=['moments','spheres','intr','stress','velocity','bstresses','jcfpm',
'cracks'])]

######################

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