On convergence of various iterative linear solvers in heterophase elastoplastic media deformation models

Simulation of the process of deformation of heterophase media is necessary for substantiating technological processes of composite forming. Representative volumes of the materials are modeled by a conglomerate of elastic and elastoplastic bodies. Numerical implementations of the models typically use the finite element method (FEM). As a part of the FEM computational procedure, it is necessary to solve a large system of simultaneous linear algebraic equations. When dealing with fine meshes, a Krylov subspace iterative solver is typically used. We present a comparison in terms of convergence of various iterative solvers in a model problem of elastoplastic deformation of a heterophase representative volume. The volume models a metal matrix composite based on the AMg6 aluminum alloy and 10 vol.% of silicone carbide reinforcement. It appears that a relatively rare method, namely QMR, shows the best results.