An orthotropic viscoelastic model for the passive myocardium: continuum basis and numerical treatment

Research output: Contribution to journalArticleResearchpeer-review

Abstract

This study deals with the viscoelastic constitutive modeling and the respective computational analysis of the human passive myocardium. We start by recapitulating the locally orthotropic inner structure of the human myocardial tissue and model the mechanical response through invariants and structure tensors associated with three orthonormal basis vectors. In accordance with recent experimental findings the ventricular myocardial tissue is assumed to be incompressible, thick-walled, orthotropic and viscoelastic. In particular, one spring element coupled with Maxwell elements in parallel endows the model with viscoelastic features such that four dashpots describe the viscous response due to matrix, fiber, sheet and fiber-sheet fragments. In order to alleviate the numerical obstacles, the strictly incompressible model is altered by decomposing the free-energy function into volumetric-isochoric elastic and isochoric-viscoelastic parts along with the multiplicative split of the deformation gradient which enables the three-field mixed finite element method. The crucial aspect of the viscoelastic formulation is linked to the rate equations of the viscous overstresses resulting from a 3-D analogy of a generalized 1-D Maxwell model. We provide algorithmic updates for second Piola-Kirchhoff stress and elasticity tensors. In the sequel, we address some numerical aspects of the constitutive model by applying it to elastic, cyclic and relaxation test data obtained from biaxial extension and triaxial shear tests whereby we assess the fitting capacity of the model. With the tissue parameters identified, we conduct (elastic and viscoelastic) finite element simulations for an ellipsoidal geometry retrieved from a human specimen.

Original languageEnglish
Pages (from-to)1647-64
Number of pages18
JournalComputer Methods in Biomechanics and Biomedical Engineering
Volume19
Issue number15
DOIs
Publication statusPublished - Nov 2016

Fingerprint

Myocardium
Tissue
Elasticity
Tensors
Fibers
Constitutive models
Free energy
Finite element method
Geometry

Keywords

  • Cardiac mechanics
  • myocardium
  • viscoelasticity
  • Finite element method
  • Orthotropy
  • collagen fiber dispersion

Cite this

An orthotropic viscoelastic model for the passive myocardium : continuum basis and numerical treatment. / Gültekin, Osman; Sommer, Gerhard; Holzapfel, Gerhard.

In: Computer Methods in Biomechanics and Biomedical Engineering, Vol. 19, No. 15, 11.2016, p. 1647-64.

Research output: Contribution to journalArticleResearchpeer-review

@article{3cd18a5a83174813a3a747ebc90e423f,
title = "An orthotropic viscoelastic model for the passive myocardium: continuum basis and numerical treatment",
abstract = "This study deals with the viscoelastic constitutive modeling and the respective computational analysis of the human passive myocardium. We start by recapitulating the locally orthotropic inner structure of the human myocardial tissue and model the mechanical response through invariants and structure tensors associated with three orthonormal basis vectors. In accordance with recent experimental findings the ventricular myocardial tissue is assumed to be incompressible, thick-walled, orthotropic and viscoelastic. In particular, one spring element coupled with Maxwell elements in parallel endows the model with viscoelastic features such that four dashpots describe the viscous response due to matrix, fiber, sheet and fiber-sheet fragments. In order to alleviate the numerical obstacles, the strictly incompressible model is altered by decomposing the free-energy function into volumetric-isochoric elastic and isochoric-viscoelastic parts along with the multiplicative split of the deformation gradient which enables the three-field mixed finite element method. The crucial aspect of the viscoelastic formulation is linked to the rate equations of the viscous overstresses resulting from a 3-D analogy of a generalized 1-D Maxwell model. We provide algorithmic updates for second Piola-Kirchhoff stress and elasticity tensors. In the sequel, we address some numerical aspects of the constitutive model by applying it to elastic, cyclic and relaxation test data obtained from biaxial extension and triaxial shear tests whereby we assess the fitting capacity of the model. With the tissue parameters identified, we conduct (elastic and viscoelastic) finite element simulations for an ellipsoidal geometry retrieved from a human specimen.",
keywords = "Cardiac mechanics, myocardium, viscoelasticity, Finite element method, Orthotropy, collagen fiber dispersion",
author = "Osman G{\"u}ltekin and Gerhard Sommer and Gerhard Holzapfel",
year = "2016",
month = "11",
doi = "10.1080/10255842.2016.1176155",
language = "English",
volume = "19",
pages = "1647--64",
journal = "Computer Methods in Biomechanics and Biomedical Engineering",
issn = "1025-5842",
publisher = "Taylor & Francis Group",
number = "15",

}

TY - JOUR

T1 - An orthotropic viscoelastic model for the passive myocardium

T2 - continuum basis and numerical treatment

AU - Gültekin, Osman

AU - Sommer, Gerhard

AU - Holzapfel, Gerhard

PY - 2016/11

Y1 - 2016/11

N2 - This study deals with the viscoelastic constitutive modeling and the respective computational analysis of the human passive myocardium. We start by recapitulating the locally orthotropic inner structure of the human myocardial tissue and model the mechanical response through invariants and structure tensors associated with three orthonormal basis vectors. In accordance with recent experimental findings the ventricular myocardial tissue is assumed to be incompressible, thick-walled, orthotropic and viscoelastic. In particular, one spring element coupled with Maxwell elements in parallel endows the model with viscoelastic features such that four dashpots describe the viscous response due to matrix, fiber, sheet and fiber-sheet fragments. In order to alleviate the numerical obstacles, the strictly incompressible model is altered by decomposing the free-energy function into volumetric-isochoric elastic and isochoric-viscoelastic parts along with the multiplicative split of the deformation gradient which enables the three-field mixed finite element method. The crucial aspect of the viscoelastic formulation is linked to the rate equations of the viscous overstresses resulting from a 3-D analogy of a generalized 1-D Maxwell model. We provide algorithmic updates for second Piola-Kirchhoff stress and elasticity tensors. In the sequel, we address some numerical aspects of the constitutive model by applying it to elastic, cyclic and relaxation test data obtained from biaxial extension and triaxial shear tests whereby we assess the fitting capacity of the model. With the tissue parameters identified, we conduct (elastic and viscoelastic) finite element simulations for an ellipsoidal geometry retrieved from a human specimen.

AB - This study deals with the viscoelastic constitutive modeling and the respective computational analysis of the human passive myocardium. We start by recapitulating the locally orthotropic inner structure of the human myocardial tissue and model the mechanical response through invariants and structure tensors associated with three orthonormal basis vectors. In accordance with recent experimental findings the ventricular myocardial tissue is assumed to be incompressible, thick-walled, orthotropic and viscoelastic. In particular, one spring element coupled with Maxwell elements in parallel endows the model with viscoelastic features such that four dashpots describe the viscous response due to matrix, fiber, sheet and fiber-sheet fragments. In order to alleviate the numerical obstacles, the strictly incompressible model is altered by decomposing the free-energy function into volumetric-isochoric elastic and isochoric-viscoelastic parts along with the multiplicative split of the deformation gradient which enables the three-field mixed finite element method. The crucial aspect of the viscoelastic formulation is linked to the rate equations of the viscous overstresses resulting from a 3-D analogy of a generalized 1-D Maxwell model. We provide algorithmic updates for second Piola-Kirchhoff stress and elasticity tensors. In the sequel, we address some numerical aspects of the constitutive model by applying it to elastic, cyclic and relaxation test data obtained from biaxial extension and triaxial shear tests whereby we assess the fitting capacity of the model. With the tissue parameters identified, we conduct (elastic and viscoelastic) finite element simulations for an ellipsoidal geometry retrieved from a human specimen.

KW - Cardiac mechanics

KW - myocardium

KW - viscoelasticity

KW - Finite element method

KW - Orthotropy

KW - collagen fiber dispersion

U2 - 10.1080/10255842.2016.1176155

DO - 10.1080/10255842.2016.1176155

M3 - Article

VL - 19

SP - 1647

EP - 1664

JO - Computer Methods in Biomechanics and Biomedical Engineering

JF - Computer Methods in Biomechanics and Biomedical Engineering

SN - 1025-5842

IS - 15

ER -