Hyperelastic Material Modeling Using First- and Second-Order Differential Operations with Complex Step Derivative Approximation

This study proposes a method for calculating hyperelastic materials’ stress and tangent modulus using a generalized complex-step derivative approximation. This method is accurate and requires only the operational code of the strain energy density function. Moreover, because it uses only standard complex operations in Fortran, no prior preparation is required, making the method convenient and practical. The effectiveness of the method was confirmed by applying it to the implementation of isotropic and anisotropic hyperelastic models.