Nestacionární odezva

Schéma výpočtu

Postup výpočtu

1. Příprava výpočtu

Řešení této úlohy musí předcházet výpočet matic tuhosti prvků, přičemž platí následující podmínky:

Použít se mohou všechny prvky jako v lineární elastostatické úloze.
Nelze využít obecnou plochu symetrie nebo periodicitu.

Na pravé straně, tj. zatěžovacím stavu, v tomto případě nezáleží, proto lze v AS 1 uvažovat jen veličiny specifikující materiálové vlastnosti a posuvové okrajové podmínky.

Program: RMD2, RPD2, SRH2 (2D úloha) nebo RMD3, RPD3, SRH3 (3D úloha)
Vstupy: name.i1, name.i2, name.i3
Protokol: name.o1, name.o2, name.o3
Výstupy: binární soubory
Detaily: Lineární elastostatika, Organizace výpočtu, Referenční příručka: Vstupy

2. Výpočet matic hmotnosti

Vstupní data se zapíší do textového souboru name.iM. Program generuje konzistentní (a tedy pozitivně definitní) matice hmotnosti $\mathbf {M}$ jednotlivých prvků.

Program: HMOT (2D i 3D úloha)
Vstupy: name.iM, binární soubory z předchozího výpočtu
Protokol: name.oM
Výstupy: binární soubory
Detaily: Organizace výpočtu, Referenční příručka: Vstupy

3. Výpočet matic tlumení

Vstupní data se zapíší do textového souboru name.iC. Program generuje Rayleighovy matice proporcionálního tlumení $\mathbf {C}$ jednotlivých prvků (o typu definitnosti matic $\mathbf{C}$ nelze předem nic prohlásit).

Program: HCRE (2D i 3D úloha)
Vstupy: name.iC, binární soubory z předchozího výpočtu
Protokol: name.oC
Výstupy: binární soubory
Detaily: Organizace výpočtu, Referenční příručka: Vstupy

Poznámka
Lze počítat i s jinými typy matic tlumení; v takovém případě se však uživatel musí postarat o jejich uložení do binárního souboru name.AMP a výpočet programem HCRE se neprovádí.

4. Faktorizace matice soustavy rovnic

Vstupní data se zapíší do textového souboru name.iR. Program faktorizuje matici $$\sum(\mathbf{K}+a_0\mathbf{M}+a_1\mathbf{C}),$$ kde symbol $\sum(\dots)$ značí, že se jedná o celkové (nikoli prvkové) matice. Konzistentní matice hmotnosti $\mathbf{M}$ je pozitivně definitní a stejnou vlastnost obvykle vykazuje i matice ($a_0\mathbf{M}+a_1\mathbf{C}$). Protože zde nesledujeme kvazistatické děje, bude pak i matice $\sum(\mathbf{K}+a_0\mathbf{M}+a_1\mathbf{C})$ pozitivně definitní, a to i v případě, že $\sum(\mathbf{K})$ je jen pozitivně semidefinitní. V protokolu name.oR je vhodné zkontrolovat, zda se nevyskytuje hlášení pivotů menších než $\mathtt{PIVOT}.$

Program: HFRO (2D i 3D úloha)
Vstupy: name.iR, binární soubory z předchozího výpočtu
Protokol: name.oR
Výstupy: binární soubory
Detaily: Organizace výpočtu, Referenční příručka: Vstupy

5. Integrace pohybových rovnic

Vstupní data se zapíší do textového souboru name.iW. Algoritmem řešení je Newmarkova metoda přímé integrace. Jsou k dispozici obecné prostředky pro zadání počátečních podmínek, silového nebo kinematického buzení. Rovněž řízení výstupu výsledků umožňuje udržet jejich množství v rozumných mezích a usnadňuje jejich následné grafické zpracování. Metodu lze použít jak pro tlumené, tak i pro netlumené děje.

Program: HNEW (2D i 3D úloha)
Vstupy: name.iW, binární soubory z předchozího výpočtu
Protokol: name.oW
Výstupy: binární soubory (řešení je v souboru name.S)
Detaily: Organizace výpočtu, Referenční příručka: Vstupy

6. 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: STR2 (2D úloha) nebo STR3 (3D úloha)
Vstupy: name.i5, binární soubory z předchozího výpočtu
Protokol: name.o5
Výstupy: name.STR (volitelně), name.STB (volitelně)
Detaily: Organizace výpočtu, Referenční příručka: Vstupy