A Generative Model of Hyperelastic Strain Energy Density Functions for Real-Time Simulation of Brain Tissue Deformation

Author(s):  
Alejandro Granados ◽  
Martin Schweiger ◽  
Vejay Vakharia ◽  
Andrew W. McEvoy ◽  
Anna Miserocchi ◽  
...  
Author(s):  
Alejandro Granados ◽  
Fernando Perez-Garcia ◽  
Martin Schweiger ◽  
Vejay Vakharia ◽  
Sjoerd B. Vos ◽  
...  

Abstract Purpose Estimation of brain deformation is crucial during neurosurgery. Whilst mechanical characterisation captures stress–strain relationships of tissue, biomechanical models are limited by experimental conditions. This results in variability reported in the literature. The aim of this work was to demonstrate a generative model of strain energy density functions can estimate the elastic properties of tissue using observed brain deformation. Methods For the generative model a Gaussian Process regression learns elastic potentials from 73 manuscripts. We evaluate the use of neo-Hookean, Mooney–Rivlin and 1-term Ogden meta-models to guarantee stability. Single and multiple tissue experiments validate the ability of our generative model to estimate tissue properties on a synthetic brain model and in eight temporal lobe resection cases where deformation is observed between pre- and post-operative images. Results Estimated parameters on a synthetic model are close to the known reference with a root-mean-square error (RMSE) of 0.1 mm and 0.2 mm between surface nodes for single and multiple tissue experiments. In clinical cases, we were able to recover brain deformation from pre- to post-operative images reducing RMSE of differences from 1.37 to 1.08 mm on the ventricle surface and from 5.89 to 4.84 mm on the resection cavity surface. Conclusion Our generative model can capture uncertainties related to mechanical characterisation of tissue. When fitting samples from elastography and linear studies, all meta-models performed similarly. The Ogden meta-model performed the best on hyperelastic studies. We were able to predict elastic parameters in a reference model on a synthetic phantom. However, deformation observed in clinical cases is only partly explained using our generative model.


2018 ◽  
Vol 24 (6) ◽  
pp. 1785-1799 ◽  
Author(s):  
Y. Q. Li ◽  
X.-L. Gao

The upper triangular decomposition has recently been proposed to multiplicatively decompose the deformation gradient tensor into a product of a rotation tensor and an upper triangular tensor called the distortion tensor, whose six components can be directly related to pure stretch and simple shear deformations, which are physically measurable. In the current paper, constitutive equations for hyperelastic materials are derived using strain energy density functions in terms of the distortion tensor, which satisfy the principle of material frame indifference and the first and second laws of thermodynamics. Being expressed directly as derivatives of the strain energy density function with respect to the components of the distortion tensor, the Cauchy stress components have simpler expressions than those based on the invariants of the right Cauchy-Green deformation tensor. To illustrate the new constitutive equations, strain energy density functions in terms of the distortion tensor are provided for unconstrained and incompressible isotropic materials, incompressible transversely isotropic composite materials, and incompressible orthotropic composite materials with two families of fibers. For each type of material, example problems are solved using the newly proposed constitutive equations and strain energy density functions, both in terms of the distortion tensor. The solutions of these problems are found to be the same as those obtained by applying the polar decomposition-based invariants approach, thereby validating and supporting the newly developed, alternative method based on the upper triangular decomposition of the deformation gradient tensor.


1999 ◽  
Vol 4 (5) ◽  
pp. 281-285 ◽  
Author(s):  
Uli Bockholt ◽  
Wolfgang Müller ◽  
Gerrit Voss ◽  
Ulrich Ecke ◽  
Ludger Klimek

Author(s):  
D. J. Bang ◽  
E. Madenci

This study concerns the development of peridynamic (PD) strain energy density functions for a Neo-Hookean type membrane under equibiaxial, planar, and uniaxial loading conditions. The material parameters for each loading case are determined by equating the PD strain energy density to that of the classical continuum mechanics. The PD equations of motion are derived based on the Neo-Hookean model under the assumption of incompressibility. Numerical results concern the deformation of a membrane with a defect in the form of a hole, a crack, and a rigid inclusion under equibiaxial, planar, and uniaxial loading conditions. The PD predictions are verified by comparison with those of finite element analysis.


1996 ◽  
Vol 63 (4) ◽  
pp. 869-876 ◽  
Author(s):  
Jiun-Shyan Chen ◽  
Cheng-Tang Wu ◽  
Chunhui Pan

In the first part of this paper a pressure projection method was presented for the nonlinear analysis of structures made of nearly incompressible hyperelastic materials. The main focus of the second part of the paper is to demonstrate the performance of the present method and to address some of the issues related to the analysis of engineering elastomers including the proper selection of strain energy density functions. The numerical procedures and the implementation to nonlinear finite element programs are presented. Mooney-Rivlin, Cubic, and Modified Cubic strain energy density functions are used in the numerical examples. Several classical finite elasticity problems as well as some practical engineering elastomer problems are analyzed. The need to account for the slight compressibility of rubber (finite bulk modulus) in the finite element formulation is demonstrated in the study of apparent Young’s modulus of bonded thin rubber units. The combined shear-bending deformation that commonly exists in rubber mounting systems is also analyzed and discussed.


1996 ◽  
Vol 63 (2) ◽  
pp. 419-423 ◽  
Author(s):  
R. F. Yu¨kseler

An approach for the derivation of two-dimensional strain energy density functions of compressible, rubber-like shells of revolution which undergo axisymmetric, arbitrarily large rotations and strains including transverse normal and transverse shear strains is proposed. The proposed method is applied to polyurethane rubber, circular plates and numerical results for the flexural buckling of polyurethane circular plates are obtained for an illustrative comparison with a method presented previously.


Sign in / Sign up

Export Citation Format

Share Document