Lineární stabilita

Schéma výpočtu

Postup výpočtu

1. Příprava výpočtu

Řešení této úlohy musí předcházet řešení lineární elastostatické úlohy, přičemž platí následující podmínky:

Použít se mohou jen následující prvky: pětistěn, šestistěn, semi-loof a příslušné spojovací prvky: hranové spojení, stěnové spojení.
Nelze využít obecnou symetrii a periodicitu.

Program	RMD3, RPD3, SRH3, FEFS
Vstupy	`name.i1`, `name.i2`, `name.i3`, `name.i4`
Protokol	`name.o1`, `name.o2`, `name.o3`, `name.o4`
Výstupy	binární soubory (řešení je v souboru `name.SOL`)
Detaily	Lineární elastostatika, Organizace výpočtu / Referenční příručka: Vstupy

2. Výpočet matic počáteční napjatosti

Vstupní data se zapíší do textového souboru name.iG. Program generuje matice počáteční napjatosti $\mathbf{G}$ jednotlivých prvků pro zvolený zatěžovací stav ze souboru name.SOL. Matice $\mathbf{G}$ jsou symetrické; mohou být regulární, singulární (s různou nulitou), definitní nebo indefinitní.

Program	GEO3
Vstupy	`name.iG`, binární soubory z předchozích kroků
Protokol	`name.oG`
Výstupy	binární soubory
Detaily	Organizace výpočtu / Referenční příručka: Vstupy

3. Řešení soustavy rovnic

Vstupní data se zapíší do textového souboru name.iE, přičemž klíč $\mathtt{KEVP}=1.$ Program vypočte metodou iterace podprostoru $\mathtt{NROOT}$ vlastních párů (vlastních vektorů a vlastních čísel) obecného vlastního problému $$\sum(\mathbf{K}-\lambda_i\mathbf{G})\mathbf{v}_i=0, \quad i=1,\dots,\mathtt{NROOT},$$ kde $\mathbf{v}_i$ a $\lambda_i$ jsou $i$-tým vlastním vektorem a vlastním číslem a $\sum(\dots)$ značí, že se jedná o celkové (nikoli prvkové) matice.

Program	HEIG
Vstupy	`name.iE`, binární soubory z předchozích kroků
Protokol	`name.oE`
Výstupy	binární soubory (řešení je v souboru `name.EIG`)
Detaily	Organizace výpočtu / Referenční příručka: Vstupy

4. Normalizace vlastních vektorů

Vstupní data se zapíší do textového souboru name.iS. Program ukládá pro každý z $\mathtt{NROOT}$ vlastních párů, vypočtených v předchozím kroku, dva záznamy do binárního souboru name.S. Liché záznamy obsahují normalizované vlastní vektory $\mathbf{v}_i,$ sudé záznamy jim příslušející vlastní čísla $\lambda_i.$

Program	STAB
Vstupy	`name.iS`, binární soubory z předchozích kroků
Protokol	`name.oS`
Výstupy	binární soubory (řešení je v souboru `name.S`)
Detaily	Organizace výpočtu / Referenční příručka: Vstupy

5. Výpočet deformací a napětí

Vstupní data se zapíší do textového souboru name.i5, přičemž klíč typu úlohy $\mathtt{KPROB}=1.$

Program	STR3
Vstupy	`name.i5`, binární soubory z předchozích kroků
Protokol	`name.o5`
Výstupy	`name.STR` (volitelně), `name.STB` (volitelně)
Detaily	Organizace výpočtu / Referenční příručka: Vstupy