Prediction of insertion torque and stiffness of a dental implant in bovine trabecular bone using explicit micro-finite element analysis

Effects of dental implant diameter and tapered body design on stress distribution in rigid polyurethane foam during insertion

Anchorage, evaluated by the maximum insertion torque (IT), refers to mechanical engagement between dental implant and host bone at the time of insertion without external loads. Sufficient anchorage has been highly recommended in the clinic. In several studies, the effects of implant diameter and taper body design under external loading have been evaluated after insertion; however, there are few studies, in which their effects on stress distribution during insertion have been investigated to understand establishment of anchorage. Therefore, the objective of this study was to investigate the effects of dental implant diameter and tapered body design on anchorage combining experiments, analytical modeling, and finite element analysis (FEA). Two implant designs (parallel-walled and tapered) with two implant diameters were inserted into rigid polyurethane (PU) foam with corresponding straight drill protocols. The IT was fit to the analytical model (R2 = 0.88-1.0). The insertion process was modeled using explicit FEA. For parallel-walled implants, normalized IT and final FEA contact ratio were not related to the implant diameter while the implant diameter affected normalized IT (R2 = 0.90, p < 0.05, β1 = 0.20 and β2 = 0.93, standardized regression coefficients for implant diameter and taper body design) and final FEA contact ratio of tapered implants. The taper design distributed the PU foam stress further away from the thread compared to parallel-walled implants, which demonstrated compression in PU foam established by the tapered body during insertion.

Comparison of simplified bone-screw interface models in materially nonlinear μFE simulations

Micro finite-element (μFE) simulations serve as a crucial research tool to assist laboratory experiments in the biomechanical assessment of screw anchorage in bone. However, accurately modelling the interface between bone and screw threads at the microscale poses a significant challenge. Currently, the gold-standard approach involves employing computationally intensive physical contact models to simulate this interface. This study compared nonlinear μFE predictions of deformations, whole-construct stiffness, maximum force and damage patterns of three different computationally efficient simplified interface approaches to the general contact interface in Abaqus Explicit, which was defined as gold-standard and reference model. The μCT images (resolution: 32.8 μm) of two human radii with varying bone volume fractions were utilized and a screw was virtually inserted up to 50% and 100% of the volar-dorsal cortex distance. Materially nonlinear μFE models were generated and loaded in tension, compression and shear. In a first step, the common simplification of using a fully-bonded interface was compared to the general contact interface, revealing overestimations of whole-construct stiffness (19% on average) and maximum force (26% on average), along with inaccurate damage pattern replications. To enhance predictions, two additional simplified interface models were compared: tensionally strained element deletion (TED) and a novel modification of TED (TED-M). TED deletes interface elements strained in tension based on a linear-elastic simulation before the actual simulation. TED-M extends the remaining contact interface of TED by incorporating neighboring elements to the contact area. Both TED and TED-M reduced the errors in whole-construct stiffness and maximum force and improved the replication of the damage distributions in comparison to the fully-bonded approach. TED was better in predicting whole-construct stiffness (average error of 1%), while TED-M showed lowest errors in maximum force (1% on average). In conclusion, both TED and TED-M offer computationally efficient alternatives to physical contact modelling, although the fully-bonded interface may deliver sufficiently accurate predictions for many applications.

Critical loss of primary implant stability in osteosynthesis locking screws under cyclic overloading

Primary implant stability, which refers to the stability of the implant during the initial healing period is a crucial factor in determining the long-term success of the implant and lays the foundation for secondary implant stability achieved through osseointegration. Factors affecting primary stability include implant design, surgical technique, and patient-specific factors like bone quality and morphology. In vivo, the cyclic nature of anatomical loading puts osteosynthesis locking screws under dynamic loads, which can lead to the formation of micro cracks and defects that slowly degrade the mechanical connection between the bone and screw, thus compromising the initial stability and secondary stability of the implant. Monotonic quasi-static loading used for testing the holding capacity of implanted screws is not well suited to capture this behavior since it cannot capture the progressive deterioration of peri‑implant bone at small displacements. In order to address this issue, this study aims to determine a critical point of loss of primary implant stability in osteosynthesis locking screws under cyclic overloading by investigating the evolution of damage, dissipated energy, and permanent deformation. A custom-made test setup was used to test implanted 2.5 mm locking screws under cyclic overloading test. For each loading cycle, maximum forces and displacement were recorded as well as initial and final cycle displacements and used to calculate damage and energy dissipation evolution. The results of this study demonstrate that for axial, shear, and mixed loading significant damage and energy dissipation can be observed at approximately 20 % of the failure force. Additionally, at this load level, permanent deformations on the screw-bone interface were found to be in the range of 50 to 150 mm which promotes osseointegration and secondary implant stability. This research can assist surgeons in making informed preoperative decisions by providing a better understanding of the critical point of loss of primary implant stability, thus improving the long-term success of the implant and overall patient satisfaction.

Investigation on primary stability of dental implants through considering peri-implant bone damage, caused by small and large deformations: A validated non-linear micro finite element study

Primary stability of a dental implant is defined as its ability to resist the applied load without showing excessive damage in peri-implant bone, which is a prerequisite for secondary stability, and consequently for implantation success. The main goal of this study was to develop a validated micro-finite element (μFE) approach to assess the primary stability of dental implants in terms of stiffness, stiffness reduction, and irreversible displacement of the bone-implant system, subjected to an increasing step-wise quasi-static compressive loading-unloading test. The μFE models were generated based on the μCT images of bone, taken from extracted bovine tibia trabecular bone samples after drilling and implantation. A tissue constitutive model was considered for trabecular bone by describing elasto-plasticity with a modified von Mises yield criterion and element deletion technique to account for trabecular bone damage behavior. Then, the obtained force-displacement curves from the simulation were compared with the in-vitro mechanical test curves to evaluate the validity of the model. The results showed that the proposed μFE model could be properly predict the bone-implant system mechanical response in terms of irreversible displacement (R2 = 0.99), stiffness (R2 = 0.77), and stiffness reduction (R2 = 0.72) of the bone-implant construct for all the applied displacements without a significant difference from the unit slope and zero intercept of the QQ-plot (p-value<0.05). Moreover, a qualitative agreement was seen between the peri-implant bone damage predicted by the μFE model and the observed from μCT images. The adopted methodology used in this study can predict the mechanical failure response of the bone-implant system, which can be employed as a representative tool to study the effects of various dental implant design parameters on the primary stability with the ultimate goal of optimizing dental implants design.

Accuracy of osseointegrated screw-bone construct stiffness and peri-implant loading predicted by homogenized FE models relative to micro-FE models

Computational predictions of stiffness and peri-implant loading of screw-bone constructs are highly relevant to investigate and improve bone fracture fixations. Homogenized finite element (hFE) models have been used for this purpose in the past, but their accuracy has been questioned given the numerous simplifications, such as neglecting screw threads and modelling the trabecular bone structure as a continuum. This study aimed to investigate the accuracy of hFE models of an osseointegrated screw-bone construct when compared to micro-FE models considering the simplified screw geometry and different trabecular bone material models. Micro-FE and hFE models were created from 15 cylindrical bone samples with a virtually inserted, osseointegrated screw (fully bonded interface). Micro-FE models were created including the screw with threads (=reference models) and without threads to quantify the error due to screw geometry simplification. In the hFE models, the screws were modelled without threads and four different trabecular bone material models were used, including orthotropic and isotropic material derived from homogenization with kinematic uniform boundary conditions (KUBC), as well as from periodicity-compatible mixed uniform boundary conditions (PMUBC). Three load cases were simulated (pullout, shear in two directions) and errors in the construct stiffness and the volume average strain energy density (SED) in the peri-implant region were evaluated relative to the micro-FE model with a threaded screw. The pooled error caused by only omitting screw threads was low (max: 8.0%) compared to the pooled error additionally including homogenized trabecular bone material (max: 92.2%). Stiffness was predicted most accurately using PMUBC-derived orthotropic material (error: -0.7 ± 8.0%) and least accurately using KUBC-derived isotropic material (error: +23.1 ± 24.4%). Peri-implant SED averages were generally well correlated (R² ≥ 0.76), but slightly over- or underestimated by the hFE models and SED distributions were qualitatively different between hFE and micro-FE models. This study suggests that osseointegrated screw-bone construct stiffness can be predicted accurately using hFE models when compared to micro-FE models and that volume average peri-implant SEDs are well correlated. However, the hFE models are highly sensitive to the choice of trabecular bone material properties. PMUBC-derived isotropic material properties represented the best trade-off between model accuracy and complexity in this study.

Comparison of linear and nonlinear stepwise μFE displacement predictions to digital volume correlation measurements of trabecular bone biopsies

Digital volume correlation (DVC) enables to evaluate the ability of μFE models in predicting experimental results on the mesoscale. In this study predicted displacement fields of three different linear and materially nonlinear μFE simulation methods were compared to DVC measured displacement fields at specific load steps in the elastic regime (Step_El) and after yield (Step_Ult). Five human trabecular bone biopsies from a previous study were compressed in several displacement steps until failure. At every compression step, μCT images (resolution: 36 μm) were recorded. A global DVC algorithm was applied to compute the displacement fields at all loading steps. The unloaded 3D images were then used to generate homogeneous, isotropic, linear and materially nonlinear μFE models. Three different μFE simulation methods were used: linear (L), nonlinear (NL), and nonlinear stepwise (NLS). Regarding L and NL, the boundary conditions were derived from the interpolated displacement fields at Step_El and Step_Ult, while for the NLS method nonlinear changes of the boundary conditions of the experiments were captured using the DVC displacement field of every available load step until Step_El and Step_Ult. The predicted displacement fields of all μFE simulation methods were in good agreement with the DVC measured displacement fields (individual specimens: R²>0.83 at Step_El and R²>0.59 at Step_Ult; pooled data: R²>0.97 at Step_El and R²>0.92 at Step_Ult). At Step_El, all three simulation methods showed similar intercepts, slopes, and coefficients of determination while the nonlinear μFE models improved the prediction of the displacement fields slightly in all Cartesian directions at Step_Ult (individual specimens: L: R²>0.59 and NL, NLS: R²>0.68; pooled data: L: R²>0.92 and NL, NLS: R²>0.94). Damaged/overstrained elements in L, NL, and NLS occurred at similar locations but the number of overstrained elements was overestimated when using the L simulation method. Considering the increased solving time of the nonlinear μFE models as well as the acceptable performance in displacement prediction of the linear μFE models, one can conclude that for similar use cases linear μFE models represent the best compromise between computational effort and accuracy of the displacement field predictions.

Measurement of bone damage caused by quasi-static compressive loading-unloading to explore dental implants stability: Simultaneous use of in-vitro tests, μ-CT images, and digital volume correlation

Primary stability of dental implants is the initial mechanical engagement of the implant with its adjacent bone. Implantation and the subsequent loading may cause mechanical damage in the peripheral bone, which ultimately reduces the stability of the implant. This study aimed at evaluating primary stability of dental implants through applying stepwise compressive displacement-controlled, loading-unloading cycles to obtain overall stiffness and dissipated energy of the bone-implant structure; and quantifying induced plastic strains in surrounding bone using digital volume correlation (DVC) method, through comparing μCT images in different loading steps. To this end, dental implants were inserted into the cylindrical trabecular bones, then the bone-implant structure was undergone step-wise loading-unloading cycles, and μCT images were taken in some particular steps, then comparison was made between undeformed and deformed configurations using DVC to quantify plastic strain within the trabecular bone. Comparing stiffness reduction and dissipated energy values in different loading steps, obtained from the force-displacement curve in each loading step, revealed that the maximum displacement of 0.16 mm can be deemed as a safe threshold above which damages in peri-implant bone started to increase considerably (p < 0.05). In addition, it was found here that peri-implant bone strain linearly increased with decreasing bone-implant stiffness (p < 0.05). Moreover, strain concentration in peri-implant bone region showed that the plastic strain in trabecular bone spread up to a distance of about 2.5 mm away from the implant surface. Research of this kind can be used to optimize the design of dental implants, with the ultimate goal of improving their stability, also to validate in-silico models, e.g., micro-finite element models, which can help gain a deeper understanding of bone-implant construct behavior.

Evaluation of the efficacy of modal analysis in predicting the pullout strength of fixation bone screws

Background

Pilot hole preparation has been shown to have an impact on the short and long-term stability of the screw fixation constructs.

Purpose

Investigation and comparison of two nondestructive modal analysis methods with conventional insertion torque (IT) and pullout tests in optimum pilot hole diameter detection.

Methods

Twenty conical core titanium screws were embedded in high-density polyethylene blocks with different pilot hole diameters. The maximum IT was recorded for each screw during implantation. Then, two modal analysis methods including accelerometer (classical modal analysis [CMA]) and acoustic modal analysis (AMA) were carried out to measure the natural frequency (NF) of the bone-screw structure. Finally, stiffness (S), pullout force (F_ult), displacement at F_ult (d_ult) and energy dissipation (ED) were obtained from the destructive pullout test.

Results

The IT increased, as the pilot hole diameter decreased. The maximum value of IT was observed in the smallest pilot hole diameter. The same trend was found for the F_ult and the first NF derived from both modal methods except for 5.5 mm pilot hole diameter. The natural NFs derived from CMA and AMA showed high correlations in different groups (R² = 0.94) and did not deviate from y = x hypothesis in linear regression analysis. The F_ult, d_ult, and ED were measured 4800 ± 172 N, 3.10 ± 0.08 mm and 14.23 ± 1.10 N.mm, respectively.

Discussion

No significant change was observed in "S" between the groups. The highest F_ult and first NF were obtained for the 5.5 mm pilot hole diameter. Both CMA and AMA were found to be reliable methods and can promote the undesirable contradiction between F_ult and IT.

Non-linear explicit micro-FE models accurately predict axial pull-out force of cortical screws in human tibial cortical bone

Screws are the most frequently used implants for treatment of bone fractures and play an essential role in determining fixation stability. Robust prediction of the bone-screw interface failure would enable development of improved fixation strategies and implant designs, ultimately reducing failure rates and improving outcomes of bone fracture treatments. This study aimed to compare the accuracy of micro-computed tomography image based bone volume measures, linear micro-finite element (FE) and non-linear micro-FE simulations in predicting pull-out force of 3.5 mm screws in human cadaveric tibial cortical bone. Axial pull-out experiments were performed in forty samples harvested from a single human tibia to measure ultimate force, which was correlated with bone volume around the screw and the predictions by both linear micro-FE and non-linear explicit micro-FE models. Correlation strength was similar for bone volume around the screw (R² = 0.866) and linear micro-FE (R² = 0.861), but the explicit non-linear micro-FE models were able to capture the experimental results more accurately (R² = 0.913) and quantitatively correctly. Therefore, this technique may have potential for future in silico studies aiming at implant design optimization.

Influences of screw design features on initial stability in immediate implant placement and restoration

BACKGROUND

Self-tapping screws have been extensively used for dental implants. Their biomechanical behavior is highly associated with their clinical success, particularly for screws used for immediate implant placement and restoration, because occlusal forces can directly affect the loading transfer at the bone-implant interface after implantation. The effect of implant design on the initial stability of self-tapping screws remains unclear. This study explored the biomechanical behaviors of implant stability in standardized implants with different design features.

METHODS

Six types of dental implants were designed using computer-aided design/computer-aided manufacturing technology, including three types of cutting flute shapes (spiral, straight, and non-self-tapping) combined with two types of screw features. Peak insertion torque values were first recorded; initial stability levels were subsequently evaluated in terms of the maximum force and resistance to lateral loads using an electrodynamic test system.

FINDINGS

The peak insertion torque values, maximum force, and resistance to lateral loads of the non-self-tapping groups were higher than those of the self-tapping groups by 17%-90% (p < 0.01). The peak insertion torque values of the Straumann implant with a spiral flute was higher than that of the original straight flute by 20% (p < 0.001). However, compared with the original spiral flute, the Nobel Biocare implant with straight flute had a 23% higher maximum force (p = 0.016) and 24.5% higher resistance (p = 0.012) under lateral loading.

INTERPRETATION

Changing the flute design would affect initial implant stability. Non-self-tapping implants exhibited superior initial stability than did self-tapping implants.

Multi-objective design optimization of dental implant geometrical parameters

In-silico investigations are becoming an integral part of the development of novel biomedical devices, including dental implants. Using computer simulations can streamline the process by tuning different geometrical and structural features, emphasizing the osseointegration of the implant design a priori, leading to the optimal designs in preparation for in-vivo trails. This research aims to elucidate the interrelationship between 12 geometrical variables that holistically define the shape of the implant. The approach to achieve optimality hinged on coupling the finite element analysis results with the fractional factorial design method. The latter was used to determine the most influential variables during the screening process, followed by the parameter optimization process using the response surface method, regarding four different objectives, namely: bone-implant contact area, volume of trabecular bone dead cells, volume of cortical bone dead cells, and axial displacement. This resulted in reducing the number of virtual experiments and substantially decreasing the computational cost without compromising the accuracy of the solution. It was found that the optimized values improved the performance significantly. The validity of all models was verified by comparing optimized responses with simulation results. A sensitivity analysis was performed on all five optimized models to address the effect of friction coefficient on the implant-bone joint interaction. It was shown that the mechanical behavior of implant-bone would be independent in higher friction coefficients. The significance of this study is demonstrated in determining the most effective and optimized values of all possible geometrical parameters considering their singular or interactive effects.

Homogenized finite element models can accurately predict screw pull-out in continuum materials, but not in porous materials

BACKGROUND AND OBJECTIVE

Bone screw fixation can be estimated with several test methods such as insertion torque, pull-out, push-in and bending tests. A basic understanding of the relationship between screw fixation and bone microstructure is still lacking. Computational models can help clarify this relationship. The objective of the paper is to evaluate homogenized finite element (hFE) models of bone screw pull-out.

METHODS

Experimental pull-out tests were performed on three materials: two polyurethane (PU) foams having a porous microstructure, and a high density polyethylene (HDPE) which is a continuum material. Forty-five titanium pedicle screws were inserted to 10, 20, and 30 mm in equally sized blocks of all three materials (N = 5/group). Pull-out characteristics i.e. stiffness (S), yield force (F_y), peak pull-out force (F_ult) and displacement at F_ult (d_ult) were measured. hFE models were created replicating the experiments. The screw was modeled as a rigid body and 5 mm axial displacement was applied to the head of the screw. Simulations were performed evaluating two different conditions at the bone-screw interface; once in which the screw fitted the pilot hole exactly ("free-stressed") and once in which interface stresses resulting from the insertion process were taken into account ("pre-stressed").

RESULTS

The simulations representing the pre-stressed condition in HDPE matched the experimental data well; S, F_y, and F_ult differed less than 11%, 2% and 0.5% from the experimental data, respectively, whereas d_ult differed less than 16%. The free-stressed simulations were less accurate, especially stiffness (158% higher than the pre-stressed condition) and d_ult (30% lower than pre-stressed condition) were affected. The simulations representing PU did not match the experiments well. For the 20 mm insertion depth, S, F_y and F_ult differed by more than 104%, 89% and 66%, respectively from the experimental values. Agreement did not improve for 10 and 30 mm insertion depths.

CONCLUSIONS

We found that hFE models can accurately quantify screw pull-out in continuum materials such as HDPE, but not in materials with a porous structure, such as PU. Pre-stresses in the bone induced by the insertion process cannot be neglected and need to be included in the hFE simulations.

High-Performance Computing Comparison of Implicit and Explicit Nonlinear Finite Element Simulations of Trabecular Bone

BACKGROUND AND OBJECTIVE

Finite element models built from micro-computed tomography scans have become a powerful tool to investigate the mechanical properties of trabecular bone. There are two types of solving algorithms in the finite element method: implicit and explicit. Both of these methods have been utilized to study the trabecular bone. However, an investigation comparing the results obtained using the implicit and explicit solvers is lacking. Thus, in this paper, we contrast implicit and explicit procedures by analyzing trabecular bone samples as a case study.

METHODS

Micro-computed tomography-based finite element analysis of trabecular bone under a direct quasi-static compression was done using implicit and explicit methods. The differences in the predictions of mechanical properties and computational time of the two methods were studied using high-performance computing.

RESULTS

Our findings indicate that the results using implicit and explicit solvers are well comparable, given that similar problem set up is carefully utilized. Also, the parallel scalability of the two methods was similar, while the explicit solver performed about five times faster than the implicit method. Along with faster performance, the explicit method utilized significantly less memory for the analysis, which shows another benefit of using an explicit solver for this case study.

CONCLUSIONS

The comparison of the implicit and explicit methods for the simulation of trabecular bone samples should be highly valuable to the bone modeling community and researchers studying complex cellular and architectured materials.

Evaluation of Bone-Implant Interface Stress and Strain Using Heterogeneous Mandibular Bone Properties Based on Different Empirical Correlations

Objectives  The purpose of this study was to compare methods used for calculating heterogeneous patient-specific bone properties used in finite element analysis (FEA), in the field of implant dentistry, with the method based on homogenous bone properties.

Materials and Methods  In this study, three-dimensional (3D) computed tomography data of an edentulous patient were processed to create a finite element model, and five identical 3D implant models were created and distributed throughout the dental arch. Based on the calculation methods used for bone material assignment, four groups—groups I to IV—were defined. Groups I to III relied on heterogeneous bone property assignment based on different equations, whereas group IV relied on homogenous bone properties. Finally, 150 N vertical and 60-degree-inclined forces were applied at the top of the implant abutments to calculate the von Mises stress and strain.

Results  Groups I and II presented the highest stress and strain values, respectively. Based on the implant location, differences were observed between the stress values of group I, II, and III compared with group IV; however, no clear order was noted. Accordingly, variable von Mises stress and strain reactions at the bone–implant interface were observed among the heterogeneous bone property groups when compared with the homogenous property group results at the same implant positions.

Conclusion  Although the use of heterogeneous bone properties as material assignments in FEA studies seem promising for patient-specific analysis, the variations between their results raise doubts about their reliability. The results were influenced by implants’ locations leading to misleading clinical simulations.

Experimental and numerical investigation of secondary screw perforation in the human proximal humerus

Surgical treatment of proximal humerus fractures remains challenging, with a reported failure rate ranging from 15% to 35%. The dominant failure mode is secondary, i.e. post-operative screw perforation through the glenohumeral joint. A better understanding and the ability to predict this complication could lead to improved fracture fixation and decreased failure rate. The aims of this study were (1) to develop an experimental model for single screw perforation in the human humeral head and (2) to evaluate the ability of densitometric measures and micro finite element (microFE) analyses to predict the experimental failure event. Screw perforation was investigated experimentally in twenty cuboidal specimens cut from four pairs of fresh-frozen human cadaveric proximal humeral heads. A centrally inserted 3.5 mm screw was pushed quasi-statically at a constant displacement rate until perforation of the articular cartilage in each specimen. Force and displacement were recorded and evaluated at both initial screw loosening and perforation events. Bone volume was calculated around and in front of the screw and tip-to-joint distance was measured on the combined pre- and post-instrumentation micro computed tomography (microCT) scans. Implicit linear and explicit non-linear microFE models were created based on the microCT scans. The strength of these densitometric, geometrical and microFE methods to predict the experimental results was evaluated via correlation analysis. The bone volume measures were optimized in a parametric analysis to maximize correlation coefficients. The strongest and quantitatively correct predictions of perforation force (R² = 0.93) and displacement (R² = 0.77) were achieved using the explicit, non-linear microFE models. Linear microFE simulations provided the strongest predictions of loosening force (R² = 0.90). Correlation strengths reached by optimized bone volume measures for predicting experimental force and by tip-to-joint distance for predicting displacement were only slightly inferior compared to the results of microFE models. The strong correlations achieved with densitometric and geometric measures indicate that monotonic perforation of single screws through the articular surface of the humeral head can be well predicted with these easily accessible measures. However, non-linear microFE models delivered even stronger correlations and quantitatively correct predictions of perforation force and displacement. This indicates that if computational resources are available, non-linear simulations may have a high potential to investigate more complex fixations and loading scenarios.