Package for Machine Design

Finite Element Analysis in Structural Mechanics

Uživatelské nástroje

Nástroje pro tento web


cs:user:problem:dyn:solver

Řídící parametry řešiče

Řídící parametry řešiče HDYN se zadávají v souboru name.iN:

IP KMET KOUT 2*0 NINT KTPR NGD RP 3*0 PENAL TSTEP BETA

Metodou řešení je metoda centrálních diferencí ($\mathtt{KMET}=1$), která je podmínečně stabilní. Velikost integračního kroku je dána parametrem $\mathtt{TSTEP}.$

Význam ostatních parametrů:

  • $\mathtt{KOUT}$ aktivuje výstup vektorů řešení do binárního souboru name.PLS, kterým pak v programu STR2/STR3 odpovídá pořadové číslo $\mathtt{ILC}$ (viz name.i5). Pokud byla úloha restartována s $\mathtt{KREST}=2$ (viz Posloupnost zatěžovacích stavů), soubor name.PLS se přepíše, tj. řešení z předchozího běhu se musí zpracovat před restartem. Číslo v závorce udává počet zapsaných vektorů řešení:
    • $=0$ … kontrola vstupních dat ($0$)
    • $=1$ … po každém zatěžovacím stavu ($\mathtt{NCYC}\cdot\mathtt{NLC}$)
    • $=2$ … po každém cyklu ($\mathtt{NCYC}$)
    • $=3$ … jen výsledné řešení ($1$)
  • $\mathtt{NINT}$ specifikuje dělení integračního intervalu a souvisí s přesností řešení. Doporučuje se zadat $\mathtt{NINT}=1.$
  • $\mathtt{NGD}=1$ aktivuje jednobodovou integraci pro isoparametrické prvky. Je-li $\mathtt{NGD}=0,$ použije se řád integrace $\mathtt{NG}$ přiřazený prvkům v souboru name.i1.
  • $\mathtt{PENAL}$ se zadává u kontaktních úloh a na jeho správné volbě závisí stabilita i přesnost řešení.
  • $\mathtt{BETA}$ aktivuje speciální tvar Rayleighovy matice tlumení $\mathbf{C}=\beta\mathbf{M}.$

Počáteční podmínky

Počáteční podmínky se zadávají v souboru name.iN jednou z následujících dávek:

Varianta 1

IC číslo sady T KQT R $[X_0]$ definuje počáteční vektor. Počet hodnot musí přesně odpovídat počtu stupňů volnosti sítě.

Varianta 2

IC číslo sady T KQT I IREC přečte počáteční vektor z $\mathtt{IREC}$-tého záznamu v souboru name.SOL.

Varianta 3

IC číslo sady T KQT R $X_0$ $Y_0$ $Z_0$ definuje konstantní počáteční vektor ve směru os globálního systému. Ve 2D úlohách se $Z_0$ nezapisuje. Klíč počátečních podmínek $\mathtt{KQT}=1$ specifikuje počáteční posunutí a $\mathtt{KQT}=2$ počáteční rychlost.

Pokud není IC dávka definována, automaticky se předpokládá, že jsou počáteční podmínky homogenní.

cs/user/problem/dyn/solver.txt · Poslední úprava: 2024-11-07 13:18 autor: Petr Pařík