A multiscale predictor/corrector scheme for efficient elastoplastic voxel finite element analysis, with application to CT-based bone strength prediction

Lam H. Nguyen, Dominik Schillinger

Research output: Research - peer-reviewArticle

Abstract

Voxel finite elements combined with plasticity have been shown to accurately predict the evolution of bone failure, but involve a prohibitive computational cost when applied to high-resolution CT scans of a complete bone. We present a simple multiscale predictor/corrector scheme that uses elasticity and the finite cell method on a coarse-scale mesh, complemented by plasticity and fine-scale voxel finite elements in regions where failure occurs. The core components of our method are top-down displacement and bottom-up stress projectors for the exchange of information between coarse and fine scales. Our choice of projectors eliminates communication of fine-scale voxel elements beyond boundaries of coarse-scale cells, which enables the solution in terms of a series of small uncoupled systems at a fraction of the computing power and memory required by the fully coupled fine-scale system. At the same time, we illustrate that the multiscale approach yields the same accuracy as the full-resolution voxel finite element method, if we appropriately balance the approximation power of coarse-scale and fine-scale meshes. We demonstrate the advantages of our method for the load capacity analysis of a patient-specific vertebra.

LanguageEnglish (US)
Pages598-628
Number of pages31
JournalComputer Methods in Applied Mechanics and Engineering
Volume330
DOIs
StatePublished - Mar 1 2018

Fingerprint

Plasticity
Bone
Finite element method
bones
predictions
Computerized tomography
Elasticity
Data storage equipment
Communication
Costs
projectors
plastic properties
mesh
cells
vertebrae
finite element method
elastic properties
communication
costs
high resolution

Keywords

  • Bone strength prediction
  • Finite cell method
  • Multiscale predictor/corrector scheme
  • Parallel computing
  • Plasticity
  • Voxel finite element method

Cite this

@article{a81ca3306ad5469e9489cec4a1b17484,
title = "A multiscale predictor/corrector scheme for efficient elastoplastic voxel finite element analysis, with application to CT-based bone strength prediction",
abstract = "Voxel finite elements combined with plasticity have been shown to accurately predict the evolution of bone failure, but involve a prohibitive computational cost when applied to high-resolution CT scans of a complete bone. We present a simple multiscale predictor/corrector scheme that uses elasticity and the finite cell method on a coarse-scale mesh, complemented by plasticity and fine-scale voxel finite elements in regions where failure occurs. The core components of our method are top-down displacement and bottom-up stress projectors for the exchange of information between coarse and fine scales. Our choice of projectors eliminates communication of fine-scale voxel elements beyond boundaries of coarse-scale cells, which enables the solution in terms of a series of small uncoupled systems at a fraction of the computing power and memory required by the fully coupled fine-scale system. At the same time, we illustrate that the multiscale approach yields the same accuracy as the full-resolution voxel finite element method, if we appropriately balance the approximation power of coarse-scale and fine-scale meshes. We demonstrate the advantages of our method for the load capacity analysis of a patient-specific vertebra.",
keywords = "Bone strength prediction, Finite cell method, Multiscale predictor/corrector scheme, Parallel computing, Plasticity, Voxel finite element method",
author = "Nguyen, {Lam H.} and Dominik Schillinger",
year = "2018",
month = "3",
doi = "10.1016/j.cma.2017.11.014",
volume = "330",
pages = "598--628",
journal = "Computer Methods in Applied Mechanics and Engineering",
issn = "0374-2830",
publisher = "Elsevier",

}

TY - JOUR

T1 - A multiscale predictor/corrector scheme for efficient elastoplastic voxel finite element analysis, with application to CT-based bone strength prediction

AU - Nguyen,Lam H.

AU - Schillinger,Dominik

PY - 2018/3/1

Y1 - 2018/3/1

N2 - Voxel finite elements combined with plasticity have been shown to accurately predict the evolution of bone failure, but involve a prohibitive computational cost when applied to high-resolution CT scans of a complete bone. We present a simple multiscale predictor/corrector scheme that uses elasticity and the finite cell method on a coarse-scale mesh, complemented by plasticity and fine-scale voxel finite elements in regions where failure occurs. The core components of our method are top-down displacement and bottom-up stress projectors for the exchange of information between coarse and fine scales. Our choice of projectors eliminates communication of fine-scale voxel elements beyond boundaries of coarse-scale cells, which enables the solution in terms of a series of small uncoupled systems at a fraction of the computing power and memory required by the fully coupled fine-scale system. At the same time, we illustrate that the multiscale approach yields the same accuracy as the full-resolution voxel finite element method, if we appropriately balance the approximation power of coarse-scale and fine-scale meshes. We demonstrate the advantages of our method for the load capacity analysis of a patient-specific vertebra.

AB - Voxel finite elements combined with plasticity have been shown to accurately predict the evolution of bone failure, but involve a prohibitive computational cost when applied to high-resolution CT scans of a complete bone. We present a simple multiscale predictor/corrector scheme that uses elasticity and the finite cell method on a coarse-scale mesh, complemented by plasticity and fine-scale voxel finite elements in regions where failure occurs. The core components of our method are top-down displacement and bottom-up stress projectors for the exchange of information between coarse and fine scales. Our choice of projectors eliminates communication of fine-scale voxel elements beyond boundaries of coarse-scale cells, which enables the solution in terms of a series of small uncoupled systems at a fraction of the computing power and memory required by the fully coupled fine-scale system. At the same time, we illustrate that the multiscale approach yields the same accuracy as the full-resolution voxel finite element method, if we appropriately balance the approximation power of coarse-scale and fine-scale meshes. We demonstrate the advantages of our method for the load capacity analysis of a patient-specific vertebra.

KW - Bone strength prediction

KW - Finite cell method

KW - Multiscale predictor/corrector scheme

KW - Parallel computing

KW - Plasticity

KW - Voxel finite element method

UR - http://www.scopus.com/inward/record.url?scp=85036468250&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85036468250&partnerID=8YFLogxK

U2 - 10.1016/j.cma.2017.11.014

DO - 10.1016/j.cma.2017.11.014

M3 - Article

VL - 330

SP - 598

EP - 628

JO - Computer Methods in Applied Mechanics and Engineering

T2 - Computer Methods in Applied Mechanics and Engineering

JF - Computer Methods in Applied Mechanics and Engineering

SN - 0374-2830

ER -