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.     

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.     

Trabecular bone remodeling around smooth and porous implants in an equine patellar model

The objective of this investigation was to examine the stress-morphology relationships for trabecular bone around implants with different surface characteristics. Stainless steel spheres with either a polished surface or a sintered-bead porous coating were implanted unilaterally into equine patellae and maintained for a 6 month period. Stereological methods were used to quantify the trabecular bone morphology and finite element analyses were performed to predict the trabecular bone stresses. In general, the remodeling response around the smooth implants was greater than that around those porous implants that exhibited bone ingrowth. In accordance with these differences, the finite element models predicted greater changes in the stresses adjacent to the smooth implants due to the nonlinear boundary conditions. However, it did not appear that the trajectorial theory, in its simplest form, was applicable to the remodeling induced by the implants. A linear relationship between the change in bone areal density and the change in von Mises effective stress provides support for the hypothesis that the architecture of trabecular bone corresponds to an optimal structure. The results also demonstrated that, under certain circumstances, small changes in the stress state may result in large changes in the principal material orientation.     

Multi-factorial analysis of variables influencing the bone loss of an implant placed in the maxilla: Prediction using FEA and SED bone remodeling algorithm

The aim of this study was to investigate the interactions of implant position, implant–abutment connection and loading condition influencing bone loss of an implant placed in the maxilla using finite element (FE) analysis and mathematical bone remodeling theory. The maxilla section contours were acquired using CT images to construct FE models containing RS (internal retaining-screw) and the TIS (taper integrated screwed-in) implants placed in SC (along the axis of occlusal force) and RA (along the axis of residual ridge) positions. The adaptive strain energy density (SED) algorithm was combined with FE approach to study the preliminary bone remodeling around implant systems under different load conditions. The simulated results showed that the implant position obviously influenced the bone loss. An implant placed in the RA position resulted in substantially increased bone loss. Implant receiving a lateral load slightly increased bone loss compared with an axial load. The implant type did not significantly influence bone loss. It was found that buccal site suffered the most bone loss around the implant, followed by distal, lingual and mesial sites. The implant position primarily influenced bone loss and it was found most obviously at the buccal site. Implant placed along the axial load direction of a proposed prosthesis could obtain less bone loss around the implant. Attaining proper occlusal adjustments to reduce the lateral occlusal force is recommended in implant–bone–prosthesis system. Abutments of internal engagement with or without taper-fit did not affect the bone loss in the surrounding bone.     

A repeated sampling bone chamber methodology for the evaluation of tissue differentiation and bone adaptation around titanium implants under controlled mechanical conditions

A repeated sampling bone chamber methodology was developed for the study of the influence of the mechanical environment on skeletal tissue differentiation and bone adaptation around titanium implants. Via perforations, bone grows into the implanted outer bone chamber, containing an inner bone chamber with a central test implant. An actuator—easily mounted on the outer bone chamber—allows a controlled mechanical stimulation of the test implant. After each experiment, the inner bone chamber—with its content—can be harvested and analysed. A new inner bone chamber with a central implant can be inserted consecutively in the outer bone chamber and a new experiment can start. Pilot studies led to a reliable surgical protocol and showed the applicability of the methodology, offering the possibility to study skeletal tissue differentiation and adaptation around implants under well-controlled mechanical conditions, and this protected from external loading. Repeated sampling of the bone chamber allows conducting several experiments within the same animal at the same site, thereby excluding subject- and site-dependent variability and reducing the amount of experimental animals.     

Prediction of bone density around orthopedic implants delivering bisphosphonate

The fixation of an orthopedic implant depends strongly upon its initial stability. Peri-implant bone may resorb shortly after the surgery. This resorption is directly followed by new bone formation and implants fixation strengthening, the so-called secondary fixation. If the initial stability is not reached, the resorption continues and the implant fixation weakens, which leads to implant loosening. Studies with rats and dogs have shown that a solution to prevent peri-implant resorption is to deliver bisphosphonate from the implant surface.The aims of the study were, first, to develop a model of bone remodeling around an implant delivering bisphosphonate, second, to predict the bisphosphonate dose that would induce the maximal peri-implant bone density, and third to verify in vivo that peri-implant bone density is maximal with the calculated dose.The model consists of a bone remodeling equation and a drug diffusion equation. The change in bone density is driven by a mechanical stimulus and a drug stimulus. The drug stimulus function and the other numerical parameters were identified from experimental data. The model predicted that a dose of 0.3 μg of zoledronate on the implant would induce a maximal bone density. Implants with 0.3 μg of zoledronate were then implanted in rat femurs for 3, 6 and 9 weeks. We measured that peri-implant bone density was 4% greater with the calculated dose compared to the dose empirically described as best.The approach presented in this paper could be used in the design and analysis processes of experiments in local delivery of drug such as bisphosphonate.     

Mandibular bone remodeling induced by dental implant

The ability to assess the effects of an implant on bone remodeling is of particular importance to prosthesis placement planning and associated treatment assurance. Prediction of on-going bone responses will enable us to improve the performance of a restoration. Although the bone remodeling for long bones had been extensively studied, there have been relatively few reports for dental scenarios despite its increasing significance with more and more dental implant placements. This paper aimed to develop a systematic protocol to assess mandibular bone remodeling induced by dental implantation, which extends the remodeling algorithms established for the long bones into dental settings. In this study, a 3D model for a segment of a human mandible was generated from in vivo CT scan images, together with a titanium implant embedded to the mandible. The results examined the changes in bone density and stiffness as a result of bone remodeling over a period of 48 months. Resonance frequency analysis was also performed to relate natural frequencies to bone remodeling. The density contours are qualitatively compared with clinical follow-up X-ray images, thereby providing validity for the bone remodeling algorithm presented in dental bone analysis.     

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.     

Influence of gait loads on implant integration in rat tibiae: Experimental and numerical analysis

Implanted rat bones play a key role in studies involving fracture healing, bone diseases or drugs delivery among other themes. In most of these studies the implants integration also depends on the animal daily activity and musculoskeletal loads, which affect the implants mechanical environment. However, the tissue adaption to the physiological loads is often filtered through control groups or not inspected. This work aims to investigate experimentally and numerically the effects of the daily activity on the integration of implants inserted in the rat tibia, and to establish a physiological loading condition to analyse the peri-implant bone stresses during gait. Two titanium implants, single and double cortex crossing, are inserted in the rat tibia. The animals are caged under standard conditions and divided in three groups undergoing progressive integration periods. The results highlight a time-dependent increase of bone samples with significant cortical bone loss. The phenomenon is analysed through specimen-specific Finite Element models involving purpose-built musculoskeletal loads. Different boundary conditions replicating the post-surgery bone–implant interaction are adopted. The effects of the gait loads on the implants integration are quantified and agree with the results of the experiments. The observed cortical bone loss can be considered as a transient state of integration due to bone disuse atrophy, initially triggered by a loss of bone–implant adhesion and subsequently by a cyclic opening of the interface.     

Peri-implant bone remodeling after total hip replacement combined with systemic alendronate treatment: a finite element analysis

In order to decrease the peri-implant bone loss during the life-time of the implant, oral use of anti-osteoporosis drugs (like bisphosphonates) has been suggested. In this study, bone remodeling parameters identified from clinical trials of alendronate were used to simulate the effect of those drugs used after total hip arthroplasty on the peri-implant bone density. Results of the simulation show that the oral administrated drugs increase bone density around the implant and decreases, at the same time, the micromovements between the implant and the surrounding bone tissue. Incorporation of drug effect in numerical studies of bone remodeling is a promising tool especially to predetermine safe bisphosphonate doses that could be used with orthopedic implants.     

Stress-shielding induced bone remodeling in cementless shoulder resurfacing arthroplasty: a finite element analysis and in vivo results

Cementless surface replacement arthroplasty (CSRA) of the shoulder was designed to preserve the individual anatomy and humeral bone stock. A matter of concern in resurfacing implants remains the stress shielding and bone remodeling processes. The bone remodeling processes of two different CSRA fixation designs, conical-crown (Epoca RH) and central-stem (Copeland), were studied by three-dimensional (3-D) finite element analysis (FEA) as well as evaluation of contact radiographs from human CSRA retrievals. FEA included one native humerus model with a normal and one with a reduced bone stock quality. Compressive strains were evaluated before and after virtual CSRA implantation and the results were then compared to the bone remodeling and stress-shielding pattern of eight human CSRA retrievals (Epoca RH n=4 and Copeland n=4). FEA revealed for both bone stock models increased compressive strains at the stem and outer implant rim for both CSRA designs indicating an increased bone formation at those locations. Unloading of the bone was seen for both designs under the central implant shell (conical-crown 50–85%, central-stem 31–93%) indicating high bone resorption. Those effects appeared more pronounced for the reduced than for the normal bone stock model. The assumptions of the FEA were confirmed in the CSRA retrieval analysis which showed bone apposition at the outer implant rim and stems with highly reduced bone stock below the central implant shell. Overall, clear signs of stress shielding were observed for both CSRAs designs in the in vitro FEA and human retrieval analysis. Especially in the central part of both implant designs the bone stock was highly resorbed. The impact of these bone remodeling processes on the clinical outcome as well as long-term stability requires further evaluation.     

Numerical analysis of an osseointegrated prosthesis fixation with reduced bone failure risk and periprosthetic bone loss

Currently available implants for direct attachment of prosthesis to the skeletal system after transfemoral amputation (OPRA system, Integrum AB, Sweden and ISP Endo/Exo prosthesis, ESKA Implants AG, Germany) show many advantages over the conventional socket fixation. However, restraining biomechanical issues such as considerable bone loss around the stem and peri-prosthetic bone fractures are present. To overcome these limiting issues a new concept of the direct intramedullary fixation was developed. We hypothesize that the new design will reduce the peri-prosthetic bone failure risk and adverse bone remodeling by restoring the natural load transfer in the femur. Generic CT-based finite element models of an intact femur and amputated bones implanted with 3 analyzed implants were created and loaded with a normal walking and a forward fall load. The strain adaptive bone remodeling theory was used to predict long-term bone changes around the implants and the periprosthetic bone failure risk was evaluated by the von Mises stress criterion. The results show that the new design provides close to physiological distribution of stresses in the bone and lower bone failure risk for the normal walking as compared to the OPRA and the ISP implants. The bone remodeling simulations did not reveal any overall bone loss around the new design, as opposed to the OPRA and the ISP implants, which induce considerable bone loss in the distal end of the femur. This positive outcome shows that the presented concept has a potential to considerably improve safety of the rehabilitation with the direct fixation implants.     

Importance of diameter-to-length ratio in selecting dental implants: a methodological finite element study

Implant dimensions greatly influence load transfer characteristics and the lifetime of a dental system. Excessive stresses at peri-implant area may result in bone failure. Finding the critical point at the implant–bone interface and evaluating the influence of implant diameter-to-length ratio on adjacent bone stresses makes it possible to select implant dimensions. For this, different cylindrical implants were numerically analysed using geometrical models generated from computed tomography images of mandible with osseointegrated implants. All materials were assumed to be linearly elastic and isotropic. Masticatory load was applied in its natural direction, oblique to occlusal plane. Maximum von Mises stresses were located around the implant neck at the critical point of its intersection with the plane of loading and were functions of implant diameter-to-length ratio. It was demonstrated that there exists a certain spectrum of diameter-to-length ratios, which will keep maximum bone stresses at a preset level chosen in accordance with patient's bone strength.     

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.     

Automated three-dimensional finite element modelling of bone: a new method

Three-dimensional finite element stress analysis of bone is a key to understanding bone remodelling, assessing fracture risk, and designing prostheses; however, the cost and complexity of predicting the stress field in bone with accuracy has precluded the routine use of this method. A new, automated method of generating patient-specific three-dimensional finite element models of bone is presented — it uses digital computed tomographic (CT) scan data to derive the geometry of the bone and to estimate its inhomogeneous material properties. Cubic elements of a user-specified size are automatically defined and then individually assigned the CT scan-derived material properties. The method is demonstrated by predicting the stress, strain, and strain energy in a human proximal femur in vivo. Three-dimensional loading conditions corresponding to the stance phase of gait were taken from the literature and applied to the model. Maximum principal compressive stresses of 8–23 MPa were computed for the medial femoral neck. Automated generation of additional finite element models with larger numbers of elements was used to verify convergence in strain energy.     

Simulated bone remodeling around tilted dental implants in the anterior maxilla

Dental implants have to be placed with the long axis in different angulations due to the change in bone morphology. The objective of this study was to investigate the different bone remodeling response induced by the tilted dental implants and to assess whether it could lead to bone loss and implant failure. In this study, bone remodeling due to palato-labially inclined dental implants placed in the anterior maxillary incisor region was simulated. CT-based finite element models of a maxillary bone with dental implants were created herein. Five dental implants were placed at \(+10^{\circ }\), \(+5^{\circ }\), \(0^{\circ }\), \(-5^{\circ }\) and \(-10^{\circ }\), respectively. The remodeling progression was recorded and compared. Model \(-10^{\circ }\) (palatal side) shows the highest bone density values, but the inclined implant at \(+10^{\circ }\) (labial side) leads to significant bone loss. From a biomechanical perspective, it is speculated that a palatally inclined implant is more likely to enhance the bone density in the maxillary anterior region, but labial inclination of implant could jeopardize its stability.     

Influence of the stiffness of bone defect implants on the mechanical conditions at the interface--a finite element analysis with contact

The study focused on the influence of the implant material stiffness on stress distribution and micromotion at the interface of bone defect implants. We hypothesized that a low-stiffness implant with a modulus closer to that of the surrounding trabecular bone would yield a more homogeneous stress distribution and less micromotion at the interface with the bony bed. To prove this hypothesis we generated a three-dimensional, non-linear, anisotropic finite element (FE) model. The FE model corresponded to a previously developed animal model in sheep. A prismatic implant filled a standardized defect in the load-bearing area of the trabecular bone beneath the tibial plateau. The interface was described by face-to-face contact elements, which allow press fits, friction, sliding, and gapping. We assumed a physiological load condition and calculated contact pressures, shear stresses, and shear movements at the interface for two implants of different stiffness (titanium: E=110GPa; composite: E=2.2GPa). The FE model showed that the stress distribution was more homogeneous for the low-stiffness implant. The maximum pressure for the composite implant (2.1 MPa) was lower than for the titanium implant (5.6 MPa). Contrary to our hypothesis, we found more micromotion for the composite (up to 6 microm) than for the titanium implant (up to 4.5 microm). However, for both implants peak stresses and micromotion were in a range that predicts adequate conditions for the osseointegration. This was confirmed by the histological results from the animal studies.     

Stress distribution in cylindrical and conical implants under rotational micromovement with different boundary conditions and bone properties: 3-D FEA

Factors related to micromovements at bone-implant interface have been studied because they are considered adverse to osseointegration. Simplifications are commonly observed in these FEA evaluations. The aim of this study was to clarify the influence of FEA parameters (boundary conditions and bone properties) on the stress distribution in peri-implant bone tissue when micromovements are simulated in implants with different geometries. Three-dimensional models of an anterior section of the jaw with cylindrical or conical titanium implants (4.1 mm in width and 11 mm in length) were created. Micromovement (50, 150, or 250 μm) was applied to the implant. The FEA parameters studied were linear vs. non-linear analyses, isotropic vs. orthogonal anisotropic bone, friction coefficient (0.3) vs. frictionless bone-implant contact. Data from von Mises, shear, maximum, and minimum principal stresses in the peri-implant bone tissue were compared. Linear analyses presented a relevant increase of the stress values, regardless of the bone properties. Frictionless contact reduced the stress values in non-linear analysis. Isotropic bone presented lower stress than orthogonal anisotropic. Conical implants behave better, in regard to compressive stresses (minimum principal), than cylindrical ones, except for nonlinear analyses when micromovement of 150 and 250 μm were simulated. The stress values raised as the micromovement amplitude increased. Non-linear analysis, presence of frictional contact and orthogonal anisotropic bone, evaluated through maximum and minimum principal stress should be used as FEA parameters for implant-micromovement studies.     

Peri-implant Bone Remodeling after Total Hip Replacement Combined with Systemic Alendronate Treatment: A Finite Element Analysis

In order to decrease the peri-implant bone loss during the life-time of the implant, oral use of anti-osteoporosis drugs (like bisphosphonates) has been suggested.

In this study, bone remodeling parameters identified from clinical trials of alendronate were used to simulate the effect of those drugs used after total hip arthroplasty on the peri-implant bone density. Results of the simulation show that the oral administrated drugs increase bone density around the implant and decreases, at the same time, the micromovements between the implant and the surrounding bone tissue.

Incorporation of drug effect in numerical studies of bone remodeling is a promising tool especially to predetermine safe bisphosphonate doses that could be used with orthopedic implants.