I use the simple isotropic elastic material in uniaxial tension as an example. First the isotropic elastic behavior is directly coded inside the UMAT subroutine and it works fine below (OK.for).