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.     

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.     

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.     

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.     

Current density in a model of a human body with a conductive implant exposed to ELF electric and magnetic fields

A numerical model of a human body with an intramedullary nail in the femur was built to evaluate the effects of the implant on the current density distribution in extremely low frequency electric and magnetic fields. The intramedullary nail was chosen because it is one of the longest high conductive implants used in the human body. As such it is expected to alter the electric and magnetic fields significantly. The exposure was a simultaneous combination of inferior to superior electric field and posterior to anterior magnetic field both alternating at 50 Hz with the values corresponding to the ICNIRP reference levels: 5000 V m^?1 for electric field and 100 µT for magnetic flux density. The calculated current density distribution inside the model was compared to the ICNIRP basic restrictions for general public (2 mA m^?2). The results show that the implant significantly increases the current density up to 9.5 mA m^?2 in the region where it is in contact with soft tissue in the model with the implant in comparison to 0.9 mA m^?2 in the model without the implant. As demonstrated the ICNIRP basic restrictions are exceeded in a limited volume of the tissue in spite of the compliance with the ICNIRP reference levels for general public, meaning that the existing safety limits do not necessarily protect implanted persons to the same extent as they protect people without implants. Bioelectromagnetics 30:591–599, 2009. © 2009 Wiley‐Liss, Inc.     

Implant–bone interface healing and adaptation in resurfacing hip replacement

Hip resurfacing demonstrates good survivorship as a treatment for young patients with osteoarthritis, but occasional implant loosening failures occur. On the femoral side there is radiographic evidence suggesting that the implant stem bears load, which is thought to lead to proximal stress shielding and adaptive bone remodelling. Previous attempts aimed at reproducing clinically observed bone adaptations in response to the implant have not recreated the full set of common radiographic changes, so a modified bone adaptation algorithm was developed in an attempt to replicate more closely the effects of the prosthesis on the host bone. The algorithm features combined implant–bone interface healing and continuum bone remodelling. It was observed that remodelling simulations that accounted for progressive gap filling at the implant–bone interface predicted the closest periprosthetic bone density changes to clinical X-rays and DEXA data. This model may contribute to improved understanding of clinical failure mechanisms with traditional hip resurfacing designs and enable more detailed pre-clinical analysis of new designs.     

Experimental and probabilistic analysis of distal femoral periprosthetic fracture: a comparison of locking plate and intramedullary nail fixation. Part A: experimental investigation

The following is a two-part study. Part A evaluates biomechanically intramedullary (IM) nails vs. locking plates for fixation of femoral fractures in osteoporotic bone. Part B of this study introduces a deterministic finite element model of each construct type and investigates the probability of periprosthetic fracture of the locking plate compared with the retrograde IM nail using Monte Carlo simulation. For Part A, an extra-articular, metaphyseal wedge fracture pattern was created in 11 osteoporotic fourth-generation composite femurs. Fixation was performed with a locking plate or a retrograde IM nail. Axial, torsion and bending cyclic loading to simulate post-operative damage accumulation were performed followed by ramped load to failure. Locking plates proved to be more stable (using stiffness as the determining factor) in osteoporotic bone as observed under low load cycle conditions. However, some of these advantages were offset by a greater incidence of sudden periprosthetic fracture observed under ramped loading conditions. Cadaveric, osteoporotic femurs included as a case study also exhibited periprosthetic fracture, but failure was accompanied by catastrophic comminution of the cortex. Periprosthetic failure at the implant end including bone comminution is difficult to salvage with revision fixation. The weakened trabecular matrix and thinned cortex of osteoporotic bone may increase the incidence of periprosthetic fracture. It is, therefore, essential for the surgeon to consider all possible loading scenarios when recommending an ideal implant for the osteoporotic patient.     

Bone Remodelling of a Proximal Femur with the Thrust Plate Prosthesis: An In Vitro Case

The key to the development of a successful implant is an understanding of the effect of bone remodelling on its long-term fixation. In this study, clinically observed patterns of bone remodelling have been compared with computer-based predictions for one particular design of prosthesis, the Thrust Plate Prosthesis (Centerpulse Orthopedics, Winterthur, Switzerland). Three-dimensional finite-element models were created using geometrical and bone density data obtained from CT scanning. Results from the bone remodelling simulation indicated that varying the relative rate of bone deposition/resorption and the interfacial conditions between the bone and the implant could produce the trend towards the two clinically observed patterns of remodelling.     

Bone remodelling of a proximal femur with the thrust plate prosthesis: an in vitro case

The key to the development of a successful implant is an understanding of the effect of bone remodelling on its long-term fixation. In this study, clinically observed patterns of bone remodelling have been compared with computer-based predictions for one particular design of prosthesis, the Thrust Plate Prosthesis (Centerpulse Orthopedics, Winterthur, Switzerland). Three-dimensional finite-element models were created using geometrical and bone density data obtained from CT scanning. Results from the bone remodelling simulation indicated that varying the relative rate of bone deposition/resorption and the interfacial conditions between the bone and the implant could produce the trend towards the two clinically observed patterns of remodelling.     

Implant-bone interface healing and adaptation in resurfacing hip replacement

Hip resurfacing demonstrates good survivorship as a treatment for young patients with osteoarthritis, but occasional implant loosening failures occur. On the femoral side there is radiographic evidence suggesting that the implant stem bears load, which is thought to lead to proximal stress shielding and adaptive bone remodelling. Previous attempts aimed at reproducing clinically observed bone adaptations in response to the implant have not recreated the full set of common radiographic changes, so a modified bone adaptation algorithm was developed in an attempt to replicate more closely the effects of the prosthesis on the host bone. The algorithm features combined implant-bone interface healing and continuum bone remodelling. It was observed that remodelling simulations that accounted for progressive gap filling at the implant-bone interface predicted the closest periprosthetic bone density changes to clinical X-rays and DEXA data. This model may contribute to improved understanding of clinical failure mechanisms with traditional hip resurfacing designs and enable more detailed pre-clinical analysis of new designs.     

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.     

Correlation between pre-operative periprosthetic bone density and post-operative bone loss in THA can be explained by strain-adaptive remodelling.

Periprosthetic adaptive bone remodelling after total hip arthroplasty can be simulated in computer models, combining bone remodelling theory with finite element analysis. Patient specific three-dimensional finite element models of retrieved bone specimens from an earlier bone densitometry (DEXA) study were constructed and bone remodelling simulations performed. Results of the simulations were analysed both qualitatively and quantitatively. Patterns of predicted bone loss corresponded very well with the DEXA measurements on the retrievals. The amount of predicted bone loss, measured quantitatively by simulating DEXA on finite element models, was found to be inversely correlated with the initial bone mineral content. It was concluded that the same clinically observed correlation can therefore be explained by mechanically induced remodelling. This finding extends the applicability of numerical pre-clinical testing to the analysis of interaction between implant design and initial state of the bone.     

Influence of ingrowth regions on bone remodelling around a cementless hip resurfacing femoral implant

Hip resurfacing arthroplasty is an alternative to traditional hip replacement that can conserve proximal bone stock and has gained popularity but bone resorption may limit implant survival and remains a clinical concern. The goal of this study was to analyze bone remodelling patterns around an uncemented resurfacing implant and the influence of ingrowth regions on resorption. A computed tomography-derived finite element model of a proximal femur with a virtually implanted resurfacing component was simulated under peak walking loads. Bone ingrowth was simulated by six interface conditions: fully bonded; fully friction; bonded cap with friction stem; a small bonded region at the stem-cup intersection with the remaining surface friction; fully frictional, except for a bonded band along the distal end of the cap and superior half of the cap bonded with the rest frictional. Interface condition had a large influence on remodelling patterns. Bone resorption was minimized when no ingrowth occurred at the bone-implant interface. Bonding only the superior half of the cap increased bone resorption slightly but allowed for a large ingrowth region to improve secondary stability.     

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.     

Development of a locking femur nail for mice

We herein report on a novel locking intramedullary nail system in a murine closed femur fracture model. The nail system consists of a modified 24-gauge injection needle and a 0.1-mm-diameter tungsten guide wire. Rotation stability was accomplished by flattening the proximal and distal end of the needle. Torsional mechanical testing of the implants in osteotomized cadaveric femora revealed a superiority of the locking nail (3.9+/-1.0 degrees rotation at a torque of 0.9 Nmm, n=10) compared to the unmodified injection needle (conventional nail; 52.4+/-3.2 degrees, n=10, p<0.05). None of the implants, however, achieved the rotation stability of unfractured femora (0.3+/-0.5 degrees, n=10). In a second step, we tested the feasibility of the in vivo application of the locking nail to stabilize a closed femoral midshaft fracture in C57BL/6 mice. Of interest, none of the 10 animals showed a dislocation of the locking nail over a 5-week period, while 3 of 4 animals with conventional nail fracture stabilization showed a significant pin dislocation within the first 3 days (p<0.05). Mechanical testing after 5-weeks stabilization with the locking nail revealed an appropriate bone healing with a torque at failure of 71.6+/-3.4% and a peak rotation before failure of 68.4+/-5.3% relative to the unfractured contralateral femur. With the advantage that closed fractures can be fixed with rotation stability, the herein introduced model may represent an ideal tool to study bone healing in transgenic and knockout mice.     

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.     

Application of an anisotropic bone-remodelling model based on a damage-repair theory to the analysis of the proximal femur before and after total hip replacement

In this work, a new model for internal anisotropic bone remodelling is applied to the study of the remodelling behaviour of the proximal femur before and after total hip replacement (THR). This model considers bone remodelling under the scope of a general damage-repair theory following the principles of continuum damage mechanics. A "damage-repair" tensor is defined in terms of the apparent density and Cowin's "fabric tensor", respectively, associated with porosity and directionality of the trabeculae. The different elements of a thermodynamically consistent damage theory are established, including resorption and apposition criteria, evolution law and rate of remodelling. All of these elements were introduced and discussed in detail in a previous paper (García, J. M., Martinez, M. A., Doblaré, M., 2001. An anisotrophic internal-external bone adaptation model based on a combination of CAO and continuum damage mechanics technologies. Computer Methods in Biomechanics and Biomedical Engineering 4(4), 355-378.), including the definition of the proposed mechanical stimulus and the qualitative properties of the model. In this paper, the fundamentals of the proposed model are briefly reviewed and the computational aspects of its implementation are discussed. This model is then applied to the analysis of the remodelling behaviour of the intact femur obtaining densities and mass principal values and directions very close to the experimental data. The second application involved the proximal femoral extremity after THR and the inclusion of an Exeter prosthesis. As a result of the simulation process, some well-known features previously detected in medical clinics were recovered, such as the stress yielding effect in the proximal part of the implant or the enlargement of the cortical layer at the distal part of the implant. With respect to the anisotropic properties, bone microstructure and local stiffness are known to tend to align with the stress principal directions. This experimental fact is mathematically proved in the framework of this remodelling model and clearly shown in the results corresponding to the intact femur. After THR the degree of anisotropy decreases tending, specifically in the proximal femur, to a more isotropic behaviour.     

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.     

Finite element analysis of an additional implant for intramedullary nailing]

In the surgical treatment of fractured femurs, the fracture is bridged by a medullary nail fixed in the bone with interlocking screws. Failure of bone substance in the region of the interlocking screws is the most common complication in the treatment of osteoporotic bone. With the aim of preventing this complication, an additional implant was developed. A finite element analysis of an ideal bone/implant system was carried out to investigate the role of the additional implant. Three defined finite element models were generated, and the associated stress situations compared. The first model is a standard fixation without the additional implant. In the second model, the additional implant is integrated within the bone/implant system. The third model uses a modified form of the additional implant. The results show that both additional implants reduce the stresses occurring, both in the bone substance and at the screws. The modified form of the additional implant proved to be the most favorable version. In the case of the original additional implant, the negative effect of the sharp edges of the thread was demonstrable.     

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.