yade-users team mailing list archive

[Luc Scholtès <question697724@xxxxxxxxxxxxxxxxxxxxx>]

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

Luc Scholtès proposed the following answer:
Hello,

You can choose either option in the jcfpm model. Just change the boolean
smoothJoint from True to False:

https://yade-
dem.org/doc/yade.wrapper.html?highlight=smoothjoint#yade.wrapper.Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM.smoothJoint

if smothJoint=False: regular contacts geometry and mechanical properties (you can modify the properties as you want, e.g., isCohesive=0, FnMax=FsMAx=0).
if smoothJoint=True: the contact geometry will be modified according to the fracture plane and the properties modified according to the joint properties of the JCFPMmat, i.e., jointNormalStiffness, jointShearStiffness, jointFrictionAngle, jointTensileStrength, jointCohesion (you might have errors/strange behaviors if you don't define those since their default values are 0).

The identification of the contacts and particles along the pre-existing
crack/fracture  is done when you type execfile('identifBis.py'). After
that, just loop over the particles and interactions and check their
status: those which are along the fracture plane will be identified as
onJoint:

for particles: https://yade-dem.org/doc/yade.wrapper.html?highlight=onjoint#yade.wrapper.JCFpmState.onJoint
for interactions: https://yade-dem.org/doc/yade.wrapper.html?highlight=onjoint#yade.wrapper.JCFpmPhys.isOnJoint

Luc

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