227x

31.8.2023

VE0048 | Jednoosý ohyb s tlakem

Popis

Konstrukce tvořená nosníky profilu I je na levém konci vetknutá a na pravém podepřená posuvnou kloubovou podporou. Konstrukce se skládá ze dvou segmentů podle následujícího {%/#sketch skici]]. Vlastní tíha se v tomto příkladu nezohledňuje. Stanovte maximální průhyb konstrukce u_z,max, ohybový moment M_y na pevném konci, natočení φ_2,y segmentu 2 a reakční sílu R_Bz pomocí geometricky lineární analýzy a teorie druhého řádu analýzy. Verifikační příklad je založen na příkladu, který představili Gensichen a Lumpe (viz odkaz).

Materiál	Ocel	Modul pružnosti	E	210000,000	MPa
Materiál	Ocel	Poissonův součinitel	ν	0,300	-
Geometrie	Konstrukce	Délka segmentu 1	L₁	6,000	m
	Konstrukce	Délka segmentu 2	L₂	1,200	m
	Průřez	Výška	h	400,000	mm
		Šířka	b	180,000	mm
		Tloušťka stěny	s	10,000	mm
		Tloušťka pásnice	t	14,000	mm
Zatížení		Normálová síla	F_x	100,000	KN
Zatížení		Příčná síla	F_z	0,500	KN

01 Jednoosý ohyb s tlakem

Analytické řešení

Teorie I. řádu (geometricky lineární výpočet)

Nejdříve se provede geometricky lineární analýza. V tomto případě se normálová síla F_x nezohledňuje. Úlohu lze pak vyřešit stejně jako konzolu délky L₁ zatíženou pouze příčnou silou F_z. Maximální průhyb u_z,max lze spočítat pomocí Mohrova' integrálu a výsledkem je:

u_{z, \max} = \frac{F_{z} L_{1}^{3}}{3 {EI}_{y}} = 0.743 mm

F<sub>y</sub>

Kvadratický moment průřezu vzhledem k ose y

Maximální průhyb - Geometricky lineární analýza

Ohybový moment na pevném konci lze vypočítat podle následujícího vzorce:

M_{y} (0) = F_{z} L_{1} = 3.000 kNm

Maximální ohybový moment - Geometricky lineární analýza

Natočení segmentu 2 φ_2,y se vypočítá z geometrické podmínky následovně:

φ_{2, y} = \arctan (\frac{u_{z, \max}}{L_{2}}) = 0.619 mrad

Natočení segmentu 2 - Geometricky lineární analýza

Reakční sílu v kluzné spáře R_Bz se zohledněním nulového účinku normálové síly F_x lze získat z diagramu volného tělesa v následujícím {%/#náčrtu volného tělesa]].

R_{B z} = - \frac{F_{x} u_{z, \max}}{L_{2}} = 0.000

Reakční síla v kluzné spáře - Geometricky lineární analýza

02 Návrh konstrukcí zdarma

Analýza druhého řádu

Vzhledem k nezanedbatelnému účinku normálové síly F_x je třeba zohlednit analýzu druhého řádu. Tím se zohlední normálová síla F_x a vytvoří další příspěvek k ohybovému momentu. Problém lze popsat pomocí diagramu volného tělesa segmentů podle {%/#freebody skici]]. Neznámé reakční síly lze získat z rovnic rovnováhy a následně zadat vzorec pro ohybový moment.

M_{y} = - R_{A z} (L_{1} - x) - R_{A x} (u_{z, \max} - u_{z} (x))

Ohybový moment v segmentu 1 - analýza druhého řádu

Řešením je Euler-Bernoulliho diferenciální rovnice.

\frac{d^{2} u_{z}}{{dx}^{2}} = - \frac{M_{y}}{{EI}_{y}}

Euler-Bernoulliho diferenciální rovnice

Při zohlednění okrajových podmínek lze najít řešení diferenciální rovnice a vypočítat maximální průhyb konstrukce.

u_{z} (x) = \frac{F_{z} L_{2} [- \cos ({αL}_{1}) αx + \cos ({αL}_{1}) \sin (αx) - \sin ({αL}_{1}) \cos (αx) + \sin ({αL}_{1})]}{F_{x} [{αL}_{1} \cos ({αL}_{1}) + {αL}_{2} \cos ({αL}_{1}) - \sin ({αL}_{1})]}

Funkce průhybu pro segment 1 - analýza druhého řádu

u_{z, \max} = \frac{F_{z} L_{2} [{αL}_{1} \cos ({αL}_{1}) - \sin ({αL}_{1})]}{F_{x} [αcos ({αL}_{1}) (L_{1} + L_{2}) - \sin ({αL}_{1})]} = 0.878 mm

Maximální průhyb konstrukce - analýza druhého řádu

Ohybový moment na pevném konci lze vypočítat podle následujícího vzorce:

M_{y} (0) = R_{A z} L_{1} + R_{A x} u_{z, \max} = 3.527 kNm

Maximální ohybový moment - analýza druhého řádu

Natočení segmentu 2 φ_2,y se vypočítá z geometrické podmínky následovně:

φ_{2, y} = \arctan (\frac{u_{z, \max}}{L_{2}}) = 0.732 mrad

Natočení segmentu 2 - analýza druhého řádu

Výsledná reakční síla v kluzné spáře R_Bz :

R_{B z} = - \frac{F_{x} u_{z, \max}}{L_{2}} = - 0.073 kN

Reakční síla v kluzné spáře - analýza druhého řádu

Je zřejmé, že vliv normálové síly F_x je značný. Celkový průhyb konstrukce při předepsaném zatížení je v případě teorie druhého řádu asi o 18 % větší než v případě geometricky lineární analýzy. Porovnání geometricky lineární analýzy a analýzy druhého řádu je znázorněno v {%/#grafu]] při zohlednění poměru zatěžovacích sil F_z = F_x/200. Je zřejmé, že rozdíl mezi těmito analýzami je výraznější, pokud je zatížení větší. Řešení analýzy druhého řádu se blíží vodorovné asymptotě. Z numerického řešení vyplývá hodnota vodorovné asymptoty F_x,cr = 650,873 kN.

03 Porovnání geometricky lineární analýzy (přerušovaná čára) a analýzy druhého řádu (plná linie)

Nastavení programů RFEM a RSTAB

Modelováno v programech RFEM 5.05 a RSTAB 8.05 a RFEM 6.01, RSTAB 9.01
Počet prvků je 2 (jeden prvek na prut)
Počet přírůstků je 5
Je použit izotropní lineárně elastický materiálový model
Konstrukce je modelována pomocí prutů
Smyková tuhost prutů se zanedbává

Výsledky

Geometricky lineární analýza	Analytické řešení	RFEM 6	Poměrná hodnota	RSTAB 9	Poměrná hodnota
u_z,max [mm]	0,743	0,743	1,000	0,743	1,000
M_y (0) [kNm]	3,000	3,000	1,000	3,000	1,000
φ_2,y [mrad]	0,619	0,619	1,000	0,619	1,000
R_Bz [kN]	0,000	0,000	-	0,000	-

Geometricky lineární analýza	Analytické řešení	RFEM 5	Poměrná hodnota	RSTAB 8	Poměrná hodnota
u_z,max [mm]	0,743	0,743	1,000	0,743	1,000
M_y (0) [kNm]	3,000	3,000	1,000	3,000	1,000
φ_2,y [mrad]	0,619	0,619	1,000	0,619	1,000
R_Bz [kN]	0,000	0,000	-	0,000	-

Analýza druhého řádu	Analytické řešení	RFEM 6	Poměrná hodnota	RSTAB 9	Poměrná hodnota
u_z,max [mm]	0,878	0,878	1,000	0,878	1,000
M_y (0) [kNm]	3,527	3,527	1,000	3,527	1,000
φ_2,y [mrad]	0,732	0,732	1,000	0,732	1,000
R_Bz [kN]	-0,073	-0,073	1,000	-0,073	1,000

Analýza druhého řádu	Analytické řešení	RFEM 5	Poměrná hodnota	RSTAB 8	Poměrná hodnota
u_z,max [mm]	0,878	0,878	1,000	0,878	1,000
M_y (0) [kNm]	3,527	3,527	1,000	3,527	1,000
φ_2,y [mrad]	0,732	0,732	1,000	0,732	1,000
R_Bz [kN]	-0,073	-0,073	1,000	-0,073	1,000

500x

16x

Verifikační příklad 000048 | 1

Reference

LUMPE, G. a GENSITEN, V. Vyhodnocení lineární a nelineární analýzy prutů v teorii a v softwaru: Testovací příklady, příčiny selhání, detailní teorie. Ernest.

;