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.     

A paradigm for the development and evaluation of novel implant topologies for bone fixation: Implant design and fabrication

The future development of bio-integrated devices will improve the functionality of robotic prosthetic limbs. A critical step in the advancement of bio-integrated prostheses will be establishing long-term, secure fixation to the remnant bone. To overcome limitations associated with contemporary bone-anchored prosthetic limbs, we established a paradigm for developing and fabricating novel orthopedic implants undergoing specified loading. A topology optimization scheme was utilized to generate optimal implant macrostructures that minimize deformations near the bone-implant interface. Variations in implant characteristics and interfacial connectivity were investigated to examine how these variables influence the layout of the optimized implant. For enhanced tissue integration, the optimally designed macroscopic geometry of a titanium (Ti)-alloy implant was further modified by introducing optimized microstructures. The complex geometries of selected implants were successfully fabricated using selective laser sintering (SLS) technology. Fabrication accuracy was assessed by comparing volumes and cross-sectional areas of fabricated implants to CAD data. The error of fabricated volume to CAD design volume was less than 8% and differences in cross sectional areas between SEM images of fabricated implants and corresponding cross sections from CAD design were on average less than 9%. We have demonstrated that this computational design method, combined with solid freeform fabrication techniques, provides a versatile way to develop novel orthopedic implants.     

Prognosis of implant longevity in terms of annual bone loss: a methodological finite element study

Dental implant failure is mainly the consequence of bone loss at peri-implant area. It usually begins in crestal bone. Due to this gradual loss, implants cannot withstand functional force without bone overload, which promotes complementary loss. As a result, implant lifetime is significantly decreased. To estimate implant success prognosis, taking into account 0.2 mm annual bone loss for successful implantation, ultimate occlusal forces for the range of commercial cylindrical implants were determined and changes of the force value for each implant due to gradual bone loss were studied. For this purpose, finite element method was applied and von Mises stresses in implant–bone interface under 118.2 N functional occlusal load were calculated. Geometrical models of mandible segment, which corresponded to Type II bone (Lekholm & Zarb classification), were generated from computed tomography images. The models were analyzed both for completely and partially osseointegrated implants (bone loss simulation). The ultimate value of occlusal load, which generated 100 MPa von Mises stresses in the critical point of adjacent bone, was calculated for each implant. To estimate longevity of implants, ultimate occlusal loads were correlated with an experimentally measured 275 N occlusal load (Mericske-Stern & Zarb). These findings generally provide prediction of dental implants success.     

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.     

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.     

Computational simulation of dental implant osseointegration through resonance frequency analysis

The aim of this study is to predict the evolution of the resonance frequency of the bone-implant interface in a dental implant by means of finite element simulation. A phenomenological interface model able to simulate the mechanical effects of the osseointegration process at the bone-implant interface is applied and compared with some experimental results in rabbits. An early stage of slow bone ingrowth, followed by a faster osseointegration phase until final stability is predicted by the simulations. The evolution of the resonance frequency of the implant and surrounding tissues along the simulation period was also obtained, observing a 3-fold increase in the first principal frequency. These findings are in quantitative agreement with the experimental measurements and suggest that the model can be useful to evaluate the influence of mechanical factors such as implant geometry or implant loading on the indirect evaluation of the process of implant osseintegration.     

Micromotion-induced peri-prosthetic fluid flow around a cementless femoral stem

Micromotion-induced interstitial fluid flow at the bone-implant interface has been proposed to play an important role in aseptic loosening of cementless implants. High fluid velocities are thought to promote aseptic loosening through activation of osteoclasts, shear stress induced control of mesenchymal stem cells differentiation, or transport of molecules. In this study, our objectives were to characterize and quantify micromotion-induced fluid flow around a cementless femoral stem using finite element modeling. With a 2D model of the bone-implant interface and full-factorial design, we first evaluated the relative influence of material properties, and bone-implant micromotion and gap on fluid velocity. Transverse sections around a femoral stem were built from computed tomography images, while boundary conditions were obtained from experimental measurements on the same femur. In a second step, a 3D model was built from the same data-set to estimate the shear stress experienced by cells hosted in the peri-implant tissues. The full-factorial design analysis showed that local micromotion had the most influence on peak fluid velocity at the interface. Remarkable variations in fluid velocity were observed in the macrostructures at the surface of the implant in the 2D transverse sections of the stem. The 3D model predicted peak fluid velocities extending up to 2.2 mm/s in the granulation tissue and to 3.9 mm/s in the trabecular bone. Peak shear stresses on the cells hosted in these tissues ranged from 0.1 to 12.5 Pa. These results offer insight into mechanical stimuli encountered at the bone-implant interface.     

A bone preservation protocol that enables evaluation of osseointegration of implants with micro- and nanotextured surfaces

Development of surface treatments has enabled secure attachment of dental implants in less than 1 month. Consequently, it is necessary to characterize accurately the osseointegration of the implant surface in the region of the bone-implant contact (BIC). We developed a method for sample preparation that preserves both bone and BIC to permit analysis of the contact interface. We prepared eight nanotextured implants and implanted them in rabbit tibias. After healing for 30 days, outcomes were analyzed using both our bone preservation protocol and routine decalcification followed by preparation of histological sections stained by hematoxylin and eosin (H & E). Pull-out tests for implant osseointegration were performed after healing. Non-implanted samples of rabbit mandible were used as a control for assessing organic and mineralized bone characteristics and bone structure. Our bone preservation protocol enabled evaluation of many of the same bone characteristics as histological sections stained with H & E. Our protocol enables analysis of implant samples, implant surfaces and osseointegration without risk of BIC damage.     

A histological and biomechanical study of bone stress and bone remodeling around immediately loaded implants

Immediate loading(IL)increases the risk of marginal bone loss.The present study investigated the biomechanical response of peri-implant bone in rabbits after IL,aiming at optimizing load management.Ninety-six implants were installed bilaterally into femurs of 48 rabbits.Test implants on the left side created the maximal initial stress of 6.9 and 13.4 MPa in peri-implant bone and unloaded implants on the contralateral side were controls.Bone morphology and bone-implant interface strength were measured with histological examination and push-out testing during a 12-week observation period.Additionally,the animal data were incorporated into finite element(FE)models to calculate the bone stress distribution at different levels of osseointegration.Results showed that the stress was concentrated in the bone margin and the bone stress gradually decreased as osseointegration proceeded.A stress of about 2.0 MPa in peri-implant bone had a positive effect on new bone formation,osseointegration and bone-implant interface strength.Bone loss was observed in some specimens with stress exceeding 4.0 MPa.Data indicate that IL significantly increases bone stress during the early postoperative period,but the load-bearing capacity of peri-implant bone increases rapidly with an increase of bone-implant contact.Favorable bone responses may be continually promoted when the stress in peri-implant bone is maintained at a definite level.Accordingly,the progressive loading mode is recommended for IL implants.     

A finite element comparison between the mechanical behaviour of rigid and resilient oral implants with respect to immediate loading

In this paper, a qualitative comparison between two types of dental implants with respect to their behaviour under immediate loading is presented. This analysis has been carried out using the finite element method. Since micromotions (and not the load) are responsible of the appearance of a fibrous interface avoiding osseointegration, the relative displacement between the bone surface and the implant has been the main variable analyzed at different loading states and for the two implant types here considered. The implants analyzed differ in their mechanical behavior: rigid or resilient. Their main difference lies in the joining between the different pieces that make up the dental system. While in the rigid implant all the pieces are screwed, in the resilient implant a relative displacement between the pieces is allowed, with the additional introduction of a silicone gasket that acts like the periodontal ligament. Both implants were considered with a similar geometry and under two different loading scenarios, one equivalent to the force of chewing applied to a molar and another which corresponds to a premolar. For the resilient implant, a hyperelastic behaviour for the silicone and contact conditions between the different mobile parts of the implant are considered. The displacements of the emerging-body in both designs are also compared with the values obtained by several authors. However, the results show that both implants fulfill the constraint of the immediate loading protocol. The micromotions of the resilient implant are lower to those of the rigid one, favouring therefore a good osseointegration process while keeping the stresses in the implant under admissible maximum values.     

Modeling the debonding process of osseointegrated implants due to coupled adhesion and friction

Cementless implants have become widely used for total hip replacement surgery. The long-term stability of these implants is achieved by bone growing around and into the rough surface of the implant, a process called osseointegration. However, debonding of the bone–implant interface can still occur due to aseptic implant loosening and insufficient osseointegration, which may have dramatic consequences. The aim of this work is to describe a new 3D finite element frictional contact formulation for the debonding of partially osseointegrated implants. The contact model is based on a modified Coulomb friction law by Immel et al. (2020), that takes into account the tangential debonding of the bone-implant interface. This model is extended in the direction normal to the bone-implant interface by considering a cohesive zone model, to account for adhesion phenomena in the normal direction and for adhesive friction of partially bonded interfaces. The model is applied to simulate the debonding of an acetabular cup implant. The influence of partial osseointegration and adhesive effects on the long-term stability of the implant is assessed. The influence of different patient- and implant-specific parameters such as the friction coefficient \(\mu _\text {b}\), the trabecular Young’s modulus \(E_\text {b}\), and the interference fit \(I\!F\) is also analyzed, in order to determine the optimal stability for different configurations. Furthermore, this work provides guidelines for future experimental and computational studies that are necessary for further parameter calibration.

     

Bone-implant interface shear modulus and ultimate stress in a transcortical rabbit model of open-pore Ti6Al4V implants

This experimental study on laser-textured implants aimed to evaluate periimplant bone elasticity and ultimate stress of the bone-implant interface in a rabbit femur model. After randomization, two cylindrical Ti6Al4V samples (3.5 mm wide, 5.5 mm long) were transcortically implanted in each femur of 15 female New Zealand White Rabbits. Polished implants had been laser-textured with 100, 200, and 300 microm diameter pores, and another corundum blasted implant was additionally textured with 200 microm pores. Twelve weeks into the experiment, a modified push-out test was performed. The median shear modulus indicating the elasticity of the periimplant bone was 41.12 MPa for the proximal implant location and 25.38 MPa for the distal, without evidence for significant differences between implant types. Taking into account the median ultimate shear stress for 200 microm implants with and without corundum blasting, no significant difference could be demonstrated. However, for blasted 200 microm implants a statistically significant (p<0.025) relative gain in ultimate shear stress of 41% and 17% was proven in comparison with 100 and 300 microm implants, respectively. Non-blasted 200 microm implants reached 48% relative gain in respect of 100 microm samples.     

Modulation of bone ingrowth of rabbit femur titanium implants by in vivo axial micromechanical loading.

Titanium implants commonly used in orthopedics and dentistry integrate into host bone by a complex and coordinated process. Despite increasingly well illustrated molecular healing processes, mechanical modulation of implant bone ingrowth is poorly understood. The objective of the present study was to determine whether micromechanical forces applied axially to titanium implants modulate bone ingrowth surrounding intraosseous titanium implants. We hypothesized that small doses of micromechanical forces delivered daily to the bone-implant interface enhance implant bone ingrowth. Small titanium implants were placed transcortically in the lateral aspect of the proximal femur in 15 New Zealand White rabbits under general anesthesia and allowed to integrate with the surrounding bone for 6 wk. Micromechanical forces at 200 mN and 1 Hz were delivered axially to the right femur implants for 10 min/day over 12 consecutive days, whereas the left femur implants served as controls. The average bone volume 1 mm from mechanically loaded implants (n = 15) was 73 +/- 12%, which was significantly greater than the average bone volume (52 +/- 21%) of the contralateral controls (n = 15) (P < 0.01). The average number of osteoblast-like cells per endocortical bone surface was 55 +/- 8 cells/mm(2) for mechanically loaded implants, which was significantly greater than the contralateral controls (35 +/- 6 cells/mm(2)) (P < 0.01). Dynamic histomorphometry showed a significant increase in mineral apposition rate and bone-formation rate of mechanically stressed implants (3.8 +/- 1.2 microm/day and 2.4 +/- 1.0 microm(3).microm(-2).day(-1), respectively) than contralateral controls (2.2 +/- 0.92 microm/day and 1.2 +/- 0.60 microm(3).microm(-2).day(-1), respectively; P < 0.01). Collectively, these data suggest that micromechanical forces delivered axially on intraosseous titanium implants may have anabolic effects on implant bone ingrowth.     

Influence of custom-made implant designs on the biomechanical performance for the case of immediate post-extraction placement in the maxillary esthetic zone: a finite element analysis

Due to the increasing adoption of immediate implantation strategies and the rapid development of the computer aided design/computer aided manufacturing technology, a therapeutic concept based on patient-specific implant dentistry has recently been reintroduced by many researchers. However, little information is available on the designs of custom-made dental implant systems, especially their biomechanical behavior. The influence of the custom-made implant designs on the biomechanical performance for both an immediate and a delayed loading protocol in the maxillary esthetic zone was evaluated by means of the finite element (FE) method. FE models of three dental implants were considered: a state of the art cylindrical implant and two custom-made implants designed by reverse engineering technology, namely a root-analogue implant and a root-analogue threaded implant. The von Mises stress distributions and micro-motions around the bone-implant interfaces were calculated using ANSYS software. In a comparison of the three implant designs for both loading protocols, a favorable biomechanical performance was observed for the use of root-analogue threaded implant which approximated the geometry of natural anterior tooth and maintained the original long-axis. The results indicated that bone-implant interfacial micro-motion was reduced and a favorable stress distribution after osseointegration was achieved.     

Load transfer along the bone–dental implant interface

In this paper the variation of normal and shear stresses along a path defined on the bone–dental implant interface is investigated. In particular, the effects of implant diameter, collar length and slope, body length, and the effects of four different types of external threads on the interfacial stress distribution are studied. The geometry of the bone is digitized from a CT scan of a mandibular incisor and the surrounding bone. The bone and the implant are assumed to be perfectly bonded. The finite element method with 2D plane strain assumption is used to compute interfacial stresses. Highest continuous interfacial stresses are encountered in the region where the implant collar engages the cortical region, and near the apex of the implant in the subcortical region. Stress concentrations in the interfacial stresses occur near the geometric discontinuities on the implant contour, and jumps in stress values occur where the elastic modulus of the bone transitions between the cortical and trabecular bone values. Among the six contour parameters, the slope and the length of the implant collar, and the implant diameter influence the interfacial stress levels the most, and the effects of changing these parameters are significantly noticed only in the cortical bone (alveolar ridge) area. External threads cause significant stress concentrations in interfacial stresses in otherwise smoothly varying regions. This work shows that the presence of external threads could cause significant variations in both normal and shear stresses along the bone–implant interface, but not reduction in shear stress as previously thought.     

The discrete nature of trabecular bone microarchitecture affects implant stability

Small endosseous implants, such as screws, are important components of modern orthopedics and dentistry. Hence they have to reliably fulfill a variety of requirements, which makes the development of such implants challenging. Finite element analysis is a widely used computational tool used to analyze and optimize implant stability in bone. For these purposes, bone is generally modeled as a continuum material. However, bone failure and bone adaptation processes are occurring at the discrete level of individual trabeculae; hence the assessment of stresses and strains at this level is relevant. Therefore, the aim of the present study was to investigate how peri-implant strain distribution and load transfer between implant and bone are affected by the continuum assumption. We performed a computational study in which cancellous screws were inserted in continuum and discrete models of trabecular bone; axial loading was simulated. We found strong differences in bone-implant stiffness between the discrete and continuum bone model. They depended on bone density and applied boundary conditions. Furthermore, load transfer from the screw to the surrounding bone differed strongly between the continuum and discrete models, especially for low-density bone. Based on our findings we conclude that continuum bone models are of limited use for finite element analysis of peri-implant mechanical loading in trabecular bone when a precise quantification of peri-implant stresses and strains is required. Therefore, for the assessment and improvement of trabecular bone implants, finite element models which accurately represent trabecular microarchitecture should be used.     

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.     

Multiscale design and multiobjective optimization of orthopedic hip implants with functionally graded cellular material

Revision surgeries of total hip arthroplasty are often caused by a deficient structural compatibility of the implant. Two main culprits, among others, are bone-implant interface instability and bone resorption. To address these issues, in this paper we propose a novel type of implant, which, in contrast to current hip replacement implants made of either a fully solid or a foam material, consists of a lattice microstructure with nonhomogeneous distribution of material properties. A methodology based on multiscale mechanics and design optimization is introduced to synthesize a graded cellular implant that can minimize concurrently bone resorption and implant interface failure. The procedure is applied to the design of a 2D left implanted femur with optimized gradients of relative density. To assess the manufacturability of the graded cellular microstructure, a proof-of-concept is fabricated by using rapid prototyping. The results from the analysis are used to compare the optimized cellular implant with a fully dense titanium implant and a homogeneous foam implant with a relative density of 50%. The bone resorption and the maximum value of interface stress of the cellular implant are found to be over 70% and 50% less than the titanium implant while being 53% and 65% less than the foam implant.     

Mechanisms of reduced implant stability in osteoporotic bone

The determining factors for the fixation of uncemented screws in bone are the bone-implant interface and the peri-implant bone. The goal of this work was to explore the role of the peri-implant bone architecture on the mechanics of the bone-implant system. In particular, the specific aims of the study were to investigate: (i) the impact of the different architectural parameters, (ii) the effects of disorder, and (iii) the deformations in the peri-implant region. A three-dimensional beam lattice model to describe trabecular bone was developed. Various microstructural features of the lattice were varied in a systematic way. Implant pull-out tests were simulated, and the stiffness and strength of the bone-implant system were computed. The results indicated that the strongest decrease in pull-out strength was obtained by trabecular thinning, whereas pull-out stiffness was mostly affected by trabecular removal. These findings could be explained by investigating the peri-implant deformation field. For small implant displacements, a large amount of trabeculae in the peri-implant region were involved in the load transfer from implant to bone. Therefore, trabecular removal in this region had a strong negative effect on pull-out stiffness. Conversely, at higher displacements, deformations mainly localized in the trabeculae in contact with the implant; hence, thinning those trabeculae produced the strongest decrease in the strength of the system. Although idealized, the current approach is helpful for a mechanical understanding of the role played by peri-implant bone.     

The influence of mechanical loading on osseointegration: an animal study

Osseointegration of implant provides a stable support for the prosthesis under functional loads. The timing of loading is a critical parameter that can govern the success of the osseointegration of implant. However, it is not clear whether the early loading can affect the success of osseointegration, or whether the no-loading healing period can be shortened. This paper presents an animal study conducted to investigate how external loads influence the osseointegration at the initial stage of healing. Titanium implants were inserted into the goat tibia laterally, and different axial loadings were applied to the implants in 4 weeks after surgery. After the 2 weeks period of early loading, animals were sacrificed and the tibia bones with the implants were cut off from the bodies. Then mechanical test was employed to find out the differences in the pull-out force, and shear strength at the bone-implant interface between the non-loaded and the loaded implants. The implant-bone interfaces were analyzed by histomorphometric method, SEM (scanning electron micrograph) and EDS (energy density spectrum). The results indicated that the bone-implant interface did not well integrate 4 weeks after surgery, and the fibrous tissue could be found at the interfaces of the specimens without loadings. While the results of loaded specimens with 10 N axial force showed that that parts of the interface were well integrated, indicating that the early mild loading may play a positive role in the process of the osseointegration. The results support that a certain range of external loading would influence the process of osseointegration, and appropriate mechanical loading can be applied to shorten the osseointegration period after surgery. Supported by the National Natural Science Foundation of China (Grant Nos. 30370376, 10529202 and 10672015).