Products: Abaqus/Standard Abaqus/Explicit
The constitutive behavior of hyperelastic materials is discussed in “Hyperelastic material behavior,” Section 4.6.1, in the context of isotropic response. However, many materials of industrial and technological interest exhibit anisotropic elastic behavior due to the presence of preferred directions in their microstructure. Examples of such materials include common engineering materials (such as fiber-reinforced composites, reinforced rubber, and wood) as well as soft biological tissues (such as those found in arterial walls and heart tissues). Under large deformations these materials exhibit highly anisotropic and nonlinear elastic behavior due to rearrangements in their microstructure, such as reorientation of the fiber directions with deformation. The simulation of these nonlinear effects requires constitutive models formulated within the framework of anisotropic hyperelasticity.
Hyperelastic materials are described in terms of a “strain energy potential,” which defines the strain energy stored in the material per unit of reference volume (volume in the initial configuration) as a function of the deformation at that point in the material. Two distinct formulations are used for the representation of the strain energy potential of anisotropic hyperelastic materials: strain-based and invariant-based.
In this case the strain energy function is expressed directly in terms of the components of a suitable strain tensor, such as the Green strain tensor (see “Strain measures,” Section 1.4.2):
is Green's strain, 
is the right Cauchy-Green strain tensor, 
is the deformation gradient, and 
is the identity matrix. Without loss of generality, the strain energy function can be written in the form where 
is the modified Green strain tensor, 
is the distortional part of the right Cauchy-Green strain, and 
is the volume change.The underlying assumption in models based on the strain-based formulation is that the preferred materials directions are initially aligned with an orthogonal coordinate system in the reference (stress-free) configuration. These directions may become nonorthogonal only after deformation. Examples of this form of strain energy function include the generalized Fung-type form (see “Generalized Fung form,” below).
From Equation 4.6.3–1 the variation of 
is given as
For a compressible material the strain variations are arbitrary, so this equation defines the stress components for such a material as
When the material response is almost incompressible, the pure displacement formulation, in which the strain invariants are computed from the kinematic variables of the finite element model, can behave poorly. One difficulty is that from a numerical point of view the stiffness matrix is almost singular because the effective bulk modulus of the material is so large compared to its effective shear modulus, thus causing difficulties with the solution of the discretized equilibrium equations. Similarly, in Abaqus/Explicit the high bulk modulus increases the dilatational wave speed, thus reducing the stable time increment substantially. To avoid such problems, Abaqus/Standard offers a “mixed” formulation for such cases (refer to “Hyperelastic material behavior,” Section 4.6.1).
Using the continuum theory of fiber-reinforced composites (Spencer, 1984), the strain energy function can be expressed directly in terms of the invariants of the deformation tensor and fiber directions. For example, consider a composite material that consists of an isotropic hyperelastic matrix reinforced with 
families of fibers. The directions of the fibers in the reference configuration are characterized by a set of unit vectors 
, (
). Assuming that the strain energy depends not only on deformation, but also on the fiber directions, the following form is postulated:
and the vectors : where 
and 
are the first and second deviatoric strain invariants; 
is the volume ratio (or third strain invariant); and 
and 
are the pseudo-invariants of 
, , and 
, defined as
are geometrical constants (independent of deformation) equal to the cosine of the angle between the directions of any two families of fibers in the reference configuration,
Unlike in the case of the strain-based formulation, in the invariant-based formulation the fiber directions need not be orthogonal in the initial configuration. An example of the invariant-based energy function is the form proposed by Holzapfel, Gasser, and Ogden (2000) for arterial walls (see “Holzapfel-Gasser-Ogden form,” below).
From Equation 4.6.3–4 the variation of is given as
Using the principle of virtual work (Equation 4.6.3–3) and after some lengthy derivations, the stress components for a compressible material are found to be given as
and 
. Several particular forms of the strain energy potential are available in Abaqus. The incompressible or almost incompressible models make up:
the generalized Fung form and
the Holzapfel-Gasser-Ogden form.
The generalized Fung strain energy potential in Abaqus is based on the two-dimensional exponential form proposed by Fung et al. (1979), suitably generalized to arbitrary three-dimensional states following Humphrey (1995). It has the following form:
is the elastic volume ratio, and 
is defined as 
is a dimensionless symmetric fourth-order tensor of anisotropic material constants that can be temperature dependent and 
are the components of the modified Green strain tensor.The elastic volume ratio, , relates the total volume ratio, J, and the thermal volume ratio, 
:
is given by 
are the principal thermal expansion strains that are obtained from the temperature and the thermal expansion coefficients. The initial deviatoric elasticity tensor, 
, and bulk modulus, 
, are given by
Abaqus supports two forms of the generalized Fung model: anisotropic and orthotropic. The number of independent components that must be specified depends on the level of anisotropy of the material: 21 for the fully anisotropic case and 9 for the orthotropic case.
The form of the strain energy potential is based on that proposed by Holzapfel, Gasser, and Ogden (2000) and Gasser, Ogden, and Holzapfel (2006) for modeling arterial layers with distributed collagen fiber orientations:
, D, 
, 
, and 
are temperature-dependent material parameters; is the number of families of fibers (
); is the first invariant of ; is the elastic volume ratio, as defined previously; and 
are pseudo-invariants of and . The model assumes that the directions of the collagen fibers within each family are dispersed (with rotational symmetry) about a mean preferred direction. The parameter (
) describes the level of dispersion in the fiber directions. If 
is the orientation density function that characterizes the distribution (it represents the normalized number of fibers with orientations in the range 
with respect to the mean direction), the parameter is defined as
the fibers are perfectly aligned (no dispersion). When 
the fibers are randomly distributed and the material becomes isotropic; this corresponds to a spherical orientation density function.The strain-like quantity 
characterizes the deformation of the family of fibers with mean direction . 
for perfectly aligned fibers (
), and 
for randomly distributed fibers (
).
The first two terms in the expression of the strain energy function represent the distortional and volumetric contributions of the noncollagenous isotropic ground material; and the third term represents the contributions from the different families of collagen fibers, taking into account the effects of dispersion. A basic assumption of the model is that collagen fibers can support only tension, as they would buckle under compressive loading. Thus, the anisotropic contribution in the strain energy function appears only when the strain of the fibers is positive or, equivalently, when 
. This condition is enforced by the term 
, where the operator 
stands for the Macauley bracket and is defined as 
.
The initial deviatoric elasticity tensor, , and bulk modulus, , are given by
is the fourth-order unit tensor and 
is the Heaviside unit step function.Abaqus also allows other forms of strain-based energy potentials to be defined via user subroutines UANISOHYPER_STRAIN in Abaqus/Standard and VUANISOHYPER_STRAIN in Abaqus/Explicit by programming the first and second derivatives of the strain energy potential with respect to the components of the modified Green strain and the elastic volume ratio, .
Abaqus also allows other forms of invariant-based energy potentials to be defined via user subroutines UANISOHYPER_INV in Abaqus/Standard and VUANISOHYPER_INV in Abaqus/Explicit by programming the first and second derivatives of the strain energy potential with respect to each invariant.
