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,
Absolutely. Here is my MWE. Thanks.
Best
Jiannan Wang
#############
from builtins import range
from yade import pack,utils
import math
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=SpherePack()
sp=pack.randomDensePack(pack.inAlignedBox(mn,mx),spheresInCell=2000,radius=1/20.,returnSpherePack=True)
O.bodies.append([sphere(s[0],s[1],color=(0.6+0.15*rand(),0.5+0.15*rand(),0.15+0.15*rand()),material='mat1')
for s in sp])
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=intR)],
[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),
]
O.engines=O.engines+[VTKRecorder(Key='identifier',label='vtk',iterPeriod=iterper,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.
_______________________________________________
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
