In this work we discuss the design and implementation of a computer program to simulate the wear on knee implants over a standardized gait cycle.
The contact problems during the gait cycle are formulated using linear elasticity with a linear non-penetration condition.
Archard's wear law is used to estimate long-term wear.
The algebraic problem after discretization is solved using the Truncated Nonsmooth Newton Multigrid (TNNMG) algorithm.
The linear correction step of the algorithm leads to indefinite systems in our application for which we describe a novel heuristic damping algorithm.
Our problem sizes allow using a direct solver for these systems and with some modifications we are able to reuse the Cholesky factorization of the system matrix over multiple iterations of the TNNMG algorithm.
For our application we are able to significantly improve over the performance of IPopt, a state-of-the-art solver using penalty-based methods.
The program is then used for in-silico simulations of wear of two different total knee implant models for which in-vitro results are also available.
Both wear volume and spatial distribution of wear that are computed compare well with the experimental results.
| Identifer | oai:union.ndltd.org:DRESDEN/oai:qucosa:de:qucosa:79216 |
| Date | 18 May 2022 |
| Creators | Burchardt, Ansgar |
| Contributors | Sander, Oliver, Gräser, Carsten, Technische Universität Dresden |
| Source Sets | Hochschulschriftenserver (HSSS) der SLUB Dresden |
| Language | English |
| Detected Language | English |
| Type | info:eu-repo/semantics/publishedVersion, doc-type:doctoralThesis, info:eu-repo/semantics/doctoralThesis, doc-type:Text |
| Rights | info:eu-repo/semantics/openAccess |
Page generated in 0.0022 seconds