I checked failel without GISSMO. If (hsv(2).gt.1.) then failel = .true. The results are even more strange. Right now “failure” (if I can call it like that) occurs when hsv(2) = 0.97 and the stress distribution becomes not uniform. And in the same time, the specimen remains solid like particles are still bonded and there is no material spalling.
You are navigating away from the AIS Discovery experience