Процесс разрушения

Здравствуйте, подскажите пожалуйста, как можно смоделировать процесс развития трещины через usermat? Возможно ли сделать так "при достижении  определённой переменной в элементе, обнулять жёсткость элемента?  Метод 0.0001*dsdeEl -> stress=dsdeEl*strainEl , приводит к ошибки, что нужно добавить или изменить? Спасибо