name.IB (XRPD)

name.IB (XRPD)

; řízení programu
IP KREST 0 KOUT INT3 NSAX NSTEPX KSU
RP TIMS ERAL EDIF TOL DTRUN PIVOT PENALTA

; nepovinný popis nezávisle proměnných
IV JIV T IV V x₁ x₂ ... x_N

; materiál
MP ISET T 1V l r×c

; přechodový odpor
MP ISET T 2V b

; popis teplot všech uzlů sítě - počáteční podmínka
GV ISET T 1 V T₁ T₂ ... T_NNOD
GV ISET T 1D 4 IREC

; objemový tepelný zdroj
VV ISET T 6V

; přestup tepla konvekcí
SV ISET T 1/11 V a T_o

; přenos tepla zářením
SV ISET T 2/12 V c T_o

; obecný přestup tepla
SV ISET T 3/13 V c₁ c₂ c₃ T_o

; tepelný tok
SV ISET T 4/14V

; přestup tepla konvekcí na hraně prvku semi-loof
LV ISET T 1/11Va T_o

; přenos tepla zářením na hraně prvku semi-loof
LV ISET T 2/12 V c T_o

; obecný přestup tepla na hraně prvku semi-loof
LV ISET T 3/13V c₁ c₂ c₃ T_o

; tepelný tok na hraně prvku semi-loof
LV ISET T 4/14V

; teplota v uzlu
NV ISET T 1/11V T_i

; koncentrovaný tepelný tok
NV ISET T 2/12 V

; řízení výpočtu (platí pouze v místě přiřazení)
AV ISET T 6 N KAPPR KAUTO KPRED

V 0 0 0 0

; První zatěžovací stav. Veličiny s KQT < 6 (viz B7) platí
; pro celý uvažovaný děj.
AS1

; přiřazení MP sad
/MISET /M ISET E seznam prvků

; řízení výpočtu (přiřazení AV sady)
/RTIMX STEP TSC
/AISET

; přiřazení GV sady
/GISET

; přiřazení VV sad
/VISET /V ISET E seznam prvků

; přiřazení SV sad
/SISET E seznam prvkůS č. plochy

; přiřazení LV sad
/L ISET Eseznam prvkůL č. hrany

; přiřazení NV sad
/NISET /N ISET N seznam uzlů

; Další zatěžovací stavy. Přiřazená veličina platí jen v daném zatěžovacím stavu a musí být KQT > 5 - viz B7.
AS2 / ....... / .........

; ukončení souboru
EN
EN

c Radiační konstanta [W/m²K⁴]

c Viz r×c = měrná tepelná kapacita [J/m³K].

c₁ Konstanta obecného přestupu [W/m²K^c2+c3].

c₂ Konstanta obecného přestupu [-].

c₃ Konstanta obecného přestupu [-].

DTRUN Elementární časový krok [s]. Délka skutečného časového kroku se zaokrouhlí na celistvý násobek DTRUN. Uvažuje se jen pro DTRUN > 10^-6 [s].

EDIF Kritérium konvergence přírůstku teplot [^oC]: ||T⁽ⁱ⁾ - T^(i-1)||_MAX < EDIF. Doporučeno 1 < EDIF < 5 ^oC. (uplatní se pouze při KAPPR = 1)

ERAL Residuální kritérium konvergence ||Res T⁽ⁱ⁾|| < ERAL*||T⁽ⁱ⁾|| .Doporučeno 0.0001 < ERAL < 0.1. (uplatní se pouze při KAPPR = 1)

INT3 Pořadové číslo integračního (časového) kroku, od kterého se má pokračovat s KREST = 3. Jestliže KREST = 1, zapíše se INT3 = 0.

IREC Pořadové číslo záznamu v binárním souboru name.IC.

ISET Rozlišovací číslo sady.

IV Identifikační číslo proměnné dle seznamu B5.

JIV Číslo IV dávky.

KAPPR Klíč postupných aproximací.

= 0 Bez použití iterační metody.

= 1 S iteracemi řízenými kriteriem ERAL a EDIF.

KAUTO Klíč automatické volby kroku.

= 0 Řízení uživatelem.

= 1 Automatické řízení.

KOUT Klíč výstupu do souboru name.OT.

= 1 Posloupnost všech aproximací.

= 2 Jen řešení.

KPRED Klíč predikce termofyzikálních vlastností.

= 0 Bez predikce.

= 1 S predikcí (může urychlit výpočet).

KQT Identifikační číslo dle seznamu B7.

KREST Klíč restartu - pouze pro nestacionární úlohu.

= 1 Nová úloha.

= 3 Pokračování v úspěšně dokončené úloze.

KSU Klíč stacionárnosti.

= 0 Nestacionární děj.

= 1 Stacionární děj.

NSTEPX Maximální počet časových kroků (pro automatické nastavování délky kroku.

NNOD Počet uzlů v síti.

NSAX Maximální počet iterací modifikovanou N-R metodou. Doporučeno 10 < NSAX < 20.

PENALTA Hodnota penalty pro spojovací prvky. Default = 10⁶.

PIVOT Minimální dovolená hodnota pivotu. Default = 10^-6.

Tepelný tok [W/m²].

Koncentrovaný tepelný tok [W].

STEP Délka kroku [s]. Při automatickém řízení délky kroku (KAUTO = 1) je STEP délka prvního kroku.

T Teplota [^oC].

T_i Řádkový vektor teplot
= T_i ... pro uzel mimo prvek semi-loof (ref A8, A9)
= T_i DT_i ... pro uzel prvku semi-loof (ref A8, A9)

T_i Teplota v i-tém uzlu výpočtové sítě [^oC]; jedná-li se o uzel prvku semi-loof (ref A8, A9), je to teplota na jeho střednicové ploše.

DT_i Rozdíl teplot na horní a spodní straně prvku semi-loof (ref A8, A9) pro i-tý uzel výpočtové sítě [^oC].

T_o Teplota okolí tělesa [^oC].

TIMS Čas od kterého začíná řešení [s]. U stacionárních úloh (KSU = 1), nebo při KREST = 3 je TIMS = 0.

TIMX Konec časového úseku [s].

TOL Tolerance chyby v jednom časovém kroku [^oC], využívaná pouze pro automatické nastavování délky kroku (při KAUTO = 1). Doporučeno 1 < TOL < 10 [^oC].

TSC Konstanta integrační metody <0,1>. TSC = 0 odpovídá explicitní metodě, TSC = 1 představuje plně implicitní schema (doporučeno).

Tepelný zdroj [W/m³].

x Globální osa [m].

x_i Diskretní hodnoty nezávisle proměnné, která patří do seznamu B5.

y Globální osa [m].

z Globální osa [m].

a Součinitel přestupu tepla [W/m²K]

b Přechodový odpor [W/m²K].

l Tepelná vodivost [W/mK].

r Hustota [kg/m³]. Viz r×c.

r×c Měrná tepelná kapacita [J/m³K].

c	Radiační konstanta [W/m²K⁴]
c	Viz r×c = měrná tepelná kapacita [J/m³K].
c₁	Konstanta obecného přestupu [W/m²K^c2+c3].
c₂	Konstanta obecného přestupu [-].
c₃	Konstanta obecného přestupu [-].
DTRUN	Elementární časový krok [s]. Délka skutečného časového kroku se zaokrouhlí na celistvý násobek DTRUN. Uvažuje se jen pro DTRUN > 10^-6 [s].
EDIF	Kritérium konvergence přírůstku teplot [^oC]: \|\|T⁽ⁱ⁾ - T^(i-1)\|\|_MAX < EDIF. Doporučeno 1 < EDIF < 5 ^oC. (uplatní se pouze při KAPPR = 1)
ERAL	Residuální kritérium konvergence \|\|Res T⁽ⁱ⁾\|\| < ERAL*\|\|T⁽ⁱ⁾\|\| .Doporučeno 0.0001 < ERAL < 0.1. (uplatní se pouze při KAPPR = 1)
INT3	Pořadové číslo integračního (časového) kroku, od kterého se má pokračovat s KREST = 3. Jestliže KREST = 1, zapíše se INT3 = 0.
IREC	Pořadové číslo záznamu v binárním souboru name.IC.
ISET	Rozlišovací číslo sady.
IV	Identifikační číslo proměnné dle seznamu B5.
JIV	Číslo IV dávky.
KAPPR	Klíč postupných aproximací.
	= 0 Bez použití iterační metody.
	= 1 S iteracemi řízenými kriteriem ERAL a EDIF.
KAUTO	Klíč automatické volby kroku.
	= 0 Řízení uživatelem.
	= 1 Automatické řízení.
KOUT	Klíč výstupu do souboru name.OT.
	= 1 Posloupnost všech aproximací.
	= 2 Jen řešení.
KPRED	Klíč predikce termofyzikálních vlastností.
	= 0 Bez predikce.
	= 1 S predikcí (může urychlit výpočet).
KQT	Identifikační číslo dle seznamu B7.
KREST	Klíč restartu - pouze pro nestacionární úlohu.
	= 1 Nová úloha.
	= 3 Pokračování v úspěšně dokončené úloze.
KSU	Klíč stacionárnosti.
	= 0 Nestacionární děj.
	= 1 Stacionární děj.
NSTEPX	Maximální počet časových kroků (pro automatické nastavování délky kroku.
NNOD	Počet uzlů v síti.
NSAX	Maximální počet iterací modifikovanou N-R metodou. Doporučeno 10 < NSAX < 20.
PENALTA	Hodnota penalty pro spojovací prvky. Default = 10⁶.
PIVOT	Minimální dovolená hodnota pivotu. Default = 10^-6.
	Tepelný tok [W/m²].
	Koncentrovaný tepelný tok [W].
STEP	Délka kroku [s]. Při automatickém řízení délky kroku (KAUTO = 1) je STEP délka prvního kroku.
T	Teplota [^oC].
T_i	Řádkový vektor teplot = T_i ... pro uzel mimo prvek semi-loof (ref A8, A9) = T_i DT_i ... pro uzel prvku semi-loof (ref A8, A9)
T_i	Teplota v i-tém uzlu výpočtové sítě [^oC]; jedná-li se o uzel prvku semi-loof (ref A8, A9), je to teplota na jeho střednicové ploše.
DT_i	Rozdíl teplot na horní a spodní straně prvku semi-loof (ref A8, A9) pro i-tý uzel výpočtové sítě [^oC].
T_o	Teplota okolí tělesa [^oC].
TIMS	Čas od kterého začíná řešení [s]. U stacionárních úloh (KSU = 1), nebo při KREST = 3 je TIMS = 0.
TIMX	Konec časového úseku [s].
TOL	Tolerance chyby v jednom časovém kroku [^oC], využívaná pouze pro automatické nastavování délky kroku (při KAUTO = 1). Doporučeno 1 < TOL < 10 [^oC].
TSC	Konstanta integrační metody <0,1>. TSC = 0 odpovídá explicitní metodě, TSC = 1 představuje plně implicitní schema (doporučeno).
	Tepelný zdroj [W/m³].
x	Globální osa [m].
x_i	Diskretní hodnoty nezávisle proměnné, která patří do seznamu B5.
y	Globální osa [m].
z	Globální osa [m].
a	Součinitel přestupu tepla [W/m²K]
b	Přechodový odpor [W/m²K].
l	Tepelná vodivost [W/mK].
r	Hustota [kg/m³]. Viz r×c.
r×c	Měrná tepelná kapacita [J/m³K].