Skip to main content

Lund University Publications

LUND UNIVERSITY LIBRARIES

Phase-field simulation of crack growth in cortical bone microstructure : parameter identification and comparison against experiments

Carlsson, Jenny LU orcid ; Karlsson, Olivia ; Isaksson, Hanna LU orcid and Gustafsson, Anna LU orcid (2025) In Biomechanics and Modeling in Mechanobiology
Abstract

Computational models are commonly used to investigate how the cortical bone microstructure affects fracture resistance; recently, phase-field models have been introduced for this purpose. However, experimentally measured material parameters for the microstructural tissues are lacking. Moreover, as no validation studies have been published, it remains unclear to what extent classical phase-field methods, assuming linear-elastic, brittle fracture, accurately represent bone. In this study, we address both these shortcomings by first applying a design-of-experiments methodology to calibrate a set of material parameters for a two-dimensional phase-field finite element model of bovine osteonal microstructure. This was achieved by comparing... (More)

Computational models are commonly used to investigate how the cortical bone microstructure affects fracture resistance; recently, phase-field models have been introduced for this purpose. However, experimentally measured material parameters for the microstructural tissues are lacking. Moreover, as no validation studies have been published, it remains unclear to what extent classical phase-field methods, assuming linear-elastic, brittle fracture, accurately represent bone. In this study, we address both these shortcomings by first applying a design-of-experiments methodology to calibrate a set of material parameters for a two-dimensional phase-field finite element model of bovine osteonal microstructure. This was achieved by comparing the outcomes from simulation to data from single-edge notched bending experiments on bovine osteonal bone and subsequent imaging of the crack path. Second, we used these parameters in new bone geometries to evaluate the parameters and the predictive performance of the model. Reasonable agreement was achieved between prediction and experiments in terms of peak load, crack initiation toughness and crack path. However, the model is unable to capture the experimentally observed gradual evolution of damage, leading to a nonlinear force response before the onset of visible crack extension. Nor does it capture the similarly observed increase in toughness with increasing crack length. These limitations are inherent to all classical phase-field methods since they originate from theories of brittle fracture, and alternative formulations are discussed. This is the first study attempting to validate classical phase-field methods in simulation of cortical bone fracture, and it highlights both potential and limitations to be addressed in future work.

(Less)
Please use this url to cite or link to this publication:
author
; ; and
organization
publishing date
type
Contribution to journal
publication status
epub
subject
keywords
Cement line, Crack propagation, Fracture toughness, Osteon, Phase-field method
in
Biomechanics and Modeling in Mechanobiology
article number
106034
pages
15 pages
publisher
Springer
external identifiers
  • scopus:85218928641
  • pmid:40025294
ISSN
1617-7959
DOI
10.1007/s10237-025-01929-8
project
Fracture mechanics of bone tissue – impact of age and disease
language
English
LU publication?
yes
additional info
Publisher Copyright: © The Author(s) 2025.
id
cabf5976-8caa-416b-8041-3ba694f80257
date added to LUP
2025-03-12 11:38:58
date last changed
2025-06-19 03:04:48
@article{cabf5976-8caa-416b-8041-3ba694f80257,
  abstract     = {{<p>Computational models are commonly used to investigate how the cortical bone microstructure affects fracture resistance; recently, phase-field models have been introduced for this purpose. However, experimentally measured material parameters for the microstructural tissues are lacking. Moreover, as no validation studies have been published, it remains unclear to what extent classical phase-field methods, assuming linear-elastic, brittle fracture, accurately represent bone. In this study, we address both these shortcomings by first applying a design-of-experiments methodology to calibrate a set of material parameters for a two-dimensional phase-field finite element model of bovine osteonal microstructure. This was achieved by comparing the outcomes from simulation to data from single-edge notched bending experiments on bovine osteonal bone and subsequent imaging of the crack path. Second, we used these parameters in new bone geometries to evaluate the parameters and the predictive performance of the model. Reasonable agreement was achieved between prediction and experiments in terms of peak load, crack initiation toughness and crack path. However, the model is unable to capture the experimentally observed gradual evolution of damage, leading to a nonlinear force response before the onset of visible crack extension. Nor does it capture the similarly observed increase in toughness with increasing crack length. These limitations are inherent to all classical phase-field methods since they originate from theories of brittle fracture, and alternative formulations are discussed. This is the first study attempting to validate classical phase-field methods in simulation of cortical bone fracture, and it highlights both potential and limitations to be addressed in future work.</p>}},
  author       = {{Carlsson, Jenny and Karlsson, Olivia and Isaksson, Hanna and Gustafsson, Anna}},
  issn         = {{1617-7959}},
  keywords     = {{Cement line; Crack propagation; Fracture toughness; Osteon; Phase-field method}},
  language     = {{eng}},
  publisher    = {{Springer}},
  series       = {{Biomechanics and Modeling in Mechanobiology}},
  title        = {{Phase-field simulation of crack growth in cortical bone microstructure : parameter identification and comparison against experiments}},
  url          = {{http://dx.doi.org/10.1007/s10237-025-01929-8}},
  doi          = {{10.1007/s10237-025-01929-8}},
  year         = {{2025}},
}