Hi Asarma,
I'm not very familiar with CZM but wanted to point out that in case you decided to pursue with Ansys LS-DYNA (as you indicated), this type of application should be handled relatively easy using the MAT_ADD_EROSION card that is also available to use via Workbench LS-DYNA.
Just a suggestion.