MATHE

Bulk Data Entry Defines material properties for nonlinear hyperelastic materials. The Polynomial form is available and various material types 3 can be defined by specifying the corresponding coefficients.

Attention: Valid for Implicit and Explicit Analysis

Format A

Generalized Mooney-Rivlin Polynomial (MOONEY), Reduced Polynomial (RPOLY), Physical Mooney-Rivlin (MOOR), Neo-Hookean (NEOH), and Yeoh model (YEOH):
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
MATHE MID Model NU RHO TEXP TREF
C10 C01 D TAB1 TAB2 TAB4 TABD
C20 C11 C02 D NA ND
C30 C21 C12 C03 D
C40 C31 C22 C13 C04 D
C50 C41 C32 C23 C14 C05 D
MODULI MTIME

Format B

Arruda-Boyce Model (Model=ABOYCE):
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
MATHE MID Model NU RHO TEXP TREF
C λmλm TAB1 TAB2 TAB4
D
MODULI MTIME

Format C

Ogden Material Model (Model=OGDEN):
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
MATHE MID Model NA NU RHO TEXP TREF
MU1 ALPHA1 D TAB1 TAB2 TAB4
MU2 ALPHA2 MU3 ALPHA3
MU4 ALPHA4 MU5 ALPHA5
MODULI MTIME

Format D

Hill Foam Material Model (Model=FOAM):
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
MATHE MID Model NA NU RHO TEXP TREF
MU1 ALPHA1 BETA1 TAB1 TAB2 TAB4
MU2 ALPHA2 BETA2 MU3 ALPHA3 BETA3
MU4 ALPHA4 BETA4 MU5 ALPHA5 BETA5
MODULI MTIME

Format E

Marlow Material Model (Model=MARLOW):
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
MATHE MID Model NU RHO TEXP TREF
D TAB1 TAB2 TAB4 TABD
MODULI MTIME

Example

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
MATHE 2 MOONEY
80 20 0.001

Definitions

Field Contents SI Unit Example
MID Unique material identification number.

No default (Integer > 0)

Model Hyperelastic material model type. See Comments.
MOONEY (Default)
Generalized Mooney-Rivlin hyperelastic model
MOOR
Physical Mooney-Rivlin model
RPOLY
Reduced Polynomial model
NEOH
Neo-Hookean model
YEOH
Yeoh model
ABOYCE
Arruda-Boyce model
OGDEN
Ogden model
FOAM
Hill foam model
MARLOW
Marlow model
blank

(Character)

NU Poisson's ratio.

Default = 0.495 for all models (except FOAM)

Default = 0.0 for FOAM (Real)

RHO Material density.

No default (Real)

TEXP Coefficient of thermal expansion.

No default (Real)

TREF Reference temperature.

No default (Real)

NA Order of the distortional strain energy polynomial function if the type of the model is generalized polynomial (MOONEY) or Reduced Polynomial (RPOLY).

It is also the Order of the Deviatoric Part of the Strain Energy Function of the OGDEN material (Format C).

Default = 2 (0 < Integer ≤ 5)

ND Order of the volumetric strain energy polynomial function. 3

Default = 1 (Integer > 0)

Cpq Material constants related to distortional deformation.

No default (Real)

Dp Material constant related to volumetric deformation (MODEL=BOYCE).

No defaults (Real ≥ 0.0)

TAB1 Table identification number of a TABLES1 entry that contains simple tension-compression data to be used in the estimation of the material constants, related to distortional deformation. The x-values in the TABLES1 entry should be the stretch ratios and y-values should be values of the engineering stress.

(Integer > 0 or blank)

TAB2 Table identification number of a TABLES1 entry that contains equi-biaxial tension data to be used in the estimation of the material constants, related to distortional deformation. The x-values in the TABLES1 entry should be the stretch ratios and y-values should be values of the engineering stress.

(Integer > 0 or blank)

TAB4 Table identification number of a TABLES1 entry that contains pure shear data to be used in the estimation of the material constants, related to distortional deformation. The x-values in the TABLES1 entry should be the stretch ratios and y-values should be values of the engineering stress.

(Integer > 0 or blank)

TABD Table identification number of a TABLES1 entry that contains volumetric part (Dp) of the data to be used in the estimation of the material constants. The x-values in the TABLES1 entry should be the volumetric ratio and y-values should be values of the pressure.

TABD can only be used to fit volumetric data for formats A, B, C, and E. Additionally, only first-order fitting is currently supported (only D value is sourced from the TABD data).

(Integer > 0 or blank)

C Initial shear modulus (Model = ABOYCE). 4

No default (Real)

λmλm Maximum locking stretch ratio.

Used to calculate the value of ββ (Model = ABOYCE). 4

No default (Real)

MUi, ALPHAi Material constants for the Ogden Material Model (Model = OGDEN) 5; or

Hill Foam Material Model (Model = FOAM. 6

BETAi Material constants for Hill Foam Material Model (Model=FOAM). 6
MODULI Continuation line flag for moduli temporal property. 10
MTIME Material temporal property. This field controls the interpretation of the input material property for viscoelasticity.
INSTANT
This material property is considered as the Instantaneous material input for viscoelasticity on the MATVE entry.
LONG (Default)
This material property is considered as the Long-term relaxed material input for viscoelasticity on the MATVE entry.

Comments

  1. If the Cpq and TAB# fields are input, the Cpq (≠ 0.0) values are overwritten with the curve fit values based on the corresponding TAB# tables. However, any Cpq values set to 0.0 are not overwritten.
  2. The Generalized polynomial form (MOONEY) of the hyperelastic material model is written as a combination of the deviatoric and volumetric strain energy of the material. The potential or strain energy density ( WW ) is written in polynomial form, as:

    Generalized polynomial form (MOONEY):

    W=N1p+q=1Cpq(ˉI13)p(ˉI23)q+N2p=11Dp(Jelas1)2pW=N1p+q=1Cpq(¯I13)p(¯I23)q+N2p=11Dp(Jelas1)2p

    Where,
    N1N1
    Order of the distortional strain energy polynomial function (NA).
    N2N2
    Order of the volumetric strain energy polynomial function (ND). Currently only first order volumetric strain energy functions are supported (ND=1).
    CpqCpq
    The material constants related to distortional deformation ( CpqCpq ).
    ˉI1¯I1 , ˉI2¯I2
    Strain invariants, calculated internally by OptiStruct.
    DpDp
    Material constants related to volumetric deformation ( DpDp ). These values define the compressibility of the material.
    JelasJelas
    Elastic volume strain, calculated internally by OptiStruct.
  3. The polynomial form can be used to model the following material types by specifying the corresponding coefficients ( CpqCpq , DpDp ) on the MATHE entry.

    Physical Mooney-Rivlin Material (MOOR):

    N1 = N2 =1

    W=C10(ˉI13)+C01(ˉI23)+1D1(Jelas1)2W=C10(¯I13)+C01(¯I23)+1D1(Jelas1)2

    Reduced Polynomial (RPOLY):

    q=0, N2 =1

    W=N1p=1Cp0(ˉI13)p+1D1(Jelas1)2W=N1p=1Cp0(¯I13)p+1D1(Jelas1)2

    Neo-Hooken Material (NEOH):

    N1= N2 =1, q=0

    W=C10(ˉI13)+1D1(Jelas1)2W=C10(¯I13)+1D1(Jelas1)2

    Yeoh Material (YEOH):

    N1 =3 N2 =1, q=0

    W=C10(ˉI13)+1D1(Jelas1)2+C20(ˉI13)2+C30(ˉI13)3W=C10(¯I13)+1D1(Jelas1)2+C20(¯I13)2+C30(¯I13)3

    Some other material models from the Generalized Mooney Rivlin model are:

    Three term Mooney-Rivlin Material:

    W=C10(ˉI13)+C01(ˉI23)+C11(ˉI13)(ˉI23)+1D1(Jelas1)2W=C10(¯I13)+C01(¯I23)+C11(¯I13)(¯I23)+1D1(Jelas1)2

    Signiorini Material:

    W=C10(ˉI13)+C01(ˉI23)+C20(ˉI13)2+1D1(Jelas1)2W=C10(¯I13)+C01(¯I23)+C20(¯I13)2+1D1(Jelas1)2

    Third Order Invariant Material:

    W=C10(ˉI13)+C01(ˉI23)+C11(ˉI13)(ˉI23)+C20(ˉI13)2+1D1(Jelas1)2W=C10(¯I13)+C01(¯I23)+C11(¯I13)(¯I23)+C20(¯I13)2+1D1(Jelas1)2

    Third Order Deformation Material (James-Green-Simpson):

    W=C10(ˉI13)+C01(ˉI23)+C11(ˉI13)(ˉI23)+C20(ˉI13)2+C30(ˉI13)3+1D1(Jelas1)2

  4. The Arruda-Boyce model (ABOYCE) is defined as:
    W=C15i=1αiβi1(ˉIi13i)+1D1(Jelas1)2

    Where,

    β=1N=1λ2m
    N
    Measure of the limiting locking stretch ratio.
    λm
    Maximum locking stretch ratio.
    D1
    Related to volumetric deformation. It defines the compressibility of the material.
    ˉI1
    First strain invariant, internally calculated by OptiStruct.
    Wherein, ˉI1=I1J23 .
    Jelas
    Elastic volume strain, internally calculated by OptiStruct.
    C1
    Initial shear modulus.

    α1=12;α2=120;α3=111050;α4=197000;α5=519673750

  5. The Ogden Material model (OGDEN) is defined as:
    W=N1i=12μiα2i(ˉλαi1+ˉλαi2+ˉλαi33)+1D1(Jelas1)2
    Where,
    ˉλ1,ˉλ2,ˉλ3
    The three deviatoric stretch ratios (deviatoric stretch ratios are related to principal stretch ratios by ˉλi=J13λi ).
    μi
    Defined by the MUi fields.
    αi
    Defined by the ALPHAi fields.
    N1
    Order of the deviatoric part of the strain energy function defined on the NA field.
  6. The Hill Foam Material model (FOAM) is defined as:
    W=N1i=12μiα2i(λαi1+λαi2+λαi33+1βi(Jαiβi1))
    Where,
    λ1,λ2,λ3
    Principle stretch ratios.
    μi
    Defined by the MUi fields.
    αi
    Defined by the ALPHAi fields.
    βi
    Defined by the BETAi fields.
    N1
    Order of the strain energy function defined on the NA field.

    Additionally, MUi/ALPHAi can instead be fitted using TAB# table data, and BETAi are user specified values.

    If the TAB# fields are input, the MUi/ALPHAi values are overwritten by the fitted values. Any user specified values of MUi/ALPHAi will be overwritten.

    If both Poisson’s ratio NU (non-zero) and TAB# are specified, BETAi values will all be determined or overwritten by:

    βi=ν12ν

    If Poisson’s ratio NU is 0.0 or not specified, then it will be ignored. For parameter fitting, only the first value BETA1 will be used and BETA2-BETA5 are not used. It is recommended to use the same value of BETAi for parameter fitting.

    The Hill FOAM material is supported for both Implicit and Explicit Nonlinear analysis.

  7. The Marlow model is a hyperelastic material model which directly defines the potential based on the experiment test data; there are no mathematical expressions based on the deformation tensors’ invariants or the deformation stretches for the potential. The isochoric deformation potential is determined by TAB1, TAB2 or TAB4. Only one test can be specified.

    A uniaxial tension test is equivalent to an equi-biaxial compression test; a uniaxial compression test is equivalent to an equi-biaxial tension test; a planar tension test is equivalent to a planar compression test. Either tension or compression test data can be specified but not at the same time.

    For Marlow, D, TABD, or Poisson’s ratio can be defined to specify the volumetric behavior. Either D or TABD can be specified, but not both.
    1. If D or TABD is specified, the volumetric behavior is determined by D or TABD.
    2. If D and TABD are not specified and Poisson’s ratio is specified, Poisson’s ratio is used to determine volumetric behavior.
    3. If D, TABD, or Poisson’s ratio are all not specified, the default Poisson’s ratio of 0.495 is used to determine volumetric behavior.
    4. If Poisson’s ratio and one of D or TABD are defined, D or TABD take precedence.
  8. If Poisson’s ratio and D or TABD are both defined, Poisson’s ratio takes precedence.
  9. The initial modulus used for linear analysis is:
    • Mooney, Neo-Hookean, Mooney-Rivlin, Yeoh, Reduced Polynomial

      G=2(C10+C01) and K=2D1

    • Ogden

      G=μ1+μ2+μ3+μ4+μ5=5t=1μt and K=2D1

    • Arruda-Boyce

      G=C1(1+35λ2m+99175λ4m+513875λ6m+4203967375λ8m) and K=2D1

    • Hill Foam

      G=μ1+μ2+μ3+μ4+μ5=5i=1μi and K=5i=12μi(13+βi)

    Additional treatments on bulk modulus K are as:
    • If Poisson's ratio ν , is non-zero, bulk modulus K is replaced with:
      K=2G(1+ν)3(12ν)
    • If K=0, K is set to be 30G
    • If K30G, K is set to be 30G

    The Young’s modulus and Poisson's ratio are given by:

    E=9KG3K+G

    and

    ν=3K2G6K+2G

    Where,
    E
    Young’s modulus.
    G
    Shear modulus.
    K
    Bulk modulus.
    ν
    Poisson's ratio.
    C10 , C01 , D1 , μi , and C1
    Material coefficients.
    λm
    Stretch ratio at which the polymer chain network is locked.
  10. MODULI continuation line is only applicable when used together with the MATVE entry. Refer to MATVE which provides additional information on how this material input is interpreted.
  11. The support information for the available material models (in Model field) is:
    Analysis Type Support Information
    Implicit Analysis All the material models are supported with:
    Explicit Dynamic Analysis All the material models are supported with:
  12. Temperature-dependent hyperelastic material data can be defined via the MATTHE entry.
  13. This card is represented as a material in HyperMesh.