Vzper a strata stability: teória za konvergenčným zlyhaním

V inom článku som opísal, ako mi nelineárny výpočet tenkej škrupiny opakovane umieral a ako som zistil, že za tým stojí strata stability, nie chyba modelu. Tento článok je krok naspäť — k teórii, ktorá tomu dáva zmysel. Bez nej je vzper len číslo zo solvera; s ňou viete, čo to číslo znamená a kedy mu môžete veriť.

Kedy je rovnováha stabilná

Predstavte si guľôčku. V jamke je stabilná — po malom postrčení sa vráti. Na vrchole kopca je nestabilná — najmenšie postrčenie ju pošle preč. Na rovnej ploche je indiferentná — ostane, kam ju dáte. Konštrukcia sa správa rovnako, len namiesto polohy guľôčky sledujeme jej deformovaný tvar a namiesto gravitácie pôsobí zaťaženie.

Strata stability je bod, v ktorom pôvodne stabilná rovnovážna poloha prestane byť stabilná. Kritické zaťaženie je hodnota, pri ktorej sa to stane. Nie je to bod, kde sa materiál poruší — je to bod, kde konštrukcia prejde do iného tvaru, hoci napätie môže byť ešte hlboko pod medzou pevnosti. Práve preto je vzper zákerný: klasická kontrola napätia ho nezachytí. [1]

Stabilita rovnováhy: guľôčka v jamke sa po vychýlení vráti, na kopci ujde, na rovine ostane. Kritické zaťaženie je bod, keď sa jamka mení na kopec.

Euler: prečo o vzpere rozhoduje tuhosť, nie pevnosť

Klasický prípad je tlačený prút. Leonhard Euler ukázal, že sa vybočí pri kritickej sile

$$P_{\mathrm{kr}} = \frac{\pi^2 E I}{L^2},$$

kde $E$ je modul pružnosti, $I$ kvadratický moment prierezu a $L$ dĺžka. Na tomto vzorci je kľúčové, čo v ňom nie je: medza pevnosti materiálu. Kritické zaťaženie vzperu závisí od tuhosti ($E$), geometrie ($I$, $L$) a spôsobu uloženia — nie od toho, aký pevný materiál je. [1]

To má priamy dôsledok, ktorý som videl aj v praxi: či je akrylát krehký alebo ťažný, kritické zaťaženie straty stability sa nezmení. Materiál rozhoduje o tom, čo sa stane po strate stability (praskne, alebo sa ohne a prežije), a o napätiach — ale samotný bod nestability ťahá tuhosť a tvar. Preto sa dá tenká konštrukcia navrhnúť tak, že vyhovie na pevnosť a napriek tomu sa vybočí pri zlomku návrhového zaťaženia.

Pre aké uloženie ten vzorec platí — a čo je K

Tvar vyššie platí pre obojstranne kĺbovo uložený prút (kĺb–kĺb): oba konce sa môžu voľne otáčať, ale nemôžu sa bočne posunúť. To je referenčný, „základný Eulerov" prípad a vybočí sa do jednej polvlny sínusu. [1]

Ostatné uloženia sa pokryjú súčiniteľom vzpernej dĺžky $K$ (nazýva sa aj súčiniteľ efektívnej dĺžky). Zavedie sa efektívna (vzperná) dĺžka $L_{\mathrm{ef}} = K L$ a vzorec sa zovšeobecní na

$$P_{\mathrm{kr}} = \frac{\pi^2 E I}{(K L)^2}.$$

$K$ vyjadruje, akú dĺžku „vidí" prút pri vzpere v porovnaní s kĺbovo uloženým. Menšie $K$ = tuhšie uloženie = vyššie kritické zaťaženie. Pozor: toto $K$ nemá nič spoločné s tuhostnou maticou $\mathbf{K}$ z ďalšej sekcie — je to len historicky rovnaké písmeno.

Dva stĺpce $K$ majú svoj dôvod. Teoretické hodnoty vychádzajú z presného riešenia pre dokonalé uloženie (dokonalý kĺb, dokonalé votknutie). Odporúčané návrhové hodnoty (AISC 360, komentár, tabuľka C-A-7.1) sú zámerne vyššie — reálny spoj nikdy nie je nekonečne tuhý, takže sa počíta konzervatívnejšie. [2]

A platí to len za predpokladov, ktoré k Eulerovi patria: ideálne priamy, centricky tlačený, prizmatický prút, lineárne elastický materiál pod medzou úmernosti, malé deformácie až po bod vybočenia, bez imperfekcií — a len pre dostatočne štíhle prúty. Krátke a tučné zlyhajú skôr plasticky a Euler by ich nadhodnotil; tam sa prechádza na inelastické krivky (napr. AISC 360, kap. E). [1][2]

Základný Eulerov prípad (kĺb–kĺb, K = 1). Súčiniteľ vzpernej dĺžky K pretvára vzorec na ostatné uloženia cez efektívnu dĺžku KL.

Dva mechanizmy: bifurkácia verzus limitný bod

Nie každá strata stability vyzerá rovnako. [1][4]

Bifurkačný vzper je Eulerov prípad. Dokonalý prút sedí na priamej rovnovážnej ceste — stláča sa, ale ostáva rovný. Pri kritickom zaťažení sa objaví druhá, vybočená rovnovážna cesta a konštrukcia na ňu „preskočí". Je to vetvenie rovnováhy — do kritického bodu jedna cesta, za ním dve.

Limitný bod (snap-through) nemá vetvenie. Zaťažovacia cesta plynule rastie, dosiahne maximum a preklopí sa — konštrukcia „presnapne" do vzdialenejšieho tvaru. Typické pre plytké oblúky a škrupiny: plochá kupola zatlačená zhora v istom bode náhle presnapne dovnútra. Na krivke sila–posun to vidno ako vrchol, za ktorým sila klesá.

Dva mechanizmy straty stability. Bifurkácia (vľavo) vetví rovnovážnu cestu v kritickom bode; snap-through (vpravo) dosiahne maximum a preklopí sa. Lineárna analýza vlastných hodnôt nájde len prvý typ.

Rozdiel nie je akademický. Bifurkáciu spoľahlivo nájde lineárna analýza vlastných hodnôt; limitný bod nie — ten treba dohnať nelineárne, lebo je to jav veľkých deformácií. A do toho vstupuje tretí rozmer: lokálna verzus globálna nestabilita. Globálna preklopí celú konštrukciu (a je to reálna medza únosnosti). Lokálna prepadne len oblasť, kým zvyšok ešte nesie — a solver na nej môže umrieť, hoci konštrukcia ako celok ešte zďaleka nekolabuje. Odlíšiť tieto dve je v praxi to najdôležitejšie rozhodnutie.

Imperfekčná citlivosť: prečo realita nedosahuje Eulera

Eulerov vzorec platí pre dokonalú konštrukciu. Reálna má odchýlky tvaru, excentricitu zaťaženia, zvyškové napätia — a tie rozhodujú viac, než by človek čakal. Prút znesie asi 90 % teoretickej hodnoty. Tenká škrupina v tlaku je iný svet: valec zaťažený osovo unesie bežne len 20–30 % toho, čo predpovedá lineárna teória, a to rozpätie riadia práve imperfekcie. Preto sa v návrhu škrupín pracuje so znižujúcimi súčiniteľmi (knock-down factors) — teoretická hodnota je horný, optimistický odhad. [4][6]

Southwellova metóda: ako z podkritických dát získať kritickú silu

S imperfekciami súvisí elegantný nástroj, ktorý som použil ako nezávislé overenie kritického zaťaženia. Southwellova metóda ho odhadne bez toho, aby ste konštrukciu museli dotlačiť až do vybočenia — pracuje len s podkritickými meraniami. [3]

Odvodenie (pre kĺbovo uložený prút)

Vezmite kĺbovo uložený prút s malou počiatočnou geometrickou imperfekciou v tvare prvého vlastného tvaru (polvlna sínusu) s amplitúdou $a_0$ v strede. Označte pomer zaťaženia $\alpha = P / P_{\mathrm{kr}}$.

Osové zaťaženie počiatočné vybočenie zosilňuje. Celkový priehyb v strede rastie ako

$$\Delta = \frac{a_0}{1 - \alpha}.$$

Prídavný priehyb (nad počiatočný, nezaťažený tvar) je teda

$$\delta = \Delta - a_0 = \frac{a_0\,\alpha}{1 - \alpha}.$$

Z tejto rovnice po vydelení $P$ a úprave vypadne lineárny vzťah — jadro celej metódy:

$$\frac{\delta}{P} = \frac{1}{P_{\mathrm{kr}}}\,\delta + \frac{a_0}{P_{\mathrm{kr}}}.$$

To znamená: keď vynesiete $\delta/P$ na zvislú os proti $\delta$ na vodorovnú os, dostanete priamku, ktorej

• sklon $= 1 / P_{\mathrm{kr}}$, teda kritické zaťaženie $P_{\mathrm{kr}} = 1 / \text{sklon}$,
• úsek na osi $= a_0 / P_{\mathrm{kr}}$, teda amplitúda imperfekcie $a_0 = \text{úsek} \times P_{\mathrm{kr}} = \text{úsek} / \text{sklon}$.

Z priamky preloženej podkritickými bodmi sa čítajú dve veci naraz: prevrátený sklon dáva kritické zaťaženie, úsek na zvislej osi amplitúdu počiatočnej imperfekcie. Zakrivené body pri nízkych zaťaženiach do preloženia nepatria.

Ako to aplikovať v praxi

1. Zaťažujte prút postupne a pri niekoľkých hladinách $P$ v podkritickej oblasti merajte bočný priehyb v strede $\delta$ (ideálne prídavný, teda voči nezaťaženému tvaru).
2. Pre každý bod spočítajte podiel $\delta/P$.
3. Vyneste body $(\delta\,;\ \delta/P)$ a preložte nimi priamku (lineárna regresia).
4. Kritické zaťaženie je prevrátená hodnota sklonu. Úsek na zvislej osi dáva amplitúdu imperfekcie $a_0$.
5. Kontrola kvality: body musia ležať na priamke. Zakrivenie pri nízkych $P$ býva chyba nuly merania priehybu — tie počiatočné body vynechajte a preložte len lineárnu časť.

Existuje aj variant s prehodenými osami — $\delta$ proti $\delta/P$ — kde je sklon rovný priamo $P_{\mathrm{kr}}$. Výsledok je identický, len sa pozeráte na prevrátenú priamku.

Poznámka k rozsahu platnosti: pôvodná Southwellova metóda je odvodená pre elastický bifurkačný vzper prizmatického prúta. Pre plasticitu, dosky a škrupiny existujú rozšírenia (napr. Donnell, 1938), ktoré menia, ktorú veličinu proti ktorej vynášať. [3]

Ako sa vzper počíta v MKP — tri úrovne

Lineárna analýza vlastných hodnôt (eigenvalue buckling). Matematicky je vzper problém vlastných hodnôt. V najbežnejšom, učebnicovom tvare sa zapisuje ako

$$(\mathbf{A} - \lambda \mathbf{I})\,\mathbf{x} = \mathbf{0}$$

a netriviálne riešenie $\mathbf{x}$ existuje len vtedy, keď je determinant matice nulový — to je charakteristická rovnica $\det(\mathbf{A} - \lambda \mathbf{I}) = 0$, ktorej korene sú vlastné hodnoty $\lambda$. Pri vzpere má úloha zovšeobecnený tvar

$$(\mathbf{K}_E + \lambda \mathbf{K}_G)\,\boldsymbol{\varphi} = \mathbf{0},$$

kde $\mathbf{K}_E$ je elastická (klasická) tuhostná matica a $\mathbf{K}_G$ geometrická tuhostná matica závislá od napätosti v konštrukcii. Podmienka netriviálneho riešenia

$$\det(\mathbf{K}_E + \lambda \mathbf{K}_G) = 0$$

dá kritické súčinitele zaťaženia $\lambda_i$ a im prislúchajúce tvary vybočenia $\boldsymbol{\varphi}_i$. Najmenšie kladné $\lambda$ je hľadaný súčiniteľ kritického zaťaženia — vynásobením ním dostanete $P_{\mathrm{kr}}$. Je to rýchle, dá horný odhad a tvary vybočenia, ale ignoruje imperfekcie, veľké deformácie aj nelineárny materiál. Preto pri škrupinách nadhodnocuje, niekedy násobne. Beriem ju ako diagnostiku a strop, nie ako výsledok. [5]

Nelineárna statika (NLGEOM) s postupným zaťažovaním. Zachytí limitný bod, lebo počíta veľké deformácie a mení tuhosť s deformáciou. Problém je, že na limitnom bode zamrzne — statická implicitná schéma nevie prejsť cez pokles sily bez pomoci.

Arc-length (Riks) alebo kvázi-statická dynamika. Toto cez limitný bod prejde. Arc-length rieši zaťaženie a posun spoločne, takže zvládne aj klesajúcu vetvu; kvázi-statický transient necháva jav prebehnúť dynamicky, ako sa deje v realite. Sem patrí aj pravidlo, že kvázi-statika platí, kým kinetická energia ostane malá voči deformačnej.

Postupnosť je jasná: eigenvalue buckling na rýchlu diagnostiku a strop, nelineárna statika na nájdenie limitného bodu, arc-length alebo dynamika na prechod cezeň a postkritickú cestu.

Kde teória narazí na prax

Presne túto reťaz som prešiel v spomínanom prípade tenkej škrupiny. Eigenvalue buckling mi dal súčinitele pod 1 a povedal „pozor, tu je vzper" — ale ako škrupinovému hornému odhadu som mu neveril ako výsledku. Southwell z podkritických dát potvrdil kritické zaťaženie nezávisle. A rozdiel medzi lokálnou a globálnou nestabilitou rozhodol o tom, či konvergenčné zlyhanie znamená kolaps, alebo len limit statickej metódy.

Teória vzperu nie je len kapitola z učebnice. Je to rámec, ktorý zo série varovných hlásení solvera spraví zrozumiteľný fyzikálny príbeh — a ochráni vás pred dvomi opačnými chybami: uveriť optimistickému lineárnemu číslu, alebo vyhlásiť kolaps tam, kde zlyhala len metóda.

Ako to vyzeralo na konkrétnej úlohe, je v článku Keď MKP výpočet stratí stabilitu.

Zdroje

1. Timoshenko, S. P., Gere, J. M. Theory of Elastic Stability, 2. vydanie, McGraw-Hill, 1961.
2. AISC 360, Specification for Structural Steel Buildings, komentár, tabuľka C-A-7.1 (teoretické a odporúčané návrhové hodnoty súčiniteľa vzpernej dĺžky).
3. Southwell, R. V. On the Analysis of Experimental Observations in Problems of Elastic Stability, Proc. Roy. Soc. London, Series A, 135, 601–616 (1932). DOI 10.1098/rspa.1932.0055. Rozšírenie: Donnell, L. H., 1938.
4. Bažant, Z. P., Cedolin, L. Stability of Structures: Elastic, Inelastic, Fracture and Damage Theories, Oxford University Press.
5. Cook, R. D., Malkus, D. S., Plesha, M. E., Witt, R. J. Concepts and Applications of Finite Element Analysis, 4. vydanie, Wiley, 2002; Bathe, K.-J. Finite Element Procedures, Prentice Hall, 1996.
6. NASA SP-8007, Buckling of Thin-Walled Circular Cylinders (znižujúce súčinitele pre škrupiny).