TY - GEN A1 - Sander, Oliver A1 - Klapproth, Corinna A1 - Youett, Jonathan A1 - Kornhuber, Ralf A1 - Deuflhard, Peter T1 - Towards an Efficient Numerical Simulation of Complex 3D Knee Joint Motion N2 - We present a time-dependent finite element model of the human knee joint of full 3D geometric complexity. Its efficient numerical simulation requires advanced numerical algorithms that have been developed just recently. Up to now, the model comprises bones, cartilage, and the major ligaments (patella and menisci are still missing). Bones (femur, tibia, and fibula) are modelled by linear elastic materials, cartilage by viscoelastic materials, ligaments by one-dimensional so-called Cosserat rods. In order to capture the dynamical contact problems correctly, we solve the full PDEs of elasticity in the presence of strict contact inequalities. For the total spatio-temporal discretization we apply a method of layers approach (first time, then space discretization). For the time discretization of the elastic and viscoelastic parts, we apply 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, and the arising large algebraic systems are solved by a recent fast multigrid solver, the truncated non-smooth Newton multigrid method. T3 - ZIB-Report - 12-06 KW - contact problem KW - contact-stabilized Newmark method KW - domain decomposition KW - energy-momentum method KW - geodesic finite elements KW - knee model Y1 - 2012 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-14512 SN - 1438-0064 ER - TY - GEN A1 - Klapproth, Corinna T1 - The Contact-Stabilized Newmark Method - Consistency Error of a Spatiotemporal Discretization N2 - The paper considers an improved variant of the contact-stabilized Newmark method by Deuflhard et al., which provides a spatiotemporal numerical integration of dynamical contact problems between viscoelastic bodies in the frame of the Signorini condition. Up no now, the question of consistency in the case of contact constraints has been discussed for time integrators in function space under the assumption of bounded total variation of the solution. Here, interest focusses on the consistency error of the Newmark scheme in physical energy norm after discretization both in time and in space. The resulting estimate for the local discretization error allows to prove global convergence of the Newmark scheme under an additional assumption on the active contact boundaries. T3 - ZIB-Report - 12-18 KW - dynamical contact problems, contact-stabilized Newmark method, consistency error Y1 - 2012 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-15198 SN - 1438-0064 ER - TY - THES A1 - Klapproth, Corinna T1 - Adaptive numerical integration for dynamical contact problems Y1 - 2012 ER -