An open source software tool to assign the material properties of bone for ABAQUS finite element simulations

Research output: Contribution to journalArticle

3 Citations (Scopus)
95 Downloads (Pure)

Abstract

A new software tool to assign the material properties of bone to an ABAQUS finite element mesh was created and compared with Bonemat, a similar tool originally designed to work with Ansys finite element models. Our software tool (py_bonemat_abaqus) was written in Python, which is the chosen scripting language for ABAQUS. The purpose of this study was to compare the software packages in terms of the material assignment calculation and processing speed. Three element types were compared (linear hexahedral (C3D8), linear tetrahedral (C3D4) and quadratic tetrahedral elements (C3D10)), both individually and as part of a mesh.

Comparisons were made using a CT scan of a hemi-pelvis as a test case. A small difference, of -0.05 kPa on average, was found between Bonemat version 3.1 (the current version) and our python package. Errors were found in the previous release of Bonemat (version 3.0 downloaded from www.biomedtown.org) during calculation of the quadratic tetrahedron Jacobian, and conversion of the apparent density to modulus when integrating over the Young’s modulus field. These issues caused up to 2 GPa error in the modulus assignment. For these reasons, we recommend users upgrade to the most recent release of Bonemat.

Processing speeds were assessed for the three different element types. Our Python package took significantly longer (110 s on average) to perform the calculations compared with the Bonemat software (10 s). Nevertheless, the workflow advantages of the package and added functionality makes ‘py_bonemat_abaqus’ a useful tool for ABAQUS users.
Original languageEnglish
Pages (from-to)3116-3121
JournalJournal of Biomechanics
Volume49
Issue number13
Early online date8 Aug 2016
DOIs
Publication statusPublished - Sep 2016

Fingerprint

ABAQUS
Boidae
Materials properties
Bone
Software
Bone and Bones
Computerized tomography
Processing
Software packages
Workflow
Elastic Modulus
Elastic moduli
Pelvis
Language
Open source software

Keywords

  • Material properties
  • Bone
  • Finite Element
  • Python

Cite this

@article{2166dd7ea0a94e41a6fe489c8586651f,
title = "An open source software tool to assign the material properties of bone for ABAQUS finite element simulations",
abstract = "A new software tool to assign the material properties of bone to an ABAQUS finite element mesh was created and compared with Bonemat, a similar tool originally designed to work with Ansys finite element models. Our software tool (py_bonemat_abaqus) was written in Python, which is the chosen scripting language for ABAQUS. The purpose of this study was to compare the software packages in terms of the material assignment calculation and processing speed. Three element types were compared (linear hexahedral (C3D8), linear tetrahedral (C3D4) and quadratic tetrahedral elements (C3D10)), both individually and as part of a mesh. Comparisons were made using a CT scan of a hemi-pelvis as a test case. A small difference, of -0.05 kPa on average, was found between Bonemat version 3.1 (the current version) and our python package. Errors were found in the previous release of Bonemat (version 3.0 downloaded from www.biomedtown.org) during calculation of the quadratic tetrahedron Jacobian, and conversion of the apparent density to modulus when integrating over the Young’s modulus field. These issues caused up to 2 GPa error in the modulus assignment. For these reasons, we recommend users upgrade to the most recent release of Bonemat. Processing speeds were assessed for the three different element types. Our Python package took significantly longer (110 s on average) to perform the calculations compared with the Bonemat software (10 s). Nevertheless, the workflow advantages of the package and added functionality makes ‘py_bonemat_abaqus’ a useful tool for ABAQUS users.",
keywords = "Material properties, Bone, Finite Element, Python",
author = "Elise Pegg and Harinderjit Gill",
year = "2016",
month = "9",
doi = "10.1016/j.jbiomech.2016.07.037",
language = "English",
volume = "49",
pages = "3116--3121",
journal = "Journal of Biomechanics",
issn = "0021-9290",
publisher = "Elsevier",
number = "13",

}

TY - JOUR

T1 - An open source software tool to assign the material properties of bone for ABAQUS finite element simulations

AU - Pegg, Elise

AU - Gill, Harinderjit

PY - 2016/9

Y1 - 2016/9

N2 - A new software tool to assign the material properties of bone to an ABAQUS finite element mesh was created and compared with Bonemat, a similar tool originally designed to work with Ansys finite element models. Our software tool (py_bonemat_abaqus) was written in Python, which is the chosen scripting language for ABAQUS. The purpose of this study was to compare the software packages in terms of the material assignment calculation and processing speed. Three element types were compared (linear hexahedral (C3D8), linear tetrahedral (C3D4) and quadratic tetrahedral elements (C3D10)), both individually and as part of a mesh. Comparisons were made using a CT scan of a hemi-pelvis as a test case. A small difference, of -0.05 kPa on average, was found between Bonemat version 3.1 (the current version) and our python package. Errors were found in the previous release of Bonemat (version 3.0 downloaded from www.biomedtown.org) during calculation of the quadratic tetrahedron Jacobian, and conversion of the apparent density to modulus when integrating over the Young’s modulus field. These issues caused up to 2 GPa error in the modulus assignment. For these reasons, we recommend users upgrade to the most recent release of Bonemat. Processing speeds were assessed for the three different element types. Our Python package took significantly longer (110 s on average) to perform the calculations compared with the Bonemat software (10 s). Nevertheless, the workflow advantages of the package and added functionality makes ‘py_bonemat_abaqus’ a useful tool for ABAQUS users.

AB - A new software tool to assign the material properties of bone to an ABAQUS finite element mesh was created and compared with Bonemat, a similar tool originally designed to work with Ansys finite element models. Our software tool (py_bonemat_abaqus) was written in Python, which is the chosen scripting language for ABAQUS. The purpose of this study was to compare the software packages in terms of the material assignment calculation and processing speed. Three element types were compared (linear hexahedral (C3D8), linear tetrahedral (C3D4) and quadratic tetrahedral elements (C3D10)), both individually and as part of a mesh. Comparisons were made using a CT scan of a hemi-pelvis as a test case. A small difference, of -0.05 kPa on average, was found between Bonemat version 3.1 (the current version) and our python package. Errors were found in the previous release of Bonemat (version 3.0 downloaded from www.biomedtown.org) during calculation of the quadratic tetrahedron Jacobian, and conversion of the apparent density to modulus when integrating over the Young’s modulus field. These issues caused up to 2 GPa error in the modulus assignment. For these reasons, we recommend users upgrade to the most recent release of Bonemat. Processing speeds were assessed for the three different element types. Our Python package took significantly longer (110 s on average) to perform the calculations compared with the Bonemat software (10 s). Nevertheless, the workflow advantages of the package and added functionality makes ‘py_bonemat_abaqus’ a useful tool for ABAQUS users.

KW - Material properties

KW - Bone

KW - Finite Element

KW - Python

UR - https://doi.org/10.1016/j.jbiomech.2016.07.037

U2 - 10.1016/j.jbiomech.2016.07.037

DO - 10.1016/j.jbiomech.2016.07.037

M3 - Article

VL - 49

SP - 3116

EP - 3121

JO - Journal of Biomechanics

JF - Journal of Biomechanics

SN - 0021-9290

IS - 13

ER -