Contact modeling from images using cut finite element solvers

• Published in: Advanced modeling and simulation in engineering sciences Vol. 8; no. 1; pp. 1 - 23
• Main Authors: Claus, S., Kerfriden, P., Moshfeghifar, F., Darkner, S., Erleben, K., Wong, C.
• Format: Journal Article
• Language: English
• Published: Cham Springer International Publishing 01-06-2021; Springer Nature B.V; Springer; SpringerOpen
• Subjects: Algorithms; Classical and Continuum Physics; Computational Science and Engineering; Computed tomography; Computer simulation; Contact; CutFEM; Digital imaging; Engineering; Femur; Ghost penalty; High-Performance Adaptive Solution Algorithms and Advanced Digital Pipelines; Image segmentation; Information flow; Joint geometry; LaTIn; Mathematical analysis; Mechanics; Medical imaging; Nonconforming finite element; Nonlinear algorithm; Physics; Recent Developments in Unfitted Finite Element Methods: Numerical Analysis; Research Article; Robustness (mathematics); Structural mechanics; Theoretical and Applied Mechanics; LaTIn; Ghost penalty; Nonconforming finite element; Nonlinear algorithm; CutFEM; Contact; LATIN; CONTACT; NONCONFORMING FINITE ELEMENT; GHOST PENALTY; CUTFEM; NONLINEAR ALGORITHM
• ISSN: 2213-7467; 2213-7467
• DOI: 10.1186/s40323-021-00197-2

This paper presents a robust digital pipeline from CT images to the simulation of contact between multiple bodies. The proposed strategy relies on a recently developed immersed finite element algorithm that is capable of simulating unilateral contact between solids without meshing (Claus and Kerfriden in Int J Numer Methods Eng 113(6):938–966, 2018). It was shown that such an approach reduces the difficulties associated with the digital flow of information from analytically defined geometries to mechanical simulations. We now propose to extend our approach to include geometries, which are not defined mathematically but instead are obtained from images, and encoded in 3D arrays of voxels. This paper introduces two novel elements. Firstly, we reformulate our contact algorithm into an extension of an augmented Lagrangian CutFEM algorithm. Secondly, we develop an efficient algorithm to convert the surface data generated by standard segmentation tools used in medical imaging into level-set functions. These two elements give rise to a robust digital pipeline with minimum user intervention. We demonstrate the capabilities of our algorithm on a hip joint geometry with contact between the femur and the hip bone.