Obsah
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) / 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 / 3D úloha) |
---|---|
Vstupy | name.iM , binární soubory z předchozích kroků |
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 / 3D úloha) |
---|---|
Vstupy | name.iC , binární soubory z předchozích kroků |
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 souboruname.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 / 3D úloha) |
---|---|
Vstupy | name.iR , binární soubory z předchozích kroků |
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 / 3D úloha) |
---|---|
Vstupy | name.iW , binární soubory z předchozích kroků |
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) / STR3 (3D úloha) |
---|---|
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 |