Package for Machine Design

Finite Element Analysis in Structural Mechanics

Uživatelské nástroje

Nástroje pro tento web


cs:user:problem:heat:nonlin

Řídící parametry řešiče pro stacionární nelineární úlohu

Řídící parametry řešiče XT2S/XT3S se zapisují do souboru name.iB:

IP 1 0 KOUT 0 NSAX NSTEPX 1 RP 0 ERAL EDIF 0 0 PIVOT PENAL

Význam parametrů

  • $\mathtt{KOUT}$ je klíč výstupu do protokolu:
    • $=1$ … posloupnost všech aproximací
    • $=2$ … jen řešení
  • $\mathtt{NSAX}$ je maximální počet postupných aproximací před spuštěním nové faktorizace matice. Doporučená hodnota je $10<\mathtt{NSAX}<20.$
  • Je-li $\mathtt{NSTEPX}<0,$ nepoužije se akcelerace line search.
  • $\mathtt{EDIF}$ je kritérium konvergence přírůstku teplot $[^\circ\text{C}].$ Doporučená hodnota je $1<\mathtt{EDIF}<10.$ Uplatní se pouze při $\mathtt{KAPPR}=1.$
  • $\mathtt{ERAL}$ je kritérium konvergence pro reziduum. Doporučená hodnota je $10^{-5}<\mathtt{ERAL}<10^{-2}.$ Uplatní se pouze při $\mathtt{KAPPR}=1.$
  • $\mathtt{PIVOT}$ je minimální dovolená hodnota pivotu při faktorizaci matice. Výchozí hodnota je $10^{-6}.$
  • $\mathtt{PENAL}$ je hodnota pokutové funkce pro spojovací prvky všech typů. Výchozí hodnota je $10^6.$

Přiřazení v AS dávce /R TIMX STEP TSC/ se nepoužije. Vhodné je však aktivovat AV sadu s klíčem $\mathtt{KAPPR}=1$:

AV 1 T 6 N 1 0 0 V 4*0

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

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

Každému zatěžovacímu stavu přísluší jeden záznam v souboru name.STR a name.TEM, který obsahuje stacionární teplotní pole.

Komentář

U nelineární úlohy je nutné vyjádřit všechny veličiny (okrajové podmínky a materiálové vlastnosti) pro výsledné uzlové teploty, které však nejsou předem známy. Z tohoto důvodu se řešení koriguje modifikovanou, příp. akcelerovanou ($\mathtt{NSTEPX}\ge0$) Newton-Raphsonovou metodou. Kritériem konvergence je malá změna teplot ve dvou po sobě následujících iteracích a zároveň velikost rezidua, tedy $$\begin{array}{lll} \max|\mathbf{T}^{(i+1)}-\mathbf{T}^{(i)}| < \mathtt{EDIF} & \land & ||\operatorname{Res}\mathbf{T}^{(i)}|| < \mathtt{ERAL}\,||\mathbf{T}^{(i)}||. \end{array}$$

Pokud je konvergence pomalá, je možné sestavit novou matici soustavy, vypočtenou pro momentální aproximaci teplot. Tento postup však vyžaduje novou faktorizaci, která je časově velmi náročná. Parametr $\mathtt{NSAX}$ určuje počet postupných aproximací, po kterých se vytvoří nová matice soustavy, nedošlo-li doposud ke konvergenci.

Důležitým parametrem je klíč $\mathtt{KAPPR},$ kterým se spouští proces korekcí řešení. Pokud není $\mathtt{KAPPR}=1$ nebo chybí AV dávka (nebyla přiřazena), kontrola konvergence se neprovádí a výpočet se ukončí po první iteraci. Následkem je stacionární lineární řešení s veličinami vyjádřenými pro výchozí teplotu.

Pomocí klíče $\mathtt{KOUT}$ lze sledovat průběh konvergence. Jestliže jsou potřebné jen výsledky, nastaví se $\mathtt{KOUT}=2.$

Při použití iterační metody (AV sada s $\mathtt{KAPPR}=1$) je nutné zadat výchozí aproximaci teplotního pole pomocí GV sady.

cs/user/problem/heat/nonlin.txt · Poslední úprava: 2024-11-07 13:16 autor: Petr Pařík