There is a field variable:

DG All components of the total deformation gradient. Available only for hyperelasticity, hyperfoam, and material models defined in user subroutine UMAT. For fully integrated first-order quadrilaterals and hexahedra, the selectively reduced integration technique is used. A modified deformation gradient is output for these elements.

As you can see, usage is limited.

You can also try doing it by user subroutine UVARM, but I have no experience here.

Thanks Stanisław Burzyński . I tried UVARM, which doesn't work for my mesh, S3R element.

Are you using user-materials? UMAT and VUMAT have access to DFGRD and defgrad, respectively.