Tracking the Crack

You can use the jcfpm argument in the vtkrecorder:

VTKRecorder(recorders=['spheres',intr','jcfpm'])

This recorder tracks jcfpm specific interaction characteristics, such as "isCohesive" which will show the bonds that are no longer cohesive in paraview.

You can also use the 'recordCracks=True' argument for the jcpfm law:

Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(recordCracks=True)

which will allow you to track the type of cracks (tensile or shear). I've had mixed luck with visualizing that vtu in paraview, so I find it easier to use python to plot the cracks.txt file that is also generated.

Hi Robert,

Thanks so much for your recommendation, can you please explain what you me exactly by:

 " I've had mixed luck with visualizing that vtu in paraview, so I find it easier to use python to plot the cracks.txt file that is also generated."

Thanks,
Seti

This question was expired because it remained in the 'Open' state without activity for the last 15 days.