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.
We present a discretization for dynamic large deformation contact problems without friction. Our model is based on Hamilton’s principle, which avoids the explicit appearance of the contact forces. The resulting differential inclusion is discretized in time using a modified midpoint rule. This modification, which concerns the evaluation of the generalized gradient, allows to achieve energy dissipativity. For the space discretization we use a dual-basis mortar method. The resulting spatial algebraic
problems are nonconvex minimization problems with nonconvex inequality constraints. These can be solved efficiently using a trust-region SQP framework with a monotone multigrid inner solver.
We present a globally convergent method for the solution of frictionless large deformation contact problems involving hyperelastic materials.
For the discretisation we apply the dual mortar method which is known to be more stable than node-to-segment approaches. The resulting non-convex constrained minimisation problems are solved using a filter–trust-region scheme.
This method combines several techniques from non-linear optimisation to achieve global convergence towards first-order optimal points. To speed up the method inexact linearisations of the non-penetration constraint are used
whenever the current iterate is far away from a critical point. A monotone multigrid method is applied for the fast solution of the constrained Newton
problems.