Implementace plné linearizace matic v semihladké Newtonově metodě pro dynamické kontaktní problémy

Abstract

Na katedře 470 Vysoké Školy Báňské -- Technické Univerzitě Ostrava jsou v současnosti kontaktní úlohy řešeny pomocí minimalizace kvadratického funkcionálu s rovnostními omezeními a nerovnostními omezeními speciálního tvaru. Problém vyvstává v případě nelinearity materiálu a počítání v přírůstcích, kdy matice soustavy není pozitivně definitní. Jedním z možných řešení je přeformulování úlohy minimalizace na řešení soustavy nelineárních rovnic a přepsání omezujících nerovnostních omezení na rovnostní omezení, ve kterých se vyskytují nehladké funkce. Tímto způsobem přeformulovanou úlohu můžeme řešit pomocí semihladké Newtonovy metody. V takovém případě je však nutné linearizovat (derivovat) matice tuhosti a kontaktní matice. Naším cílem bude implementace semihladké Newtonovy metody do knihovny MatSol, která je vyvíjená na katedře 470, a srovnání výsledků s existujícími benchmarky.

At the department 470 of VŠB -- Technical University of Ostrava, a contact problems are solved by minimalization of quadratic functional with equality constraints and inequality constraints in a special form. The problem arises at the case when the material is nonlinear and when the matrix of solved system isn't positive definite in increment counting. One of the possible ways is to reformulate minimalization problem into problem of solving nonlinearity system and rewrite inequality constraints to equality constraints containing non-smooth functions. Thus rewritten problem could be solved by semismooth Newton method. In this case, linearizations (derivations) of stiffness matrix and contact matrices are required. Our goal is to implementat semismooth Newton method to MatSol library, which is being developed at the department 470 and to compare results with the existing benchmarks.