arrow
Return

Matrix-free higher-order finite element methods for hyperelasticity

delete2025-02-01
delete0
delete
OA
AI
R
Richard Schussnig *
N
Niklas Fehn
P
Peter Münch
M
Martin Kronbichler
DOI:10.1016/j.cma.2024.117600delete
deleteOriginal
deleteShare
deleteSave
View PDF
Abstract

Abstract

En 中文
This work presents a matrix-free finite element solver for finite-strain elasticity adopting an hpmultigrid preconditioner. Compared to classical algorithms relying on a global sparse matrix, matrix-free solution strategies significantly reduce memory traffic by repeated evaluation of the finite element integrals. Following this approach in the context of finite-strain elasticity, the precise statement of the final weak form is crucial for performance, and it is not clear a priori whether to choose problem formulations in the material or spatial domain. With a focus on hyperelastic solids in biomechanics, the arithmetic costs to evaluate the material law at each quadrature point might favor an evaluation strategy where some quantities are precomputed in each Newton iteration and reused in the Krylov solver for the linearized problem. Hence, we discuss storage strategies to balance the compute load against memory access in compressible and incompressible neo-Hookean models and an anisotropic tissue model. Additionally, numerical stability becomes increasingly important using lower/mixed-precision ingredients and approximate preconditioners to better utilize modern hardware architectures. Application of the presented method to a patient-specific geometry of an iliac bifurcation shows significant speed-ups, especially for higher polynomial degrees, when compared to alternative approaches with matrix-based geometric or black-box algebraic multigrid preconditioners.
Keywords:
Finite-strain problem
Matrix-free
Finite-element method
Hyperelasticity
Geometric multigrid
AI Summary

AI Summary

Key information extracted from the uploaded paper, including a brief overview, abstract, background, key highlights, visual analysis, and future outlook.

Journal

Computer Methods in Applied Mechanics and Engineering cover
Computer Methods in Applied Mechanics and Engineering
IF:
7.3
Papers:
1.3W
Citations:
5.6W

Organization

U
University of Augsburg
Scholars:
3.6K
Papers: 3.0K
Citations: 5.2K
U
uppsala university
Scholars:
3.7W
Papers: 3.4W
Citations: 47
R
ruhr university bochum
Scholars:
2.3W
Papers: 1.9W
Citations: 14
researcher View more organizations