Skip to content

Material Data

Elastic Static Analysis, Linear Dynamic Analysis and Eigenvalue Analysis

In the elastic static analysis and eigenvalue analysis of FrontISTR, it is necessary to use isotropic elastic material and define the Young's modulus, Poisson's ratio density and the coefficient of linear expansion for each element. The property value of these materials is defined in the header !SECTION and !MATERIAL of the data. An example of the definition is in the following.

(Example)

!! Definition of SECTION
!SECTION,TYPE=SOLID,EGRP=ALL,MATERIAL=M1

In the above !SECTION, the material data name of the element belonging to the group name=ALL in the solid type element is defined as M1.

Next, the definition method of the material data is shown in the following.

(Example)

!! Intention of defining three types of property values in the material name M1
!MATERIAL,NAME=M1,ITEM=3

!! Young's modulus and Poisson's ratio are defined in !ITEM=1 (mandatory)
!ITEM=1,SUBITEM=2
 4000.0,0.3

!! The mass density must be defined in !ITEM=2 (mandatory in the case of ITEM=3)
!ITEM=2
 8.0102E-10

!! The coefficient of linear expansion must be defined in !ITEM=3
!ITEM=3
 1.0E-5

As long as each item number and the type of physical property are compatible, and the ITEM number is correct, it can be defined in any order. However, ITEM=1 must be defined in the order of the Young's modulus and Poisson's ratio.

Heat Conduction Analysis

In the heat conduction analysis of FrontISTR, isotropic material taking into consideration the temperature dependency can be used. Regarding the property values to be defined, the density, specific heat and thermal conductivity are defined for the link, plane surface, solid and shell element, and the gap heat transfer coefficient and the gap radiation factor are defined in the interface element. An example of the definition method of these property values is shown in the following.

Defined in the !SECTION and !MATERIAL header.

(Example)

!! Definition of section
!SECTION,TYPE=SOLID,EGRP=ALL,MATERIAL=M1

In the above !SECTION, the material data name of the element belonging to the group name=ALL in the solid type element is defined as M1. The following shows the definition method of the material data.

(Example)

!! Intention of defining three types of property values
!! in the material of material name M1
!MATERIAL,NAME=M1,ITEM=3

!! The density and temperature are defined in !ITEM=1 (mandatory)
!ITEM=1,SUBITEM=1
 7850., 300.
 7790., 500.
 7700., 800.

!! The specific heat and temperature are defined in !ITEM=2 (mandatory)
!ITEM=2
 0.465, 300.
 0.528, 500.
 0.622, 800.

!! The thermal conductivity and temperature are defined in !ITEM=3 (mandatory)
!ITEM=3
 43., 300.
 38.6, 500.
 27.7, 800.

Since the type of physical properties are identified by each item number in FrontISTR, it can be defined in any order as long as both are consistent.

(2) In the case of an interface element

Defined in the !SECTION header. (Material data is not required)

(Example)

!! Definition of section
!SECTION,TYPE=INTERFACE,EGRP=GAP
 1.0, 20.15, 8.99835E-9, 8.99835E-9

In the above !SECTION, the gap parameter of the element belonging to the group name=GAP in the interface element is defined.

  • 1st parameter : Gap width
  • 2nd parameter : Gap heat transfer coefficient
  • 3rd parameter : Gap radiation factor 1
  • 4th parameter : Gap radiation factor 2

(3) In the case of a shell element

Defined in the !SECTION and MATERIAL header.

(Example)

!! Definition of section
!SECTION, TYPE=SHELL, EGRP=SH, MATERIAL=M2
 10.0, 5

In the above !SECTION, the shell characteristics of the element belonging to the group name=SH in the shell type element is defined.

  • 1st parameter : Shell thickness
  • 2nd parameter : Integration points in the thickness direction

This also means that the material data name of the material physical properties of the element which belong to this group is defined as M2.

For single layered isotropic shell, the definition method of the material physical properties is the same as in the case of the solid element. Refer to the description of the solid element.

In addition, another definition method stated below is available for multi-layered and/or orthotropic shell.

(Example) Isotropic, single-layered shell

!MATERIAL, NAME=M1, ITEM=1
!ITEM=1,SUBITEM=4
 0, 200000, 0.3, 2.0
  • 1st parameter : 0=isotropic
  • 2nd parameter : Young’s modulus
  • 3rd parameter : Poisson’s ratio
  • 4th parameter : Shell thickness

(Example) Isotropic, double-layered shell

!MATERIAL, NAME=M1, ITEM=1
!ITEM=1,SUBITEM=7
 0, 200000, 0.3, 2.0, 200000, 0.3, 2.0
  • 1st parameter : 0=isotropic
  • 2nd parameter : Young’s modulus (1st layer)
  • 3rd parameter : Poisson’s ratio (1st layer)
  • 4th parameter : Shell thickness (1st layer)
  • 5th parameter : : Young’s modulus (2nd layer)
  • 6th parameter : Poisson’s ratio (2nd layer)
  • 7th parameter : Shell thickness (2nd layer)

By specifying as many material constants as needed for the number of layers, the material is automatically treated as multi-layered shell. The total thickness of the shell is the sum of the thickness of all the layers. The layers are ordered from the top surface of the shell with respect to the normal direction.

(Example) Orthotropic, single-layered shell

!MATERIAL, NAME=M1, ITEM=1
!ITEM=1, SUBITEM=9
 1, 28600., 0.15, 32.3 , 28600., 12434., 12434., 12434., 0.0
  • 1st parameter : 1=orthotropic
  • 2nd parameter : Young’s modulus E1
  • 3rd parameter : Poisson’s ratio n12
  • 4th parameter : Shell thickness
  • 5th parameter : : Young’s modulus E2
  • 6th parameter : Shear modulus G12
  • 7th parameter : Shear modulus G23
  • 8th parameter : Shear modulus G13
  • 9th parameter : Orthotropy angle θ

Orthotropy angle is specified in θ [degree].

(Example) Orthotropic, double-layered shell

!MATERIAL, NAME=M1, ITEM=1
!ITEM=1, SUBITEM=17
 1, 28600., 28600., 0.15, 32.3 , 12434., 12434., 28600., 12434., 12434., 12434., 0.0
  • 1st parameter : 1=orthotropic
  • 2nd parameter : Young’s modulus E1 (1st layer)
  • 3rd parameter : Poisson’s ratio n12 (1st layer)
  • 4th parameter : Shell thickness (1st layer)
  • 5th parameter : : Young’s modulus E2 (1st layer)
  • 6th parameter : Shear modulus G12 (1st layer)
  • 7th parameter : Shear modulus G23 (1st layer)
  • 8th parameter : Shear modulus G13 (1st layer)
  • 9th parameter : Orthotropy angle q (1st layer)
  • 10th parameter : Young’s modulus E1 (2nd layer)
  • 11th parameter : Poisson’s ratio n12 (2nd layer)
  • 12th parameter : Shell thickness (2nd layer)
  • 13th parameter : Young’s modulus E2 (2nd layer)
  • 14th parameter : Shear modulus G12 (2nd layer)
  • 15th parameter : Shear modulus G23 (2nd layer)
  • 16th parameter : Shear modulus G13 (2nd layer)
  • 17th parameter : Orthotropy angle θ (2nd layer)

Nonlinear Static Analysis

In the nonlinear static analysis of FrontISTR, in addition to the method to define the !SECTION and !MATERIAL shown in item 4.2.1, !ELASTIC, !HYPERELASTIC, !PLASTIC and etc. in the analysis control data can also be defined. An example of the definition is shown in the following.

(Example) Definition of hyperelastic material

!MATERIAL

!! Definition of the Neo Hooke hyperelastic material
!HYPERELASTIC, TYPE=NEOHOOKE

!! C_10 and D are defined (mandatory)
 1000.0, 0.00005

(Example) Definition of the elastoplastic material

!MATERIAL

!! Definition of the isotropic elastic material
!ELASTIC, TYPE=ISOTROPIC

!! The Young's modulus and the Poisson's ratio are defined (mandatory)
 21000.0, 0.3

!! Definition of the Drucker-Prager plastic material
!PLASTIC, TYPE=DRUCKER-PRAGER

!! Viscosity, angle of friction and hardening factor are defined (mandatory)
 500.0, 4.0, 10.0

Viscoelastic Material

In the viscoelastic analysis of FrontISTR, a generalized Maxwell model with Prony series is used. Temperature dependence can be considered through Thermorheological Simplicity (TRS).

Shift Factor Formulation

The DEFINITION parameter of the !TRS keyword selects one of the following two shift factor types.

WLF type (DEFINITION=WLF)

The shift factor \(A\) is defined as:

\[\log_{10}(A) = \frac{C_1(\theta - \theta_0)}{C_2 + \theta - \theta_0}\]
  • \(\theta_0\): Reference temperature
  • \(C_1\), \(C_2\): WLF material constants

Arrhenius type (DEFINITION=ARRHENIUS)

The shift factor \(A\) is defined as:

\[\ln(A) = \frac{E_0}{R}\left(\frac{1}{\theta - \theta^Z} - \frac{1}{\theta_0 - \theta^Z}\right)\]
  • \(\theta_0\): Reference temperature
  • \(C_1 = E_0/R\): ratio of activation energy \(E_0\) to gas constant \(R\)
  • \(C_2 = \theta^Z\): absolute zero offset (e.g., 0 for Kelvin, -273.15 for Celsius)