Keď MKP výpočet stratí stabilitu: ako odlíšiť fyziku od numeriky

Nelineárny statický výpočet sa zastaví. Solver hlási nekonvergenciu, časový krok sa rúti k nule, hlásenia o penetrácii sa množia. Výpočtársky reflex je: zjemniť sieť, pridať tlmenie, zmenšiť krok. Niekedy je to správne. Inokedy je to najdrahšia chyba vo výpočte — pretože zlyhaná konvergencia nie je porucha modelu, ale správa od konštrukcie, že stráca stabilitu. [1][2] Kto ju prekryje umelým tlmením, dostane číslo, ktoré nevie obhájiť.

Nedávno som riešil presne takýto prípad: tenkostennú priehľadnú polymérovú škrupinu z krehkého akrylátu, lokálne zaťažovanú cez kontakt s druhým telesom. Výpočet mi opakovane umieral. Toto je blog o tom, ako nad podobným problémom premýšľať a ako postupovať až k fyzikálne správnemu riešeniu.

Úloha: tenká škrupina lokálne zaťažovaná cez kontakt s iným telesom. Zaťaženie vstupuje do škrupiny cez toto teleso, nie priamo.

Dve tváre straty stability

Prvý a najdôležitejší krok vôbec nie je technický. Je to otázka: akú stratu stability vlastne vidím?

Globálny limitný bod (snap-through celej konštrukcie). Zaťaženie dosiahne maximum — konštrukcia už viac neunesie a krivka sila–posun sa „preklopí". Toto je reálne vyčerpanie únosnosti. Ak solver umrie tu, je to medzný stav celej konštrukcie. [2]

Lokálna nestabilita. Odizolovaná oblasť sa prepadne alebo zvlní (kúsok tenkej škrupiny sa zaboril), zatiaľ čo globálne zaťaženie stále rastie. Statická implicitná schéma cez tento skok neprejde — ale to nie je dôkaz kolapsu. Je to limit metódy, nie limit konštrukcie.

Globálna nestabilita preklopí celú konštrukciu a je to medzný stav; lokálna prepadne len oblasť, kým zvyšok nesie ďalej. Rozlíšiť ich je kľúčové rozhodnutie.

Rozdiel je zásadný, lebo si žiada úplne inú reakciu. Globálny limit = mám únosnosť. Lokálna nestabilita = potrebujem metódu, ktorá cez ňu prejde, inak si nechávam ujsť zvyšok zaťaženia. (Teóriu vzperu, bifurkácie a snap-through rozoberám v samostatnom článku Vzper a strata stability.)

Ako som to rozlíšil — čítanie signálov, nie hádanie

Zlyhaná konvergencia sama o sebe nič nehovorí. Diagnostika áno.

Reziduá po uzloch. Pozrel som sa na silové reziduá v jednotlivých uzloch cez iterácie. Ukázalo sa, že reziduum je hyperlokalizované — prevažnú časť L2 normy niesla hŕstka uzlov z niekoľkých desiatok tisíc — a drvivo v normálovom smere. To mi naraz vylúčilo dve hypotézy: globálnu nestabilitu (tá by dala rozptýlené reziduum po celej ploche) aj stick-slip trenie (to je tangenciálne, nie normálové). [4]

Reziduum sústredené do hŕstky uzlov a v normálovom smere — podpis lokálnej kontaktnej nestability, nie globálneho javu ani stick-slip trenia.

Podpis chatteringu. Keď som porovnal „hot spot" medzi iteráciami, presúval sa o jeden rad uzlov. To je typický podpis kontaktného preklápania, nie postupujúceho porušenia.

Skok kontaktu. Počet aktívnych kontaktných alebo cieľových elementov skočil z ustálenej hodnoty takmer na dvojnásobok v jedinej iterácii. Celá plocha zabrala naraz — pod pritláčaným telesom sa škrupina lokálne zaborila a obopla ho, takže do kontaktu vstúpila väčšia plocha. A riadený stupeň voľnosti (predpísané posunutie alebo zaťaženie) tento vnútorný stupeň voľnosti neriadi, takže ho solver nedokáže staticky preklenúť.

Náhly skok počtu kontaktných elementov v jedinom prírastku je podpis lokálneho snapu — celá plocha zabrala naraz. Pri globálnej nestabilite počet kontaktov v tom istom okamihu naopak prepadne: konštrukcia sa odklopí a kontakt stráca.

Tvar krivky tuhosti. Prírastková tuhosť klesla o desiatky percent a potom sa ustálila — nesmerovala k nule. Globálny limitný bod sa ohlasuje tak, že tuhosť ide k nule a sila dosiahne maximum. [2] U mňa sila rástla monotónne, bez náznaku maxima. Záver: lokálna nestabilita, nie globálny limit.

Krížová kontrola: dve nezávislé cesty k rovnakému číslu

Dôvera v nelineárny výsledok neprichádza z jedného behu. Prichádza z toho, že dve nezávislé metódy ukážu na to isté.

Eigenvalue buckling ako diagnostika potvrdil, že konštrukcia je v tomto režime náchylná na vzper — súčinitele zaťaženia mi vyšli hlboko pod 1. A Southwellov odhad z úplne iného behu pristál zhruba na polovici návrhového zaťaženia, kým nelineárny solver umieral prakticky na tej istej úrovni. [6] Keď sa dva nezávislé prístupy zhodnú, že „okolo tejto úrovne sa deje niečo reálne", z podozrenia sa stáva zistenie.

Fyzikálny obraz: tenká škrupina v ťahovom poli

Ešte jeden signál, ktorý sa ľahko označí za chybu a pritom je čistou fyzikou: minimálne hlavné napätie mi kleslo na veľkých plochách takmer na nulu. Tenká škrupina neunesie tlak — namiesto toho sa zvlní, a tlakové pole skolabuje do ťahového. [3] Priehyb dosiahol rádovo desaťnásobok hrúbky panela; som hlboko v geometrickej nelinearite, kde dominuje membránové pôsobenie. Kto tento stav nerozpozná, počíta na modeli, ktorý neplatí.

Nastavenia kontaktu, ktoré beh odblokovali

Aby som vôbec dosiahol skutočnú nestabilitu, musel model najprv prežiť numerické artefakty kontaktu. Zmeny, ktoré ho odblokovali:

Účinok: chybové hlásenia aj penetrácie zmizli úplne, substep, ktorý predtým beh zabíjal, prešiel za pár iterácií, a časový krok namiesto kolapsu rástol.

Stabilizácia — a prečo v tomto prípade stačila

Nelineárna stabilizácia je barlička — umelá disipovaná energia, ktorú je potrebné logicky obhájiť a ktorá môže skresliť výsledok. Poctivejšie cesty pre snap-through existujú: arc-length metóda, ktorá rieši zaťaženie a posun spoločne a prejde aj cez klesajúcu vetvu, alebo kvázi-statický transient, kde jav prebehne dynamicky. [1]

Napokon som ani jednu z nich nepotreboval. Stačilo si strážiť prírastok stabilizačnej energie voči deformačnej ($\Delta \text{STEN} / \Delta \text{SENE}$) cez celú zaťažovaciu dráhu. Kým tento pomer ostáva malý, umelá energia je zanedbateľná oproti tomu, čo v konštrukcii reálne pracuje — a stabilizované riešenie sa dá brať ako regulárna rovnováha, nie ako artefakt tlmenia. Znížením tlmenia na úroveň, kde prírastkový pomer ostal na jednotkách percent, som dostal výsledok, ktorý drží aj bez arc-length či transientu — a kritické zaťaženie mi navyše nezávisle potvrdili eigenvalue buckling aj Southwell.

Poznámka: nízky prírastkový pomer (jednotky percent) je štandardné akceptačné kritérium pre stabilizáciu [5] a spolu s nezávislou zhodou bucklingu a Southwella dáva istotu, že kritické zaťaženie nie je artefakt tlmenia. Arc-length by pri samotnom kritickom zaťažení dal prakticky to isté — rozdiel by ostal v rámci diskretizačnej neistoty. Spresnil by hlavne postkritickú vetvu za limitným bodom, ktorú v tejto úlohe nepotrebujem.

Ako k takým úlohám pristupujem

• Konvergenčné zlyhanie je diagnostika, nie len prekážka — čítam reziduá skôr, než siahnem po tlmení.
• Rozlišujem lokálnu nestabilitu od globálneho limitného bodu: sledujem, či zaťažovacia cesta má maximum a či tuhosť smeruje k nule.
• Dôveru v nelineárny výsledok staviam na nezávislom potvrdení (buckling, Southwell), nie na jednom behu.
• Nula v minimálnom hlavnom napätí na tenkej škrupine je ťahové pole, nie chyba.
• Stabilizácia je prijateľná, ak prírastok stabilizačnej energie ostane malý voči deformačnej — vtedy nemusím siahať po arc-length ani transiente.

Keď mi výpočet nechce skonvergovať, najprv si položím otázku, či mi model nehovorí niečo o konštrukcii. Riešenie často nie je viac tlmenia — je to lepšia diagnostika.

---

Ste konštruktér alebo výpočtár a máte na stole úlohu, s ktorou si neviete rady? Nelineárny model, ktorý nekonverguje, kontakt, čo padá, alebo výsledok, ktorému neveríte — presne takéto veci riešim. Napíšte mi na ales@makewise.sk alebo cez kontaktný formulár a pozrieme sa na to spolu.

Zdroje

1. Crisfield, M. A. Non-Linear Finite Element Analysis of Solids and Structures, Wiley, 1991.
2. Timoshenko, S. P., Gere, J. M. Theory of Elastic Stability, 2. vyd., McGraw-Hill, 1961.
3. Timoshenko, S. P., Woinowsky-Krieger, S. Theory of Plates and Shells, 2. vyd., McGraw-Hill, 1959.
4. Wriggers, P. Computational Contact Mechanics, 2. vyd., Springer, 2006.
5. Ansys, Inc. Mechanical User's Guide / Theory Reference — nelineárna stabilizácia (kritérium stabilizačnej vs. deformačnej energie) a kontaktná diagnostika NLDIAG.
6. Southwell, R. V. On the Analysis of Experimental Observations in Problems of Elastic Stability, Proc. Roy. Soc. London A, 135, 601–616 (1932). DOI 10.1098/rspa.1932.0055.