BLUF: I want to define an abaqus load step that I fully expect to cause plastic deformation in a model, but to abort it as soon as the plastic energy (output, ElPD) is measurably greater than zero, and report that load
I am attempting to verify a complex user-generated vumat on a single element (also, in parallel trying to apply the loads to a cylinder, but I can't get that sucker to run). After plastically deforming the material, the constitutive model is expected to take into account preferential orientation of the metallic grains. I would like, after a significant plastic deformation, to sketch the shape of the yield surface on the yield plane by probing it in several directions. I cannot just willy-nilly ramp to huge values in those directions and extract the yield in post-processing, because large plastic strains will change the behavior of the material. Aside from running a series of increasingly larger load steps, and requesting restart file after each, is there a more elegant solution?