, kde 
F=F-Fcr.) Přímý tvar prutu se stává nestabilním stavem. Tento typ vidličkové bifurkace má přívlastek superkritická.1Ing. Petr Frantík, Ph.D., Vysoké učení technické v Brně, Fakulta stavební, Ústav stavební mechaniky, Veveří 331/95, 602 00 Brno, e-mail: frantik.p at fce.vutbr dot cz
2doc. RNDr. Jiří Macur, CSc., Vysoké učení technické v Brně, Fakulta stavební, Ústav automatizace inženýrských úloh a informatiky, Veveří 331/95, 602 00 Brno, e-mail: macur.j at fce.vutbr dot cz
Úvod
V minulém století nově vzniklé obory, opírající se o studium silně nelineárních systémů, jako jsou teorie katastrof, teorie bifurkací a teorie chaosu, přinesly nejen do mechaniky nové způsoby náhledu na numerické modelování. Zásluhy na jejich vzniku má zejména Henri Poincaré, jenž již koncem 19. století odhaluje složitý pohyb tří gravitačně vázaných těles. Ukázalo se, že díky univerzálnosti projevů nelineárních systémů, jejím tzv. generickým vlastnostem, lze výstižnost modelu reálných soustav a procesů posuzovat odlišnými kritérii, než bylo doposud zvykem. Důležitým atributem modelu přestává být vysoká úroveň detailu či značný počet stupňů volnosti. Naopak, v návaznosti na obtížnost a časovou náročnost řešení silně nelineárních systémů je zvýhodněn model s minimálním počtem stupňů volnosti. Podstatným se stává kvalitativní a nikoliv kvantitativní hledisko vystižení reálného systému či procesu. Jinak řečeno model má vykazovat stejné nelineární projevy jako reálný systém.
Typickými nelineárními projevy - generickými vlastnostmi - jsou výskyty katastrofických přechodů, bifurkací, chaotického chování. Souhrnně jsou tyto projevy doprovázeny ztrátou stability jak reálného systému, tak i jeho numerického modelu. Vznikají zvláštní případy, kdy není jednoduché rozhodnout o správnosti získaných výsledků. Bylo a je třeba vytvořit netradiční přístupy pro porozumění takovým jevům.
Jednou z částí mechaniky, kde se setkáváme s potřebou užití nových metod, jsou pružné štíhlé konstrukce, které vlivem zatížení výrazně mění svou geometrii. Dochází k velkým deformacím částí konstrukce, aniž by materiál z něhož je složena, ztrácel schopnost vracet se do původního stavu. Ačkoliv aplikace v tomto odvětví vyžadují nelineární matematický aparát, vykazují výbornou shodu s experimenty. Proto hrají důležitou roli v celé problematice nelineárních systémů.
V mechanice konstrukcí je velmi známým nelineárním projevem tzv. vidličková bifurkace. Dochází k ní mimo jiné při zatížení dokonale symetrického prutu dvěma stejně velkými silami působícími na jeho koncích proti sobě, viz obr. 1.
Obr. 1: Vzpěr přímého prutu
Dosáhne-li velikost síly F kritické hodnoty označované Fcr, přestává být přímý tvar prutu stabilní a dochází k prudkému příčnému vybočení prutu na jednu ze dvou stran (v rovinné úloze). Vhodné znázornění této závislosti vypadá jako vidlička, viz obr. 2. Odtud její název vidličková bifurkace.
Obr. 2: Superkritická vidličková bifurkace
Na obr. 2 je zobrazen schematický graf příčného vybočení středního průřezu prutu (vzepětí), označeného symbolem w, v závislosti na velikosti vzpěrné síly F. Vidíme, že dosažením kritické síly dochází k nekonečně rychlému vzdalování dvou možných stabilních stavů. (Vzdálenost těchto stavů je úměrná , kde F=F-Fcr.) Přímý tvar prutu se stává nestabilním stavem. Tento typ vidličkové bifurkace má přívlastek superkritická.
Úloha vzpěru prutu a různá její řešení jsou dnes velmi dobře známy. Již v roce 1744 publikoval Leonhard Euler první aproximaci ztráty stability prutu při dosažení kritické síly. Jeho řešení odvozené užitím tzv. teorie druhého řádu (V teorii druhého řádu uvažujeme vliv změny geometrie konstrukce na podmínky rovnováhy. Posunutí způsobující změnu geometrie prutu je předpokládáno pouze v příčném směru a je vztahováno k původnímu nedeformovanému tvaru.)
nazýváme klasické. Předpokládá štíhlý prut tvořený lineárně se chovajícím materiálem a pro hodnotu kritické síly Fcr dává přesný vztah:
 | (1) |
kde L je délka prutu a EI jeho ohybová tuhost. Klasické řešení nám dává dobrou představu pouze pro velikost síly F menší než hodnota kritická Fcr. Díky němu víme, že je před dosažením kritické síly přímý tvar prutu stabilní. Obraz o ztrátě stability prutu je ovšem mlhavý. Řešení pouze ukazuje možné deformační tvary, kterými prut může náhle vybočit.
Přesné řešení pokritického vzpěru štíhlého prutu odvozuje Alexej Rufovič Ržanicyn v polovině dvacátého století užitím eliptických integrálů.
V knize [Ržanicyn 1955] publikuje explicitní vyjádření vzepětí prutu v závislosti na velikosti vzpěrné síly.
Imperfektované systémy
Výzkum doprovázející vznik teorie katastrof ukázal, že studované systémy lze rozdělit do dvou skupin na tzv. typické (též regulární), které nejsou citlivé na malou změnu svých parametrů a na systémy netypické (též degenerované), jejichž malá změna způsobí, že se z nich stanou systémy typické, viz [Arnold 1983]. Netypické systémy, jako je např. vzpěr dokonale symetrického prutu, se reálně nevyskytují. Jejich přechod na systém typický je doprovázen ztrátou symetrie, jako je tomu právě při vzpěru prutu. Již Henri Poincaré poukázal, že při studiu systémů s proměnlivými parametry se nelze vyhnout degenerovaným systémům, jelikož tyto rozdělují oblasti parametrů odpovídajících kvalitativně odlišným systémům typickým. Degenerované systémy tvoří neodstranitelné množiny a je třeba je brát do úvahy právě tehdy, když se studuje systém s proměnlivými parametry.
Prvek, jenž vnáší do systému nesymetrii, se nazývá imperfekce (odchylka od dokonalosti). Přidáním imperfekce do systému s vidličkovou bifurkací se výrazně mění chování systému. Při změně řídícího parametru (V případě vzpěru prutu je řídícím parametrem velikost síly F.)
se již nesetkáváme s výskytem bifurkace na obr. 2. Stav systému se plynule vyvíjí, zatížený prut vybočuje jen na jednu stranu, viz obr. 3. Na opačné straně se v určitý okamžik zrodí "ze vzduchu" dvojice větví, z nichž vnější je stabilní a vnitřní je nestabilní.
Obr. 3: Imperfektovaný systém se superkritickou vidličkovou bifurkací
Jelikož ve skutečnosti existují pouze systémy regulární, nabízí se otázka určování tak význačné hodnoty, jakou je kritická síla. Vhodným bodem pro měření bifurkační hodnoty Fcr je poloha inflexe závislosti w(F). Označme velikost síly v inflexním bodě Fin (viz obr. 3). V publikaci [Macur & Frantík 2004] je dokázáno, že poloha inflexe Fin je v případě kubického modelu totožná s hodnotou bifurkační Fcr a nezávislá na velikosti imperfekce. U složitějších systémů s vidličkovou bifurkací lze předpokládat, že závislost hodnoty Fin na velikosti imperfekce bude buď konstantní, nebo alespoň monotónní. Jak bude ukázáno dále, je-li imperfekce "malá", je poloha inflexního bodu od bifurkační hodnoty jen velmi málo odlišná.
Poznamenejme, že v inflexním bodě závislosti w(F) se stav systému nejrychleji mění (první derivace má lokální extrém), což dobře odpovídá jedné z vlastností kritické síly.
Druhým význačným bodem, který se limitně ztotožní s bifurkační hodnotou řídícího parametru, je místo zrodu dvojice větví. Označme hodnotu síly pro tento bod For, viz obr.3. Nachází-li se systém ve stavu odpovídajícímu stabilní větvi takto zrozené, pak při poklesu řídícího parametru dojde při dosažení místa zrodu ke katastrofě - tzv. tvrdé ztrátě stability. Systém přejde na stabilní větev "na straně imperfekce" a bez vnějšího zásahu se již do stavu na opačné straně nevrátí.
Jak bude ukázáno dále, hodnota For je na imperfekci velmi citlivá. Proto ji lze využít jako vhodné měřítko velikosti imperfekce.
Numerický model prutu
Pro ověření popsaných teoretických závěrů byl vytvořen diskrétní model sloužící k výpočtům velkých deformací štíhlých prutů, poprvé popsaný v publikaci
[Frantík 2004]. Prut je v modelu reprezentován zvoleným počtem tuhých dílců s vnitřními liniovými pružinami. Dílce jsou vzájemně spojeny klouby s rotačními pružinami, viz obr. 4. Užitou metodu lze nazvat fyzikální diskretizací. (Fyzikální diskretizací se myslí vytvoření takového modelu konstrukce, který je svou podstatou diskrétní, ovšem lze jej chápat jako možnou reálnou strukturu. Protějškem je diskretizace matematická, kdy se spojitý problém diskretizuje parametrizovanými funkcemi (tj. superpozicí funkcí s konečným počtem parametrů).)
Zařazení modelu do metod mechaniky přehledně uvádí studie [Henrych 1985].
Obr. 4: Imperfektovaný systém se superkritickou vidličkovou bifurkací
Oba typy pružin uvažujme jako lineární. Potenciální energii Ep, která se akumuluje v těchto pružinách lze zapsat ve tvaru:
 | (2) |
kde nl je počet liniových pružin, nf je počet rotačních pružin, kl je tuhost liniových pružin, kf je tuhost rotačních pružin, dli je protažení i-tého dílce a d
i je natočení rotační pružiny (vzájemné pootočení dílců, které pružina spojuje). Předpokládejme, že deformační stav modelu prutu jednoznačně určují polohy všech jeho kloubů. Nechť je poloha každého kloubu dána dvojicí souřadnic (xi,yi), kde i je index kloubu. Pro protažení liniové pružiny dli lze pak psát:
 | (3) |
kde l je původní délka dílce (dílec bez napětí) a li je délka dílce po deformaci. Pro pootočení rotační pružiny di lze psát:
| (3) |
kde i je pootočení i-tého dílce. Ze vztahů (3) a (4) je patrné, že je užita geometricky přesná formulace modelu. Lze dokázat, že odezva modelu konverguje k přesnému řešení [Ržanicyn 1955], pokud zvyšujeme počet dílců a tuhost liniových pružin.
Model lze řešit například Newtonovou iterační metodou s užitím numerické derivace potenciální energie Ep. Dodejme, že obohacení modelu o zatížení a imperfekci (Imperfekce je v použitém modelu reprezentována pomocí síly působící uprostřed prutu ve směru osy y, viz obr. 4.)
lze snadno provést rozšířením vztahu (2) o jejich energetickou bilanci.
Na obr. 5 jsou vidět vypočtené závislosti vzepětí prutu w na vzpěrné síle F včetně vlivu imperfekce. Imperfekce je určena pomocí vzepětí prutu w0, který není zatížen vzpěrnou silou F. Jeho hodnota je v grafu uváděna bezrozměrně. Ržanicynovu přesnému řešení odpovídá funkce s nulovou velikostí imperfekce. Inflexní body funkcí, respektive extrémy první derivace funkcí, jsou zvýrazněny puntíky.
Obr. 5: Graf vzepětí prutu s imperfekcí v bezrozměrném tvaru. Velikost imperfekce prutu je dána poměrem 
(v promile), kde w0 je vzepětí nezatíženého prutu
Nakolik je poloha inflexního bodu Fin a poloha zrodu dvojice větví For závislá na velikosti imperfekce, ukazuje obr. 6. Je vidět, že pro malé imperfekce je hodnota Fin prakticky shodná s hodnotou kritické síly. Z grafu je patrné, že se poloha inflexního bodu Fin jen málo liší od bifurkační hodnoty Fcr i při imperfekci odpovídající prohnutí prutu více než jeden a půl centimetru na délku jeden metr. Jinak je tomu u hodnoty For popisující polohu zrodu. Funkce For(w0) zřejmě nemá pro nulovou imperfekci derivaci.
Obr. 6: Graf závislosti polohy inflexe a zrodu na velikosti imperfekce
Závěr
V příspěvku bylo pojednáno o nelineárních systémech s vidličkovou bifurkací a vlivu imperfekce na jejich chování. Pomocí diskrétního nelineárního modelu vzpíraného prutu bylo ukázáno, že pro měření kritické síly imperfektovaného systému lze užít polohu inflexního bodu závislosti vzepětí prutu na velikosti vzpěrné síly. Velikost imperfekce systému lze stanovit pomocí polohy zrodu stabilní a nestabilní větve.
Poděkování
Tento výsledek byl získán za finančního přispění MŠMT, projekt 1M6840770001, v rámci činnosti výzkumného centra CIDEAS.
Literatura
[Arnold 1983] Arnold, V. I., 1983: Teória katastrof (orig. Teorija katastrof, vydavatelstvo Moskevské univerzity 1983), vydavateľstvo Alfa, Bratislava
[Bažant & Cedolin 1991] Bažant Z. P., Cedolin L., 1991: Stability of Structures, Elastic, Inelastic, Fracture, and Damage Theories, Oxford University Press, New York
[Frantík 2004] Frantík, P., 2004: Stability Study of the Elastic Loop, 5th International PhD Symposium in Civil Engineering, Delft, Netherlands
[Henrych 1985] Henrych, J., 1985: Úplná soustava finitních metod mechaniky a možnosti dalšího rozvoje, studie ČSAV 6.85, nakladatelství Akademia, Praha
[Macur & Frantík 2004] Macur, J., Frantík, P., 2004: Stanovení kritické síly vzepjatého nosníku, konference Pravděpodobnost porušování konstrukcí, VUT FAST v Brně
[Ržanicyn 1955] Ržanicyn, A., R., 1955: Ustojčivost' ravnovesija uprugich sistem, Gosudarstvennoje izdatel'stvo techniko-teoretičeskoj literatury, Moskva