A direct Jacobian total Lagrangian explicit dynamics finite element algorithm for real-time simulation of hyperelastic materials

by   Jinao Zhang, et al.

This paper presents a novel direct Jacobian total Lagrangian explicit dynamics (DJ-TLED) finite element algorithm for real-time nonlinear mechanics simulation. The nodal force contributions are expressed using only the Jacobian operator, instead of the deformation gradient tensor and finite deformation tensor, for fewer computational operations at run-time. Owing to this proposed Jacobian formulation, novel expressions are developed for strain invariants and constant components, which are also based on the Jacobian operator. Results show that the proposed DJ-TLED consumed between 0.70x and 0.88x CPU solution times compared to state-of-the-art TLED and achieved up to 121.72x and 94.26x speed improvements in tetrahedral and hexahedral meshes, respectively, using GPU acceleration. Compared to TLED, the most notable difference is that the notions of stress and strain are not explicitly visible in the proposed DJ-TLED but embedded implicitly in the formulation of nodal forces. Such a force formulation can be beneficial for fast deformation computation and can be particularly useful if the displacement field is of primary interest, which is demonstrated using a neurosurgical simulation of brain deformations for image-guided neurosurgery. The present work contributes towards a comprehensive DJ-TLED algorithm concerning isotropic and anisotropic hyperelastic constitutive models and GPU implementation.


page 16

page 17

page 19

page 23

page 24

page 25


Biomechanical modeling and computer simulation of the brain during neurosurgery

Computational biomechanics of the brain for neurosurgery is an emerging ...

A discrete elasticity complex on three-dimensional Alfeld splits

We construct conforming finite element elasticity complexes on the Alfel...

Suite of Meshless Algorithms for Accurate Computation of Soft Tissue Deformation for Surgical Simulation

The ability to predict patient-specific soft tissue deformations is key ...

An isogeometric finite element formulation for geometrically exact Timoshenko beams with extensible directors

An isogeometric finite element formulation for geometrically and materia...

AMPS: A Real-time Mesh Cutting Algorithm for Surgical Simulations

We present the AMPS algorithm, a finite element solution method that com...

Please sign up or login with your details

Forgot password? Click here to reset