Plastic hexahedral FEM for surgical simulation.

purposeSoft-tissue manipulations, such as collecting, stretching or tearing tissue, are a common component of surgery. When too much force is applied, these manipulations result in a residual plastic deformation that surgeons should be aware of and that should be modeled by surgical simulation.

methodsMany tissues, vessels and organs can be modeled as offsets of curved simple shapes with primary directions, e.g., radial and axial for cylinders yield a rectangular mesh whose normal offset naturally yields a hexahedral mesh that can serve as a thick shell. Other organs are easy to embed into and deform following a hex mesh. We extend existing code for the volumetric finite element method (FEM) to model tissue plasticity as hexahedral thick shells or embedded organs. Specifically, the work extends the open source Simulation Open Framework Architecture and its newest hyperelastic deformation addition, Caribou, with focus on surgical simulation. The extension factors deformation gradients into (corotational or hyperelastic) elastic factors and plastic factors and enforces volume preservation. Limits on per-element twist, twist torque, material hardening and bounds on plasticity where elements invert avoid the need for re-meshing.

resultsOur hexahedral FEM avoids the biased outcomes of asymmetric coarsely-partitioned tetrahedral FEM. Caribou's hyperelastic FEM is extended to hex-FEM stretching plasticity. Our high-order accurate blended-vertex deformation enables coarse hex meshes to model large plastic rotational and stretch deformations without re-meshing. We compare a vertex-blended to a cell-centered piecewise constant approach; contrast plasticity based on corotational FEM and hyperelastic FEM; and test the computation under mesh refinement. The volume is preserved also for large deformations.

conclusionOn-the-fly generated hexahedral meshes can directly be used as finite element domains for plastic deformation based on corotational or hyperelastic elasticity. The outcome is suitable for surgical simulation.