ANSYS Anisotropic Hyperelastic: Exponential

Constitutive Theory

Some materials are inherently anisotropic due to their microstructure or added reinforcements. Examples of this include fiber reinforced rubbers and thermoplastics, oriented thermoplastics, and biomaterials. The direction-dependent behaviors of these materials can be quite important for their overall mechanical response. If you have not modeled these materials before, the good news is that it is almost as easy to model an anisotropic material as it is to model an isotropic material. The main difference is the amount of experimental testing that is need to extract the response and guide the material model calibration.

ANSYS Mechanical supports the following exponential-function-based anisotropic hyperelastic material model energy function:

\[ W = \sum_{i=1}^3 a_i (\overline{I}_1 – 3)^i + \sum_{j=1}^3 b_j (\overline{I}_2 – 3)^j  + \frac{c_1}{2c_2}\left(\exp\left[c_2(\overline{I}_4-1)^2\right]-1\right) \\ + \frac{e_1}{2e_2}\left(\exp\left[e_2(\overline{I}_6-1)^2\right]-1\right) + \frac{1}{d} (J-1)^2 \]

When working with anisotropic polymers it is often difficult to separate out the \(I_2\) dependence from the fiber response, so I prefer to simply remove the second term giving the following simplified equation for the Helmholtz free energy. This equation is implemented in MCalibration:

\[ W = \sum_{i=1}^3 a_i (\overline{I}_1 – 3)^i + \frac{c_1}{2c_2}\left(\exp\left[c_2(\overline{I}_4-1)^2\right]-1\right) \\ + \frac{e_1}{2e_2}\left(\exp\left[e_2(\overline{I}_6-1)^2\right]-1\right) + \frac{1}{d} (J-1)^2 \]

This model consists of a Yeoh isotropic hyperelastic matrix material with independent fibers oriented in two user-defined orientations. The \(c_i\) and \(e_i\) parameters specify the fiber stiffnesses.


The following screenshot shows the ANSYS anisotropic hyperelastic polynomial material model evaluated in MCalibration.

Here are the final material models in ANSYS ADPL dat-file format:

TB, AHYPER, matid, 1, 10, EXP
TBDATA,  1, 1 ! A1
TBDATA,  2, 0 ! A2
TBDATA,  3, 0 ! A3
TBDATA,  4, 0.0 ! B1
TBDATA,  5, 0.0 ! B2
TBDATA,  6, 0.0 ! B3
TBDATA,  7, 1 ! C1
TBDATA,  8, 1 ! C2
TBDATA,  9, 0.1 ! E1
TBDATA, 10, 0 ! E2

TB, AHYPER, matid, 1, 1, PVOL
TBDATA, 1, 0.01 ! D

TB, AHYPER, matid, 1, 3, AVEC
TBDATA, 1, 1
TBDATA, 2, 0
TBDATA, 3, 0

TB, AHYPER, matid, 1, 3, BVEC
TBDATA, 1, 0
TBDATA, 2, 1
TBDATA, 3, 0
MP, DENS, matid, 1e-09

You can, of course, also calibrate the ANSYS anisotropic material model to experimental data. Since the material model has been implemented as a built-in feature in MCalibration, the parameter calibration will run extremely fast and will not use ANSYS or any ANSYS license tokens.

The following article shows how you can use the calibrated material model in ANSYS Workbench.

Share on facebook
Share on twitter
Share on linkedin

More to explore

What Is A Material Model?

I often get asked to define the terms “Material Model” and “Material Parameter”. This short article tries to explain what these terms mean.

Leave a Comment