Isotropic hyperelasticity in principal stretches: explicit elasticity tensors and numerical implementation

Research output: Contribution to journalArticlepeer-review

3 Citations (Scopus)
124 Downloads (Pure)


Elasticity tensors for isotropic hyperelasticity in principal stretches are formulated and implemented for the Finite Element Method. Hyperelastic constitutive models defined by this strain measure are known to accurately model the response of rubber, and similar materials. These models may not be available in the library of a Finite Element Analysis software, but a numerical implementation of the constitutive model may be provided by a programmed subroutine. The implementation proposed here is robust and accurate, with straightforward user input. It is presented in multiple configurations with novel features, including efficient definition of isochoric stress and elasticity coefficients, symmetric dyadic products of the principal directions, and development of a stable and accurate algorithm for equal and similar principal stretches. The proposed implementation is validated, for unique, equal and similar principal stretches. Further validation in the Finite Element Method demonstrates the developed implementation requires lower computational effort compared with an alternative, well-known implementation.
Original languageEnglish
Pages (from-to)1273-1288
Number of pages16
JournalComputational Mechanics
Issue number5
Early online date3 May 2019
Publication statusPublished - 30 Nov 2019


  • hyperelasticity
  • principal stretches
  • finite element method
  • numerical implementation
  • elasticity tensors

Cite this