!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” ලෙස ගනී.

උෂ්ණත්ව පරායත්තතාව භාවිත කරන විට අඩු උෂ්ණත්වයේ සිට වැඩි උෂ්ණත්වයට අනුපිළිවෙළින් නිර්වචනය කළ යුතුය.

උෂ්ණත්ව පරායත්තතාව භාවිත කරන විට එකම උෂ්ණත්වය දෙවරක් හෝ වැඩි වාර ගණනක් භාවිත කළහොත් දෝෂයක් වේ.

භාවිත උදාහරණය

!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>-------පොයිසන් අනුපාතය (අනිවාර්ය) <Density>--------------ස්කන්ධ ඝනත්වය (ITEMnum=2 විට අනිවාර්ය) < Expansion_coeff >--රේඛීය තාප ප්‍රසාරණ ගුණකය (ITEMnum=3 විට)

(උදාහරණය)

!MATERIAL, NAME=M1,
ITEM=3                   --- M1 ද්‍රව්‍ය සඳහා ද්‍රව්‍ය ගුණ වර්ග 3 ක් නිර්වචනය කරන බවයි
!ITEM=1, SUBITEM=2       --- !ITEM=1 හි Young මාපාංකය සහ පොයිසන් අනුපාතය නිර්වචනය කරයි (අනිවාර්ය)
4000.0, 0.3
!ITEM=2                  --- !ITEM=2 හි ස්කන්ධ ඝනත්වය නිර්වචනය කරයි (ITEM=2 විට අනිවාර්ය)
8.0102E-0
!ITEM=3                  --- !ITEM=3 හි රේඛීය තාප ප්‍රසාරණ ගුණකය නිර්වචනය කරයි
1.0E-5

සටහන:

මෝඩල් විශ්ලේෂණයේදී උෂ්ණත්ව-පරායත්ත ද්‍රව්‍ය භාවිත කළ නොහැක.

තාප සන්නයන විශ්ලේෂණය¶

සම්බන්ධක, තල, ඝන සහ ත්‍රිමාන පත්‍ර මූලද්‍රව්‍ය සඳහා

!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   --- කොටස නිර්වචනය කිරීම
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