Young patients receiving metallic bone implants after surgical resection of bone cancer require implants that last into adulthood, and ideally life-long. Porous implants with similar stiffness to bone can promote bone ingrowth and thus beneficial clinical outcomes. A mechanical remodelling stimulus, strain energy density (SED), is thought to be the primary control variable of the process of bone growth into porous implants. The sequential process of bone growth needs to be taken into account to develop an accurate and validated bone remodelling algorithm, which can be employed to improve porous implant design and achieve better clinical outcomes. A bone remodelling algorithm was developed, incorporating the concept of bone connectivity (sequential growth of bone from existing bone) to make the algorithm more physiologically relevant. The algorithm includes adaptive elastic modulus based on apparent bone density, using a node-based model to simulate local remodelling variations while alleviating numerical checkerboard problems. Strain energy density (SED) incorporating stress and strain effects in all directions was used as the primary stimulus for bone remodelling. The simulations were developed to run in MATLAB interfacing with the commercial FEA software ABAQUS and Python. The algorithm was applied to predict bone ingrowth into a porous implant for comparison against data from a sheep model.Abstract
Objectives
Methods
Current therapies for osteoporosis are limited to generalised antiresorptive or anabolic interventions, which do not target specific regions to improve skeletal health. Moreover, the adaptive changes of separate and combined pharmacological and biomechanical treatments in the ovariectomised (OVX) mouse tibia has not been studied yet. Therefore, this study combines micro- computed tomography (micro-CT) imaging and computational modelling to evaluate the efficacies of treatments in reducing bone loss. In vivo micro-CT (10.4µm/voxel) images of the right tibiae of N=18 female OVX C57BL/6 mice were acquired at weeks 14, 16, 18, 20 and 22 of age for 3 groups: mechanical loading (ML), parathyroid hormone (PTH) or combined therapies (PTHML). All mice received either injection of PTH (100μg/kg/day, 5days/week) or vehicle from week 18. The right tibiae were mechanically loaded in vivo at week 19 and 21 with a 12N peak load, 40 cycles/day and 3 days/week. Bone adaptation was quantified through spatial changes in bone mineral density (BMD) and strain distribution was obtained from micro-CT-based finite element models.Abstract
Objectives
Methodology
Prediction of bone adaptation in response to mechanical loading is useful in the clinical management of osteoporosis. However, few studies have investigated the effect of repeated mechanical loading in the mouse tibia. Therefore, this study uses a combined experimental and computational approach to evaluate the effect of mechanical loading on bone adaptation in a mouse model of osteoporosis. Six female C57BL/6 mice were ovariectomised (OVX) at week 14 and scanned using in vivo micro computed tomography (10.4µm/voxel) at week 14, 16, 18, 20 and 22. The right tibiae were mechanically loaded in vivo at week 19 and 21 with a 12N peak load, 40 cycles/day, 3 days/week. Linear isotropic homogeneous finite element (microFE) models were created from the tissue mineral density calibrated microCT images. Changes in bone adaptation, densitometric and spatial analyses were measured by comparing the longitudinal images after image registration.Abstract
Objectives
Methods
Long-term survival of massive prostheses used to treat bone cancers is associated with extra-cortical bone growth and osteointegration into a grooved hydroxyapatite coated collar positioned adjacent to the transection site on the implant shaft [1]. The survivorship at 10 years reduces from 98% to 75% where osteointegration of the shaft does not occur. Although current finite element (FE) methods successfully model bone adaption, optimisation of adventitious new bone growth and osteointegration is difficult to predict. There is thus a need to improve existing FE models by including biological processes of osteoconduction and osteoinduction. The principal bone adaptation criteria is based on the standard strain-energy remodeling algorithm, where the rate of remodeling is controlled by the difference in the stimulus against the reference value [3]. The additional concept of bone connectivity was introduced, to limit bone growth to neighbouring elements (cells) adjoining existing bone elements. The algorithm was developed on a cylindrical model before it was used on an ovine model. The geometry and material properties from two ovine tibiae were obtained from computed tomography (CT) scans and used to develop FE models of the tibiae implanted with a grooved collar. The bones were assigned inhomogeneous material properties based on the CT grey values and typical ovine walking load conditions were applied. The FE results show a region of bone tissue growth below the implanted collar and a small amount of osteointegration with the implant, which is in good agreement to clinical results. Some histological results suggest that further bone growth is possible and potential improvements to the model will be discussed. In summary, by including an algorithm that describes osteoconduction, adventitious bone growth can be predicted.
