Adaptive Numerical Integration of Dynamical Contact Problems

The present thesis deals with the e¿cient adaptive numerical integration of dynamical contact problems. This is the core problem in the fast and robust simulation of stresses arising in a real patient¿s knee joint for di¿erent kinds of loading situations. The topic is of high interest in the ¿eld of computer-assisted therapy planning, which aims at the generation of a ¿virtual patient¿. This tool allows the design of e¿ective treatment options and precise surgery strategies within a clinical environment. Potential tasks of patient-speci¿c techniques in orthopedics are osteotomic interventions and the construction and selection of implants or ¿xation devices. For realistic predictions of therapeutical manipulations, numerical simulation and optimization are applied on a detailed three-dimensional geometry of the individual patient¿s knee (obtained from anatomical CT or MRT image data). In view of a reasonable clinical application, the necessary computations have to be performed on local workstations in clinics within short time frames. Moreover, the solutions have to be resilient enough to serve as a basis for responsible medical decisions. With regard to these conditions, highest level requirements have to be set on the e¿ciency and accuracy of the applied numerical techniques. The appropriate approach to cope with this mathematical challenge is the construction of an adaptive numerical integrator for the dynamical contact problem. For this purpose, dynamical contact problems have to be analyzed precisely from both the analytical and the numerical point of view.