Obsah
name.iB
Program
XRPD
Formát
; řízení výpočtu IP KREST 0 KOUT INT3 NSAX NSTEPX KSU KSOL RP TIMS ERAL EDIF TOL DTRUN PIVAL PENAL
; nepovinný popis nezávisle proměnných IV JIV T IV V $x_1$ $x_2$ $\dots$ $x_N$
; materiálové veličiny MP ISET T 1 V $\lambda$ $\rho c$
; tepelný přechodový odpor MP ISET T 2 V $\beta$
; teploty všech uzlů sítě GV ISET T 1/6 V $[T]_0$ GV ISET T 1/6 D 4 IREC
; objemový tepelný zdroj VV ISET T 6 V $\dot w$
; přestup tepla SV ISET T 1/11 V $\alpha$ $T_o$ ; konvekcí SV ISET T 2/12 V $c$ $T_o$ ; zářením SV ISET T 3/13 V $c_1$ $c_2$ $c_3$ $T_o$ ; obecný
; tepelný tok SV ISET T 4/14 V $\dot q$
; přestup tepla na hraně prvku semi-loof LV ISET T 1/11 V $\alpha$ $T_o$ ; konvekcí LV ISET T 2/12 V $c$ $T_o$ ; zářením LV ISET T 3/13 V $c_1$ $c_2$ $c_3$ $T_o$ ; obecný
; tepelný tok na hraně prvku semi-loof LV ISET T 4/14 V $\dot q$
; teplota uzlu NV ISET T 1/11 V $T$ ; 1 složka NV ISET T 1/11 V $T$ $\Delta T$ ; 2 složky
; koncentrovaný tepelný tok NV ISET T 2/12 V $\dot q$
; řízení výpočtu AV ISET T 6 N KAPPR KAUTO KPRED V 4*0
; první zatěžovací stav (viz poznámka) AS 1
; přiřazení MP sad ␣␣/M ISET ; povinné přiřazení výchozího materiálu všem prvkům ␣␣/M ISET E $[$IE$]$
; přiřazení GV sad ␣␣/G ISET
; přiřazení VV sad ␣␣/V ISET E $[$IE$]$
; přiřazení SV sad ␣␣/S ISET E $[$IE$]$ S IS
; přiřazení LV sad ␣␣/L ISET E $[$IE$]$ L IH
; přiřazení NV sad ␣␣/N ISET N $[$IN$]$
; přiřazení AV sady ␣␣/A ISET
; řízení výpočtu ␣␣/R TIMX STEP TSC
; nepovinné další zatěžovací stavy (viz poznámka) AS 2 /$\dots$ /$\dots$ $\vdots$
; konec vstupních dat EN EN
Poznámka
V prvním zatěžovacím stavu mohou být přiřazeny všechny veličiny. Veličiny s $\mathtt{KQT}\le5$ platí ve všech zatěžovacích stavech. Ve druhém a dalších zatěžovacích stavech lze přiřadit pouze veličiny s $\mathtt{KQT}>5.$ Tyto veličiny platí pouze v daném zatěžovacím stavu.
Jsou-li uvedeny dvě hodnoty $\mathtt{KQT}$ oddělené lomítkem, první hodnota se použije pro předpis veličiny pro celý uvažovaný děj a druhá hodnota pro předpis veličiny pouze pro jeden zatěžovací stav.
Vysvětlivky
$\mathtt{KREST}$ | Klíč restartu. | |
---|---|---|
$=1$ | nové řešení | |
$=3$ | pokračování v úspěšně dokončeném řešení | |
$\mathtt{KOUT}$ | Klíč výstupu do protokolu. | |
$=1$ | posloupnost všech aproximací | |
$=2$ | jen řešení | |
$\mathtt{INT3}$ | Pořadové číslo integračního (časového) kroku, od kterého se má pokračovat při $\mathtt{KREST}=3.$ Při $\mathtt{KREST}=1$ se zapíše $\mathtt{INT3}=0.$ | |
$\mathtt{NSAX}$ | Maximální počet postupných aproximací před spuštěním nové faktorizace matice. Doporučeno $10<\mathtt{NSAX}<20.$ | |
$\mathtt{NSTEPX}$ | Maximální počet časových kroků (pro automatické nastavování délky kroku). | |
$\mathtt{KSU}$ | Klíč typu úlohy. | |
$=0$ | nestacionární děj | |
$=1$ | stacionární děj | |
$\mathtt{KSOL}$ | Klíč metody lineárního řešení. | |
$=1$ | přímý frontální řešič (výchozí) | |
$=2$ | paralelní přímý řídký řešič | |
$\mathtt{TIMS}$ | Čas, od kterého začíná řešení $[\text{s}].$ Při $\mathtt{KSU}=1$ (stacionární děj) nebo při $\mathtt{KREST}=3$ je $\mathtt{TIMS}=0.$ | |
$\mathtt{ERAL}$ | Kritérium konvergence pro reziduum, $||\operatorname{Res}\mathbf{T}^{(i)}||<\mathtt{ERAL}\cdot||\mathbf{T}^{(i)}||.$ Doporučeno $10^{-5}<\mathtt{ERAL}<10^{-2}.$ Uplatní se pouze při $\mathtt{KAPPR}=1.$ | |
$\mathtt{EDIF}$ | Kritérium konvergence přírůstku teplot $[^\circ\text{C}],$ $||\mathbf{T}^{(i)}-\mathbf{T}^{(i-1)}||_\text{MAX}<\mathtt{EDIF}.$ Doporučeno $1<\mathtt{EDIF}<5.$ Uplatní se pouze při $\mathtt{KAPPR}=1.$ | |
$\mathtt{TOL}$ | Tolerance chyby v jednom časovém kroku $[^\circ\text{C}],$ využívána pouze pro automatické nastavování délky kroku při $\mathtt{KAUTO}=1.$ Doporučeno $1<\mathtt{TOL}<10.$ | |
$\mathtt{DTRUN}$ | Elementární časový krok $[\text{s}].$ Délka skutečného časového kroku se zaokrouhlí na celistvý násobek $\mathtt{DTRUN}.$ Uvažuje se jen pro $\mathtt{DTRUN}>10^{-6}.$ | |
$\mathtt{PIVAL}$ | Minimální dovolená hodnota pivotu při faktorizaci matice. Výchozí hodnota je $10^{-6}.$ | |
$\mathtt{PENAL}$ | Hodnota pokutové funkce pro spojovací prvky všech typů. Výchozí hodnota je $10^6.$ | |
$c$ | Měrná tepelná kapacita $[\text{J}/\text{kgK}],$ zadává se jako součin $\rho c.$ | |
$c$ | Radiační konstanta $[\text{W}/\text{m}^2\text{K}^4]$ pro výpočet tepelného toku $\dot q=c(T^4-T_o^4).$ | |
$c_1$ | Konstanta obecného přestupu $[\text{W}/\text{m}^2\text{K}^{c_2+c_3}]$ pro výpočet tepelného toku $\dot q=c_1\left(T^{c_2}-T_o^{c_2}\right)^{c_3}.$ | |
$c_2$ | Konstanta obecného přestupu $[1].$ | |
$c_3$ | Konstanta obecného přestupu $[1].$ | |
$[\mathtt{IE}]$ | Seznam čísel prvků. | |
$\mathtt{IH}$ | Lokální číslo hrany prvku. | |
$[\mathtt{IN}]$ | Seznam čísel uzlů. | |
$\mathtt{IREC}$ | Pořadové číslo záznamu v binárním souboru name.TEM přejmenovaném na name.TIC . |
|
$\mathtt{IS}$ | Lokální číslo stěny prvku. | |
$\mathtt{ISET}$ | Rozlišovací číslo sady. | |
$\mathtt{IV}$ | Identifikační číslo proměnné. | |
$\mathtt{JIV}$ | Číslo IV dávky. |
|
$\mathtt{KAPPR}$ | Klíč postupných aproximací. | |
$=0$ | bez použití iterační metody | |
$=1$ | s iteracemi řízenými kriterii $\mathtt{ERAL}$ a $\mathtt{EDIF}$ (doporučeno) | |
$\mathtt{KAUTO}$ | Klíč automatické volby kroku. | |
$=0$ | řízení uživatelem | |
$=1$ | automatické řízení (doporučeno) | |
$\mathtt{KPRED}$ | Klíč predikce termofyzikálních vlastností. | |
$=0$ | bez predikce | |
$=1$ | s predikcí (může urychlit výpočet, doporučeno při $\mathtt{KAUTO}=1$) | |
$\mathtt{KQT}$ | Identifikační číslo veličiny. | |
$\dot q$ | Tepelný tok $[\text{W}/\text{m}^2].$ | |
$\dot q$ | Koncentrovaný tepelný tok $[\text{W}].$ | |
$\mathtt{STEP}$ | Délka integračního kroku $[\text{s}].$ Při automatickém řízení délky kroku ($\mathtt{KAUTO}=1$) je $\mathtt{STEP}$ délka prvního kroku. | |
$T$ | Uzlová teplota $[^\circ\text{C}].$ | |
$\Delta T$ | Rozdíl teplot na horní a spodní straně prvku semi-loof $[^\circ\text{C}].$ | |
$[T]_0$ | Globálního teplotní pole. Počet složek (délka vektoru) musí přesně odpovídat počtu stupňů volnosti sítě. Jsou-li v síti pouze uzly s jedním stupněm volnosti, je počet stupňů volnosti sítě roven počtu uzlů v síti. | |
$T_o$ | Teplota okolí $[^\circ\text{C}].$ | |
$\mathtt{TIMX}$ | Konec časového úseku $[\text{s}].$ | |
$\mathtt{TSC}$ | Konstanta integrační metody, $0\le\mathtt{TSC}\le1.$ $\mathtt{TSC}=0$ odpovídá explicitní metodě, $\mathtt{TSC}=1$ představuje plně implicitní schéma (doporučeno). | |
$\dot w$ | Výkon tepelného zdroje $[\text{W}/\text{m}^3].$ | |
$x_i$ | Diskrétní hodnoty nezávisle proměnné. | |
$\alpha$ | Součinitel přestupu tepla $[\text{W}/\text{m}^2\text{K}]$ pro výpočet tepelného toku $\dot q=\alpha(T-T_o).$ | |
$\beta$ | Tepelný přechodový odpor $[\text{W}/\text{m}^2\text{K}]$ pro výpočet tepelného toku $\dot q=\beta\Delta T.$ | |
$\lambda$ | Tepelná vodivost $[\text{W}/\text{mK}].$ | |
$\rho$ | Hustota $[\text{kg}/\text{m}^3],$ zadává se jako součin $\rho c.$ | |
$\rho c$ | Měrná tepelná kapacita na jednotku objemu $[\text{J}/\text{m}^3\text{K}].$ |