Phase-field simulation of crack growth in cortical bone microstructure : parameter identification and comparison against experiments
(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)
- author
- Carlsson, Jenny
LU
; Karlsson, Olivia ; Isaksson, Hanna LU
and Gustafsson, Anna LU
- organization
- publishing date
- 2025
- 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}}, }