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.     

Finite element modelling of implant designs and cortical bone thickness on stress distribution in maxillary type IV bone

The aims of this study were to examine the effect of implant neck design and cortical bone thickness using 3D finite element analysis and to analyse the stability of clinical evidence based on micromotion and principal stress. Four commercial dental implants for a type IV bone and maxillary segments were created. Various parameters were considered, including the osseointegration condition, loading direction and cortical bone thickness. Micromotion and principal stresses were used to evaluate the failure of osseointegration and bone overloading, respectively. It was found that the maximum stress of the peri-implant bone decreased as cortical bone thickness increased. The micromotion level in full osseointegration is less than that in non-osseointegration and it also decreases as cortical bone thickness increases. The cortical bone thickness should be measured before surgery to help select a proper implant. In the early stage of implantation, the horizontal loading component induces stress concentration in bone around the implant neck more easily than does the vertical loading component, and this may result in crestal bone loss.     

Use of micro-CT-based finite element analysis to accurately quantify peri-implant bone strains: a validation in rat tibiae

Although research has been addressed at investigating the effect of specific loading regimes on bone response around the implant, a precise quantitative understanding of the local mechanical response close to the implant site is still lacking. This study was aimed at validating micro-CT-based finite element (μFE) models to assess tissue strains after implant placement in a rat tibia. Small implants were inserted at the medio-proximal site of 8 rat tibiae. The limbs were subjected to axial compression loading; strain close to the implant was measured by means of strain gauges. Specimen-specific μFE models were created and analyzed. For each specimen, 4 different models were created corresponding to different representations of the bone–implant interface: bone and implant were assumed fully osseointegrated (A); a low stiffness interface zone was assumed with thickness of 40 μm (B), 80 μm (C), and 160 μm (D). In all cases, measured and computational strains correlated highly (R ² = 0.95, 0.92, 0.93, and 0.95 in A, B, C, and D, respectively). The averaged calculated strains were 1.69, 1.34, and 1.15 times higher than the measured strains for A, B, and C, respectively, and lower than the experimental strains for D (factor = 0.91). In conclusion, we demonstrated that specimen-specific FE analyses provide accurate estimates of peri-implant bone strains in the rat tibia loading model. Further investigations of the bone-implant interface are needed to quantify implant osseointegration.     

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.     

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.     

Are plasticity models required to predict relative risk of lag screw cut-out in finite element models of trochanteric fracture fixation?

Using a finite element model of unstable trochanteric fracture stabilized with a sliding hip screw, the benefits of two plasticity-based formulations, Drucker–Prager and crushable foam, were evaluated and compared to the commonly used linear elastic model of trabecular bone in order to predict the relative risk of lag screw cut-out for five distinct load cases. The crushable foam plasticity formulation leads to a much greater strain localization, in comparison to the other two models, with large plastic strains in a localized region. The plastic zone predicted with Drucker–Prager is relatively more diffuse. Linear elasticity associated with a minimum principal strain criterion provides the smallest volume of elements susceptible to yielding for all loading modes. The region likely to undergo plastic deformation, as predicted by the linear elastic model, is similar to that obtained from plasticity-based formulations, which indicates that this simple criterion provides an adequate estimate of the risk of cut-out.     

Effect of including damage at the tissue level in the nonlinear homogenisation of trabecular bone

Being able to predict bone fracture or implant stability needs a proper constitutive model of trabecular bone at the macroscale in multiaxial, non-monotonic loading modes. Its macroscopic damage behaviour has been investigated experimentally in the past, mostly with the restriction of uniaxial cyclic loading experiments for different samples, which does not allow for the investigation of several load cases in the same sample as damage in one direction may affect the behaviour in other directions. Homogenised finite element models of whole bones have the potential to assess complicated scenarios and thus improve clinical predictions. The aim of this study is to use a homogenisation-based multiscale procedure to upscale the damage behaviour of bone from an assumed solid phase constitutive law and investigate its multiaxial behaviour for the first time. Twelve cubic specimens were each submitted to nine proportional strain histories by using a parallel code developed in-house. Evolution of post-elastic properties for trabecular bone was assessed for a small range of macroscopic plastic strains in these nine load cases. Damage evolution was found to be non-isotropic, and both damage and hardening were found to depend on the loading mode (tensile, compression or shear); both were characterised by linear laws with relatively high coefficients of determination. It is expected that the knowledge of the macroscopic behaviour of trabecular bone gained in this study will help in creating more precise continuum FE models of whole bones that improve clinical predictions.     

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.     

Enhancement of implant osseointegration by high-frequency low-magnitude loading

Background

Mechanical loading is known to play an important role in bone remodelling. This study aimed to evaluate the effect of high- and low-frequency axial loading, applied directly to the implant, on peri-implant bone healing and implant osseointegration.

Methodology

Titanium implants were bilaterally installed in rat tibiae. For every animal, one implant was loaded (test) while the other one was not (control). The test implants were randomly divided into 8 groups according to 4 loading regimes and 2 experimental periods (1 and 4 weeks). The loaded implants were subject to an axial displacement. Within the high- (HF, 40 Hz) or low-frequency (LF, 8 Hz) loading category, the displacements varied 2-fold and were ranked as low- or high-magnitude (LM, HM), respectively. The strain rate amplitudes were kept constant between the two frequency groups. This resulted in the following 4 loading regimes: 1) HF-LM, 40 Hz-8 µm; 2) HF-HM, 40 Hz-16 µm; 3) LF-LM, 8 Hz-41 µm; 4) LF-HM, 8 Hz-82 µm. The tissue samples were processed for resin embedding and subjected to histological and histomorphometrical analyses. Data were analyzed statistically with the significance set at p<0.05.

Principal Findings

After loading for 4 weeks, HF-LM loading (40 Hz-8 µm) induced more bone-to-implant contact (BIC) at the level of the cortex compared to its unloaded control. No significant effect of the four loading regimes on the peri-implant bone fraction (BF) was found in the 2 experimental periods.

Conclusions

The stimulatory effect of immediate implant loading on bone-to-implant contact was only observed in case of high-frequency (40 Hz) low-magnitude (8 µm) loading. The applied load regimes failed to influence the peri-implant bone mass.     

Accuracy of finite element predictions in sideways load configurations for the proximal human femur

Subject-specific finite element models have been used to predict stress-state and fracture risk in individual patients. While many studies analysed quasi-axial loading configurations, only few works simulated sideways load configurations, such as those arising in a fall. The majority among these latter directly predicted bone strength, without assessing elastic strain prediction accuracy. The aim of the present work was to evaluate if a subject-specific finite element modelling technique from CT data that accurately predicted strains in quasi-axial loading configurations is suitable to accurately predict strains also when applying low magnitude loads in sideways configurations. To this aim, a combined numerical-experimental study was performed to compare finite element predicted strains with strain-gauge measurements from three cadaver proximal femurs instrumented with sixteen strain rosettes and tested non-destructively under twelve loading configurations, spanning a wide cone (0-30° for both adduction and internal rotation angles) of sideways fall scenarios. The results of the present study evidenced a satisfactory agreement between experimentally measured and predicted strains (R(2) greater than 0.9, RMSE% lower than 10%) and displacements. The achieved strain prediction accuracy is comparable to those obtained in state of the art studies in quasi-axial loading configurations. Still, the presence of the highest strain prediction errors (around 30%) in the lateral neck aspect would deserve attention in future studies targeting bone failure.     

Development of residual strains in human vertebral trabecular bone after prolonged static and cyclic loading at low load levels

Development of irreversible residual strains in trabecular bone may be a mechanism by which age-related non-traumatic vertebral fractures occur. To investigate this concept, static and cyclic loading tests were conducted at low loading levels for cylindrical cores of cadaveric vertebral trabecular bone. Stresses were applied equivalent to elastic strains of either 750 or 1,500 microstrain. Creep strains were measured during the tests, which lasted for 125,000 seconds (about 35 h), and for an additional 125,000 seconds after complete unloading. Emphasis was placed on the residual strains that developed, defined as the strain remaining at the end of the unloading phase. The results indicated that appreciable residual strains did develop, and were similar for static and cyclic loading. Irrespective of the applied load levels and loading modes, the residual strains that remained after the unloading phase were similar in magnitude to the originally applied elastic strain. Extrapolation of the observed residual strains to full recovery indicated that the time that would be required for full recovery was over 20 times longer than the duration of the applied loads. These results indicate that human vertebral trabecular bone does not creep in a linear viscoelastic fashion at low stress levels, and that creep mechanisms dominate the residual strains regardless of the loading mode. Taken together, these findings support the concept that non-traumatic vertebral fractures may be related to long-term creep effects because the trabecular bone does not have sufficient time to recover mechanically from creep deformations accumulated by prolonged static or cyclic loading.     

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.     

Influence of clinically relevant factors on the immediate biomechanical surrounding for a series of dental implant designs

The objective of the present study was to assess the influence of various clinically relevant scenarios on the strain distribution in the biomechanical surrounding of five different dental implant macrogeometries. The biomechanical environment surrounding an implant, i.e., the cortical and trabecular bone, was modeled along with the implant. These models included two different values of the study parameters including loading conditions, trabecular bone elastic modulus, cortical/trabecular bone thickness ratio, and bone loss for five implant designs. Finite element analysis was conducted on the models and strain in the bones surrounding the implant was calculated. Bone volumes having strains in four different windows of 0-200?με, 200-1000?με, 1000-3000?με, and > 3000 με were measured and the effect of each biomechanical variable and their two-way interactions were statistically analyzed using the analysis of variance method. This study showed that all the parameters included in this study had an effect on the volume of bones in all strain windows, except the implant design, which affected only the 0-200?με and >3000?με windows. The two-way interaction results showed that interactions existed between implant design and bone loss, and loading condition, bone loss in the 200-1000?με window, and between implant design and loading condition in the 0-200 με window. Within the limitations of the present methodology, it can be concluded that although some unfavorable clinical scenarios demonstrated a higher volume of bone in deleterious strain levels, a tendency toward the biomechanical equilibrium was evidenced regardless of the implant design.     

An analytical and numerical study of the stability of bone remodelling theories: dependence on microstructural stimulus.

The origin of unstable bone remodelling simulations using strain-energy-based remodelling rules was studied mathematically in order to assess whether the unstable behavior was due to the mathematical rules proposed to characterize the processes, or to the numerical approximations used to exercise the mathematical predictions. A condition which is necessary for the stability of a strain-energy-based remodelling theory was derived analytically using the calculus of variation. The analytical result was derived using a simple elastic model which consists of a long beam loaded by an axial force and a bending moment. This loading situation mimics the coupling between local density and global density distributions seen in vivo. A condition necessary for a stable remodelling scheme is arrived at, but the conditions necessary to guarantee a stable remodelling scheme are not. In this remodelling scheme, the elastic modulus is proportional to volumetric density raised to an exponent n, and the microstructural stimulus is taken as the strain energy density divided by volumetric density raised to an exponent m. In order for a remodelling scheme to be stable in this loading situation, m must be greater than n. Finite-difference time-stepping is used to verify the predictions of the analytical study. These numerical studies appear to confirm the analytical studies. Physiologic interpretation of the behavior found with n greater than m indicates that this type of unstable behavior is unlikely to be observed in vivo. Since numerical approximations are not made in deriving this stability condition, we conclude that the mathematical rules proposed to characterize bone remodelling based on strain energy density should meet this condition to be relevant to physiologic bone remodelling.     

3-dimensional registration of micro-movements of vertebral implants by a specially adapted measuring system using magnetic transducers]

The aim of this study was to develop a biomechanical experimental set-up to quantify motion of ventrally inserted spinal implants at the implant/bone interface. The model we used was the vertebral column of the calf. Lumbar vertebrae L2 to L4 were "instrumented" with a screw-rod system. The adjacent vertebrae L1 and L5 were connected to a servohydraulic testing machine and axial compression applied. Shortening of the specimen and three-dimensional movement of the most cranial implant relative to the bone was recorded using 3 electromagnetic transducers. 100,000 cycles of axial loading varying between -0.5 kN and -1 kN were applied. Static shortening of the specimen of 8.5 mm and an elastic movement of 180 to 280 microns were measured. The greatest amplitude of single-plane motion was recorded in the sagittal plane in both static and elastic modes. Implant motion within each cycle was recorded accurately as load-displacement curves within a range of 1.35 to 30 microns. With this test set-up, primary stability of different spinal implant systems can be compared. The use of electromagnetic transducers permits three-dimensional implant motion analysis even when only a mono-axial servohydraulic testing machine is available.     

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.     

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.     

Primary stability of uncemented femoral resurfacing implants for varying interface parameters and material formulations during walking and stair climbing

Primary stability of uncemented resurfacing prosthesis is provided by an interference fit between the undersized implant and the reamed bone. Dependent on the magnitude of interference, the implantation process causes high shear forces and large strains which can exceed the elastic limit of cancellous bone. Plastification of the bone causes reduced stiffness and could lead to bone damage and implant loosening. The purpose in this study was to determine press-fit conditions which allow implantation without excessive plastic bone deformation and sufficient primary stability to achieve bone ingrowth. In particular, the influence of interference, bone quality and friction on the micromotion during walking and stair-climbing was investigated. Therefore elastic and plastic finite element (FE) models of the proximal femur were developed. Implantation was realized by displacing the prosthesis onto the femur while monitoring the contact pressure, plastic bone deformation as well as implantation forces. Subsequently a physiologic gait and stair-climbing cycle was simulated calculating the micromotion at the bone-implant interface. Results indicate that plastic deformation starts at an interference of 30 μm and the amount of plastified bone at the interface increases up to 90% at 150 μm interference. This effect did not reduce the contact pressure if interference was below 80 μm. The micromotion during walking was similar for the elastic and plastic FE models. A stable situation allowing bony ingrowth was achieved for both constitutive laws (elastic, plastic) for walking and stair climbing with at least 60 μm press-fit, which is feasible with clinically used implantation forces of 4 kN.     

Sensitivity Analysis and Optimal Shape Design for Bone-Prosthesis Interfaces in a Femoral Head Surface Replacement

A numerical optimization procedure has been applied for the shape optimal design of a femoral head surface replacement. The failure modes of the prosthesis that were considered in the formulation of the objective functions concerned the interface stress magnitude and the bone remodelling activity beneath the implant. In order to find a compromising solution between different requirements demanded by the two objective functions, a two step optimization procedure has been developed. Through step 1 the minimization of interface stress was achieved, through step 2 the minimization of bone remodelling was achieved with constraints on interface stresses. The results obtained provided an optimal design that generates limited bone remodelling activity with controlled interface stress distribution. The computational procedure was based on the application of the finite element method, linked to a mathematical programming package and a design sensitivity analysis package.     

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.