Package for Machine Design

Finite Element Analysis in Structural Mechanics

Uživatelské nástroje

Nástroje pro tento web


cs:ref:name:i2

name.i2

Program

RPD2 RPD3

Formát

; řízení programu IP KREST

; nepovinný popis nezávisle proměnných IV JIV T IV V $x_1$ $x_2$ $\dots$ $x_N$

; materiálové veličiny MP ISET T 1 V $E$ $\alpha$ $\nu$ $\rho$ $\sigma_Y$ $Q_Y$ ${\dot\varepsilon}_c$ $\Phi$

; posunutí všech uzlů sítě GV ISET T 1 V $[u]_1$ $\dots$ $[u]_\mathtt{NNOD}$ GV ISET T 1 D 12 IREC

; teploty všech uzlů sítě GV ISET T 6 V $[T]_1$ $\dots$ $[T]_\mathtt{NNOD}$ GV ISET T 6 D 4 IREC

; objemová síla VV ISET T 6 V $Q_x$ $Q_y$ $Q_z$

; Winklerův podklad SV ISET T 2 V $K_n$ $K_t$ ; v lokálním systému stěny isoparametrického prvku SV ISET T 2 V $K_\xi$ $K_\eta$ $K_\zeta$ ; v lokálním systému stěny prvku semi-loof SV ISET T 3 V $K_x$ $K_y$ $K_z$ ; v globálním systému

; plošné zatížení SV ISET T 6 V $q_n$ ; ve směru normály stěny prvku SV ISET T 9 V $q_x$ $q_y$ $q_z$ ; v globálním systému

; kontaktní plochy (viz poznámka 2) SV ISET T 10 V 0 ; plocha A SV ISET T 11 V 0 ; plocha B

; pružné podepření hrany LV ISET T 2 V $K_{x_h}$ $K_{y_h}$ $K_{z_h}$ $C_{y_h}$ ; v lokálním systému hrany prvku semi-loof

; liniové zatížení LV ISET T 6 V $l_{x_h}$ $l_{y_h}$ $l_{z_h}$ $m_{y_h}$ ; v lokálním systému hrany prvku semi-loof LV ISET T 9 V $l_x$ $l_y$ $l_z$ ; v globálním systému

; posunutí uzlu NV ISET T 1 V $u$ $v$ $w$ ; 3 složky NV ISET T 1 V $u$ $v$ $w$ $\alpha$ $\beta$ ; 5 složek NV ISET T 1 V $u$ $v$ $w$ $\varphi_x$ $\varphi_y$ $\varphi_z$ ; 6 složek NV ISET T 1 C $[$IC$]$ V $[u]$ ; pouze vybrané složky

; pružina NV ISET T 2 N IN V $k_n$ $k_t$ $\cos(n,x)$ $\cos(n,y)$ $\cos(n,z)$ ; v obecném směru NV ISET T 3 N IN V $k_x$ $k_y$ $k_z$ ; v globálním systému NV ISET T 4 N IN V $k_{11}$ $k_{12}$ $k_{22}$ $k_{13}$ $k_{23}$ $k_{33}$ $\dots$ $k_{1m}$ $k_{2m}$ $\dots$ $k_{mm}$ ; symetrická matice tuhosti v globálním systému

; osamělá síla NV ISET T 6 V $F_x$ $F_y$ $F_z$ ; 3 složky NV ISET T 6 V $F_x$ $F_y$ $F_z$ $M_\alpha$ $M_\beta$ ; 5 složek NV ISET T 6 V $F_x$ $F_y$ $F_z$ $M_x$ $M_y$ $M_z$ ; 6 složek

; první zatěžovací stav (viz poznámka 1) AS 1

; přiřazení MP sad ␣␣/M ISET ; povinné přiřazení výchozího materiálu všem prvkům ␣␣/M ISET E $[$IE$]$

; předpis nulových posunutí ␣␣/B 0 N $[$IN$]$ ; pro všechny složky ␣␣/B 0 C $[$IC$]$ N $[$IN$]$ ; pouze pro vybrané složky

; přiřazení GV sad ␣␣/G ISET

; přiřazení VV sad ␣␣/V ISET E $[$IE$]$

; přiřazení SV sad ␣␣/S ISET E $[$IE$]$ S IS

; přiřazení LV sad ␣␣/L ISET E $[$IE$]$ L IH

; přiřazení NV sad ␣␣/N ISET N $[$IN$]$ ␣␣/N ISET E $[$IE$]$ ; pouze pro přiřazení uzlových pružin

; předpis konstant ␣␣/R $R_m$ $T_o$ $T_w$ $\varepsilon_{z0}$

; nepovinné další zatěžovací stavy (viz poznámka 1) AS 2 /$\dots$ /$\dots$ $\vdots$

; konec vstupních dat EN EN

Poznámka 1
V prvním zatěžovacím stavu mohou být přiřazeny všechny veličiny. Veličiny s $\mathtt{KQT}\le5$ platí ve všech zatěžovacích stavech. Ve druhém a dalších zatěžovacích stavech lze přiřadit pouze veličiny s $\mathtt{KQT}>5.$ Tyto veličiny platí pouze v daném zatěžovacím stavu.
Poznámka 2
Kontaktní dvojice (může jich být více) se zadávají nejprve sudou (plocha A) a pak lichou (plocha B) hodnotou $\mathtt{KQT}$ (např. 10 a 11, 12 a 13, 102 a 103). Vnější normály ploch v kontaktní dvojici musí být opačně orientované!

Vysvětlivky

$\mathtt{KREST}$Klíč restartu.
$=1$nová úloha
$=2$dodatečné zadání zatěžovacích stavů
$\cos(n,x)$ $\cos(n,y)$ $\cos(n,z)$směrové kosiny osy pružiny
$C_{y_h}$Momentová tuhost podepření hrany prvku semi-loof $[\text{Nm}/\text{rad}\!\cdot\!\text{m}].$
$E$Modul pružnosti $[\text{Pa}].$
$F_x,F_y,F_z$Složky síly ve směru globálních os.
$[\mathtt{IC}]$Seznam lokálních čísel složek zobecněných posunutí v uzlu. Možné hodnoty jsou: $1\equiv u,$ $2\equiv v,$ $3\equiv w,$ $4\equiv\alpha\equiv\varphi_x,$ $5\equiv\beta\equiv\varphi_y$ a $6\equiv\varphi_z.$
$[\mathtt{IE}]$Seznam čísel prvků.
$\mathtt{IH}$Lokální číslo hrany prvku.
$\mathtt{IN}$Číslo uzlu.
$[\mathtt{IN}]$Seznam čísel uzlů.
$\mathtt{IREC}$Pořadové číslo záznamu v binárním souboru name.SOL (pole posunutí) nebo name.TEM (pole teplot).
$\mathtt{IS}$Lokální číslo stěny prvku.
$\mathtt{ISET}$Rozlišovací číslo sady.
$\mathtt{IV}$Identifikační číslo proměnné.
$\mathtt{JIV}$Číslo IV dávky.
$k_{11},\dots,k_{mm}$Prvky symetrické matice tuhosti pružiny (horní trojúhelník zapsaný po sloupcích). Počet prvků musí být roven $m(m+1)/2,$ kde $m$ je počet stupňů volnosti uzlu, ke kterému je pružina připojena.
$k_n,k_t$Osová a příčná tuhost pružiny $[\text{N}/\text{m}].$
$K_n,K_t$Normálová a tangenciální tuhost Winklerova podkladu $[\text{Pa}/\text{m}].$
$\mathtt{KQT}$Identifikační číslo veličiny.
$k_x,k_y,k_z$Tuhosti pružiny ve směru globálních os $[\text{N}/\text{m}].$
$K_x,K_y,K_z$Tuhosti Winklerova podkladu ve směru globálních os $[\text{Pa}/\text{m}].$
$K_{x_h},K_{y_h},K_{z_h}$Tuhosti podepření hrany prvku semi-loof ve směru lokálních os hrany $[\text{N}/\text{m}^2].$
$K_\xi,K_\eta,K_\zeta$Tuhosti Winklerova podkladu ve směru lokálních os prvku semi-loof $[\text{Pa}/\text{m}].$
$l_x,l_y,l_z$Liniové zatížení ve směru globálních os $[\text{N}/\text{m}].$
$l_{x_h},l_{y_h},l_{z_h}$Liniové zatížení hrany prvku semi-loof ve směru lokálních os hrany $[\text{N}/\text{m}].$
$M_x,M_y,M_z$Složky momentu ve směru globálních os v uzlu nosníkového prvku $[\text{Nm}].$
$m_{y_h}$Liniový moment na hranu prvku semi-loof ve směru lokální osy hrany $y_h$ $[\text{Nm}/\text{m}].$
$M_\alpha,M_\beta$Složky momentu na hranu prvku semi-loof $[\text{Nm}].$
$q_n$Plošné zatížení ve směru vnější normály stěny prvku $[\text{Pa}].$ Tlak se zadává záporně.
$q_x,q_y,q_z$Složky plošného zatížení ve směru globálních os $[\text{Pa}].$
$Q_x,Q_y,Q_z$Složky objemové síly ve směru globálních os $[\text{N}/\text{m}^3].$
$Q_Y$Kinematická složka zpevnění $[\text{Pa}].$
$R_m$Otáčky $[1/\text{min}].$ Osou rotace je osa $z.$
$[T]_1,\dots,[T]_\mathtt{NNOD}$Globální teplotní pole. Počet složek (délka vektoru) musí přesně odpovídat počtu teplotních stupňů volnosti sítě. Jsou-li v síti pouze uzly s jedním stupněm volnosti, je počet stupňů volnosti sítě roven počtu uzlů v síti.
$T_o$Výchozí teplota (teplota okolí) $[^\circ\text{C}].$
$T_w$Výsledná teplota (pracovní teplota) $[^\circ\text{C}].$
$u,v,w$Posunutí uzlu ve směru globálních os $[\text{m}].$
$[u]$Předepsané složky zobecněných posunutí v uzlu. Pořadí hodnot složek musí odpovídat pořadí čísel složek v příslušném seznamu $[\mathtt{IC}].$
$[u]_1,\dots,[u]_\mathtt{NNOD}$Globální pole zobecněných posunutí. Počet složek (délka vektoru) musí přesně odpovídat počtu stupňů volnosti sítě. Jsou-li v síti pouze uzly s translačními stupni volnosti, je počet stupňů volnosti sítě roven počtu uzlů v síti $\times$ dimenze úlohy.
$x_1,\dots,x_N$Diskrétní hodnoty nezávisle proměnné.
$\alpha$Integrální součinitel teplotní roztažnosti $[1/\text{K}].$
$\alpha,\beta$Úhly natočení hrany prvku semi-loof $[\text{rad}].$
${\dot\varepsilon}_c$Rychlost efektivní creepové deformace $[1/\text{h}].$
$\varepsilon_{z0}$Nenulová rovinná deformace $[1].$
$\nu$Poissonovo číslo $[1].$
$\rho$Hustota $[\text{kg}/\text{m}^3].$
$\sigma_Y$Mez kluzu $[\text{Pa}].$
$\varphi_x,\varphi_y,\varphi_z$Úhly natočení uzlu nosníkového prvku kolem globálních os $[\text{rad}].$
$\Phi$Dilatační faktor $[1].$
cs/ref/name/i2.txt · Poslední úprava: 2023-09-07 08:53 autor: Petr Pařík