A method of quantification of stress shielding in the proximal femur using hierarchical computational modeling

Stress shielding is a biomechanical phenomenon causing adaptive changes in bone strength and stiffness around metallic implants, which potentially lead to implant loosening. Accordingly, there is a need for standard, objective engineering measures of the “stress shielding” performances of an implant that can be employed in the process of computer-aided implant design. To provide and test such measures, we developed hierarchical computational models of adaptation of the trabecular microarchitecture at different sites in the proximal femur, in response to insertion of orthopaedic screws and in response to hypothetical reductions in hip joint and gluteal muscle forces. By identifying similar bone adaptation outcomes from the two scenarios, we were able to quantify the stress shielding caused by screws in terms of analogous hypothetical reductions in hip joint and gluteal muscle forces. Specifically, we developed planar lattice models of trabecular microstructures at five regions of interest (ROI) in the proximal femur. The homeostatic and abnormal loading conditions for the lattices were determined from a finite element model of the femur at the continuum scale and fed to an iterative algorithm simulating the adaptation of each lattice to these loads. When screws were inserted to the femur model, maximal simulated bone loss (17% decrease in apparent density, 10% decrease in thickness of trabeculae) was at the greater trochanter and this effect was equivalent to the effect of 50% reduction in gluteal force and normal hip joint force. We conclude that stress shielding performances can be quantified for different screw designs using model-predicted hypothetical musculoskeletal load fractions that would cause a similar pattern and extent of bone loss to that caused by the implants.     

A method of quantification of stress shielding in the proximal femur using hierarchical computational modeling

Stress shielding is a biomechanical phenomenon causing adaptive changes in bone strength and stiffness around metallic implants, which potentially lead to implant loosening. Accordingly, there is a need for standard, objective engineering measures of the "stress shielding" performances of an implant that can be employed in the process of computer-aided implant design. To provide and test such measures, we developed hierarchical computational models of adaptation of the trabecular microarchitecture at different sites in the proximal femur, in response to insertion of orthopaedic screws and in response to hypothetical reductions in hip joint and gluteal muscle forces. By identifying similar bone adaptation outcomes from the two scenarios, we were able to quantify the stress shielding caused by screws in terms of analogous hypothetical reductions in hip joint and gluteal muscle forces. Specifically, we developed planar lattice models of trabecular microstructures at five regions of interest (ROI) in the proximal femur. The homeostatic and abnormal loading conditions for the lattices were determined from a finite element model of the femur at the continuum scale and fed to an iterative algorithm simulating the adaptation of each lattice to these loads. When screws were inserted to the femur model, maximal simulated bone loss (17% decrease in apparent density, 10% decrease in thickness of trabeculae) was at the greater trochanter and this effect was equivalent to the effect of 50% reduction in gluteal force and normal hip joint force. We conclude that stress shielding performances can be quantified for different screw designs using model-predicted hypothetical musculoskeletal load fractions that would cause a similar pattern and extent of bone loss to that caused by the implants.     

Finite element analysis of peri-implant bone volume affected by stresses around Morse taper implants: effects of implant positioning to the bone crest

Abstract

Objectives: The purpose of the present study was to evaluate the distribution and magnitude of stresses through the bone tissue surrounding Morse taper dental implants at different positioning relative to the bone crest. Materials and Methods: A mandibular bone model was obtained from a computed tomography scan. A three-dimensional (3D) model of Morse taper implant-abutment systems placed at the bone crest (equicrestal) and 2?mm bellow the bone crest (subcrestal) were assessed by finite element analysis (FEA). FEA was carried out on axial and oblique (45°) loading at 150 N relatively to the central axis of the implant. The von Mises stresses were analysed considering magnitude and volume of affected peri-implant bone. Results: On vertical loading, maximum von Mises stresses were recorded at 6-7?MPa for trabecular bone while values ranging from 73 up to 118?MPa were recorded for cortical bone. On oblique loading at the equiquestral or subcrestal positioning, the maximum von Mises stresses ranged from 15 to 21?MPa for trabecular bone while values at 150?MPa were recorded for the cortical bone. On vertical loading, >99.9vol.% cortical bone volume was subjected to a maximum of 2?MPa while von Mises stress values at 15?MPa were recorded for trabecular bone. On oblique loading, >99.9vol.% trabecular bone volume was subjected to maximum stress values at 5?MPa, while von Mises stress values at 35?MPa were recorded for >99.4vol.% cortical bone. Conclusions: Bone volume-based stress analysis revealed that most of the bone volume (>99% by vol) was subjected to significantly lower stress values around Morse taper implants placed at equicrestal or subcrestal positioning. Such analysis is commentary to the ordinary biomechanical assessment of dental implants concerning the stress distribution through peri-implant sites.     

Predicting bone remodeling around tissue- and bone-level dental implants used in reduced bone width

The objective of this study was to predict time-dependent bone remodeling around tissue- and bone-level dental implants used in patients with reduced bone width. The remodeling of bone around titanium tissue-level, and titanium and titanium–zirconium alloy bone-level implants was studied under 100 N oblique load for one month by implementing the Stanford theory into three-dimensional finite element models. Maximum principal stress, minimum principal stress, and strain energy density in peri-implant bone and displacement in x- and y- axes of the implant were evaluated. Maximum and minimum principal stresses around tissue-level implant were higher than bone-level implants and both bone-level implants experienced comparable stresses. Total strain energy density in bone around titanium implants slightly decreased during the first two weeks of loading followed by a recovery, and the titanium–zirconium implant showed minor changes in the axial plane. Total strain energy density changes in the loading and contralateral sides were higher in tissue-level implant than other implants in the cortical bone at the horizontal plane. The displacement values of the implants were almost constant over time. Tissue-level implants were associated with higher stresses than bone-level implants. The time-dependent biomechanical outcome of titanium–zirconium alloy bone-level implant was comparable to the titanium implant.     

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.     

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.     

Trabecular deformations during screw pull-out: a micro-CT study of lapine bone

The mechanical fixation of endosseous implants, such as screws, in trabecular bone is challenging because of the complex porous microstructure. Development of new screw designs to improve fracture fixation, especially in high-porosity osteoporotic bone, requires a profound understanding of how the structural system implant/trabeculae interacts when it is subjected to mechanical load. In this study, pull-out tests of screw implants were performed. Screws were first inserted into the trabecular bone of rabbit femurs and then pulled out from the bone inside a computational tomography scanner. The tests were interrupted at certain load steps to acquire 3D images. The images were then analysed with a digital volume correlation technique to estimate deformation and strain fields inside the bone during the tests. The results indicate that the highest shear strains are concentrated between the inner and outer thread diameter, whereas compressive strains are found at larger distances from the screw. Tensile strains were somewhat smaller. Strain concentrations and the location of trabecular failures provide experimental information that could be used in the development of new screw designs and/or to validate numerical simulations.     

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.     

Simulation of bone strain by orthodontic implants using the finite element method]

Load direction of applied forces, implant geometry and other biomechanical parameters lead to varying reactions in the surrounding bone structure. Three types of endosseous implant measuring 9 mm in length and 3.3 mm in diameter with and without superperiosteal step, and a threaded surface were investigated with the aid of a finite element method using the COSMOS/M 2.5 program. The load on the implant was investigated under vertical, horizontal, and diagonal forces of between 0.01 N and 100 N. Vertical loading of simple implants caused bone deformation of more than 600 mu eps. The application of the superperiosteal step clearly reduced the deformation. The largest deformations under vertical loading were observed in the trabecular bone with all 3 implant geometries. On horizontal loading the deformation shifted from the trabecular to the cortical bone and was particularly marked at the transition between the two. The smallest deformations, less than 300 mu eps, were measured at implants with a superperiosteal step under diagonal loading. The thread did not improve loading capacity. Implants with a superperiosteal step are recommended since they contribute to more rapid healing and strengthening of the bone.     

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.     

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

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

Finite element investigation of implant-supported fixed partial prosthesis in the premaxilla in immediately loaded and osseointegrated states

The aim of this study was to gain insight into the behaviour of the stresses and strains at the bone–implant interface of an implant-supported fixed partial prosthesis (FPP) in the premaxilla under immediate loading and osseointegrated conditions. Finite element models of a four-unit FPP were generated. An extreme condition was simulated, using only two immediately loaded implants in order to derive recommendations for possible clinical application. Straight and 20°-angled abutments and bonded or sliding contact between the bridge and abutment were simulated. In addition, two models were generated with two completely osseointegrated implants. A 150 N load to the prosthesis at a 45° angle to the long axis of each implant was applied. Minor differences were observed in implant displacements, stress and strain distributions of the two abutment designs. However, bone loading exceeded the physiological limits, including a risk of bone atrophy. A considerable decrease in implant displacements and bone loading was observed in the osseointegrated cases. An FPP supported by only two implants cannot be recommended for immediate loading.     

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.     

High-resolution finite element models with tissue strength asymmetry accurately predict failure of trabecular bone

The ability to predict trabecular failure using microstructure-based computational models would greatly facilitate study of trabecular structure–function relations, multiaxial strength, and tissue remodeling. We hypothesized that high-resolution finite element models of trabecular bone that include cortical-like strength asymmetry at the tissue level, could predict apparent level failure of trabecular bone for multiple loading modes. A bilinear constitutive model with asymmetric tissue yield strains in tension and compression was applied to simulate failure in high-resolution finite element models of seven bovine tibial specimens. Tissue modulus was reduced by 95% when tissue principal strains exceeded the tissue yield strains. Linear models were first calibrated for effective tissue modulus against specimen-specific experimental measures of apparent modulus, producing effective tissue moduli of (mean±S.D.) 18.7±3.4 GPa. Next, a parameter study was performed on a single specimen to estimate the tissue level tensile and compressive yield strains. These values, 0.60% strain in tension and 1.01% strain in compression, were then used in non-linear analyses of all seven specimens to predict failure for apparent tensile, compressive, and shear loading. When compared to apparent yield properties previously measured for the same type of bone, the model predictions of both the stresses and strains at failure were not statistically different for any loading case (p>0.15). Use of symmetric tissue strengths could not match the experimental data. These findings establish that, once effective tissue modulus is calibrated and uniform but asymmetric tissue failure strains are used, the resulting models can capture the apparent strength behavior to an outstanding level of accuracy. As such, these computational models have reached a level of fidelity that qualifies them as surrogates for destructive mechanical testing of real specimens.     

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.     

The influence of implant number and abutment design on the biomechanical behaviour of bone for an implant-supported fixed prosthesis: a finite element study in the upper anterior region

It is always recommended to use more implants for supporting a prosthesis in the immediate loading condition than in the classical two-stage treatment procedure. By means of the finite element (FE) method, the influence of the number of implants used in immediately loaded fixed partial prosthesis (FPP) on the load distribution was investigated, considering the abutment geometry. Two 3D FE models were studied employing four implants to support a FPP in the premaxilla. One model was designed with straight abutments and the other with 20°-angled abutments. The results concerning implant displacements, stresses and strains were compared with those of two implant-supported FPPs, obtained in a previous study. A noticeable reduction in the determined biomechanical bone loading was observed with the use of more implants in supporting an immediately loaded prosthesis. This study confirms that the use of additional numbers of implants in an immediately loaded prosthesis is highly recommended.     

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.     

The influence of implant number and abutment design on the biomechanical behaviour of bone for an implant-supported fixed prosthesis: a finite element study in the upper anterior region

It is always recommended to use more implants for supporting a prosthesis in the immediate loading condition than in the classical two-stage treatment procedure. By means of the finite element (FE) method, the influence of the number of implants used in immediately loaded fixed partial prosthesis (FPP) on the load distribution was investigated, considering the abutment geometry. Two 3D FE models were studied employing four implants to support a FPP in the premaxilla. One model was designed with straight abutments and the other with 20°-angled abutments. The results concerning implant displacements, stresses and strains were compared with those of two implant-supported FPPs, obtained in a previous study. A noticeable reduction in the determined biomechanical bone loading was observed with the use of more implants in supporting an immediately loaded prosthesis. This study confirms that the use of additional numbers of implants in an immediately loaded prosthesis is highly recommended.     

Influence of bicortical techniques in internal connection placed in premaxillary area by 3D finite element analysis

The aim of study was to evaluate the stress distribution in implant-supported prostheses and peri-implant bone using internal hexagon (IH) implants in the premaxillary area, varying surgical techniques (conventional, bicortical and bicortical in association with nasal floor elevation), and loading directions (0°, 30° and 60°) by three-dimensional (3D) finite element analysis. Three models were designed with Invesalius, Rhinoceros 3D and Solidworks software. Each model contained a bone block of the premaxillary area including an implant (IH, Ø4 × 10 mm) supporting a metal-ceramic crown. 178 N was applied in different inclinations (0°, 30°, 60°). The results were analyzed by von Mises, maximum principal stress, microstrain and displacement maps including ANOVA statistical test for some situations. Von Mises maps of implant, screws and abutment showed increase of stress concentration as increased loading inclination. Bicortical techniques showed reduction in implant apical area and in the head of fixation screws. Bicortical techniques showed slight increase stress in cortical bone in the maximum principal stress and microstrain maps under 60° loading. No differences in bone tissue regarding surgical techniques were observed. As conclusion, non-axial loads increased stress concentration in all maps. Bicortical techniques showed lower stress for implant and screw; however, there was slightly higher stress on cortical bone only under loads of higher inclinations (60°).     

Incorporating tissue anisotropy and heterogeneity in finite element models of trabecular bone altered predicted local stress distributions

Trabecular bone is composed of organized mineralized collagen fibrils, which results in heterogeneous and anisotropic mechanical properties at the tissue level. Recently, biomechanical models computing stresses and strains in trabecular bone have indicated a significant effect of tissue heterogeneity on predicted stresses and strains. However, the effect of the tissue-level mechanical anisotropy on the trabecular bone biomechanical response is unknown. Here, a computational method was established to automatically impose physiologically relevant orientation inherent in trabecular bone tissue on a trabecular bone microscale finite element model. Spatially varying tissue-level anisotropic elastic properties were then applied according to the bone mineral density and the local tissue orientation. The model was used to test the hypothesis that anisotropy in both homogeneous and heterogeneous models alters the predicted distribution of stress invariants. Linear elastic finite element computations were performed on a 3 mm cube model isolated from a microcomputed tomography scan of human trabecular bone from the distal femur. Hydrostatic stress and von Mises equivalent stress were recorded at every element, and the distributions of these values were analyzed. Anisotropy reduced the range of hydrostatic stress in both tension and compression more strongly than the associated increase in von Mises equivalent stress. The effect of anisotropy was independent of the spatial redistribution high compressive stresses due to tissue elastic heterogeneity. Tissue anisotropy and heterogeneity are likely important mechanisms to protect bone from failure and should be included for stress analyses in trabecular bone.