A smoothed finite element method for analysis of anisotropic large deformation of passive rabbit ventricles in diastole

Smoothed finite element methods in simulation of active contraction of myocardial tissue samples

In modelling and simulation of cardiac mechanics, tetrahedral meshes are often used due to the easy availability of efficient meshing algorithms. This is beneficial in particular when complex geometries such as cardiac structures are considered. The gold standard in simulating the cardiac cycle is to solve the mechanical balance equations with the finite element method (FEM). However, using linear shape functions in the FEM in combination with nearly-incompressible material models is known to produce overly stiff approximations, whereas higher order elements are computationally more expensive. To overcome these problems, smoothed finite element methods (S-FEMs) have been proposed by Liu and co-workers. So far, S-FEMs in 3D have been utilised only in simulations of passive mechanics. In the present work, different S-FEMs are for the first time used for simulation of an active cardiac contraction on three-dimensional myocardial tissue samples. Further, node-based S-FEM (NS-FEM), face-based S-FEM (FS-FEM) and selective FS/NS-FEM are for the first time implemented as user subroutine in the commercial software Abaqus. Our results confirm that all S-FEMs perform softer than linear FEM and volumetric locking is reduced. The FS/NS-FEM produces solutions with the relative error in maximum displacement and rotation being less than 5% with respect to the reference solution obtained by the quadratic FEM for all considered mesh sizes, although linear shape functions are used. We therefore conclude that in particular FS/NS-FEM is an efficient and accurate numerical method in the simulation of an active cardiac muscle contraction.

Towards the simulation of active cardiac mechanics using a smoothed finite element method

In the last decades, various computational models have been developed to simulate cardiac electromechanics. The most common numerical tool is the finite element method (FEM). However, this method crucially depends on the mesh quality. For complex geometries such as cardiac structures, it is convenient to use tetrahedral discretisations which can be generated automatically. On the other hand, such automatic meshing with tetrahedrons together with large deformations often lead to elements distortion and volumetric locking. To overcome these difficulties, different smoothed finite element methods (S-FEMs) have been proposed in the recent years. They are known to be volumetric locking free, less sensitive to mesh distortion and so far have been used e.g. in simulation of passive cardiac mechanics. In this work, we extend for the first time node-based S-FEM (NS-FEM) towards active cardiac mechanics. Firstly, the sensitivity to mesh distortion is tested and compared to that of FEM. Secondly, an active contraction in circumferentially aligned fibre direction is modelled in the healthy and the infarcted case. We show, that the proposed method is more robust with respect to mesh distortion and computationally more efficient than standard FEM. Being furthermore free of volumetric locking problems makes S-FEM a promising alternative in modelling of active cardiac mechanics, respectively electromechanics.