Predictive structural models resulting in a trabecular bone topology closely resembling real bone would be a step toward 3D printing of sympathetic prosthetics. This study modifies an established trabecular bone structural adaptation approach, with the objective of achieving an improved adapted topology, specifically connectivity, compared to CT imaging studies; whilst retaining continuum level mechanical properties consistent with those reported in experimental studies. Strain driven structural adaptation models successfully identify trabecular trajectories, although tend to overpredict connectivity and skew trabecular radii distribution towards the smallest radius included in the adaptation. Radius adaptation of each trabecula is driven by a mechanostat approach with a target strain (1250 µɛ) below which radius is decreased (resorption), and above which radius is increased (apposition). Simulations include a lazy zone, in which neither resorption nor apposition takes place (1000 to 1500 µɛ); and a dead zone (<250 µɛ) in which complete resorption of trabeculae with the smallest included radius takes place. This study assesses the impact of increasing the dead zone threshold from <250 µɛ to <1000 µɛ, the lower limit of the lazy zone. In-silico structural models with an initial connectivity (number of trabeculae connecting at each joint) of 14 were generated using a nearest neighbour approach applied to a random cloud of points. Trabeculae were modelled using circular beams whose radii were adapted in response to normal strains caused by the axial force and bending moments due to a vertical pressure of 1 MPa applied to the top of the lattice, with the bottom of the lattice fixed in the vertical direction. Lattices in which nodes are either able (rigid jointed) or unable (pin jointed) to transmit bending moments were considered. Five virtual samples of each lattice type were used, and each simulation repeated twice: with a dead zone of either <250 µɛ or <1000 µɛ.Abstract
1.0 Objectives
2.0 Methods
