!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	តម្លៃលក្ខណៈសម្បត្តិ (មិនពឹងផ្អែកលើសីតុណ្ហភាព)

ចំណាំ:

ការប្រើឈ្មោះសម្ភារៈដដែលច្រើនដងនឹងបង្កកំហុស។

ប្រសិនបើ MATERIAL ដែលត្រូវបានយោងដោយជម្រើស !SECTION មិនបានកំណត់ នឹងបង្កកំហុស។

ប្រសិនបើលក្ខណៈសម្បត្តិត្រូវបានបញ្ចូលសម្រាប់ធាតុនីមួយៗដោយប្រើប៉ារ៉ាម៉ែត្រ MATITEM ក្នុងជម្រើស !ELEMENT តម្លៃនោះមានអាទិភាព។ ក្នុងករណីនេះ តម្លៃលក្ខណៈសម្បត្តិដែលបញ្ចូលដោយជម្រើស !MATERIAL មិនត្រូវបានប្រើទេ។

ប្រសិនបើចំនួនជម្រើសរង !ITEM=m មិនត្រូវគ្នានឹងតម្លៃប៉ារ៉ាម៉ែត្រ ITEM ឬមានជម្រើសរងដែលមិនបានកំណត់ នឹងបង្កកំហុស។

ជម្រើសរង !ITEM=m មិនចាំបាច់រៀបតាមលំដាប់ m ពីតូចទៅធំទេ។

នៅពេលប្រើជម្រើសរង SUBITEM=k ឬការពឹងផ្អែកលើសីតុណ្ហភាព តម្លៃដែលរំលងមានតម្លៃ 0.0។

នៅពេលប្រើការពឹងផ្អែកលើសីតុណ្ហភាព ត្រូវកំណត់តាមលំដាប់សីតុណ្ហភាពពីទាបទៅខ្ពស់។

នៅពេលប្រើការពឹងផ្អែកលើសីតុណ្ហភាព ការប្រើសីតុណ្ហភាពដូចគ្នាពីរដងឬច្រើនជាងនេះនឹងបង្កកំហុស។

ឧទាហរណ៍ការប្រើប្រាស់

!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>	ឈ្មោះសម្ភារៈ ដែលត្រូវគ្នានឹង MATERIAL ក្នុង !SECTION
ITEM	<ITEMnum>	ចំនួនធាតុលក្ខណៈសម្បត្តិកំណត់ដោយអ្នកប្រើ (1 ឬច្រើនជាង) <Young_modulus>-----ម៉ូឌុល Young (ចាំបាច់) <Poisson_ratio>-------សមាមាត្រ Poisson (ចាំបាច់) <Density>--------------ដង់ស៊ីតេម៉ាស (ចាំបាច់នៅពេល ITEMnum=2) < Expansion_coeff >--មេគុណពង្រីកលីនេអ៊ែរ (នៅពេល ITEMnum=3)

(ឧទាហរណ៍)

!MATERIAL, NAME=M1,
ITEM=3                   --- មានន័យថាកំណត់តម្លៃលក្ខណៈសម្បត្តិ 3 ប្រភេទសម្រាប់សម្ភារៈឈ្មោះ M1
!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

Note:

មិនអាចប្រើសម្ភារៈដែលមានលក្ខណៈសម្បត្តិពឹងផ្អែកលើសីតុណ្ហភាពក្នុងការវិភាគម៉ូដបានទេ។

ការវិភាគចម្លងកម្តៅ¶

សម្រាប់ធាតុ link, plane, សូលីដ និង plate 3D

!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>	ឈ្មោះសម្ភារៈ ដែលត្រូវគ្នានឹង MATERIAL ក្នុង !SECTION
`ITEM`	<ITEMnum>	ចំនួនធាតុលក្ខណៈសម្បត្តិកំណត់ដោយអ្នកប្រើ (តែងតែ 3) <Density>----------ដង់ស៊ីតេ <Specific_heat>----កម្តៅជាក់លាក់ <Conductivity>-----សមត្ថភាពដឹកនាំកម្តៅ <Temperature>------សីតុណ្ហភាព

(ឧទាហរណ៍)

!MATERIAL, NAME=M1,
ITEM=3                   ---　មានន័យថាកំណត់តម្លៃលក្ខណៈសម្បត្តិ 3 ប្រភេទសម្រាប់សម្ភារៈឈ្មោះ M1

!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   ---　កំណត់ ស៊ិចស្យុង
1.0, 20.15, 8.99835E-9, 8.99835E-9

ក្នុង !SECTION ខាងលើ ប៉ារ៉ាម៉ែត្រ គម្លាត ត្រូវបានកំណត់សម្រាប់ធាតុ ចំណុចប្រទាក់ ដែលស្ថិតក្នុងក្រុមឈ្មោះ GAP។

ប៉ារ៉ាម៉ែត្រទី 1: ទទឹង គម្លាត
ប៉ារ៉ាម៉ែត្រទី 2: មេគុណផ្ទេរកម្តៅ គម្លាត
ប៉ារ៉ាម៉ែត្រទី 3: មេគុណវិទ្យុសកម្ម គម្លាត 1
ប៉ារ៉ាម៉ែត្រទី 4: មេគុណវិទ្យុសកម្ម គម្លាត 2

យោង¶

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