Tymošenkova teorie ohybu nosníku

Z testwiki
Skočit na navigaci Skočit na vyhledávání

Tymošenkova teorie ohybu nosníku (přepisem z ruštiny Timošenkova teorie, též Bresseho-Tymošenkova teorie[1][2]) je teorie umožňující přibližně popsat ohyb nosníku. Teorie vychází ze dvou předpokladů o deformaci nosníku:

  • Příčné řezy nosníku, které byly před deformací nosníku kolmé na osu nosníku, zůstávají při deformování nosníku rovinné, ale ne nutně kolmé na osu nosníku.
  • Rozměry příčných řezů nosníku se při deformaci nosníku nemění.

Na základě těchto předpokladů lze deformaci nosníku zcela popsat pomocí relativního posunutí a natočení průřezů nosníku. Od klasické Bernoulliho–Navierovy hypotézy ohybu nosníku[3][4][5] se Tymošenkova teorie liší odstraněním předpokladu o zachování kolmosti průřezů nosníku na osu nosníku. Díky tomu poskytuje Tymošenkova teorie přesnější popis deformace a namáhání nosníku, což se projevuje především v případech, kdy příčné rozměry nosníku nejsou zanedbatelně malé vůči rozměru podélnému nebo když frekvence příčného kmitání nosníku je natolik vysoká, že příčné rozměry nosníku již nejsou zanedbatelně malé vůči délce příčné vlny. Nevýhodou Tymošenkovy teorie, oproti teorii Eulerově-Bernoulliho, je vyšší počet deformačních charakteristik a tím i komplikovanější charakter řešení.

Přesnější popis ohybu nosníku, než poskytuje Tymošenkova teorie, nabízí teorie ohybu nosníku s méně omezujícími předpoklady o deformaci nosníku, např. Levinsonova teorie[6] či teorie Reddyho[7], popřípadě postupy založené na popisu deformace dvoj- či trojrozměrného kontinua[8][9]. Ekvivalentem Tymošenkovy teorie pro popis ohybu desek je teorie Reissnerova-Mindlinova.[10]

Teorie je pojmenována podle Stepana Prokopovyče Tymošenka, ukrajinského inženýra a profesora inženýrské mechaniky na Stanfordově univerzitě,[11] který je autorem knihy Kurs teorii uprugosti: Steržni i plastinki z roku 1916, kde je tato teorie popsána[12].

Rovinný ohyb

V případě, že se nosník prohýbá v jedné rovině a nedochází ke kroucení nosníku, což platí například u izotropního nosníku se symetrickým průřezem, který je vystavený zatížení působícímu v rovině symetrie, lze dle předpokladů teorie popsat posuv bodů nosníku pomocí vztahů

u(x,z,t)=u0(x,t)+zφ(x,t)
v(x,z,t)=0
w(x,z,t)=w(x,t),

kde u je posuv ve směru osy nosníku, w je průhyb nosníku a v je posuv ve směru kolmém na rovinu průhybu a osu nosníku, u0 je posuv ve směru osy nosníku na ose nosníku x, z je kolmá vzdálenost od osy nosníku měřená v rovině průhybu a φ je natočení průřezu nosníku kolem osy y.

Řídicí diferenciální rovnice

Diferenciální rovnice popisující chování nosníku a správný tvar okrajových podmínek lze odvodit pomocí d'Alembertova principu[13] pro deformovatelné těleso[5]. Podle něj pro těleso bez jednostranných vazeb v kterémkoli okamžiku platí, že součet práce vnitřních sil, vnějších sil a práce setrvačných sil vykonané na jakémkoli geometricky přípustném poli posunutí, jež nezávisí na čase a vztahuje se ke stavu tělesa v daném okamžiku, se rovná nule:

δWvnitrˇnı´+δWvneˇjsˇı´+δWsetrvacˇne´=0,

kde W označuje práci a δ označuje variaci, tj. hodnotu související s jakoukoli geometricky přípustnou deformací nosníku, popsanou pomocí libovolných geometricky přípustných posuvů δu0 a δw a natočení δφ.


V případě nosníku se spojitým příčným zatížením platí (zatížení na konci nosníku nemění tvar výsledných diferenciálních rovnic) Šablona:Vzorec kde σxx je normálové napětí působící v průřezu nosníku, τxz jsou smykové napětí působící v průřezu nosníku, εxx je poměrné délkové přetvoření, γxz je poměrné zkosení, q je spojité příčné zatížení nosníku, ρ je hustota materiálu, A je plocha průřezu a L je délka nosníku.

Pro poměrné délkové přetvoření (normálové inženýrské přetvoření) a poměrné zkosení (smykové inženýrské přetvoření) platí:Šablona:Poznámka

εxx=u(x,z,t)x=u0(x,t)x+zφ(x,t)x

Šablona:Vzorec Velikost smykového přetvoření podle Tymošenkovy teorie tedy nezávisí na poloze po průřezu nosníku.

Z rovnice (Šablona:Odkaz na vzorec) lze po vyjádření přetvoření pomocí posuvů a integrace per partes odvodit Šablona:Vzorec

Znaménková konvence – část nosníku vymezená dvěma nekonečně blízkými průřezyŠablona:Poznámka a) výchozí stav, b) deformovaný stav.

kde

N=AσxxdA

je osová síla,

M=AzσxxdA

je ohybový moment kolem osy y,

Q=AτxzdA

je posouvající (smyková) síla,

η=AρdA

je hmotnost nosníku připadající na jednotku délky,

J=Aρz2dA

je moment setrvačnosti průřezu kolem osy y, a platí, že osa nosníku prochází těžišti průřezů nosníku, tj. že

L(2u0t2δφAzρdA)dx=L(2φt2δu0AzρdA)dx=0.

Z rovnice (Šablona:Odkaz na vzorec) , která musí platit pro libovolné geometricky přípustné posuvy δu0 a δw a natočení δφ, pak vyplývá soustava diferenciálních rovnic popisujících kmitání nosníku Šablona:Vzorec

Šablona:Vzorec

Šablona:Vzorec

a okrajové podmínky, které musí být specifikovány na obou koncích nosníku: N, nebo u0, M, nebo φ a Q, nebo w.

Rovnice (Šablona:Odkaz na vzorec) a (Šablona:Odkaz na vzorec) lze převést na jedinou parciální diferenciální rovnici Šablona:Vzorec

Homogenní izotropní nosník

V případě homogenního nosníku z lineárně elastického izotropního materiálu lze s pomocí Hookova zákona vyjádřit osovou sílu, ohybový moment a posouvající sílu následujícím způsobem

N=AσxxdA=AEεxxdA=AEu0xdA+AzEφxdA=EAu0x
M=AzσxxdA=AzEεxxdA=AzEu0xdA+Az2EφxdA=EIφx
Q=AτxzdA=κAGγxzdA=κAG(φ+wx)dA=κAG(φ+wx),

kde E je Youngův modul (modul pružnosti v tahu), G je modul pružnosti ve smyku,

I=Az2dA

je kvadratický modul průřezu a κ je opravný smykový součinitel upravující vztah mezi smykovým přetvořením γxz, které má dle předpokladů Tymošenkovy teorie po průřezu nosníku neměnnou hodnotu (viz (Šablona:Odkaz na vzorec)), a posouvající silou Q. Korekční součinitel byl zaveden, protože skutečné rozložení smykového napětí po průřezu nosníku je nerovnoměrné, a tudíž velikost posouvající síly určená pomocí Hookova zákona z velikosti smykového přetvoření, u kterého je předpokládáno, že je po průřezu konstantní, by neodpovídala skutečnosti.

Kmitání nosníku

Řídicí diferenciální rovnice (Šablona:Odkaz na vzorec)-(Šablona:Odkaz na vzorec) pak mohou být vyjádřeny ve tvaru Šablona:Vzorec Šablona:Vzorec Šablona:Vzorec

Je zřejmé, že rovnici (Šablona:Odkaz na vzorec) lze řešit nezávisle na rovnicích (Šablona:Odkaz na vzorec) a (Šablona:Odkaz na vzorec).

U nosníku s konstantním průřezem lze rovnice (Šablona:Odkaz na vzorec) a (Šablona:Odkaz na vzorec), za předpokladu dostatečné hladkosti všech jejich členů, převést na jedinou rovnici

EI4wx4(J+EIηκAG)4wx2t2+JηκAG4wt4+η2wt2=q+JκAG2qt2EIκAG2qx2,

respektive Šablona:Vzorec kde

cL=Eρ

je rychlost šíření podélné vlny[14].

Rovnici (Šablona:Odkaz na vzorec) lze rovněž vyjádřit ve tvaru Šablona:Vzorec kde

cS=Gρ

je rychlost šíření smykové vlny[14].

Za stejných předpokladů jako při odvození rovnice (Šablona:Odkaz na vzorec), lze rovnice (Šablona:Odkaz na vzorec) a (Šablona:Odkaz na vzorec) převést na rovnici

4φx41cL2(1+EκG)4φx2t2+ρcL2κG4φt4+AcL2I2φt2=1EIqx,

která je pro q=0 formálně totožná s (Šablona:Odkaz na vzorec).

Disperze vlnění

Postupující harmonické vlnění lze popsat vztahem Šablona:Vzorec kde D je amplituda kmitání, c je rychlost šíření vlnění (fázová rychlost), ω je úhlová frekvence a k je vlnové číslo, pro něž platí

k=2πλ=ωc,

kde λ označuje vlnovou délku.

Dosadí-li se vztah (Šablona:Odkaz na vzorec) do rovnice (Šablona:Odkaz na vzorec) a předpokládá-li se, že q=0 a D0, pak musí platit

(k2+ω2cL2)(k2+ω2κcS2)Aω2IcL2=(k2+k2c2cL2)(k2+k2c2κcS2)Ak2c2IcL2=0.

Odtud lze odvodit vztah pro rychlost šíření vlnění ve tvaru

c1,2=(κcS2+cL2+κcS2Ak2I)±(κcS2+cL2+κcS2Ak2I)24κcS2cL22.

Rychlost šíření vlnění je tedy závislá na vlnové délce – dochází k disperzi vlnění[14]. Pro vlnové délku blížící se v limitě nule je rychlost šíření vlnění konečná, čímž se Timošenkova teorie liší od teorie Eulerovy-Bernoulliho. To platí i pro grupovou rychlost[14]. Z výše uvedeného vztahu je také patrné, že existují dva typy vlnění, které se šiří různou rychlostí[14].

Statický ohyb

V případě statického ohybu lze rovnice (Šablona:Odkaz na vzorec) a (Šablona:Odkaz na vzorec) převést na rovnice vhodnější pro řešení Šablona:Vzorec Šablona:Vzorec

U nosníku s konstantním průřezem lze pak psát

4wx4=qEI1κAG2qx2.
Příklad – vetknutý nosník s konstantním spojitým zatížením
Vetknutý nosník s konstantním spojitým zatížením.

Pro vetknutý nosník (konzolu) s konstantním průřezem, který je vystavený spojitému konstantnímu svislému zatížení, lze integrací rovnic (Šablona:Odkaz na vzorec) a (Šablona:Odkaz na vzorec) při užití okrajových podmínek

w(L)=φ(L)=M(0)=Q(0)=0

získat řešení

φ=q6EI(L3x3),
w=q2κAG(L2x2)+qx24EI(x34L3)+qL48EI.

Maximální průhyb lze vyjádřit ve tvaru

wmax=qL22κAG+qL48EI=qL48EI(1+4EIκAGL2)=wmax,EB(1+4EIκAGL2),

kde člen v závorce představuje součinitel zpřesňující hodnotu průhybu získanou pomocí klasické Eulerovy-Bernoulliho teorie.

Pro nosník s obdélníkovým průřezem pak lze psát

wmax=3qL42Ebh3(1+13κEGh2L2),

kde b je šířka a h výška průřezu. Obdobný vztah platí i pro nosník z ortotropního materiálu, jehož hlavní materiálová osa je rovnoběžná s osou nosníku[15]. Z uvedeného vztahu je zřejmé, že vliv smykové deformace se výrazněji projeví u krátkých nosníků či nosníků z materiálů, jejichž modul pružnosti ve smyku je výrazně menší než modul pružnosti v tahu.

Opravný smykový součinitel

Opravný smykový součinitel formálně udává poměr mezi rovnoměrnou poměrnou příčnou smykovou deformací nosníku, dle předpokladů Tymošenkovy teorie, a velikostí posouvající síly

κ=Q(φ+wx)AGdA=Qγxz,TymosˇenkoAGdA,

přičemž pro κ= Tymošenkova teorie odpovídá Rayleighově teorii kmitání nosníku[16]. Pokud se nejedná o rovinný ohyb nosníku, je zapotřebí znát alespoň dva obecně různé korekční součinitele[17].

Někdy se pojem opravný smykový součinitel užívá pro jeho převrácenou hodnotu, tj.

κ*=1κ,

viz např. [14].

Opravný smykový součinitel lze chápat jako nástroj pro úpravu předpovědí chování nosníků dle Tymošenkovy teorie tak, aby co nejlépe odpovídaly výsledkům experimentů, či předpovědím dle přesnějších teorií[18][19]. Velikost součinitele obecně závisí na tvaru[12][20] a velikosti[21] průřezu, materiálových vlastnostech[22], zatížení a okrajových podmínkách nosníku a frekvenci kmitání nosníku[23].

Výpočtové postupy jsou založeny na různých principech, např.:

  • Velikost součinitele je určena z Žuravského vztahu pro smykové napětí na ose nosníku, tedy vztahu mezi posouvající silou a smykovým napětím, a Hookova zákona, pomocí něhož lze smykové napětí převést na smykové přetvoření[12][20][24].
  • Velikost součinitele je určena z podmínky ekvivalence deformační energie určené pomocí součinu posouvající síly a smykového přetvoření a deformační energie určené integrací po průřezu nosníku součinu smykového napětí podle Žuravského vztahu a jemu odpovídajícího smykového přetvoření[25][26][27].
  • Velikost součinitele vyplývá ze vztahu mezi vlastní frekvencí ohybového kmitání nosníku (dle Tymošenkovy teorie) a vlastní frekvence kmitání nosníku dle teorie pro popis kmitání dvojrozměrného kontinua (exaktní řešení – podmínka rovinného napětí, či přetvoření) či trojrozměrného kontinua (kruhový průřez nosníku – přibližné řešení)[28].
  • Velikost součinitele vyplývá ze srovnání hodnoty nejnižší vlastní frekvence kmitání nosníku při které je průhyb nosníku w nulový (tzv. smykový mód kmitání) s vlastní frekvencí kmitání dvojrozměrného kontinua, u něhož je povolen pouze posuv ve směry osy nosníku a platí, že povrch nosníku je nezatížený[29].
  • Velikost součinitele vyplývá ze vztahu mezi průměrným smykovým přetvořením po průřezu nosníku a posouvající silou, přičemž hodnota přetvoření je určena z exaktního analytického řešení pro napětí po průřezu vetknutého nosníku zatíženého osamělou silou či konstantním spojitým zatížením. Pro daný tvar průřezu je potřeba znát harmonickou funkci popisující borcení průřezu při kroucení nosníku[18].
  • Velikost součinitele je určena z podmínky ekvivalence deformační energie určené pomocí součinu posouvající síly a smykového přetvoření a deformační energie odpovídající smykovému napětí, které se určí ze známých napětí dle předpokladů Eulerovy-Bernoulliho teorie ohybu nosníku a podmínek rovnováhy nekonečně malého elementu nosníku. Z toho plyne potřeba najít řešení Poissonovy rovnice stejného typu jako u Saint-Venantovy teorie krutu[17][22].
Opravné smykové součinitele pro homogenní izotropní nosník dle Cowpera[18] – výběr.
Průřez Součinitel
Obdélník 10(1+ν)12+11ν
Kruh 6(1+ν)7+6ν
Elipsa 12(1+ν)a2(3a2+b2)(40+37ν)a4+(16+10ν)a2b2+νb4
Půlkruh 1+ν1.305+1.273ν
Plnostěnná trubka 6(1+ν)(1+m2)2(7+6ν)(1+m2)2+(20+12ν)m2, kde m=ba
Tenkostěnná trubka 2(1+ν)4+3ν
Tenkostěnný uzavřený profil 20(1+ν)48+39ν
Uzavřený profil 10(1+ν)(1+3m)2(12+72m+150m2+90m3)+ν(11+66m+135m2+90m3)+10n2((3+ν)m+3m2), kde m=bt1ht2 a n=bh
Profil I 10(1+ν)(1+3m)2(12+72m+150m2+90m3)+ν(11+66m+135m2+90m3)+30n2(m+m2)+5νn2(8m+9m2), kde m=2bt1ht2 a n=bh
Profil T 10(1+ν)(1+4m)2(12+96m+276m2+192m3)+ν(11+88m+248m2+216m3)+30n2(m+m2)+10νn2(4m+5m2+m3), kde m=bt1ht2 a n=bh
Korekční součinitele pro homogenní izotropní nosník s obdélníkovým a kruhovým průřezem dle různých autorů.
Reference Beamsection rectangle Kruh
Tymošenko (1916)[12] 230.667 34Šablona:Poznámka
Tymošenko (1922)[28]Šablona:Poznámka, Kaneko (1975)[19] 10(1+ν)12+11νŠablona:Poznámka 6(1+ν)27+12ν+4ν2
Goens(1931)[26] podle Föppla[25] 560.833 910
Mindlin (1953)[29] π2120.822 0.847
Renton (1991)[30] 165+(ν1+ν)2m=0n=1144(ba)4π6(2m+1)2n2[(2m+1)2(b2a)2+n2] 6(1+ν)27+14ν+8ν2
Pai (1999)[31] (1+ν)2a436((1+ν)2a430+ν2b4180ν2b52π5an=11n5tanh(nπab)) 6(1+ν)27+14ν+8ν2
Hutchinson (2001)[21] κ=2(1+ν)94a5bC+ν(1b2a2), kde C=445a3b(12a215νa2+5νb2)+n=116ν2b5(nπabtanh(nπab))(nπ)5(1+ν) 6(1+ν)27+12ν+4ν2

Historie

Teorii kmitání nosníku, jež zohledňuje smykovou deformaci a moment setrvačnosti průřezů nosníku, publikoval v roce 1859 Jacques Antoine Charles Bresse[32], profesor na École nationale des ponts et chaussées. Bresse odvodil soustavu diferenciálních rovnice popisující podélné a příčné kmitání slabě zakřiveného nosníku, ale tuto teorii ve své knize nepoužil pro řešení konkrétních problémů. Bresse nezavedl pojem korekční smykový součinitel, ale pracoval s příčnou tuhostí vyjádřenou součinem Youngova modulu pružnosti, E, a blíže nespecifikovaného součinitele k.

Tymošenko publikoval svoji teorii v ruštině v roce 1916[12] a v angličtině v roce 1920[20], nicméně obvykle je citována anglická verze z roku 1921[33]. Tymošenko se zabýval příčným kmitáním prizmatických nosníků, přičemž doplnil existující Rayleighovu teorii kmitání nosníku z roku 1877[16] zahrnující vliv rotační setrvačnosti průřezu nosníku o korekci na smykovou deformaci nosníku. Na příkladu vlastní frekvence kmitání nosníku s obdélníkovým průřezem a s podepřenými konci ukázal, že vliv smykové deformace na velikost vlastní frekvence je výraznější než vliv rotační setrvačnosti, přičemž oba vlivy klesají se vzrůstající vlnovou délkou a vzrůstajícím poměrem délky a výšky nosníku. Opravný smykový součinitel, označený k[12], κ[20] či λ[33], Tymošenko stanovil na základě předpokladu, že smyková deformace po průřezu nosníku odpovídá smykové deformaci na ose nosníku stanovené pomocí přibližné Žuravského teorie. Tento postup je však sporný, neboť závisí na volbě vztažné osy nosníku[25] a neposkytuje dostatečně přesné hodnoty smykové součinitele[19].

Ve svém článku Tymošenko Bresseho teorii nezmiňuje. Zmínku o Bresseho knize, avšak ne o Bresseho teorii, lze u Tymošenka najít v jeho pojednání o historii nauky o pružnosti a pevnosti[3] z roku 1953.

Ve svém článku z roku 1922 se Tymošenko zabýval kmitáním nosníku s obdélníkovým průřezem, jehož tloušťka je velmi malá (stav rovinné napjatosti) nebo velmi velká (stav rovinné deformace)[28], přičemž k nosníku přistupoval jako k dvojrozměrnému kontinuu. Odvodil vztah pro vlastní frekvenci kmitání a srovnal jej se vztahem dle teorie ze svého předcházejícího článku. Poté ještě srovnal vztah pro vlastní frekvenci kmitání válcového nosníku dle Pochhammera[34] se vztahem pro válcový nosník dle teorie ze svého dřívějšího článku. Pochhammerův vztah je aproximací exaktního řešení problému trojrozměrného kmitání nekonečně dlouhého nosníku s kruhovým průřezem.

K popularizaci své teorie Tymošenko později přispěl svými knihami[12][24][35].

V anglicky psané literatuře upozornil na existenci Bresseho teorie nejpozději v roce 1893 Karl Pearson v knize A history of the theory of elasticity and of the strength of materials from Galilei to the present time[36]. V souvislosti s Tymošenkovou teorií pak na Bresseho teorii upozornili Mindlin a Deresiewicz v roce 1953[29].

Korekce Eulerovy-Bernoulliho teorie zohledňující vliv smykové deformace na statický ohyb nosníku, která je identická s korekcí dle Tymošenka, byla publikována už v roce 1858 profesorem stavebního inženýrství a mechaniky na University of Glasgow Williamem Rankinem[37].

Odkazy

Poznámky

Šablona:Poznámky

Reference

Externí odkazy

Šablona:Autoritní data

Šablona:Portály