Package for Machine Design

Finite Element Analysis in Structural Mechanics

Uživatelské nástroje

Nástroje pro tento web


cs:ref:name:ib

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}].$
cs/ref/name/ib.txt · Poslední úprava: 2023-10-27 09:44 autor: Petr Pařík