An efficient algorithm for biomechanical problems based on a fully implicit nested Newton solver
-
Published:2022
Issue:2
Volume:49
Page:183-221
-
ISSN:1450-5584
-
Container-title:Theoretical and Applied Mechanics
-
language:en
-
Short-container-title:Theor appl mech (Belgr)
Author:
Knodel Markus1, di Stefano2, Nägel Arne1, Grillo Alfio3
Affiliation:
1. Goethe Center for Scientific Computing (GCSC), Universität Frankfurt, Frankfurt am Main, Germany 2. Dipartimento di Ingegneria Civile, Ambientale, del Territorio, Edile e di Chimica, Politecnico di Bari, Bari, Italy 3. Dipartimento di Scienze Matematiche “G.L. Lagrange” (DISMA), Politecnico di Torino, Torino, Italy
Abstract
Numerical simulations of the dynamics of soft biological tissues are highly non-trivial because tissues generally exhibit complex biological response to external and internal actions, including large deformations and remodeling. Combining the advantages of globally implicit approach (GIA) solvers with the general applicability of the semi-implicit General Plasticity Algorithm (GPA), introduced by some of us some years ago, we present a new, efficient plasticity algorithm, which we call Bio Mechanics Basis Plasticity Algorithm (BMBPA). This is fully implicit, based on a nested Newton solver, and naturally suited for massively parallel computations. The Bilby?Kr?ner?Lee (BKL) multiplicative decomposition of the deformation gradient tensor is employed to introduce the unknowns of our model. We distinguish between global and local unknowns, associated with local and global equations, which are connected by means of a resolution function. The BMBPA asks for very few conditions to be applied and thus can be easily employed to solve several types of biological and biomechanical problems. We demonstrate the efficacy of BMBPA by performing two numerical experiments of a monophasic model of fiber-reinforced tissues. In one case, we consider the shear-compression test of a cubic specimen of tissue, while, in the other case, we focus on the unconfined compression test of a cylinder. The BMBPA is capable of solving the deformation and the remodeling of anisotropic biological tissues by employing a computation time of hours, while the GPA, applied to the same problems as the BMBPA, needs a substantially longer amount of time. All computations were performed in parallel and, within all tests, the performance of the BMBPA displayed substantially higher than the one of the GPA. The results of our simulations permit to study the overall mechanical behavior of the considered tissue and enable further investigations in the field of tissue biomechanics.
Publisher
National Library of Serbia
Subject
Applied Mathematics,Mechanical Engineering,Computational Mechanics
Reference92 articles.
1. J. C. Simo, Numerical analysis and simulation of plasticity, in: P. G. Ciarlet, J. L. Lions (eds.), Handb. Numer. Anal. 6 (1998), 179-181, Elsevier Science, Amsterdam. 2. J. Alberty, C. Carstensen, D. Zarrabi, Adaptive numerical analysis in primal elastoplasticity with hardening., Comput. Methods Appl. Mech. Eng. 171 (1999), 175-204. 3. D. Ambrosi, M. Ben Amar, C. J. Cyron, A. De Simone, A. Goriely, J. D. Humphrey, E. Kuhl, Growth and remodelling of living tissues: perspectives, challenges and opportunities, Journal of The Royal Society Interface 16(157) (2019), 20190233. 4. D. Ambrosi, G. A. Ateshian, E. M. Arruda, and et al., Perspectives on biological growth and remodeling, J. Mech. Phys. Solids 59(4) (2011), 863-883. 5. D. Ambrosi, F. Guana, Stress-modulated growth, Math. Mech. Solids 12 (2007), 319-342.
|
|