!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 | Материал қасиетінің мәні (температураға тәуелсіз) |
Note:
- Материал атаулары қайталанса, қате пайда болады.
!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】
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>-----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
Note:
- Модальдық талдауда температураға тәуелді материалдарды қолдануға болмайды.
Жылу өткізгіштік талдауы¶
Сызықтық, жазық, көлемдік және 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