z-logo
Premium
A novel algorithm to resolve lack of convergence and checkerboard instability in bone adaptation simulations using non‐local averaging
Author(s) -
CalvoGallego José Luis,
Pivonka Peter,
GarcíaAznar José Manuel,
MartínezReina Javier
Publication year - 2021
Publication title -
international journal for numerical methods in biomedical engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.741
H-Index - 63
eISSN - 2040-7947
pISSN - 2040-7939
DOI - 10.1002/cnm.3419
Subject(s) - checkerboard , instability , finite element method , smoothing , algorithm , torsion (gastropod) , computer science , mathematics , mechanics , mathematical optimization , topology (electrical circuits) , structural engineering , geometry , physics , engineering , medicine , surgery , combinatorics , computer vision
Checkerboard is a typical instability in finite element (FE) simulations of bone adaptation and topology optimization in general. It consists in a patchwork pattern with elements of alternating stiffness, producing lack of convergence and instabilities in the predicted bone density. Averaging techniques have been proposed to solve this problem. One of the most acknowledged techniques (node based formulation) has severe drawbacks such as: high sensitivity to mesh density and type of element integration (full vs reduced) and, more importantly, oscillatory solutions also leading to lack of convergence. We propose a new solution consisting in a non‐local smoothing technique. It defines, as the mechanical stimulus governing bone adaptation in a certain integration point of the mesh, the average of the stimuli obtained in the neighbour integration points. That average is weighted with a decay function of the distance to the centre of the neighbourhood. The new technique has been shown to overcome all the referred problems and perform in a robust way. It was tested on a hollow cylinder, resembling the diaphysis of a long bone, subjected to bending or torsion. Checkerboard instability was eliminated and local convergence of bone adaptation was achieved rapidly, in contrast to the other averaging technique and to the model without control of checkerboard instability. The new algorithm was also tested with good results on the same geometry but in a model containing a void, which produces a stress concentration that usually leads to checkerboard instability, like in other applications such as simulations of bone‐implant interfaces.

This content is not available in your region!

Continue researching here.

Having issues? You can contact us here