I use an open source DEM platform LIGGGHTS-flexible fibres, which includes elastic bonds and strength based failure criteria. Is it possible to include Elastoplastic bond model with cohesive zone failure.

One way is to modify the source C++ code, but a simpler method will be highly appreciated

Thanks