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,

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.