Feb MAR APR » 5 2010 2011 2014 » 2 captures 5 Mar 11 - 5 Apr 14

Close Help

B1 definice veličin v systému PMD
B2 konstantní veličina
B3 polynomiální závislost
B4 závislost daná tabulkou
B5 seznam proměnných
B6 seznam veličin
B7 seznam veličin pro výpočty vedení tepla

veličina: Pojmem veličina se v systému PMD rozumí množina fyzikálních dat, která se zadává formou vektoru [v₁,v₂, ...,v_N].
Příkladem mohou být materiálové vlastnosti [E,a,n,r,s_Y,Q_Y,_c,F].

závislost: Veličiny v systému PMD mohou záviset na souřadnicích, času, teplotě a dalších proměnných. Každá složka v_n přitom může být nezávislou funkcí až čtyř proměnných, tj. v_n = v_n(x₁,x₂,x₃,x₄). Seznam proměnných, charakterizovaných identifikačním číslem IV, je uveden v příloze B5.

sada: Veličina a její případná závislost se popisují sadou. Sada je uvedena klíčovými písmeny XX a číslem KQT, jež specifikují typ veličiny dle seznamu B6 nebo B7. Každá sada je dále označena rozlišovacím číslem ISET, což umožňuje zavést několik popisů téže fyzikální veličiny, např. různých materiálových vlastností.

přiřazení: Veličiny se definují nezávisle na výpočtové síti. Teprve poté se přiřazují jednotlivým objektům, např. elementům, plochám atd. V přiřazení je nutno provést odkaz na typ veličiny XX (uvádí se jen první písmeno X) a číslo sady ISET. Dovolené tvary přiřazení jednotlivých veličin jsou zahrnuty v přílohách B6 a B7.

Konstantní veličina [v₁,v₂, ...,v_N] se zadává sadou:
XX ISET T KQT V v₁  v₂  ...  v_N

XX           klíčová písmena dle seznamu B6 nebo B7
ISET        rozlišovací číslo sady
KQT        identifikační číslo dle seznamu B6 nebo B7

Uvažujme veličinu [v₁,v₂, ...,v_N], kde každá složka je polynomem v₁=P₁(x₁,x₂,x₃,x₄), v₂=P₂(x₁,x₂,x₃,x₄), ... , v_N=P_N(x₁,x₂,x₃,x₄).

Přípustné tvary polynomů jsou:

P(x₁) = a₁ + a₂x₁ + a₃x₁² + a₄x₁³

P(x₁,x₂) = a₁ + a₂x₁ + a₃x₂ + a₄x₁² + a₅x₂² + a₆x₁x₂ +
+ a₇x₁³ + a₈x₂³ + a₉x₁²x₂ + a₁₀x₁x₂²

P(x₁,x₂,x₃) = a₁ + a₂x₁ + a₃x₂ + a₄x₃ + a₅x₁² + a₆x₂² +
+ a₇x₃² + a₈x₁x₂ + a₉x₁x₃ + a₁₀x₂x₃ +
+ a₁₁x₁³ + a₁₂x₂³ + a₁₃x₃³ + a₁₄x₁²x₂ +
+ a₁₅x₁²x₃ + a₁₆x₂²x₁ + a₁₇x₂²x₃

P(x₁,x₂,x₃,x₄) = a₁ + a₂x₁ + a₃x₂ + a₄x₃ + a₅x₄ +
+ a₆x₁² + a₇x₂² + a₈x₃² + a₉x₄² + a₁₀x₁x₂ +
+ a₁₁x₁x₃ + a₁₂x₁x₄ + a₁₃x₂x₃ +
+ a₁₄x₂x₄ + a₁₅x₃x₄ + a₁₆x₁³ + a₁₇x₂³ +
+ a₁₈x₃³ + a₁₉x₄³

Každé složce v_n tak přísluší vektor koeficientů [a]_n.
Fyzikální význam proměnných x₁,x₂,x₃,x₄ se specifikuje vektorem identifikačních čísel [IV] = [IV₁,IV₂,IV₃,IV₄] dle B5.

Veličina [v₁,v₂, ...,v_N] se zadává sadou:
XX ISET  T -KQT I  [IV]  V  [a]₁ ....   V  [a]_N

XX       klíčová písmena dle seznamu B6 nebo B7
ISET     rozlišovací číslo sady
KQT     identifikační číslo dle seznamu B6 nebo B7 (zapsáno záporně)

Uvažujme veličinu [v₁,v₂, ...,v_N], kde každá složka je funkcí v₁=v₁(x₁,x₂,x₃,x₄), v₂=v₂(x₁,x₂,x₃,x₄), ... , v_N=v_N(x₁,x₂,x₃,x₄).
Předpokládejme, že pro každou nezávisle proměnnou je zadán vektor diskretních hodnot [x]₁, [x]₂, [x]₃, [x]₄ a v těchto bodech jsou známy funkční hodnoty v_n^ijkl = v_n(x₁ⁱ,x₂^j,x₃^k,x₄^l).
Funkční hodnoty v_n^ijkl lze sestavit do vektorů
[v]_n = [v_n¹¹¹¹, v_n¹¹¹², v_n¹¹¹³, ..... , v_n¹¹²¹, v_n¹¹²², v_n¹¹²³, ...].

Nejprve se popíší hodnoty nezávisle proměnných [x]_IV dávkami
IV  JIV  T IV V  [x]_IV

JIV        číslo IV dávky
IV          identifikační číslo proměnné dle B5

Veličina [v₁,v₂, ...,v_N] se zadává sadou:
XX  ISET  T KQT   I  seznam  JIV  V [v]₁  ....  V [v]_N

XX       klíčová písmena dle seznamu B6 nebo B7
ISET     rozlišovací číslo sady
KQT     identifikační číslo dle seznamu B6 nebo B7

poznámka: Seznam čísel JIV specifikuje typ a pořadí proměnných x₁,x₂,x₃,x₄ odkazem na IV dávky.
Závisle proměnná pro argument mimo rozsah zadaný IV-dávkou je nahrazena hodnotou, jíž nabývá pro nejblíže definovaný argument (extrapolace konstantou). V případě závislosti na více argumentech platí obdobné pravidlo.

IV = 1      x [m] globální souřadnice 'x'
IV = 2      y [m] globální souřadnice 'y'
IV = 3      z [m] globální souřadnice 'z'
IV = 4      t [s] čas
IV = 5      T [^oC] teplota
IV = 6      DT [^oC] teplotní rozdíl po tloušťce prvku semi-loof
IV = 7 [-] efektivní plastická deformace
IV = 8      W_p [J/m³] plastická práce
IV = 9     s_e [Pa] efektivní napětí HMH
IV = 10 [-] efektivní creepová deformace

poznámka: Pokud byl použit zobecněný konstitutivní model, nahradí se první dvě proměnné následujícím způsobem:
IV = 1     s_m [Pa] střední napětí
IV = 2     cos(3Q) [-] podobnostní parametr (Lode)

MP  /M  ISET  /M ISET   E  seznam prvků/
        KQT=1 Materiálové vlastnosti v jednotkách SI.
        [E,a,n,r,s_Y,Q_Y,,F] má rozměr [1/h]

GV  /G  ISET /
       KQT=1 Posunutí všech uzlů sítě.
       [u₁,v₁,w₁, ... ,u_NNOD,v_NNOD,w_NNOD] [m, rad]
       KQT=6 Teploty všech uzlů sítě.
       [T₁, ... ,T_NNOD] [^oC]

VV  /V  ISET  /V ISET   E  seznam prvků /
       KQT=6 Objemové síly.
       [Q_x,Q_y,Q_z] [N/m³]

SV  /S  ISET  E  seznam prvků  S  č. plochy/
       KQT=2     Winklerův podklad v lokálním systému plochy prvku.
                         [K_n,K_t] [Pa/m] masivní prvky
                         [K_x,K_h,K_z] [Pa/m] prvky semi-loof
       KQT=3     Winklerův podklad v globálním systému.
                         [K_x,K_y,K_z] [Pa/m]
       KQT=6     Plošné zatížení ve směru normály ploch prvků.
                        [q_n] [Pa]
       KQT=9     Plošné zatížení (napětí) v globálním systému.
                        [q_x,q_y,q_z] [Pa]

LV  /L  ISET  E  seznam prvků   Lč. hrany/
       KQT=2    Pružné podepření hrany prvku semi-loof v lokálním systému A10.
                        [K_xh,K_yh,K_zh,C_yh] [N/m², Nm/rad×m]
       KQT=6    Liniové zatížení hrany prvku semi-loof v lokálním systému A10.
                        [l_xh,l_yh,l_zh,m_yh] [N/m, Nm/m]
       KQT=9    Liniové zatížení hrany v globálním systému.
                        [l_x,l_y,l_z] [N/m]
 

NV  /N  ISET /N  ISET   E  seznam prvků/  N  ISET   N seznam uzlů/
       KQT=1    Předpis složek posunutí uzlu.
                        [u,v,w] [m] DOF=3
                        [u,v,w,a,b] [m, rad] DOF=5
                        [u,v,w,a,b,g] [m, rad] DOF=6
       KQT=2    Pružina ve směru zadaném cosiny.
                        [k_n,k_t,cos(x),cos(y),cos(z)] [N/m]
       KQT=3    Pružina v globálním systému.
                        [k_x,k_y,k_z] [N/m]
       KQT=6    Osamělá síla v globálním systému.
                        [F_x,F_y,F_z] [N] DOF=3
                        [F_x,F_y,F_z,M_a,M_b] [N, Nm] DOF=5
                        [F_x,F_y,F_z,M_x,M_y,M_z] [N, Nm] DOF=6
 

poznámka: Pokud již existuje soubor name.SOL nebo soubor name.TEM, může mít GV sada alternativní tvar:
GV   ISET  T KQT   D  IFILE   IREC
   ISET            rozlišovací číslo sady
   KQT =1/6    pro popis posuvů/teplot
   IFILE          identifikační číslo souboru
                       =12/4 pro čtení posuvů/teplot
   IREC           číslo zatěžovacího stavu v souboru
                      name.SOL pro čtení posuvů
                      name.TEM pro čtení teplot

MP  /M  ISET  /M ISET  E  seznam prvků/
        KQT=1     Materiálové vlastnosti v jednotkách SI.
                          [l,r×c]  [W/mK, J/m³K]
        KQT=2     Přechodový odpor = bDT.
                          [b] [W/m²K]

GV  /G  ISET /
       KQT=1      Teploty všech uzlů sítě - počáteční podmínka.
                          [T₁, ... ,T_NNOD] [^oC]
       KQT=6      Teploty všech uzlů sítě - výchozí aproximace.
                          [T₁, ... ,T_NNOD] [^oC]

VV  /V ISET   /V ISET  E  seznam prvků/
       KQT=6  Objemový tepelný zdroj.
[W/m³]

SV  /S  ISET  E  seznam prvků  S č. plochy/
      KQT=1/11    Přestup tepla konvekcí.
                             [a,T_o] [W/m²K, ^oC]
      KQT=2/12     Přenos tepla zářením .
                             [c,T_o] [W/m²K⁴, ^oC]
      KQT=3/13     Obecný přestup tepla .
                            [c₁,c₂,c₃,T_o] [W/m²K^c2+c3, -, -, ^oC]
      KQT=4/14     Tepelný tok.
[W/m²]

LV  /L  ISET E   seznam prvků  Lč. hrany/
       KQT=1/11     Přestup tepla konvekcí.
                              [a,T_o] [W/m²K, ^oC]
       KQT=2/12     Přenos tepla zářením .
                              [c,T_o] [W/m²K⁴, ^oC]
       KQT=3/13     Obecný přestup tepla .
                              [c₁,c₂,c₃,T_o] [W/m²K^c2+c3, -, -, ^oC]
       KQT=4/14     Tepelný tok.
[W/m²]

NV  /N  ISET  /N ISET  E  seznam prvků  /N ISET  N seznam uzlů/
       KQT=1/11    Teplota v uzlu.
                            [T] [^oC] DOF=1
                            [T,DT] [^oC] DOF=2
       KQT=2/12   Koncentrovaný tepelný tok.
[W]

AV   /AISET/
        KQT=6 Řízení výpočtu.
        Za klíčovým písmenem N     [KAPPR,KAUTO,KPRED] [0/1]
        Za klíčovým písmenem V     4*0

poznámka: Pokud již existuje soubor name.TEM, lze tento soubor přejmenovat na name.IC a GV sada může mít alternativní tvar:

GV  ISET  T KQT  D4  IREC
   ISET           rozlišovací číslo sady
   KQT =1/6   pro popis počátečních podmínek/výchozí aproximace.
   IREC          číslo zatěžovacího stavu v souboru name.IC.

poznámka: Okrajové podmínky je možno zadat dvěma hodnotami KQT. Pokud KQT < 6, podmínka platí pro celý uvažovaný děj a je nutno ji přiřadit v AS 1. Je-li KQT > 6, platí podmínka jen v tom zatěžovacím stavu, ve kterém byla přiřazena.

poznámka:AV sada platí v rámci zatěžovacího stavu, ve kterém byla přiřazena. Přepínače KAPPR, KAUTO a KPRED spouštějí postupné aproximace, automatické řízení délky kroku a predikci fyzikálních parametrů. Tyto tři přepínače platí v celém zatěžovacím stavu.