Understanding the deformation gradient in Abaqus and key guidelines for anisotropic hyperelastic user material subroutines (UMATs)

This tutorial paper provides a step-by-step guide to developing a comprehensive understanding of the different forms of the deformation gradient used in Abaqus, and outlines a number of key issues that must be considered when developing an Abaqus user defined material subroutine (UMAT) in which the Cauchy stress is computed from the deformation gradient. Firstly, we examine the “classical” forms of global and local deformation gradients. We then show that Abaqus/Standard does not use the classical form of the local deformation gradient when continuum elements are used, and we highlight the important implications for UMAT development. We outline the key steps that must be implemented in developing an anisotropic fibre-reinforced hyperelastic UMAT for use with continuum elements and local orientation systems. We also demonstrate that a classical local deformation gradient is provided by Abaqus/Standard if structural (shell and membrane) elements are used, and by Abaqus/Explicit for all element types. We emphasise, however, that the majority of biomechanical simulations rely on the use of continuum elements with a local coordinate system in Abaqus/Standard, and therefore the development of a hyperelastic UMAT requires an in-depth and precise understanding of the form of the non-classical deformation gradient provided as input by Abaqus. Several worked examples and case studies are provided for each section, so that the details and implications of the form of the deformation gradient can be fully understood. For each worked example in this tutorial paper the source files and code (Abaqus input files, UMATs, and Matlab script files) are provided, allowing the reader to efficiently explore the implications of the form of the deformation gradient in the development of a UMAT.