Computational simulation of internal bone remodelling around dental implants: a sensitivity analysis

This study aimed to predict the distribution of bone trabeculae, as a density change per unit time, around a dental implant based on applying a selected mathematical remodelling model. The apparent bone density change as a function of the mechanical stimulus was the base of the applied remodelling model that describes disuse and overload bone resorption. The simulation was tested in a finite element model of a screw-shaped dental implant in an idealised bone segment. The sensitivity of the simulation to different mechanical parameters was investigated; these included element edge length, boundary conditions, as well as direction and magnitude of the implant loads. The alteration in the mechanical parameters had a significant influence on density distribution and model stability, in particular at the cortical bone region. The remodelling model could succeed to achieve trabeculae-like structure around osseointegrated dental implants. The validation of this model to a real clinical case is required.     

Numerical simulation of the remodelling process of trabecular architecture around dental implants

Dental implants may alter the mechanical environment in the jawbone, thereby causing remodelling and adaptation of the surrounding trabecular bone tissues. To improve the efficacy of dental implant systems, it is necessary to consider the effect of bone remodelling on the performance of the prosthetic systems. In this study, finite element simulations were implemented to predict the evolution of microarchitecture around four implant systems using a previously developed model that combines both adaptive and microdamage-based mechano-sensory mechanisms in bone remodelling process. Changes in the trabecular architecture around dental implants were mainly focused. The simulation results indicate that the orientational and ladder-like architecture around the implants predicted herein is in good agreement with those observed in animal experiments and clinical observations. The proposed algorithms were shown to be effective in simulating the remodelling process of trabecular architecture around dental implant systems. In addition, the architectural features around four typical dental implant systems in alveolar bone were evaluated comparatively.     

Numerical investigation of bone remodelling around immediately loaded dental implants using sika deer (Cervus nippon) antlers as implant bed

This study combines finite element method and animal studies, aiming to investigate tissue remodelling processes around dental implants inserted into sika deer antler and to develop an alternative animal consuming model for studying bone remodelling around implants. Implants were inserted in the antlers and loaded immediately via a self-developed loading device. After 3, 4, 5 and 6 weeks, implants and surrounding tissue were taken out. Specimens were scanned by μCT scanner and finite element models were generated. Immediate loading and osseointegration conditions were simulated at the implant-tissue interface. A vertical force of 10 N was applied on the implant. During the healing time, density and Young’s modulus of antler tissue around the implant increased significantly. For each time point, the values of displacement, stresses and strains in the osseointegration model were lower than those of the immediate loading model. As the healing time increased, the displacement of implants was reduced. The 3-week immediate loading model (9878 ± 1965 μstrain) illustrated the highest strains in the antler tissue. Antler tissue showed similar biomechanical properties as human bone in investigating the bone remodelling around implants, therefore the use of sika deer antler model is a promising alternative in implant biomechanical studies.     

Bone remodelling around the tibia due to total ankle replacement: effects of implant material and implant–bone interfacial conditions

Abstract

One of the major causes of implant loosening is due to excessive bone resorption surrounding the implant due to bone remodelling. The objective of the study is to investigate the effects of implant material and implant–bone interface conditions on bone remodelling around tibia bone due to total ankle replacement. Finite element models of intact and implanted ankles were developed using CT scan data sets. Bone remodelling algorithm was used in combination with FE analysis to predict the bone density changes around the ankle joint. Dorsiflexion, neutral, and plantar flexion positions were considered, along with muscle force and ligaments. Implant–bone interfacial conditions were assumed as debonded and bonded to represent non-osseointegration and fully osseointegration at the porous coated surface of the implant. To investigate the effect of implant material, three finite element models having different material combinations of the implant were developed. For model 1, tibial and talar components were made of Co–Cr–Mo, and meniscal bearing was made of UHMWPE. For model 2, tibial and talar components were made of ceramic and meniscal bearing was made of UHMWPE. For model 3, tibial and talar components were made of ceramic and meniscal bearing was made of CFR-PEEK. Changes in implant material showed no significant changes in bone density due to bone remodelling. Therefore, ceramic appears to be a viable alternative to metal and CFR-PEEK can be used in place of UHMWPE. This study also indicates that proper bonding between implant and bone is essential for long-term survival of the prosthetic components.     

To reduce the maximum stress and the stress shielding effect around a dental implant–bone interface using radial functionally graded biomaterials

In a dental implant system, the value of stress and its distribution plays a pivotal role on the strength, durability and life of the implant–bone system. A typical implant consists of a Titanium core and a thin layer of biocompatible material such as the hydroxyapatite. This coating has a wide range of clinical applications in orthopedics and dentistry due to its biocompatibility and bioactivity characteristics. Low bonding strength and sudden variation of mechanical properties between the coating and the metallic layers are the main disadvantages of such common implants. To overcome these problems, a radial distributed functionally graded biomaterial (FGBM) was proposed in this paper and the effect of material property on the stress distribution around the dental implant–bone interface was studied. A three-dimensional finite element simulation was used to illustrate how the use of radial FGBM dental implant can reduce the maximum von Mises stress and, also the stress shielding effect in both the cortical and cancellous bones. The results, of course, give anybody an idea about optimized behaviors that can be achieved using such materials. The finite element solver was validated by familiar methods and the results were compared to previous works in the literature.     

Numerical estimation of bone density and elastic constants distribution in a human mandible

In this paper, we try to predict the distribution of bone density and elastic constants in a human mandible, based on the stress level produced by mastication loads using a mathematical model of bone remodelling. These magnitudes are needed to build finite element models for the simulation of the mandible mechanical behavior. Such a model is intended for use in future studies of the stability of implant-supported dental prostheses. Various models of internal bone remodelling, both phenomenological and more recently mechanobiological, have been developed to determine the relation between bone density and the stress level that bone supports. Among the phenomenological models, there are only a few that are also able to reproduce the level of anisotropy. These latter have been successfully applied to long bones, primarily the femur. One of these models is here applied to the human mandible, whose corpus behaves as a long bone. The results of bone density distribution and level of anisotropy in different parts of the mandible have been compared with various clinical studies, with a reasonable level of agreement.     

Three-dimensional design optimisation of patient-specific femoral plates as a means of bone remodelling reduction

Previous investigations into the optimisation of internal plates have mostly focused on the material properties of the implant. In this work, we optimise the shape, size and placement of the plate for successfully minimising bone remodelling around the implant. A design optimisation algorithm based on strain energy density criterion, combined with the finite element analysis, has been used in this study. The main optimisation goal was to reduce this change and keep it close to the conditions of an intact femur. The results suggest that the anterolateral side of the bone would be the optimum location for the plate, as for the geometry, the optimum moves towards having a thick, wide and short plate. These important results could be directly applicable to orthopaedic surgeons treating a femur fracture with internal plates. Since the optimisation algorithm remains the same for any patient, this advancement provides the surgeon with a tool to minimise the post surgery remodelling by trying to maintain the natural structure of the bone.     

Evaluation of bone remodeling around single dental implants of different lengths: a mechanobiological numerical simulation and validation using clinical data

Algorithmic models have been proposed to explain adaptive behavior of bone to loading; however, these models have not been applied to explain the biomechanics of short dental implants. Purpose of present study was to simulate bone remodeling around single implants of different lengths using mechanoregulatory tissue differentiation model derived from the Stanford theory, using finite elements analysis (FEA) and to validate the theoretical prediction with the clinical findings of crestal bone loss. Loading cycles were applied on 7-, 10-, or 13-mm-long dental implants to simulate daily mastication and bone remodeling was assessed by changes in the strain energy density of bone after a 3, 6, and 12 months of function. Moreover, clinical findings of marginal bone loss in 45 patients rehabilitated with same implant designs used in the simulation (n = 15) were computed to validate the theoretical results. FEA analysis showed that although the bone density values reduced over time in the cortical bone for all groups, bone remodeling was independent of implant length. Clinical data showed a similar pattern of bone resorption compared with the data generated from mathematical analyses, independent of implant length. The results of this study showed that the mechanoregulatory tissue model could be employed in monitoring the morphological changes in bone that is subjected to biomechanical loads. In addition, the implant length did not influence the bone remodeling around single dental implants during the first year of loading.     

Comparative analysis of bone remodelling models with respect to computerised tomography-based finite element models of bone

Subject-specific finite element models are an extensively used tool for the numerical analysis of the biomechanical behaviour of human bones. However, bone modelling is not an easy task due to the complex behaviour of bone tissue, involving non-homogeneous and anisotropic mechanical properties. Moreover, bone is a living tissue and therefore its microstructure and mechanical properties evolve with time in a known process called bone remodelling. This phenomenon has been widely studied, many being the numerical models that have been formulated to predict density distribution and its evolution in several bones. The aim of the present study is to assess the capability of a bone remodelling model to predict the bone density distribution of different types of human bone (femur, tibia and mandible) comparing the obtained results with the bone density estimated by means of computerised tomography. Good accuracy was observed for the bone remodelling predictions including the thickness of the cortical layer.     

The effect of muscle loading on the simulation of bone remodelling in the proximal femur

A large number of finite element analyses of the proximal femur rely on a simplified set of muscle and joint contact loads to represent the boundary conditions of the model. In the context of bone remodelling analysis around hip implants, muscle loading affects directly the spatial distribution of the remodelling signal. In the present study we performed a sensitivity analysis on the effect of different muscle loading configurations on the outcome of the bone remodelling simulation. An anatomical model of the femur with the implanted stem in place was constructed using the CT data of the Visible Human Project dataset of the National Institute of Health. The model was loaded with three muscle force configurations with increasing level of complexity. A strain adaptive remodelling rule was employed to simulate the post-operative bone changes around the implant stem and the results of the simulation were assessed quantitatively in terms of the bone mineral content changes in 18 periprosthetic regions of interest. The results showed considerable differences in the amount of bone loss predicted between the three cases. The simplified models generally predicted more pronounced bone loss. Although the overall remodelling patterns observed were similar, the bone conserving effect of additional muscle forces in the vicinity of their areas of attachment was clear. The results of this study suggest that the loading configuration of the FE model does play an important role in the outcome of the remodelling simulation.     

Bone remodelling analysis of a bovine femur for a veterinary implant design

The response of bovine bone to the presence of an implant is analysed with the aim of simulating bone remodelling in a developing model of a polymeric intramedullary interlocking nail for veterinary use. A 3-D finite element model of the femur diaphysis is built based on computed tomography images and using a CAD-based modelling pipeline. The bone remodelling process after the surgery is analysed and compared with the healthy bone. The remodelling law assumes that bone adapts to the mechanical environment. For the analyses a consistent set of loads is determined for the bovine walk cycle. The remodelling results reproduce the morphologic features of bone and provide evidence of the difference on the bone behaviour when comparing metallic and polymeric nails. Our findings indicate that an intramedullary polymeric nail has the advantage over the metallic one of improving long-term bone healing and possibly avoiding the need of the implant removal.     

Bone remodelling around a cementless glenoid component

Post-operative change in the mechanical loading of bone may trigger its (mechanically induced) adaptation and hamper the mechanical stability of prostheses. This is especially important in cementless components, where the final fixation is achieved by the bone itself. The aim of this study is, first, to gain insight into the bone remodelling process around a cementless glenoid component, and second, to compare the possible bone adaptation when the implant is assumed to be fully bonded (best case scenario) or completely loose (worst case scenario). 3D finite element models of a scapula with and without a cementless glenoid component were created. 3D geometry of the scapula, material properties, and several physiological loading conditions were acquired from or estimated for a specific cadaver. Update of the bone density after implantation was done according to a node-based bone remodelling scheme. Strain energy density for different loading conditions was evaluated, weighted according to their frequencies in activities of daily life and used as a mechanical stimulus for bone adaptation. The average bone density in the glenoid increased after implantation. However, local bone resorption was significant in some regions next to the bone-implant interface, regardless of the interface condition (bonded or loose). The amount of bone resorption was determined by the condition imposed to the interface, being slightly larger when the interface was loose. An ideal screw, e.g. in which material fatigue was not considered, was enough to keep the interface micromotions small and constant during the entire bone adaptation simulation.     

Activity intensity,assistive devices and joint replacement influence predicted remodelling in the proximal femur

Bone morphology and density changes are commonly observed following joint replacement, may contribute to the risks of implant loosening and periprosthetic fracture and reduce the available bone stock for revision surgery. This study was presented in the ‘Bone and Cartilage Mechanobiology across the scales’ WCCM symposium to review the development of remodelling prediction methods and to demonstrate simulation of adaptive bone remodelling around hip replacement femoral components, incorporating intrinsic (prosthesis) and extrinsic (activity and loading) factors. An iterative bone remodelling process was applied to finite element models of a femur implanted with a cementless total hip replacement (THR) and a hip resurfacing implant. Previously developed for a cemented THR implant, this modified process enabled the influence of pre- to post-operative changes in patient activity and joint loading to be evaluated. A control algorithm used identical pre- and post-operative conditions, and the predicted extents and temporal trends of remodelling were measured by generating virtual X-rays and DXA scans. The modified process improved qualitative and quantitative remodelling predictions for both the cementless THR and resurfacing implants, but demonstrated the sensitivity to DXA scan region definition and appropriate implant–bone position and sizing. Predicted remodelling in the intact femur in response to changed activity and loading demonstrated that in this simplified model, although the influence of the extrinsic effects were important, the mechanics of implantation were dominant. This study supports the application of predictive bone remodelling as one element in the range of physical and computational studies, which should be conducted in the preclinical evaluation of new prostheses.     

Predictions of bone remodeling around dental implant systems

This study presents the implementation of a mathematical bone remodeling algorithm to bone adaptation in the premolar area of the mandible around various dental implant systems, and thus sheds a new perspective to the complex interactions in dental implant mechanics. A two-dimensional, plane strain model of the bone was built from a CT-scan. The effect of implant contour on internal bone remodeling was investigated by considering four dental implant systems with contours similar to commercially available ones and another four with cylindrical and conical cross-sections. The remodeling algorithm predicts non-homogeneous density/elastic modulus distribution; and, implant contour has some effect on how this is distributed. Bone density is predicted to increase on the tips of the threads of the implants, but to decrease inside the grooves. Threadless implants favor to develop a softer bone around their periphery, compared to implant systems that have threads. The overall contour (dimensions and the shape) of an implant affect the bone density redistribution, but the differences between different implant systems are relatively small.     

Simulation of bone remodelling in orthodontic treatment

Orthodontic treatments not only displace irregular teeth but also induce responses in surrounding bone tissues. Bone remodelling is regarded as the regulatory mechanism triggered by mechanical loading. This study was aimed at investigating the effect of orthodontic loading on both tooth movement and neighbouring bone density distribution. A set of computational algorithms incorporating both external and internal remodelling mechanisms was implemented into a patient-specific 3D finite element (FE) model to investigate and analyse orthodontic treatment under four typical modes of orthodontic loading. The consequence of orthodontic treatment was reproduced numerically by using this FE-based technique. The results indicated that the diverse modes of orthodontic loading would result in different magnitudes of tooth movement and particular morphology of bone density distribution. It is illuminated that the newly developed algorithms may replicate the clinical situation more closely compared with the previous proposed method.     

Numerical simulation of load-induced bone structural remodelling using stress-limit criterion

A simple and efficient numerical method for predicting the remodelling of adaptive materials and structures under applied loading was presented and implemented within a finite element framework. The model uses the trajectorial architecture theory of optimisation to predict the remodelling of material microstructure and structural organisation under mechanical loading. We used the proposed model to calculate the density distribution of proximal femur in the frontal plane. The loading considered was the hip joint contact forces and muscular forces at the attachment sites of the muscles to the bone. These forces were estimated from a separate finite element calculation using a heterogeneous three-dimensional model of the proximal femur. The density distributions obtained by this procedure has a qualitative similarity with in vivo observations. Solutions displayed the characteristic high-density channels that are evident in the Dual X-ray Absorptiometry scan. There is also evidence of the intramedullary canal, as well as low-density regions in the femoral neck. Several parametric studies were carried out to highlight the advantages of the proposed method, which includes fast convergence and low-computational cost. The potential applications of the proposed method in predicting bone structural remodelling in cancer are also briefly discussed.     

Subject-specific bone remodelling of the scapula

Finite element analyses, with increasing levels of detail and complexity, are becoming effective tools to evaluate the performance of joint replacement prostheses and to predict the behaviour of bone. As a first step towards the study of the complications of shoulder arthroplasty, the aim of this work was the development and validation of a 3D finite element model of an intact scapula for the prediction of the bone remodelling process based on a previously published model that attempts to follow Wolff's law. The boundary conditions applied include full muscle and joint loads taken from a multibody system of the upper limb based on the same subject whose scapula was here analysed. To validate the bone remodelling simulations, qualitative and quantitative comparisons between the predicted and the specimen's bone density distribution were performed. The results showed that the bone remodelling model was able to successfully reproduce the actual bone density distribution of the analysed scapula.     

Influence of mastication and edentulism on mandibular bone density

The aim of this study was to demonstrate that external loading due to daily activities, including mastication, speech and involuntary open–close cycles of the jaw contributes to the internal architecture of the mandible. A bone remodelling algorithm that regulates the bone density as a function of stress and loading cycles is incorporated into finite element analysis. A three-dimensional computational model is constructed on the basis of computerised tomography (CT) images of a human mandible. Masticatory muscle activation involved during clenching is modelled by static analysis using linear optimisation. Other loading conditions are approximated by imposing mandibular flexure. The simulations predict that mandibular bone density distribution results in a tubular structure similar to what is observed in the CT images. Such bone architecture is known to provide the bone optimum strength to resist bending and torsion during mastication while reducing the bone mass. The remodelling algorithm is used to simulate the influence of edentulism on mandibular bone loss. It is shown that depending on the location and number of missing teeth, up to one-third of the mandibular bone mass can be lost due to lack of adequate mechanical stimulation.     

A Comparative Analysis of Different Treatments for Distal Femur Fractures using the Finite Element Method

The main objective of this work is the evaluation, by means of the finite element method (FEM) of the mechanical stability and long-term microstructural modifications in bone induced to three different kinds of fractures of the distal femur by three types of implants: the Condyle Plate, the less invasive stabilization system plate (LISS) and the distal femur nail (DFN). The displacement and the stress distributions both in bone and implants and the internal bone remodelling process after fracture and fixation are obtained and analysed by computational simulation. The main conclusions of this work are that distal femoral fractures can be treated correctly with the Condyle Plate, the LISS plate and the DFN. The stresses both in LISS and DFN implant are high especially around the screws. When respect to remodelling, the LISS produces an important resorption in the fractured region, while the other two implants do not strongly modify bone tissue microstructure.     

Trabecular bone strains around a dental implant and associated micromotions—A micro-CT-based three-dimensional finite element study

The first objective of this computational study was to assess the strain magnitude and distribution within the three-dimensional (3D) trabecular bone structure around an osseointegrated dental implant loaded axially. The second objective was to investigate the relative micromotions between the implant and the surrounding bone. The work hypothesis adopted was that these virtual measurements would be a useful indicator of bone adaptation (resorption, homeostasis, formation).In order to reach these objectives, a μCT-based finite element model of an oral implant implanted into a Berkshire pig mandible was developed along with a robust software methodology. The finite element mesh of the 3D trabecular bone architecture was generated from the segmentation of μCT scans. The implant was meshed independently from its CAD file obtained from the manufacturer. The meshes of the implant and the bone sample were registered together in an integrated software environment. A series of non-linear contact finite element (FE) analyses considering an axial load applied to the top of the implant in combination with three sets of mechanical properties for the trabecular bone tissue was devised. Complex strain distribution patterns are reported and discussed. It was found that considering the Young’s modulus of the trabecular bone tissue to be 5, 10 and 15 GPa resulted in maximum peri-implant bone microstrains of about 3000, 2100 and 1400. These results indicate that, for the three sets of mechanical properties considered, the magnitude of maximum strain lies within an homeostatic range known to be sufficient to maintain/form bone. The corresponding micro-motions of the implant with respect to the bone microstructure were shown to be sufficiently low to prevent fibrous tissue formation and to favour long-term osseointegration.