A novel algorithm to resolve lack of convergence and checkerboard instability in bone adaptation simulations using non-local averaging.

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.