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.     

Stress analysis of a centrally fractured rib fixated by an intramedullary screw

The stress on an intramedullary screw rib fixation device holding together a centrally fractured human rib under in vivo force loadings was studied using finite element analysis (FEA). Validation of the FEA modelling using pullout from porcine ribs proved FEA to be suitable for assessing the structural integrity of screw/bone systems such as rib fixated by a screw. In the human rib fixation investigation, it was found that intramedullary bioresorbable Bioretec screws can fixate centrally fractured human ribs under normal breathing conditions. However, under coughing conditions, simulation showed Bioretec fixating screws to bend substantially. High stresses in the screw are mainly the result of flexion induced by the force loading, and are restricted to thin regions on the outside of the screw shaft. Stiffer screws result in less locally intense stress concentrations in bone, indicating that bone failure in the bone/screw contact regions can be averted with improvements in screw stiffness.     

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 customized fixation plate with novel structure designed by topological optimization for mandibular angle fracture based on finite element analysis

Background

The purpose of this study was to design a customized fixation plate for mandibular angle fracture using topological optimization based on the biomechanical properties of the two conventional fixation systems, and compare the results of stress, strain and displacement distributions calculated by finite element analysis (FEA).

Methods

A three-dimensional (3D) virtual mandible was reconstructed from CT images with a mimic angle fracture and a 1 mm gap between two bone segments, and then a FEA model, including volume mesh with inhomogeneous bone material properties, three loading conditions and constraints (muscles and condyles), was created to design a customized plate using topological optimization method, then the shape of the plate was referenced from the stress concentrated area on an initial part created from thickened bone surface for optimal calculation, and then the plate was formulated as “V” pattern according to dimensions of standard mini-plate finally. To compare the biomechanical behavior of the “V” plate and other conventional mini-plates for angle fracture fixation, two conventional fixation systems were used: type A, one standard mini-plate, and type B, two standard mini-plates, and the stress, strain and displacement distributions within the three fixation systems were compared and discussed.

Results

The stress, strain and displacement distributions to the angle fractured mandible with three different fixation modalities were collected, respectively, and the maximum stress for each model emerged at the mandibular ramus or screw holes. Under the same loading conditions, the maximum stress on the customized fixation system decreased 74.3, 75.6 and 70.6% compared to type A, and 34.9, 34.1, and 39.6% compared to type B. All maximum von Mises stresses of mandible were well below the allowable stress of human bone, as well as maximum principal strain. And the displacement diagram of bony segments indicated the effect of treatment with different fixation systems.

Conclusions

The customized fixation system with topological optimized structure has good biomechanical behavior for mandibular angle fracture because the stress, strain and displacement within the plate could be reduced significantly comparing to conventional “one mini-plate” or “two mini-plates” systems. The design methodology for customized fixation system could be used for other fractures in mandible or other bones to acquire better mechanical behavior of the system and improve stable environment for bone healing. And together with SLM, the customized plate with optimal structure could be designed and fabricated rapidly to satisfy the urgent time requirements for treatment.

     

Trabecular shear stress in human vertebral cancellous bone: intra- and inter-individual variations

Correlation of the mean and standard deviation of trabecular stresses has been proposed as a mechanism by which a strong relationship between the apparent strength and stiffness of cancellous bone can be achieved. The current study examined whether the relationship between the mean and standard deviation of trabecular von Mises stresses can be generalized for any group of cancellous bone. Cylindrical human vertebral cancellous bone specimens were cut in the infero-superior direction from T12 of 23 individuals (inter-individual group). Thirty nine additional specimens were prepared similarly from the T4-T12 and L2-L5 vertebrae of a 63 year old male (intra-individual group). The specimens were scanned by micro-computed tomography (microCT) and trabecular von Mises stresses were calculated using finite element modeling. The expected value, standard deviation and coefficient of variation of the von Mises stress were calculated form a three-parameter Weibull function fitted to von Mises stress data from each specimen. It was found that the average and standard deviation of trabecular von Mises shear stress were: (i) correlated with each other, supporting the idea that high correlation between the apparent strength and stiffness of cancellous bone can be achieved through controlling the trabecular level shear stress variations, (ii) dependent on anatomical site and sample group, suggesting that the variation of stresses are correlated to the mean stress to different degrees between vertebrae and individuals, and (iii) dependent on bone volume fraction, consistent with the idea that shear stress is less well controlled in bones with low BV/TV. The conversion of infero-superior loading into trabecular von Mises stresses was maximum for the tissue at the junction of the thoracic and lumbar spine (T12-L1) consistent with this junction being a common site of vertebral fracture.     

Stress changes in the femoral head due to porous ingrowth surface replacement arthroplasty

Finite element stress analyses were conducted of the canine femoral head before and after implantation of various surface replacement-type components. The femoral head was replaced by four implant geometries; (a) shell, (b) shell with peg, (c) shell with rod, and (d) a new epiphyseal replacement design. All implants were modelled to simulate bony ingrowth along the underside of the shell and along the surfaces of the peg and rod. The results indicated that in the normal femur the forces are transferred from the articular surface through the femoral head cancellous bone to the inferior cortical shell of the femoral neck. After shell-type surface replacement, forces were transferred more distally at the rim of the shell and at the end of the peg or rod, thereby reducing the stresses in the proximal head cancellous bone. Computer simulation of bone remodelling due to proximal bone stress reduction was shown to accentuate the abnormality of the stress fields. Surface replacement with a lower modulus material created a less abnormal redistribution of bone stresses. The new epiphyseal replacement design resulted in stress distributions similar to those in the normal femoral head and minimal shear stresses at the implant/bone interface. These findings suggest that the epiphyseal replacement concept may provide better initial mechanical integrity and create a more benign milieu for adaptive bone remodelling than conventional, shell-type surface replacement components.     

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.     

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.     

Influence of cartilage cap modelling on FEM simulation of femoral head stress]

The present study reports on the finite element analysis (FEA) of the femoral head in a process of preparation for a program for the realistic simulation of correctional osteotomies of the proximal femur. While the material properties have been studied extensively, only few publications consider the influence of the cartilage layer geometry on FE stimulation of the hip joint. Various models of the femoral head with and without the cartilage layer are generated and analysed. On looking at the maximum surface stresses, we found a strong influence of the cartilage layer and the subchondral osseous layer on the magnitude of the von Mises equivalent stress. The model with an anatomically realistic cartilage layer and compact bone shows stresses of between 4 and 5.5 MPa, depending on the position of the joint, while the model with a concentric cartilage layer has a maximum von Mises stress of 0.8 MPa. Only on simulation of a "realistic" cartilage layer, with a maximum thickness at the "pole" and minimum thickness at the "equator" do the changes in stress distribution--determined by changes in the position of the femoral head--become visible. Owing to major artefacts and the inability to create a realistic cartilage layer, voxel-based models of the femur are not suitable for the simulation of the femoral head surface.     

A Comparative Analysis of Different Treatments for Distal Femur Fractures using the Finite Element Method

The main objective of this work is the evaluation, by means of the finite element method (FEM) of the mechanical stability and long-term microstructural modifications in bone induced to three different kinds of fractures of the distal femur by three types of implants: the Condyle Plate, the less invasive stabilization system plate (LISS) and the distal femur nail (DFN). The displacement and the stress distributions both in bone and implants and the internal bone remodelling process after fracture and fixation are obtained and analysed by computational simulation. The main conclusions of this work are that distal femoral fractures can be treated correctly with the Condyle Plate, the LISS plate and the DFN. The stresses both in LISS and DFN implant are high especially around the screws. When respect to remodelling, the LISS produces an important resorption in the fractured region, while the other two implants do not strongly modify bone tissue microstructure.     

A comparative analysis of different treatments for distal femur fractures using the finite element method

The main objective of this work is the evaluation, by means of the finite element method (FEM) of the mechanical stability and long-term microstructural modifications in bone induced to three different kinds of fractures of the distal femur by three types of implants: the Condyle Plate, the less invasive stabilization system plate (LISS) and the distal femur nail (DFN). The displacement and the stress distributions both in bone and implants and the internal bone remodelling process after fracture and fixation are obtained and analysed by computational simulation. The main conclusions of this work are that distal femoral fractures can be treated correctly with the Condyle Plate, the LISS plate and the DFN. The stresses both in LISS and DFN implant are high especially around the screws. When respect to remodelling, the LISS produces an important resorption in the fractured region, while the other two implants do not strongly modify bone tissue microstructure.     

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.     

Morphological and numerical studies of the stress compatibility of patella implants]

Failure of total knee arthroplasty is relatively often caused by problems of the patellofemoral replacement. The purpose of this study was to analyze the distribution of stresses within an anatomical patella and the changes in stress distribution after patellar resurfacing with a Miller-Galante I patellar implant using two- and three dimensional finite element models (FEM). To assess validity, FEM results were compared with morphological findings from contact radiographs and densitographs. Internal orientation of bone trabeculae is in good agreement with the orientation of theoretically calculated principal stresses. Almost unchanged principal tensile stresses after implantation, together with the lack of extreme stress peaks within the cancellous bone ensure stress compatibility of the implant. In the case of a firmly seated implant with good bone ingrowth, increased von Mises stresses are found near the fixation peg/plate junction. Their relevance for improved bone ingrowth near this part of the interface is emphasized. At the same time, material failure at the peg/plate junction can be better understood. An analysis of the early postoperative period assuming nonlinear interface conditions failed to demonstrate an uniform distribution of normal and tangential interface forces.     

Age-related changes in human trabecular bone: Relationship between microstructural stress and strain and damage morphology

Accumulation of microdamage in aging and disease can cause skeletal fragility and is one of several factors contributing to osteoporotic fractures. To better understand the role of microdamage in fragility fracture, the mechanisms of bone failure must be elucidated on a tissue-level scale where interactions between bone matrix properties, the local biomechanical environment, and bone architecture are concurrently examined for their contributions to microdamage formation. A technique combining histological damage assessment of individual trabeculae with linear finite element solutions of trabecular von Mises and principal stress and strain was used to compare the damage initiation threshold between pre-menopausal (32-37 years, n=3 donors) and post-menopausal (71-80 years, n=3 donors) femoral cadaveric bone. Strong associations between damage morphology and stress and strain parameters were observed in both groups, and an age-related decrease in undamaged trabecular von Mises stress was detected. In trabeculae from younger donors, the 95% CI for von Mises stress on undamaged regions ranged from 50.7-67.9MPa, whereas in trabeculae from older donors, stresses were significantly lower (38.7-50.2, p<0.01). Local microarchitectural analysis indicated that thinner, rod-like trabeculae oriented along the loading axis are more susceptible to severe microdamage formation in older individuals, while only rod-like architecture was associated with severe damage in younger individuals. This study therefore provides insight into how damage initiation and morphology relate to local trabecular microstructure and the associated stresses and strains under loading. Furthermore, by comparison of samples from pre- and post-menopausal women, the results suggest that trabeculae from younger individuals can sustain higher stresses prior to microdamage initiation.     

Uniform stress state in bone structure with residual stress.

Residual stress and strain in living tissues have been investigated from the viewpoint of mechanical optimality maintained by adaptive remodeling. In this study, the residual stresses in the cortical-cancellous bone complex of bovine coccygeal vertebrae were examined. Biaxial strain gages were bonded onto the cortical surface, so that the gage axes were aligned in the cephalocaudal and circumferential directions. Strains induced by removal of the end-plate and the cancellous bone were recorded sequentially. The results revealed the existence of compressive residual stress in the cortical bone and tensile residual stress in the cancellous bone in both the cephalocaudal and the circumferential direction. The observed strains were examined on the basis of the uniform stress hypothesis using a three-bar model for the cephalocaudal direction and a three-layered cylinder model for the circumferential direction. In this model study, the magnitude of effective stresses, which is defined as the macroscopic stress divided by the area fraction of bone material, was found not to differ significantly between cephalocaudal and circumferential directions, although they were evaluated using independent models. These results demonstrate that the uniform stress state of the cortical-cancellous bone structure is consistent with results obtained in the cutting experiment when the existence of residual stress is taken into account.     

Pullout performance of self-tapping medical screws

This study investigates the effect of the pilot hole size, implant depth, synthetic bone density, and screw size on the pullout strength of the self-tapping screw using analytical, finite element, and experimental methodologies. Stress distribution and failure propagation mode around the implant thread zone are also investigated. Based on the finite element analysis (FEA) results, an analytical model for the pullout strength of the self-tapping screw is constructed in terms of the (synthetic) bone mechanical properties, screw size, and the implant depth. The pullout performance of self-tapping screws is discussed. Results from the analytical and finite element models are experimentally validated.     

A comparison of the peri-implant bone stress generated by the preload with screw tightening between ‘bonded’ and ‘contact’ model

A number of finite element analyses (FEAs) for the dental implant were performed without regard for preload and with all interfaces ‘fixed-bonded’. The purpose of this study was comparing the stress distributions between the conventional FEA model with all contacting interfaces ‘fixed-bonded’ (bonded model) and the model with the interfaces of the components in ‘contact’ with friction simulated as a preloaded implant (contact model). We further verified the accuracy of the result of the FEA using model experiment. In the contact model, the stress was more widely distributed than in the bonded model. From the model study, the preload induced by screw tightening generated strain at the peri-implant bone, even before the application of external force. As a result, the bonded model could not reproduce the mechanical phenomena, whereas the contact model is considered to be appropriate for analysing mechanical problems.     

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°).     

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.