Scalable algorithms for solving elasto-plastic problems

Abstract

V této disertační práci představíme algoritmy pro efektivní paralelní implementaci elastoplastických problémů se zpevněním, založených na metodě rozložení oblastí zvané TFETI (Total Finite Element Tearing and Interconnecting). Uvažujeme tři odlišné asociativní elastoplastické modely: von Misesův model s izotropním zpevněním, von Misesův model s kinematickým zpevněním a Drucker-Prager perfektně plastický model. Tyto modely jsou diskretizovány v čase pomocí implicitní Eulerovy metody a každý časový krok elastoplastického modelu je diskretizován v prostoru pomocí metody konečných prvků. Tímto dostaneme pro každý časový krok systém nelineárních rovnic pro rovnostní omezení nebo systém nelineárních nerovnic se silně nehladkým a silně monotónním operátorem. Nehladká Newtonova metoda se aplikuje na tento nelineární systém za účelem jeho linearizace. Odpovídající linearizované problémy vyskytující se v Newtonových iteracích jsou paralelně nebo sekvenčně řešeny pomocí metody TFETI. Zmíněné algoritmy založené na TFETI byly implementovány v Matlabovském paralením prostředí a jejich funkčnost byla ilustrovaná na 2D a 3D elastoplastických příkladech. Numerické výsledky prezentujeme pro různé časové diskretizace a různé konečno-prvkové sítě, pozorujeme kvadratickou konvergenci nehladké Newtonovy metody. Také demonstrujeme paralelní a numerickou škálovatelnost těchto algoritmů na řešení úloh elastoplasticity.

In the thesis we propose an algorithm for the efficient parallel implementation of elastoplastic problems with hardening based on the so-called TFETI (Total Finite Element Tearing and Interconnecting) domain decomposition method. We consider three different associated elasto-plastic models: the von Mises model with isotropic hardening, the von Mises model with kinematic hardening and the Drucker-Prager perfectly plastic model. Such models are discretized by the implicit Euler method in time and the consequent one time step elasto-plastic problem by the finite element method in space. The latter results in a system of nonlinear equations for equality constraints or system of nonlinear inequations with a strongly semismooth and strongly monotone operator. The semismooth Newton method is applied to solve this nonlinear system. Corresponding linearized problems arising in the Newton iterations are solved in parallel or sequentially by the above mentioned TFETI domain decomposition method. The proposed TFETI based algorithm was implemented in Matlab parallel environment and its performance is illustrated on 2D and 3D elasto-plastic benchmarks. Numerical results for different time discretizations and mesh levels are presented and discussed and a local quadratic convergence of the semismooth Newton method is observed. We also demonstrate parallel and numerical scalability of the proposed algorithms for solving elasto-plasticity.