Translation status

This page has not yet been translated into Français, or its translation is out of date. The current English version is shown as a fallback.

!MATERIAL¶

Definition of material physical properties

If a physical property is temperature-dependent, table input is possible for each corresponding temperature.

1st line

!MATERIAL, NAME=<name> [, optional parameter]

Parameter
`NAME`	Material name (required)
`ITEM`	Number of physical property items (optional; becomes "1" if omitted)
`INPUT`	External file name (optional)

Parameter Name	Parameter Value	Description
`NAME`	<name>	Material name
`ITEM`	<ITEMnum>	Number of user-defined physical property items
`INPUT`	<filename>	External file name (optional); may also be used together with the 2nd and subsequent lines

2nd and subsequent lines

(2nd line) !ITEM=1, SUBITEM=<k>
(3rd line) VAL1-1-1, VAL1-1-2, … VAL1-1-k, TEMP1-1
(4th line) VAL1-2-1，VAL1-2-2, … VAL1-2-k, TEMP1-2
...

(Line L+2) VAL1-L-1，VAL1-L-2, … VAL1-L-k, TEMP1-L

Repeat the definition through !ITEM=<ITEMnum>.

Subparameter (for `!ITEM`)
`SUBITEM`	Number of sub-property items defined for each physical property item
	(optional; becomes "1" if omitted)

Subparameter Name	Parameter Value	Description
`SUBITEM`	<subITEMnum>	Number of user-defined sub-property items

【When the m-th physical property is temperature-dependent】

If the number of entries in the temperature-dependent table is N, enter them as follows:

!ITEM=m, SUBITEM=k
VALm1-1, ..., VALm1-k, TEMPm1
VALm2-1, ..., VALm2-k, TEMPm2
...
VALmN-1, ..., VALmN-k, TEMPm-N

Parameter Name	Attribute	Description
VALmn-k	R	Physical property value (temperature-dependent)
TEMPmn	R	Corresponding temperature

TEMPm1 < TEMPm2 < ... < TEMPmN must hold.

When the temperature is at or below TEMPm1, VALm1 is used; when it is at or above TEMPmN, VALmN is used.

【When the m-th physical property is not temperature-dependent】

!ITEM=m, SUBITEM=k
VALm1-1, ..., VALm1-k
VALm2-1, ..., VALm2-k
...
VALmN-1, ..., VALmN-k

Parameter Name	Attribute	Description
VALmn-k	R	Physical property value (not temperature-dependent)

Note:

An error occurs if material names are duplicated.

An error occurs if the MATERIAL referenced by the !SECTION option is not defined.

If physical properties are entered for each element using the MATITEM parameter of the !ELEMENT option, those values take precedence. In this case, the physical property values entered using the !MATERIAL option are not used.

An error occurs if the number of !ITEM=m suboptions is inconsistent with the number specified by the ITEM parameter, or if any suboption is undefined.

The !ITEM=m suboptions do not need to be ordered by increasing m.

When the !SUBITEM=k suboption and temperature dependence are used, omitted values become "0.0".

When temperature dependence is used, definitions must be given in ascending order of temperature.

When temperature dependence is used, an error occurs if the same temperature is used two or more times.

Example

!MATERIAL, NAME= STEEL, ITEM= 2
!ITEM=1 No temperature dependence
35.0
!ITEM=2
40.0, 0.0
45.0, 100.0
50.0, 200.0

!MATERIAL, NAME= CUPPER         Number of items = 1 (default value)
!ITEM=1 No temperature dependence
80.0

Incorrect examples

Example 1【The number specified by the ITEM parameter and the number of !ITEM=m suboptions are inconsistent -1】

!MATERIAL, NAME= STEEL, ITEM= 2
!ITEM=3
20.0
!ITEM=1
35.0
!ITEM= 2
40.0

Example 2【The number specified by the ITEM parameter and the number of !ITEM=m suboptions are inconsistent -2】

!MATERIAL, NAME= STEEL, ITEM= 3
!ITEM=3
20.0
!ITEM= 2
40.0
!MATERIAL, NAME= CUPPER
...

!MATERIAL, NAME=<name>, ITEM=<ITEMnum>
!ITEM=1, SUBITEM=2
<Young_modulus>, <Poisson_ratio>
!ITEM=2
<Density>
!ITEM=3
<Expansion_coeff>

Parameter Name	Parameter Value	Description
NAME	<name>	Material name; corresponds to MATELIAL in !SECTION
ITEM	<ITEMnum>	Number of user-defined physical property items (1 or more) <Young_modulus>-----Young's modulus (required) <Poisson_ratio>-------Poisson's ratio (required) <Density>--------------Mass density (required when ITEMnum=2) < Expansion_coeff >--Coefficient of linear thermal expansion (when ITEMnum=3)

(Example)

!MATERIAL, NAME=M1,
ITEM=3                   --- Defines three types of property values for material M1
!ITEM=1, SUBITEM=2       --- Defines Young's modulus and Poisson's ratio in !ITEM=1 (required)
4000.0, 0.3
!ITEM=2                  --- Defines mass density in !ITEM=2 (required when ITEM=2)
8.0102E-0
!ITEM=3                  --- Defines the coefficient of linear thermal expansion in !ITEM=3
1.0E-5

Note:

Temperature-dependent materials cannot be used in modal analysis.

Heat Conduction Analysis¶

For link, plane, solid, and three-dimensional plate elements

!MATERIAL, NAME=<name>, ITEM=3
!ITEM=1, SUBITEM=2
<Density>, <Temperature>

!ITEM=2, SUBITEM=2
<Specific_heat>, <Temperature>

!ITEM=3, SUBITEM=2
<Conductivity>, <Temperature>

Parameter Name	Parameter Value	Description
`NAME`	<name>	Material name; corresponds to MATELIAL in !SECTION
`ITEM`	<ITEMnum>	Number of user-defined physical property items (always 3) <Density>----------Density <Specific_heat>----Specific heat <Conductivity>-----Thermal conductivity <Temperature>------Temperature

(Example)

!MATERIAL, NAME=M1,
ITEM=3　　　             ---　Defines three types of property values for material M1

!ITEM=1, SUBITEM=１　 　 ---　Defines density and temperature in !ITEM=1 (required)
7850., 300.
7790., 500.
7700., 800.

!ITEM=2, SUBITEM=1       ---　Defines specific heat and temperature in !ITEM=2 (required)
0.465, 300.
0.528, 500.
0.622, 800.

!ITEM=3                  ---　Defines thermal conductivity and temperature in !ITEM=3 (required)
43., 300.
38.6, 500.
27.7, 800.

For interface elements

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

(Example)

!SECTION, TYPE=INTERFACE, EGRP=GAP　 ---　Section definition
1.0, 20.15, 8.99835E-9, 8.99835E-9

The above !SECTION defines the gap parameters for elements belonging to the group named GAP in an interface element.

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

Reference¶

program TEST
use hecmw
implicit REAL*8 (A-H,O-Z)
type (hecmwT_local_mesh) :: hecMESH

!C
!C    !MATERIAL, NAME=SUS304, ITEM=3
!C    !ITEM=1, SUBITEM= 3
!C      100.0, 200.0, 300.0, 0.00
!C      101.0, 210.0, 301.0, 1.00
!C      102.0, 220.0, 302.0, 2.00
!C      103.0, 230.0, 303.0, 3.00
!C    !ITEM=3, SUBITEM= 2
!C      1000.0, , 0.00
!C      1001.0, 1., 1.00
!C      1002.0, 2., 2.00
!C      1003.0, 3., 3.00
!C    !ITEM=2
!C      5000.0
!C
!C    !MATERIAL, NAME=FEC, ITEM=2
!C    !ITEM=1, SUBITEM= 3
!C      2100.0, 2200.0, 2300.0, 0.00
!C      2101.0, 2210.0, 2301.0, 1.00
!C      2102.0, 2220.0, 2302.0, 2.00
!C      2103.0, 2230.0, 2303.0, 3.00
!C      3103.0, 3230.0, 2304.0, 4.00
!C    !ITEM=2
!C      6000.0, 10.0
!C      6500.0, 30.0
!C

hecMESH%material%n_mat = 2

nn= hecMESH%material%n_mat
allocate (hecMESH%material%mat_name(nn))

hecMESH%material%mat_name(1)= 'SUS304'
hecMESH%material%mat_name(2)= 'FEC'

nn= hecMESH%material%n_mat
allocate (hecMESH%material%mat_ITEM_index(0:nn))
hecMESH%material%mat_ITEM_index(0)= 0
hecMESH%material%mat_ITEM_index(1)= 3
hecMESH%material%mat_ITEM_index(2)= hecMESH%material%mat_ITEM_index(1) + 2

hecMESH%material%n_mat_ITEM= hecMESH%material%mat_ITEM_index(hecMESH%material%n_mat)

nn= hecMESH%material%n_mat_ITEM
allocate (hecMESH%material%mat_subITEM_index(0:nn))

hecMESH%material%mat_subITEM_index(0)= 0
hecMESH%material%mat_subITEM_index(1)= 3
hecMESH%material%mat_subITEM_index(2)= hecMESH%material%mat_subITEM_index(1) + 1
hecMESH%material%mat_subITEM_index(3)= hecMESH%material%mat_subITEM_index(2) + 2
hecMESH%material%mat_subITEM_index(4)= hecMESH%material%mat_subITEM_index(3) + 3
hecMESH%material%mat_subITEM_index(5)= hecMESH%material%mat_subITEM_index(4) + 1

hecMESH%material%n_mat_subITEM=
&        hecMESH%material%mat_subITEM_index(hecMESH%material%n_mat_ITEM)

nn= hecMESH%material%n_mat_subITEM
allocate (hecMESH%material%mat_TABLE_index(0:nn))
hecMESH%material%mat_TABLE_index( 0)= 0
hecMESH%material%mat_TABLE_index( 1)= 4
hecMESH%material%mat_TABLE_index( 2)= hecMESH%material%mat_TABLE_index( 1) + 4
hecMESH%material%mat_TABLE_index( 3)= hecMESH%material%mat_TABLE_index( 2) + 4
hecMESH%material%mat_TABLE_index( 4)= hecMESH%material%mat_TABLE_index( 3) + 1
hecMESH%material%mat_TABLE_index( 5)= hecMESH%material%mat_TABLE_index( 4) + 4
hecMESH%material%mat_TABLE_index( 6)= hecMESH%material%mat_TABLE_index( 5) + 4
hecMESH%material%mat_TABLE_index( 7)= hecMESH%material%mat_TABLE_index( 6) + 5
hecMESH%material%mat_TABLE_index( 8)= hecMESH%material%mat_TABLE_index( 7) + 5
hecMESH%material%mat_TABLE_index( 9)= hecMESH%material%mat_TABLE_index( 8) + 5
hecMESH%material%mat_TABLE_index(10)= hecMESH%material%mat_TABLE_index( 9) + 2

hecMESH%material%n_mat_TABLE=
&        hecMESH%material%mat_TABLE_index(hecMESH%material%n_mat_subITEM)

nn= hecMESH%material%n_mat_TABLE
allocate (hecMESH%material%mat_VAL (nn))
allocate (hecMESH%material%mat_TEMP(nn))

hecMESH%material%mat_VAL = 0.d0
hecMESH%material%mat_TEMP= 0.d0

hecMESH%material%mat_VAL ( 1)= 100.0d0
hecMESH%material%mat_TEMP( 1)=   0.0d0
hecMESH%material%mat_VAL ( 2)= 101.0d0
hecMESH%material%mat_TEMP( 2)=   1.0d0
hecMESH%material%mat_VAL ( 3)= 102.0d0
hecMESH%material%mat_TEMP( 3)=   2.0d0
hecMESH%material%mat_VAL ( 4)= 103.0d0
hecMESH%material%mat_TEMP( 4)=   3.0d0

hecMESH%material%mat_VAL ( 5)= 200.0d0
hecMESH%material%mat_TEMP( 5)=   0.0d0

hecMESH%material%mat_VAL (13)= 5000.0d0

hecMESH%material%mat_VAL (14)= 1000.0d0
hecMESH%material%mat_TEMP (14)=   0.0d0
hecMESH%material%mat_VAL (15)= 1001.0d0
hecMESH%material%mat_TEMP (15)=   1.0d0
hecMESH%material%mat_VAL (16)= 1002.0d0
hecMESH%material%mat_TEMP (16)=   2.0d0
hecMESH%material%mat_VAL (17)= 1003.0d0
hecMESH%material%mat_TEMP (17)=   3.0d0

hecMESH%material%mat_VAL (18)=  0.0d0
hecMESH%material%mat_TEMP (18)= 0.0d0
hecMESH%material%mat_VAL (19)=  1.0d0
hecMESH%material%mat_TEMP (19)= 1.0d0
hecMESH%material%mat_VAL (20)=  2.0d0
hecMESH%material%mat_TEMP (20)= 2.0d0
hecMESH%material%mat_VAL (21)=  3.0d0
hecMESH%material%mat_TEMP (21)= 3.0d0

hecMESH%material%mat_VAL (22)= 2100.0d0
hecMESH%material%mat_TEMP (22)=   0.0d0
hecMESH%material%mat_VAL (23)= 2101.0d0
hecMESH%material%mat_TEMP (23)=   1.0d0
hecMESH%material%mat_VAL (24)= 2102.0d0
hecMESH%material%mat_TEMP (24)=   2.0d0
hecMESH%material%mat_VAL (25)= 2103.0d0
hecMESH%material%mat_TEMP (25)=   3.0d0
hecMESH%material%mat_VAL (26)= 3103.0d0
hecMESH%material%mat_TEMP (26)=   4.0d0

write (*,'(a,i10)') '%n_mat_ITEM ', hecMESH%material%n_mat_ITEM
write (*,'(a,i10)') '%n_mat_subITEM', hecMESH%material%n_mat_subITEM
write (*,'(a,i10)') '%n_mat_TABLE ', hecMESH%material%n_mat_TABLE

end program TEST

!MATERIAL¶

Elastic Static Analysis and Modal Analysis¶

Heat Conduction Analysis¶

Reference¶