Skip to main content

Lund University Publications

LUND UNIVERSITY LIBRARIES

Phase field models of interface failure for bone application - evaluation of open-source implementations

Gustafsson, Anna LU orcid and Isaksson, Hanna LU orcid (2022) In Theoretical and Applied Fracture Mechanics 121.
Abstract

Bone is a composite material where tubular osteon structures surrounded by thin interfaces reinforce the tissue at the microscale. These interfaces provide alternative weak paths for crack growth and give rise to potent extrinsic toughening mechanisms. Numerical models can help to gain a deeper understanding of fracture resistance in bone, however, capturing the interface mechanics is key. For this purpose, the phase field method for fracture is appealing due to its capability of capturing complex cracking phenomena. Still, it is far from well-established in the field of biomechanics. In this study, we evaluated a selection of recent open-source implementations of the phase field method to find the best approach for simulating crack... (More)

Bone is a composite material where tubular osteon structures surrounded by thin interfaces reinforce the tissue at the microscale. These interfaces provide alternative weak paths for crack growth and give rise to potent extrinsic toughening mechanisms. Numerical models can help to gain a deeper understanding of fracture resistance in bone, however, capturing the interface mechanics is key. For this purpose, the phase field method for fracture is appealing due to its capability of capturing complex cracking phenomena. Still, it is far from well-established in the field of biomechanics. In this study, we evaluated a selection of recent open-source implementations of the phase field method to find the best approach for simulating crack growth in bone tissue. We also proposed a new method for correcting the interface fracture toughness in a phase field model. The selected implementations were compared using benchmark tests, including single edge notched tensile specimens with and without interfaces, as well as models representing typical bone geometries. We found the quasi-Newton monolithic solvers to be superior both in terms of computational speed and robustness compared to the evaluated staggered solvers. We also showed that correcting the fracture toughness of the interface is key for simulating crack patterns that are consistent with analytical predictions from linear elastic fracture mechanics.

(Less)
Please use this url to cite or link to this publication:
author
and
organization
publishing date
type
Contribution to journal
publication status
published
subject
keywords
Biomechanics, Cement line, Crack propagation, FEA, Fracture, Osteon
in
Theoretical and Applied Fracture Mechanics
volume
121
article number
103432
publisher
Elsevier
external identifiers
  • scopus:85131797040
ISSN
0167-8442
DOI
10.1016/j.tafmec.2022.103432
project
Fracture mechanics of bone tissue – impact of age and disease
language
English
LU publication?
yes
id
24452896-7966-467d-98cf-d44dd287cba4
date added to LUP
2022-08-26 14:48:38
date last changed
2023-10-23 14:00:55
@article{24452896-7966-467d-98cf-d44dd287cba4,
  abstract     = {{<p>Bone is a composite material where tubular osteon structures surrounded by thin interfaces reinforce the tissue at the microscale. These interfaces provide alternative weak paths for crack growth and give rise to potent extrinsic toughening mechanisms. Numerical models can help to gain a deeper understanding of fracture resistance in bone, however, capturing the interface mechanics is key. For this purpose, the phase field method for fracture is appealing due to its capability of capturing complex cracking phenomena. Still, it is far from well-established in the field of biomechanics. In this study, we evaluated a selection of recent open-source implementations of the phase field method to find the best approach for simulating crack growth in bone tissue. We also proposed a new method for correcting the interface fracture toughness in a phase field model. The selected implementations were compared using benchmark tests, including single edge notched tensile specimens with and without interfaces, as well as models representing typical bone geometries. We found the quasi-Newton monolithic solvers to be superior both in terms of computational speed and robustness compared to the evaluated staggered solvers. We also showed that correcting the fracture toughness of the interface is key for simulating crack patterns that are consistent with analytical predictions from linear elastic fracture mechanics.</p>}},
  author       = {{Gustafsson, Anna and Isaksson, Hanna}},
  issn         = {{0167-8442}},
  keywords     = {{Biomechanics; Cement line; Crack propagation; FEA; Fracture; Osteon}},
  language     = {{eng}},
  publisher    = {{Elsevier}},
  series       = {{Theoretical and Applied Fracture Mechanics}},
  title        = {{Phase field models of interface failure for bone application - evaluation of open-source implementations}},
  url          = {{http://dx.doi.org/10.1016/j.tafmec.2022.103432}},
  doi          = {{10.1016/j.tafmec.2022.103432}},
  volume       = {{121}},
  year         = {{2022}},
}