Continuum Constitutive Modeling for Isotropic Hyperelastic Materials

The partial differential equation for isotropic hyperelastic constitutive models has been postulated and derived from the balance between stored energy and stress work done. The partial differential equation as a function of three invariants has then been solved by Lie group methods. With geometric meanings of deformations, the general solution boils down to a particular threeterm solution. The particular solution has been applied for several isotropic hyperelastic materials. For incompressible materials, vulcanized rubber containing 8% sulfur and Entec Enflex S4035A thermoplastic elastomer, three coefficients have been determined from uniaxial tension data and applied to predict the pure shear and equibiaxial tension modes. For a slightly compressible rubber material, the coefficients have also been extracted from the confined volumetric test data.


Introduction
Many materials with wide-ranging applications, including rubbers, thermoplastic elastomers, and biological tissues, can be characterized as isotropic hyperelastic materials.Constitutive models play an important role in design and analyses.The established models have been classified into two main categories: micro-mechanical and phenomenological.Micro-mechanical models include 3-chain, 8-chain, Unit sphere, Full-network, Flory-Erman, Tube, and Extended tube.Phenomenological models can be subdivided into two groups: stretch-based and invariant-based.Stretch-based phenomenological models cover Ogden, Shariff, and Attard models.Invariant-based phenomenological models include Neo-Hooke, Mooney-Rivlin, Isihara, Gent-Thomas, Swanson, Yeoh, Arruda-Boyce, Gent, Yeoh-Fleming, Carroll, Hart-Smith, Alexander, van der Waals, Pucci-Saccomandi, and Lopez-Pamies models.The established models have been thoroughly reviewed and evaluated by Boyce and Arruda (2000) [1], Steinmann, Hossain, and Possart (2012) [2], and Hossain and Steinmann (2013) [3].Current multiscale approaches for modeling rubber-like materials and soft tissues have been briefly reviewed and compared with classical macroscopic phenomenological models by Puglisi and Saccomandi (2016) [4].
For constitutive models, great progress has been made in continuum mechanics.The thermodynamic foundations of constitutive equations at finite strains are archived in the classic treatise on continuum mechanics by Truesdell and Noll (1965) [5].Advanced mathematical foundations of elasticity are presented in monographs by Marsden and Hughes (1983) [6], Ogden (1984) [7], and Holzapfel (2000) [8].
For isothermal and adiabatic processes, the stored energy function, Ψ, is a function of the right Cauchy-Green tensor, C, as ( ), Ψ = Ψ C (1) in which T = C F F , and F is the deformation gradient tensor.With Clausius-Duhem inequality, the second Piola-Kirchhoff stress tensor, S, takes the form of ( ) When hyperelastic materials under deformation are isotropic, the stored energy function depends on the invariants of the right Cauchy-Green tensor, which was established by the classical work of Rivlin (1948) [9] ( ) , , , in which the three invariants, I 1 , I 2 , and I 3 , of right Cauchy-Green tensor are given by ( ) ( ) where ( ) are principal stretches in three mutually orthogonal directions.Substituting (3) into (2) and using the chain rule of differentiation yields ( ) The derivatives of the three invariants with respect to C are given by the following three equations, respectively , Substituting the derivatives (8), (9), and (10) into constitutive Equation (7) produces the most general form of a stress as a function of the three invariants for isotropic hyperelastic materials under finite strains For a tangent operator, a fourth-order elasticity tensor,  , is usually introduced to relate changes in stress dS and changes in the right Cauchy-Green tensor dC through where the fourth-order elasticity tensor is written as ( ) .
The most general fourth-order elasticity tensor in terms of the three principal invariants of the right Cauchy-Green tensor is documented by Holzapfel (2000) [8] ( ) ) where the fourth-order tensors ∂ ∂ C are given by , , ( ) and where  and  are fourth-order zero and unit tensors, and the parameters 1 2

8
, , ,  Although the second Piola-Kirchhoff stress tensor (11) and the fourth-order elasticity tensor (13) are the fundamental constitutive equations in the theory of finite hyperelasticity, the stored energy function, Ψ, is the key to determining all of them.Thus, the main objective is to establish, solve, and apply the stored energy function to construct a general constitutive model for isotropic hyperelastic materials.

Constitutive Modeling for Isotropic Hyperelastic Materials
For a general isotropic hyperelastic constitutive model based on the theory of continuum mechanics, the partial differential equation of a stored energy function has been postulated and derived.Once the partial differential equation is solved and the function is determined, the constitutive model can be obtained under several geometric conditions.

Stored Energy Partial Differential Equation
A stored energy at a finite deformation is postulated to be equal to its stress work done.The invariant of the right Cauchy-Green tensor-based stored energy function can be formulated by : . 2 Substituting ( 11) into (21) generates Eliminating the right Cauchy-Green tensor related terms with (4) 1 , (5) 1 , noticing , and rearranging, yields the following partial differential equation Now, the partial differential equation for the stored energy function (24) has been simplified into a function of three scalar invariants.The stored energy partial differential equation can generally be solved by Lie group methods.

General Solution for Stored Energy Function
Based on Lie group methods documented by Ibragimov (2009) [10], the characteristic system for the stored energy partial differential Equation (24) reads Taking its three first-integrals

I I ψ =
, and 3 1 I ψ = Ψ , the general solution to (24) can be solved and written as ( ) ( ) where f and g are two arbitrary functions.

Model for Isotropic Hyperelastic Materials
The general solution (26) has two arbitrary functions to be determined for practical applications.For isotropic hyperelastic materials, translational, rotational, and ellipsoidal deformations are indispensable.Thus, the curvatures of the three types of deformations are used as geometric conditions to select the two arbitrary functions.

A Particular Solution for Stored Energy Function
For the accurate representation of a translational deformation and a rotational deformation, the first arbitrary function, f, is selected as In order to accurately describe an ellipsoidal deformation, a simple inverse function is required.Thus, the second arbitrary function, g, is chosen as ( ) ( ) Substituting ( 27) and ( 28) into (26) gives the particular solution of the stored energy function where the three coefficients, c 1 , c 2 , and c 3 , are unknown constants to be determined by experimental tests.

Stress and Elasticity Tensors
Stress and elasticity tensors are indispensable elements in constitutive models for numerical solutions of boun-dary value problems by iterative approaches such as Newton-Raphson methods.
Applying the general stress Equation (11) to the particular solution (29) gives the specific second Piola-Kirchhoff stress tensor The eight parameters for the elasticity tensor have been worked out by substituting the particular solution (29) into (17) through (20) Seven out of eight parameters are non-zeros and therefore the particular solution is relatively general.

Geometric Interpretations of Particular Solution
Geometric interpretations of the three terms in the particular solution can also be revealed through the second Piola-Kirchhoff stress tensor and the elasticity tensor.The term with the coefficient c 1 only appears in the second Piola-Kirchhoff stress tensor Equation (30), describing a constant modulus but does not show up in the parameters for the elasticity tensor, meaning that the curvature of a deformation is zero due to principal translational deformations.The term related to the coefficient c 2 models the curvatures of rotational deformations.The term involving the coefficient c 3 captures the curvature of an ellipsoidal deformation, coupling a volumetric or dilatational deformation and an isochoric or distortional deformation.Specific equations for the mean curvature of a general spacial surface and the Gaussian curvature of an ellipsoid can be found in the classic monograph of differential geometry by Carmo (1976) [11].Nevertheless, ellipsoidal deformations preserve numerical stabilities.
The particular solution as a constitutive model is applicable to both incompressible and compressible isotropic hyperelastic materials.The current model will be used to fit the experimental data of vulcanized rubber containing 8% sulfur by Treloar (1944) [12], Entec Enflex S4035A thermoplastic elastomer (TPE), and a slightly compressible rubber material.

Model Validation
Several experimental test methods can be used to characterize hyperelastic materials.Most commonly used methods are uniaxial tension, uniaxial compression, equibiaxial tension, pure shear, and confined volumetric tests in which detailed reviews are given by Charlton, Yang, and Teh (1994) [13].In experiments, nominal stresses and stretches are usually recorded.The three principal nominal stresses in terms of a stored energy function and principal stretches are given by , 1, 2, 3.

Modeling for Different Deformation Modes
In order to validate the continuum constitutive model (29), the nominal stress, the first Piola-Kirchhoff stress, as the function of principal stretches of uniaxial tension, pure shear, equibiaxial tension, and confined volumetric tests will, respectively, be derived for the continuum model.

Uniaxial Tension Test
In a uniaxial tension test, we have 1  .

Fitting Continuum Model to Experimental Data
First, Treloar's uniaxial tension experimental data has been chosen to fit the uniaxial tension model (36) by a linear least square method.The total squared error between the model and experiment is given by ( ) ( ) where e u P is experimental engineering stress data in uniaxial tension tests and n is the total number of pairs of experimental data.The error,  , can be minimized with the following conditions 0, 1, 2, 3.

Discussion
The most important concept to notice about hyperelastic materials is that its deformation is not directly proportional to the applied load in which it exhibits a nonlinear behavior.For hyperelastic materials, independent load paths are usually assumed.This behavior, however, is only approximately observed in many hyperelastic materials in Belytschko, Liu, and Moran (2000) [14].The independent load path is only applicable to linear elastic materials.In the continuum constitutive modeling of isotropic hyperelastic materials (21), independent load paths are irrelevant.
The invariant of stain-based stored energy function can be postulated and formulated by a double contraction of the second Piola-Kirchhoff stress tensor and the Green-Lagrange strain tensor as ( ) ( ) ( ) Equation ( 49) has also been solved by Lie group methods, which results in where q and r are also two arbitrary functions.When the three invariants of Green-Lagrange strain tensor, ( ) .
The five strain measures were unified in terms of stretches with ( ) , in which i = 1, 2, 3, by Seth (1964) [15].In the cases of n = −2, −1, 0, 1, and 2, it gives Almansi, Swainger, Hencky, Cauchy, and Green-Lagrange strains, respectively.A stretch is ubiquitous in the strain measures.Compared to strains, not only do stretches capture the characteristics of a deformation but they also avoid singularities and imaginary numbers in constitutive modeling.Thus, the invariant of stretch-based solution is more appropriate in establishing a constitutive model for isotropic hyperelastic materials.
Deformations of hyperelastic materials in experimental tests could be coupled with boundary, geometric, and material nonlinearities.In uniaxial compression tests, experimental results are usually overestimated due to inhomogeneous deformations and hydrostatic pressure effects [16].The inhomogeneous deformations and hydrostatic pressure effects are mainly resulted from boundary nonlinearities in which a friction force is exerted on contact surfaces of samples.In equibiaxial tension tests, albeit no friction, inhomogeneous deformations and hydrostatic pressure effects under a limited amount of deformation do exist, causing overestimated experimental results such as the equibiaxial tension test of the S4035A TPE shown in Figure 6.In general, equibiaxial tension tests cannot match the accuracy and the amount of a deformation in uniaxial tension tests.
In the general solution (26) of stored energy function, Ψ, there are two functions, f and g, that look relatively arbitrary.What is not arbitrary, however, is the framework of the solution.It is worth noticing that all terms are in the order of λ 2 even though they are already dimensionless.The continuum constitutive model (29) has been derived from the general solution of the stored energy partial differential equation, which is suitable for both incompressible and compressible isotropic hyperelastic materials.The continuum model fits various modes of deformation using the data obtained from a uniaxial tension test only for the tested incompressible hyperelastic materials.This leads to a reduced number of requirements on material testing.

Conclusions
The partial differential equation for isotropic hyperelastic constitutive models has been postulated and derived from the balance between stored energy and stress work done.The stored energy as a function of three invariants of the right Cauchy-Green tensor has then been solved by Lie group methods.

1 I
λ λ λ = = .Thus, the first two invariants for the uniaxial tension deformation mode are Thus, the first two invariants for the pure shear deformation mode are the same λ = = .Thus, the first two invariants for equibiaxial tension deformation mode are

Figure 4 .
Figure 4. Comparison between model and uniaxial tension test of S4035A TPE.

Figure 5 .
Figure 5. Model prediction for pure shear test of S4035A TPE.

Figure 7 .
Figure 7.Comparison between model and confined volumetric test.

c 2 =
32856.28089MPa, and c 3 = 126.40171MPa have been obtained for the confined volumetric test model (42) and the most complete form of the stored energy function (29).

3 I
E , are used, Equation (50) takes the same form of Equation (26) in