Package for Machine Design

Finite Element Analysis in Structural Mechanics

User Tools

Site Tools


en:user:problem:heat:trans

Solver control parameters for transient analysis

The control parameters of the XT2T/XT3T solver are written to the name.iB file:

IP KREST 0 KOUT INT3 NSAX NSTEPX 0 RP TIMS ERAL EDIF TOL DTRUN PIVOT PENAL

$\mathtt{KREST}$ is the restart key:
$=1$ … new computation
$=3$ … continuation of a successfully finished computation

$\mathtt{KOUT}$ is the report output key:
$=1$ … sequence of all approximations
$=2$ … solution only

$\mathtt{INT3}$ is the index of the (time) integration step from which to resume when $\mathtt{KREST}=3$. When $\mathtt{KREST}=1$, $\mathtt{INT3}=0$.

$\mathtt{NSAX}$ is the maximum number of successive approximations before initiating a new matrix factorization. The recommended value is $10<\mathtt{NSAX}<20$.

$\mathtt{NSTEPX}$ is the maximum allowed number of time steps for the entire calculation (when using automatic time step selection).

$\mathtt{TIMS}$ $[\text{s}]$ is the solution start time.1) When $\mathtt{KREST}=3$, $\mathtt{TIMS}=0$.

$\mathtt{ERAL}$ is the convergence criterion for the residual. The recommended value is $10^{-5}<\mathtt{ERAL}<10^{-2}$. It applies only when $\mathtt{KAPPR}=1$.

$\mathtt{EDIF}$ $[^\circ\text{C}]$ is the convergence criterion for the temperature increment. The recommended value is $1<\mathtt{EDIF}<10$. It applies only when $\mathtt{KAPPR}=1$.

$\mathtt{TOL}$ $[^\circ\text{C}]$ is the error tolerance within a single time step, used only for automatic time step selection when $\mathtt{KAUTO}=1$. The recommended value is $1<\mathtt{TOL}<10$.

$\mathtt{DTRUN}$ $[\text{s}]$ is the elementary time step. The actual time step length is rounded to an integer multiple of $\mathtt{DTRUN}$. It is considered only for $\mathtt{DTRUN}>10^{-6}$.

$\mathtt{PIVOT}$ is the minimum allowed pivot value during matrix factorization. The default value is $10^{-6}$.

$\mathtt{PENAL}$ is the penalty function value for connector elements of all types. The default value is $10^6$.

Pro každý zatěžovací stav, představující časový interval řešení, se musí zadat v rámci AS dávky:

AS číslo zatěžovacího stavu /$\dots$ /R TIMX STEP TSC /$\dots$

přičemž

  • $\mathtt{TIMX}$ je konec časového úseku $[\text{s}].$
  • $\mathtt{STEP}$ je délka integračního kroku $[\text{s}].$ Při automatickém řízení délky kroku ($\mathtt{KAUTO}=1$) je $\mathtt{STEP}$ výchozí délka prvního kroku.
  • $\mathtt{TSC}$ je konstanta integrační metody, $0\le\mathtt{TSC}\le1.$ $\mathtt{TSC}=0$ odpovídá explicitní metodě, $\mathtt{TSC}=1$ (doporučeno) představuje plně implicitní schéma.

Vhodné je též aktivovat AV sadu:

AV číslo sady T 6 N KAPPR KAUTO KPRED V 4*0

přiřazením v AS dávce

AS číslo zatěžovacího stavu /$\dots$ /A číslo sady /$\dots$

přičemž

  • $\mathtt{KAPPR}$ je 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}$ je klíč automatické volby integračního kroku:
    • $=0$ … řízení uživatelem
    • $=1$ … automatické řízení (doporučeno)
  • $\mathtt{KPRED}$ je klíč predikce termofyzikálních vlastností:
    • $=0$ … bez predikce
    • $=1$ … s predikcí (doporučeno při $\mathtt{KAUTO}=1$)

In the first load case, it is necessary to specify initial approximation of the temperature field using the GV set.

Each load case corresponds to as many records in the name.STR and name.TEM files as the number of integration steps performed.

FIXME

1)
Parametr $\mathtt{TIMS}$ je jen formální veličinou. Posun časové osy slouží ke snadnějšímu popisu časových závislostí.
en/user/problem/heat/trans.txt · Last modified: by Petr Pařík