!MATERIAL¶
የቁሳቁስ ባህሪያት ትርጉም
የቁሳቁስ ባህሪው በሙቀት ላይ የተመሰረተ ከሆነ፣ ለእያንዳንዱ ተዛማጅ ሙቀት በሰንጠረዥ መልክ ግቤት መስጠት ይቻላል።
መስመር 1
| ፓራሜትር | |
|---|---|
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ኛው የቁሳቁስ ባህሪ በሙቀት ላይ ያልተመሰረተ ከሆነ】
| የተለዋዋጭ ስም | አይነት | ይዘት |
|---|---|---|
| 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】
ምሳሌ 2【የITEM ፓራሜትር ብዛት ከ!ITEM=m ንዑስ አማራጮች ብዛት ጋር አይመሳሰልም-2】
የላስቲክ ስታቲክ ትንተና እና ሞዳል ትንተና¶
!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>-----የያንግ ሞዱለስ (አስፈላጊ) <Poisson_ratio>-------የፖይሰን ሬሾ (አስፈላጊ) <Density>--------------የጅምላ ጥግግት (ITEMnum=2 ሲሆን አስፈላጊ) < Expansion_coeff >--የመስመራዊ ሙቀት መስፋፋት ኮኤፊሺየንት (ITEMnum=3 ሲሆን) |
(ምሳሌ)
!MATERIAL, NAME=M1,
ITEM=3 --- ለቁሳቁስ ስም M1 3 ዓይነት የባህሪ እሴቶች እንደሚገለጹ ያመለክታል
!ITEM=1, SUBITEM=2 --- በ!ITEM=1 የያንግ ሞዱለስ እና የፖይሰን ሬሾ ይገለጻሉ (አስፈላጊ)
4000.0, 0.3
!ITEM=2 --- በ!ITEM=2 የጅምላ ጥግግትን ይግለጹ (ITEM=2 ሲሆን አስፈላጊ)
8.0102E-0
!ITEM=3 --- በ!ITEM=3 የመስመራዊ ሙቀት መስፋፋት ኮኤፊሺየንትን ይግለጹ
1.0E-5
ማስታወሻ:
- በሞዳል ትንተና ውስጥ በሙቀት ላይ የተመሰረቱ ቁሳቁሶችን መጠቀም አይቻልም።
የሙቀት ማስተላለፊያ ትንተና¶
ለሊንክ፣ ፕሌን፣ ሶሊድ እና 3-ልኬት ፕሌት ኤለመንቶች
!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 ውስጥ፣ ለኢንተርፌስ ኤለመንቶች የቡድን ስም=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