A Hexahedral Multigrid Approach for Simulating Cuts in Deformable Objects

Abstract

We present a hexahedral finite element method for simulating cuts in deformable bodies using the corotational formulation of strain at high computational efficiency. Key to our approach is a novel embedding of adaptive element refinements and topological changes of the simulation grid into a geometric multigrid solver. Starting with a coarse hexahedral simulation grid, this grid is adaptively refined at the surface of a cutting tool until a finest resolution level, and the cut is modeled by separating elements along the cell faces at this level. To represent the induced discontinuities on successive multigrid levels, the affected coarse grid cells are duplicated and the resulting connectivity components are distributed to either side of the cut. Drawing upon recent work on octree and multigrid schemes for the numerical solution of partial differential equations, we develop efficient algorithms for updating the systems of equations of the adaptive finite element discretization and the multigrid hierarchy. To construct a surface that accurately aligns with the cuts, we adapt the splitting cubes algorithm to the specific linked voxel representation of the simulation domain we use. The paper is completed by a convergence analysis of the finite element solver and a performance comparison to alternative numerical solution methods. These investigations show that our approach offers high computational efficiency and physical accuracy, and that it enables cutting of deformable bodies at very high resolutions.

Development of three-dimensional patient face model that enables real-time collision detection and cutting operation for a dental simulator

“…The GPGPU is available on the personal computer instead of supercomputer and easy to implement with the C language. In addition, the parallel computing algorithm with the GPGPU has already been implemented in many applications including octree algorithms 29,30) . We are working on implementing GPGPU architecture for octree algorithms in a developed VR simulator.…”

Interactive Cutting of Thin Deformable Objects

“…Based on the octree structure, Dick et al [24] connected coarse hexahedral mesh to fine hexahedral mesh. The cutting was simulated by disconnecting the links between fine hexahedral elements.…”

Surgical Augmented Reality with Topological Changes

“…However, the detection method we propose is independent of the separation algorithm, other efficient algorithms like [7] could be used as well.…”

“…Handling such mechanical models and mesh operations implies elevated computational costs and makes it challenging to maintain real-time performance, that is required by augmented reality applications. The approach presented in [7] addresses these issues. The method is based on the composite finite element method, that embeds a fine grid into a coarse uniform hexahedral grid.…”