!MATERIAL¶

Материалын шинж чанарын тодорхойлолт

Материалын шинж чанар температураас хамааралтай бол харгалзах температур бүрээр хүснэгтээр оруулж болно.

1-р мөр

!MATERIAL, NAME=<name> [, заавал биш параметр]

Параметр
`NAME`	Материалын нэр (заавал)
`ITEM`	Материалын шинж чанарын зүйлийн тоо (заавал биш; орхивол 「1」 болно)
`INPUT`	Гадаад файлын нэр (заавал биш)

Параметрийн нэр	Параметрийн утга	Тайлбар
`NAME`	<name>	Материалын нэр
`ITEM`	<ITEMnum>	Хэрэглэгчийн тодорхойлсон материалын шинж чанарын зүйлийн тоо
`INPUT`	<filename>	Гадаад файлын нэр (заавал биш); 2-р болон дараагийн мөрүүдтэй хамт ашиглаж болно

2-р болон дараагийн мөрүүд

(2-р мөр) !ITEM=1, SUBITEM=<k>
(3-р мөр) VAL1-1-1, VAL1-1-2, … VAL1-1-k, TEMP1-1
(4-р мөр) VAL1-2-1, VAL1-2-2, … VAL1-2-k, TEMP1-2
・・

(L+2-р мөр)VAL1-L-1, VAL1-L-2, … VAL1-L-k, TEMP1-L

「!ITEM=<ITEMnum>」 хүртэл тодорхойлолтыг давтана

Дэд параметр (`!ITEM`-д хамаарах)
`SUBITEM`	Материалын шинж чанарын зүйл бүрд тодорхойлох дэд шинж чанарын зүйлийн тоо
	(заавал биш; орхивол 「1」 болно.)

Дэд параметрийн нэр	Параметрийн утга	Тайлбар
`SUBITEM`	<subITEMnum>	Хэрэглэгчийн тодорхойлсон дэд шинж чанарын зүйлийн тоо

【m-р материалын шинж чанар температураас хамааралтай үед】

Температурын хамаарлын хүснэгт N зүйлтэй бол дараах байдлаар оруулна:

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

Хувьсагчийн нэр	Атрибут	Тайлбар
VALmn-k	R	Материалын шинж чанарын утга (температураас хамаарна)
TEMPmn	R	Харгалзах температур

TEMPm1 < TEMPm2 < ... < TEMPmN байх ёстой.

Температур TEMPm1-ээс бага буюу тэнцүү бол VALm1, TEMPmN-ээс их буюу тэнцүү бол VALmN-ийг ашиглана.

【m-р материалын шинж чанар температураас хамааралгүй үед】

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

Хувьсагчийн нэр	Атрибут	Тайлбар
VALmn-k	R	Материалын шинж чанарын утга (температураас хамааралгүй)

Тэмдэглэл:

Материалын нэр давхардвал алдаа гарна.

!SECTION сонголтод заасан MATERIAL тодорхойлогдоогүй бол алдаа гарна.

!ELEMENT сонголтод MATITEM параметр ашиглан элемент бүрийн материалын шинж чанарыг оруулсан бол тэдгээр утга давуу эрхтэй. Энэ тохиолдолд !MATERIAL сонголтоор оруулсан материалын шинж чанарын утгыг ашиглахгүй.

!ITEM=m дэд сонголтын тоо ITEM параметрийн тоотой нийцэхгүй эсвэл тодорхойлогдоогүй дэд сонголт байвал алдаа гарна.

!ITEM=m дэд сонголтуудыг m-ийн өсөх дарааллаар байрлуулах шаардлагагүй.

SUBITEM=k дэд сонголт ба температурын хамаарлыг ашиглах үед орхисон утгыг 「0.0」 гэж авна.

Температурын хамаарал ашиглах үед тодорхойлолтуудыг температурын өсөх дарааллаар өгөх ёстой.

Температурын хамаарал ашиглах үед ижил температурыг 2 ба түүнээс олон удаа хэрэглэвэл алдаа гарна.

Жишээ

!MATERIAL, NAME= STEEL, ITEM= 2
!ITEM=1 Температурын хамааралгүй
35.0
!ITEM=2
40.0, 0.0
45.0, 100.0
50.0, 200.0

!MATERIAL, NAME= COPPER         Зүйлийн тоо＝1 (өгөгдмөл утга)
!ITEM=1 Температурын хамааралгүй
80.0

Буруу ашиглах жишээнүүд

Жишээ 1【ITEM параметр ба !ITEM=m дэд сонголтын тоо нийцэхгүй -1】

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

Жишээ 2【ITEM параметр ба !ITEM=m дэд сонголтын тоо нийцэхгүй -2】

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

Уян харимхайн статик шинжилгээ ба модаль шинжилгээ¶

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

Параметрийн нэр	Параметрийн утга	Тайлбар
NAME	<name>	Материалын нэр, !SECTION-ийн MATERIAL-тай харгалзана
ITEM	<ITEMnum>	Хэрэглэгчийн тодорхойлсон материалын шинж чанарын зүйлийн тоо (1 буюу түүнээс их) <Young_modulus>-----Young-ийн модуль (заавал) <Poisson_ratio>-------Poisson-ийн харьцаа (заавал) <Density>--------------Массын нягт (ITEMnum=2 үед заавал) < Expansion_coeff >--Шугаман тэлэлтийн коэффициент (ITEMnum=3 үед)

(Жишээ)

!MATERIAL, NAME=M1,
ITEM=3                   --- M1 нэртэй материалд 3 төрлийн шинж чанарын утга тодорхойлно
!ITEM=1, SUBITEM=2       --- !ITEM=1-д Young-ийн модуль ба Poisson-ийн харьцааг тодорхойлно (заавал)
4000.0, 0.3
!ITEM=2                  --- !ITEM=2-д массын нягтыг тодорхойлно (ITEM=2 үед заавал)
8.0102E-0
!ITEM=3                  --- !ITEM=3-д шугаман тэлэлтийн коэффициентийг тодорхойлно
1.0E-5

Тэмдэглэл:

Модаль шинжилгээнд температураас хамааралтай материал ашиглах боломжгүй.

Дулаан дамжилтын шинжилгээ¶

Link, plane, solid болон 3 хэмжээст plate элементүүдийн үед

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

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

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

Параметрийн нэр	Параметрийн утга	Тайлбар
`NAME`	<name>	Материалын нэр, !SECTION-ийн MATERIAL-тай харгалзана
`ITEM`	<ITEMnum>	Хэрэглэгчийн тодорхойлсон материалын шинж чанарын зүйлийн тоо (үргэлж 3) <Density>----------Нягт <Specific_heat>----Хувийн дулаан багтаамж <Conductivity>-----Дулаан дамжилтын коэффициент <Temperature>------Температур

(Жишээ)

!MATERIAL, NAME=M1,
ITEM=3                   ---　M1 нэртэй материалд 3 төрлийн шинж чанарын утга тодорхойлно

!ITEM=1, SUBITEM=1       ---　!ITEM=1-д нягт ба температурыг тодорхойлно (заавал)
7850., 300.
7790., 500.
7700., 800.

!ITEM=2, SUBITEM=1       ---　!ITEM=2-д хувийн дулаан багтаамж ба температурыг тодорхойлно (заавал)
0.465, 300.
0.528, 500.
0.622, 800.

!ITEM=3                  ---　!ITEM=3-д дулаан дамжилтын коэффициент ба температурыг тодорхойлно (заавал)
43., 300.
38.6, 500.
27.7, 800.

Интерфейс элементийн үед

!SECTION толгой хэсэгт тодорхойлно. (Материалын өгөгдөл шаардлагагүй)

(Жишээ)

!SECTION, TYPE=INTERFACE, EGRP=GAP   ---　Section-ийн тодорхойлолт
1.0, 20.15, 8.99835E-9, 8.99835E-9

Дээрх !SECTION нь интерфейс элементүүдийн GAP нэртэй бүлэгт харьяалагдах элементүүдийн завсрын параметрүүдийг тодорхойлно.

1-р параметр: Завсрын өргөн
2-р параметр: Завсрын дулаан дамжуулалтын коэффициент
3-р параметр: Завсрын цацрагийн коэффициент １
4-р параметр: Завсрын цацрагийн коэффициент ２

Лавлагаа¶

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