A comparison of numerical methods used for finite element modelling of soft tissue deformation

Handle URI:
http://hdl.handle.net/10754/597237
Title:
A comparison of numerical methods used for finite element modelling of soft tissue deformation
Authors:
Pathmanathan, P; Gavaghan, D; Whiteley, J
Abstract:
Soft tissue deformation is often modelled using incompressible non-linear elasticity, with solutions computed using the finite element method. There are a range of options available when using the finite element method, in particular the polynomial degree of the basis functions used for interpolating position and pressure, and the type of element making up the mesh. The effect of these choices on the accuracy of the computed solution is investigated, using a selection of model problems motivated by typical deformations seen in soft tissue modelling. Model problems are set up with discontinuous material properties (as is the case for the breast), steeply changing gradients in the body force (as found in contracting cardiac tissue), and discontinuous first derivatives in the solution at the boundary, caused by a discontinuous applied force (as in the breast during mammography). It was found that the choice of pressure basis functions is vital in the presence of a material interface, higher-order schemes do not perform as well as may be expected when there are sharp gradients, and in general it is important to take the expected regularity of the solution into account when choosing a numerical scheme. © IMechE 2009.
Citation:
Pathmanathan P, Gavaghan D, Whiteley J (2009) A comparison of numerical methods used for finite element modelling of soft tissue deformation. The Journal of Strain Analysis for Engineering Design 44: 391–406. Available: http://dx.doi.org/10.1243/03093247JSA482.
Publisher:
SAGE Publications
Journal:
The Journal of Strain Analysis for Engineering Design
KAUST Grant Number:
KUK-CI-013-04
Issue Date:
1-May-2009
DOI:
10.1243/03093247JSA482
Type:
Article
ISSN:
0309-3247; 2041-3130
Sponsors:
PP is supported through the EPSRC-funded OXMOS project New frontiers in the mathematics of solids (Grant EP/D048400/1). JPW is partially funded by the EPSRC under Grant EP/D503035/1. This publication is based on work supported by Award KUK-CI-013-04, made by King Abdullah University of Science and Technology (KAUST).
Appears in Collections:
Publications Acknowledging KAUST Support

Full metadata record

DC FieldValue Language
dc.contributor.authorPathmanathan, Pen
dc.contributor.authorGavaghan, Den
dc.contributor.authorWhiteley, Jen
dc.date.accessioned2016-02-25T12:28:40Zen
dc.date.available2016-02-25T12:28:40Zen
dc.date.issued2009-05-01en
dc.identifier.citationPathmanathan P, Gavaghan D, Whiteley J (2009) A comparison of numerical methods used for finite element modelling of soft tissue deformation. The Journal of Strain Analysis for Engineering Design 44: 391–406. Available: http://dx.doi.org/10.1243/03093247JSA482.en
dc.identifier.issn0309-3247en
dc.identifier.issn2041-3130en
dc.identifier.doi10.1243/03093247JSA482en
dc.identifier.urihttp://hdl.handle.net/10754/597237en
dc.description.abstractSoft tissue deformation is often modelled using incompressible non-linear elasticity, with solutions computed using the finite element method. There are a range of options available when using the finite element method, in particular the polynomial degree of the basis functions used for interpolating position and pressure, and the type of element making up the mesh. The effect of these choices on the accuracy of the computed solution is investigated, using a selection of model problems motivated by typical deformations seen in soft tissue modelling. Model problems are set up with discontinuous material properties (as is the case for the breast), steeply changing gradients in the body force (as found in contracting cardiac tissue), and discontinuous first derivatives in the solution at the boundary, caused by a discontinuous applied force (as in the breast during mammography). It was found that the choice of pressure basis functions is vital in the presence of a material interface, higher-order schemes do not perform as well as may be expected when there are sharp gradients, and in general it is important to take the expected regularity of the solution into account when choosing a numerical scheme. © IMechE 2009.en
dc.description.sponsorshipPP is supported through the EPSRC-funded OXMOS project New frontiers in the mathematics of solids (Grant EP/D048400/1). JPW is partially funded by the EPSRC under Grant EP/D503035/1. This publication is based on work supported by Award KUK-CI-013-04, made by King Abdullah University of Science and Technology (KAUST).en
dc.publisherSAGE Publicationsen
dc.subjectBasis functionsen
dc.subjectBreasten
dc.subjectFinite element methoden
dc.subjectHearten
dc.subjectNon-linear elasticityen
dc.titleA comparison of numerical methods used for finite element modelling of soft tissue deformationen
dc.typeArticleen
dc.identifier.journalThe Journal of Strain Analysis for Engineering Designen
dc.contributor.institutionUniversity of Oxford, Oxford, United Kingdomen
kaust.grant.numberKUK-CI-013-04en
All Items in KAUST are protected by copyright, with all rights reserved, unless otherwise indicated.