Abstract:Computational biomechanics plays an important role in biomedical engineering: using modeling to understand pathophysiology, treatment and device design. While experimental evidence indicates that the mechanical response of most tissues is viscoelastic, current biomechanical models in the computational community often assume hyperelastic material models. Fractional viscoelastic constitutive models have been successfully used in literature to capture viscoelastic material response; however, the translation of these models into computational platforms remains limited. Many experimentally derived viscoelastic constitutive models are not suitable for three-dimensional simulations. Furthermore, the use of fractional derivatives can be computationally prohibitive, with a number of current numerical approximations having a computational cost that is <math>O(NT2)</math> and a storage cost that is <math>O(NT)</math>(<math>NT</math> denotes the number of time steps). In this paper, we present a novel numerical approximation to the Caputo derivative which exploits a recurrence relation similar to those used to discretize classic temporal derivatives, giving a computational cost that is <math>O(NT)</math> and a storage cost that is fixed over time. The approximation is optimized for numerical applications, and an error estimate is presented to demonstrate the efficacy of the method. The method, integrated into a finite element solid mechanics framework, is shown to be unconditionally stable in the linear viscoelastic case. It was then integrated into a computational biomechanical framework, with several numerical examples verifying the accuracy and computational efficiency of the method, including in an analytic test, in an analytic fractional differential equation, as well as in a computational biomechanical model problem.

Intrarenal vascular changes in adult patients with recurrent haematuria and loin pain--a clinical, histological and angiographic study.

Fast porous visco-hyperelastic soft tissue model for surgery simulation: application to liver surgery

An energy-consistent discretization of hyper-viscoelastic contact models for soft tissues

A Galerkin least-square stabilisation technique for hyperelastic biphasic soft tissue

Towards Real-time Simulation of Hyperelastic Materials

Multi-compartment poroelastic models of perfused biological soft tissues: implementation in FEniCSx

Volume Preserving Simulation of Soft Tissue with Skin

Stable Neo-Hookean Flesh Simulation

A selective smoothed finite element method with visco‐hyperelastic constitutive model for analysis of biomechanical responses of brain tissues

An accurate, robust, and efficient finite element framework for anisotropic, nearly and fully incompressible elasticity

Multifrequency 3D Elasticity Reconstruction withStructured Sparsity and ADMM

A Non-Linear Viscoelastic Tensor-Mass Model for Surgery Simulation

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

An Anisotropic Micro-Sphere Approach Applied to the Modelling of Soft Biological Tissues

Three-dimensional numerical simulation of soft-tissue wound healing using constrained-mixture anisotropic hyperelasticity and gradient-enhanced damage mechanics

Stabler Neo-Hookean Simulation: Absolute Eigenvalue Filtering for Projected Newton

Numerical Modelling of the Brain Poromechanics by High-Order Discontinuous Galerkin Methods

Towards real-time finite-strain anisotropic thermo-visco-elastodynamic analysis of soft tissues for thermal ablative therapy

An efficient multi-field dynamic model for 3D wave propagation in saturated anisotropic porous media

Hyperelastic Energy Densities for Soft Biological Tissues: A Review

An efficient and accurate method for modeling nonlinear fractional viscoelastic biomaterials