!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“ პარამეტრი მიმართავს, განსაზღვრული არ არის.

თუ „!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>-----იუნგის მოდული (სავალდებულო) <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

შენიშვნა:

მოდალურ ანალიზში ტემპერატურაზე დამოკიდებული მასალის გამოყენება შეუძლებელია.

თბოგამტარობის ანალიზი¶

ბმულის, სიბრტყის, მოცულობითი და სამგანზომილებიანი ფირფიტის ელემენტებისთვის

!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