Higher-order and higher floating-point precision numerical approximations of finite strain elasticity moduli

Stephen John Connolly, Donald MacKenzie, Yevgen Gorash

Research output: Contribution to journalArticlepeer-review

3 Citations (Scopus)
111 Downloads (Pure)


Two real‐domain numerical approximation methods for accurate computation of finite strain elasticity moduli are developed and their accuracy and computational efficiency investigated, with reference to hyperelastic constitutive models with known analytical solutions. The methods are higher‐order and higher floating‐point precision numerical approximation, the latter being novel in this context. A general formula for higher‐order approximation finite difference schemes is derived and a new procedure is proposed to implement increased floating‐point precision. The accuracy of the approximated elasticity moduli is investigated numerically using higher‐order approximations in standard double precision and increased quadruple precision. It is found that as the order of the approximation increases, the elasticity moduli tend towards the analytical solution. Using higher floating‐point precision, the approximated elasticity moduli for all orders of approximation are found to be more accurate than the standard double precision evaluation of the analytical moduli. Application of the techniques to a finite element problem shows that the numerically approximated methods obtain convergence equivalent to the analytical method but require greater computational effort. It is concluded that numerical approximation of elasticity moduli is a powerful and effective means of implementing advanced constitutive models in the finite element method without prior derivation of difficult analytical solutions.
Original languageEnglish
Pages (from-to)1184-1201
Number of pages18
JournalInternational Journal for Numerical Methods in Engineering
Issue number10
Early online date25 Jul 2019
Publication statusPublished - 7 Dec 2019


  • elasticity moduli
  • numerical differentiation
  • higher-order approximation
  • higher floating-point precision
  • hyperelasticity
  • nonlinear finite element method


Dive into the research topics of 'Higher-order and higher floating-point precision numerical approximations of finite strain elasticity moduli'. Together they form a unique fingerprint.

Cite this