Towards an Efficient Numerical Simulation of Complex 3D Knee Joint Motion

Please always quote using this URN:urn:nbn:de:0296-matheon-11385
  • We present a time-dependent finite element model of the human knee joint of full 3D geometric complexity together with advanced numerical algorithms needed for its simulation. The model comprises bones, cartilage and the major ligaments, while patella and menisci are still missing. Bones are modeled by linear elastic materials, cartilage by linear viscoelastic materials, and ligaments by one-dimensional nonlinear Cosserat rods. In order to capture the dynamical contact problems correctly, we solve the full PDEs of elasticity with strict contact inequalities. The spatio--temporal discretization follows a time layers approach (first time, then space discretization). For the time discretization of the elastic and viscoelastic parts we use a new contact-stabilized Newmark method, while for the Cosserat rods we choose an energy--momentum method. For the space discretization, we use linear finite elements for the elastic and viscoelastic parts and novel geodesic finite elements for the Cosserat rods. The coupled system is solved by a Dirichlet--Neumann method. The large algebraic systems of the bone--cartilage contact problems are solved efficiently by the truncated non-smooth Newton multigrid method.

Download full text files

Export metadata

Additional Services

Share in Twitter Search Google Scholar
Metadaten
Author:Oliver Sander, Corinna Klapproth, Jonathan Youett, Ralf Kornhuber, Peter Deuflhard
URN:urn:nbn:de:0296-matheon-11385
Referee:Christof Schütte
Document Type:Preprint, Research Center Matheon
Language:English
Date of first Publication:2012/05/18
Release Date:2012/05/18
Tag:contact problem; contact-stabilized Newmark method; domain decomposition; energy-momentum method; geodesic finite elements; knee model
Institute:Freie Universität Berlin
Zuse Institute Berlin (ZIB)
MSC-Classification:65-XX NUMERICAL ANALYSIS / 65Mxx Partial differential equations, initial value and time-dependent initial- boundary value problems / 65M55 Multigrid methods; domain decomposition
Preprint Number:962
Verstanden ✔
Diese Webseite verwendet technisch erforderliche Session-Cookies. Durch die weitere Nutzung der Webseite stimmen Sie diesem zu. Unsere Datenschutzerklärung finden Sie hier.