The discrete nature of trabecular bone microarchitecture affects implant stability

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.

Micro finite element analysis of continuously loaded mini-implants - A micro-CT study in the rat tail model

Implant migration has been described as a minor displacement of orthodontic mini-implants (OMIs) when subjected to constant forces. Aim of this study was to evaluate the impact of local stresses on implant migration and bone remodelling around constantly loaded OMIs. Two mini-implants were placed in one caudal vertebra of 61 rats, connected by a nickel‑titanium contraction spring, and loaded with different forces (0.0, 0.5, 1.0, 1.5 N). In vivo micro-CT scans were taken immediately and 1, 2 (n = 61), 4, 6 and 8 (n = 31) weeks post-op. Nine volumes of interest (VOIs) around each implant were defined. To analyse stress values, micro-finite element models were created. Bone remodelling was analysed by calculating the bone volume change between scans performed at consecutive time points. Statistical analysis was performed using a linear mixed model and likelihood-ratio-tests, followed by Tuckey post hoc tests when indicated. The highest stresses were observed in the proximal top VOI. In all VOIs, stress values tended to reach their maximum after two weeks and decreased thereafter. Bone remodelling analysis revealed initial bone loss within the first two weeks and bone gain up to week eight, which was noted especially in the highest loading group. The magnitude of local stresses influenced bone remodelling and it can be speculated that the stress related bone resorption favoured implant migration. After a first healing phase with a high degree of bone resorption, net bone gain representing consolidation was observed.

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.

Finite element analysis in implant dentistry: State of the art and future directions

OBJECTIVE

To discuss the state of the art of Finite Element (FE) modeling in implant dentistry, to highlight the principal features and the current limitations, and giving recommendations to pave the way for future studies.

METHODS

The articles' search was performed through PubMed, Web of Science, Scopus, Science Direct, and Google Scholar using specific keywords. The articles were selected based on the inclusion and exclusion criteria, after title, abstract and full-text evaluation. A total of 147 studies were included in this review.

RESULTS

To date, the FE analysis of the bone-dental implant system has been investigated by analyzing several types of implants; modeling only a portion of bone considered as isotropic material, despite its anisotropic behavior; assuming in most cases complete osseointegration; considering compressive or oblique forces acting on the implant; neglecting muscle forces and the bone remodeling process. Finally, there is no standardized approach for FE modeling in the dentistry field.

SIGNIFICANCE

FE modeling is an effective computational tool to investigate the long-term stability of implants. The ultimate aim is to transfer such technology into clinical practice to help dentists in the diagnostic and therapeutic phases. To do this, future research should deeply investigate the loading influence on the bone-implant complex at a microscale level. This is a key factor still not adequately studied. Thus, a multiscale model could be useful, allowing to account for this information through multiple length scales. It could help to obtain information about the relationship among implant design, distribution of bone stress, and bone growth. Finally, the adoption of a standardized approach will be necessary, in order to make FE modeling highly predictive of the implant's long-term stability.

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.

The effect of cement augmentation on pedicle screw fixation under various load cases : results from a combined experimental, micro-CT, and micro-finite element analysis

Aims

Anchorage of pedicle screw rod instrumentation in the elderly spine with poor bone quality remains challenging. Our study aims to evaluate how the screw bone anchorage is affected by screw design, bone quality, loading conditions, and cementing techniques.

Methods

Micro-finite element (µFE) models were created from micro-CT (μCT) scans of vertebrae implanted with two types of pedicle screws (L: Ennovate and R: S⁴). Simulations were conducted for a 10 mm radius region of interest (ROI) around each screw and for a full vertebra (FV) where different cementing scenarios were simulated around the screw tips. Stiffness was calculated in pull-out and anterior bending loads.

Results

Experimental pull-out strengths were excellently correlated to the µFE pull-out stiffness of the ROI (R² > 0.87) and FV (R² > 0.84) models. No significant difference due to screw design was observed. Cement augmentation increased pull-out stiffness by up to 94% and 48% for L and R screws, respectively, but only increased bending stiffness by up to 6.9% and 1.5%, respectively. Cementing involving only one screw tip resulted in lower stiffness increases in all tested screw designs and loading cases. The stiffening effect of cement augmentation on pull-out and bending stiffness was strongly and negatively correlated to local bone density around the screw (correlation coefficient (R) = -0.95).

Conclusion

This combined experimental, µCT and µFE study showed that regional analyses may be sufficient to predict fixation strength in pull-out and that full analyses could show that cement augmentation around pedicle screws increased fixation stiffness in both pull-out and bending, especially for low-density bone.

Cite this article: Bone Joint Res 2021;10(12):797–806.

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.

Three-Dimensional Finite Element Investigation into Effects of Implant Thread Design and Loading Rate on Stress Distribution in Dental Implants and Anisotropic Bone

Variations in the implant thread shape and occlusal load behavior may result in significant changes in the biological and mechanical properties of dental implants and surrounding bone tissue. Most previous studies consider a single implant thread design, an isotropic bone structure, and a static occlusal load. However, the effects of different thread designs, bone material properties, and loading conditions are important concerns in clinical practice. Accordingly, the present study performs Finite Element Analysis (FEA) simulations to investigate the static, quasi-static and dynamic response of the implant and implanted bone material under various thread designs and occlusal loading directions (buccal-lingual, mesiodistal and apical). The simulations focus specifically on the von Mises stress, displacement, shear stress, compressive stress, and tensile stress within the implant and the surrounding bone. The results show that the thread design and occlusal loading rate have a significant effect on the stress distribution and deformation of the implant and bone structure during clinical applications. Overall, the results provide a useful insight into the design of enhanced dental implants for an improved load transfer efficiency and success rate.

Finite Element Analysis of Fracture Fixation

PURPOSE OF REVIEW

Fracture fixation aims to provide stability and promote healing, but remains challenging in unstable and osteoporotic fractures with increased risk of construct failure and nonunion. The first part of this article reviews the clinical motivation behind finite element analysis of fracture fixation, its strengths and weaknesses, how models are developed and validated, and how outputs are typically interpreted. The second part reviews recent modeling studies of the femur and proximal humerus, areas with particular relevance to fragility fractures.

RECENT FINDINGS

There is some consensus in the literature around how certain modeling aspects are pragmatically formulated, including bone and implant geometries, meshing, material properties, interactions, and loads and boundary conditions. Studies most often focus on predicted implant stress, bone strain surrounding screws, or interfragmentary displacements. However, most models are not rigorously validated. With refined modeling methods, improved validation efforts, and large-scale systematic analyses, finite element analysis is poised to advance the understanding of fracture fixation failure, enable optimization of implant designs, and improve surgical guidance.

On the limits of finite element models created from (micro)CT datasets and used in studies of bone-implant-related biomechanical problems

Patient-specific approach is gaining a wide popularity in computational simulations of biomechanical systems. Simulations (most often based on the finite element method) are to date routinely created using data from imaging devices such as computed tomography which makes the models seemingly very complex and sophisticated. However, using a computed tomography in finite element calculations does not necessarily enhance the quality or even credibility of the models as these depend on the quality of the input images. Low-resolution (medical-)CT datasets do not always offer detailed representation of trabecular bone in FE models and thus might lead to incorrect calculation of mechanical response to external loading. The effect of image resolution on mechanical simulations of bone-implant interaction has not been thoroughly studied yet. In this study, the effect of image resolution on the modeling procedure and resulting mechanical strains in bone was analyzed on the example of cranial implant. For this purpose, several finite element models of bone interacting with fixation-screws were generated using seven computed tomography datasets of a bone specimen but with different image resolutions (ranging from micro-CT resolution of 25 μm to medical-CT resolution of 1250 μm). The comparative analysis revealed that FE models created from images of low resolution (obtained from medical computed tomography) can produce biased results. There are two main reasons: 1. Medical computed tomography images do not allow generating models with complex trabecular architecture which leads to substituting of the intertrabecular pores with a fictitious mass; 2. Image gray value distribution can be distorted resulting in incorrect mechanical properties of the bone and thus in unrealistic or even completely fictitious mechanical strains. The biased results of calculated mechanical strains can lead to incorrect conclusion, especially when bone-implant interaction is investigated. The image resolution was observed not to significantly affect stresses in the fixation screw itself; however, selection of bone material representation might result in significantly different stresses in the screw.

Biomechanical comparison of sacral and transarticular sacroiliac screw fixation

STUDY DESIGN

A detailed finite element analysis of screw fixation in the sacrum and pelvis.

OBJECTIVE

To biomechanically assess and compare the fixation performance of sacral and transarticular sacroiliac screws. Instrumentation constructs are used to achieve fixation and stabilization for the treatment of spinopelvic pathologies. The optimal screw trajectory and type of bone engagement to caudally anchor long fusion constructs are not yet known.

METHODS

A detailed finite element model of the sacroiliac articulation with two different bone densities was developed. Two sacral and one transarticular sacroiliac screw trajectories were modeled with different diameters (5.5 and 6.5 mm) and lengths (uni-cortical, bi-cortical and quad-cortical purchase). Axial pullout and flexion/extension toggle forces were applied on the screws representing intra and post-operative loads. The force-displacement results and von Mises stresses were used to characterize the failure pattern.

RESULTS

Overall, sacroiliac screws provided forces to failure 2.75 times higher than sacral fixation screws. On the contrary, the initial stiffness was approximately half as much for sacroiliac screws. High stresses were located at screw tips for the sacral trajectories and near the cortical bone screw entry points for the sacroiliac trajectory. Overall, the diameter and length of the screws had significant effects on the screw fixation (33% increase in force to failure; 5% increase in initial stiffness). A 20% drop in bone mineral density (lower bone quality) decreased the initial stiffness by 25% and the force to failure by 5-10%. High stresses and failure occurred at the screw tip for uni- and tri-cortical screws and were close to trabecular-cortical bone interface for bi-cortical and quad-cortical screws.

CONCLUSIONS

Sacroiliac fixation provided better anchorage than sacral fixation. The transarticular purchase of the sacroiliac trajectory resulted in differences in failure pattern and fixation performance.

Time-dependent behaviour of bone accentuates loosening in the fixation of fractures using bone-screw systems

Aims

Loosening is a well-known complication in the fixation of fractures using devices such as locking plates or unilateral fixators. It is believed that high strains in the bone at the bone-screw interface can initiate loosening, which can result in infection, and further loosening. Here, we present a new theory of loosening of implants. The time-dependent response of bone subjected to loads results in interfacial deformations in the bone which accumulate with cyclical loading and thus accentuates loosening.

Methods

We used an 'ideal' bone-screw system, in which the screw is subjected to cyclical lateral loads and trabecular bone is modelled as non-linear viscoelastic and non-linear viscoelastic-viscoplastic material, based on recent experiments, which we conducted.

Results

We found that the interfacial deformation in the bone increases with the number of cycles, and the use of a non-linear viscoelastic-viscoplastic model results in larger deformations, some of which are irrecoverable. There is an apparent trend in which interfacial deformations increase with increasing porosity of bone.

Conclusion

The developed time-dependent model of the mechanical behaviour of bone permits prediction of loosening due to cyclical loads, which has not been possible previously. Application of this model shows that implant loosening will be accentuated by cyclical loading due to physiological activities, and the risks of loosening are greater in osteoporotic patients.Cite this article: S. Xie, K. Manda, P. Pankaj. Time-dependent behaviour of bone accentuates loosening in the fixation of fractures using bone-screw systems. Bone Joint Res 2018;7:580-586. DOI: 10.1302/2046-3758.710.BJR-2018-0085.R1.

Mechanics of linear microcracking in trabecular bone

Microcracking in trabecular bone is responsible both for the mechanical degradation and remodeling of the trabecular bone tissue. Recent results on trabecular bone mechanics have demonstrated that bone tissue microarchitecture, tissue elastic heterogeneity and tissue-level mechanical anisotropy all should be considered to obtain detailed information on the mechanical stress state. The present study investigated the influence of tissue microarchitecture, tissue heterogeneity in elasticity and material separation properties and tissue-level anisotropy on the microcrack formation process. Microscale bone models were executed with the extended finite element method. It was demonstrated that anisotropy and heterogeneity of the bone tissue contribute significantly to bone tissue toughness and the resistance of trabecular bone to microcrack formation. The compressive strain to microcrack initiation was computed to increase by a factor of four from an assumed homogeneous isotropic tissue to an assumed anisotropic heterogenous tissue.

Multi-level hp-finite cell method for embedded interface problems with application in biomechanics

This work presents a numerical discretization technique for solving 3-dimensional material interface problems involving complex geometry without conforming mesh generation. The finite cell method (FCM), which is a high-order fictitious domain approach, is used for the numerical approximation of the solution without a boundary-conforming mesh. Weak discontinuities at material interfaces are resolved by using separate FCM meshes for each material sub-domain and weakly enforcing the interface conditions between the different meshes. Additionally, a recently developed hierarchical hp-refinement scheme is used to locally refine the FCM meshes to resolve singularities and local solution features at the interfaces. Thereby, higher convergence rates are achievable for nonsmooth problems. A series of numerical experiments with 2- and 3-dimensional benchmark problems is presented, showing that the proposed hp-refinement scheme in conjunction with the weak enforcement of the interface conditions leads to a significant improvement of the convergence rates, even in the presence of singularities. Finally, the proposed technique is applied to simulate a vertebra-implant model. The application showcases the method's potential as an accurate simulation tool for biomechanical problems involving complex geometry, and it demonstrates its flexibility in dealing with different types of geometric description.

Development and initial validation of a novel smoothed-particle hydrodynamics-based simulation model of trabecular bone penetration by metallic implants

A novel computational model of implant migration in trabecular bone was developed using smoothed-particle hydrodynamics (SPH), and an initial validation was performed via correlation with experimental data. Six fresh-frozen human cadaveric specimens measuring 10 × 10 × 20 mm were extracted from the proximal femurs of female donors (mean age of 82 years, range 75-90, BV/TV ratios between 17.88% and 30.49%). These specimens were then penetrated under axial loading to a depth of 10 mm with 5 mm diameter cylindrical indenters bearing either flat or sharp/conical tip designs similar to blunt and self-tapping cancellous screws, assigned in a random manner. SPH models were constructed based on microCT scans (17.33 µm) of the cadaveric specimens. Two initial specimens were used for calibration of material model parameters. The remaining four specimens were then simulated in silico using identical material model parameters. Peak forces varied between 92.0 and 365.0 N in the experiments, and 115.5-352.2 N in the SPH simulations. The concordance correlation coefficient between experimental and simulated pairs was 0.888, with a 95%CI of 0.8832-0.8926, a Pearson ρ (precision) value of 0.9396, and a bias correction factor Cb (accuracy) value of 0.945. Patterns of bone compaction were qualitatively similar; both experimental and simulated flat-tipped indenters produced dense regions of compacted material adjacent to the advancing face of the indenter, while sharp-tipped indenters deposited compacted material along their peripheries. Simulations based on SPH can produce accurate predictions of trabecular bone penetration that are useful for characterizing implant performance under high-strain loading conditions. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:1114-1123, 2018.

Micro finite element analysis of dental implants under different loading conditions

Osseointegration is paramount for the longevity of dental implants and is significantly influenced by biomechanical stimuli. The aim of the present study was to assess the micro-strain and displacement induced by loaded dental implants at different stages of osseointegration using finite element analysis (FEA). Computational models of two mandible segments with different trabecular densities were constructed using microCT data. Three different implant loading directions and two osseointegration stages were considered in the stress-strain analysis of the bone-implant assembly. The bony segments were analyzed using two approaches. The first approach was based on Mechanostat strain intervals and the second approach was based on tensile/compression yield strains. The results of this study revealed that bone surrounding dental implants is critically strained in cases when only a partial osseointegration is present and when an implant is loaded by buccolingual forces. In such cases, implants also encounter high stresses. Displacements of partially-osseointegrated implant are significantly larger than those of fully-osseointegrated implants. It can be concluded that the partial osseointegration is a potential risk in terms of implant longevity.

Micro-CT and micro-FE analysis of pedicle screw fixation under different loading conditions

Anchorage of pedicle screw instrumentation in the elderly spine with poor bone quality remains challenging. In this study, micro finite element (µFE) models were used to assess the specific influence of screw design and the relative contribution of local bone density to fixation mechanics. These were created from micro computer tomography (µCT) scans of vertebras implanted with two types of pedicle screws, including a full region-or-interest of 10 mm radius around each screw, as well as submodels for the pedicle and inner trabecular bone of the vertebral body. The local bone volume fraction (BV/TV) calculated from the µCT scans around different regions of the screw (pedicle, inner trabecular region of the vertebral body) were then related to the predicted stiffness in simulated pull-out tests as well as to the experimental pull-out and torsional fixation properties mechanically measured on the corresponding specimens. Results show that predicted stiffness correlated excellently with experimental pull-out strength (R² > 0.92, p < .043), better than regional BV/TV alone (R2 = 0.79, p = .003). They also show that correlations between fixation properties and BV/TV were increased when accounting only for the pedicle zone (R² = 0.66-0.94, p ≤  .032), but with weaker correlations for torsional loads (R² < 0.10). Our analyses highlight the role of local density in the pedicle zone on the fixation stiffness and strength of pedicle screws when pull-out loads are involved, but that local apparent bone density alone may not be sufficient to explain resistance in torsion.

A novel in silico method to quantify primary stability of screws in trabecular bone

Insufficient primary stability of screws in bone leads to screw loosening and failure. Unlike conventional continuum finite-element models, micro-CT based finite-element analysis (micro-FE) is capable of capturing the patient-specific bone micro-architecture, providing accurate estimates of bone stiffness. However, such in silico models for screws in bone highly overestimate the apparent stiffness. We hypothesized that a more accurate prediction of primary implant stability of screws in bone is possible by considering insertion-related bone damage. We assessed two different screw types and loading scenarios in 20 trabecular bone specimens extracted from 12 cadaveric human femoral heads (N = 5 for each case). In the micro-FE model, we predicted specimen-specific Young's moduli of the peri-implant bone damage region based on morphometric parameters such that the apparent stiffness of each in silico model matched the experimentally measured stiffness of the corresponding in vitro specimen as closely as possible. The standard micro-FE models assuming perfectly intact peri-implant bone overestimated the stiffness by over 330%. The consideration of insertion related damaged peri-implant bone corrected the mean absolute percentage error down to 11.4% for both loading scenarios and screw types. Cross-validation revealed a mean absolute percentage error of 14.2%. We present the validation of a novel micro-FE modeling technique to quantify the apparent stiffness of screws in trabecular bone. While the standard micro-FE model overestimated the bone-implant stiffness, the consideration of insertion-related bone damage was crucial for an accurate stiffness prediction. This approach provides an important step toward more accurate specimen-specific micro-FE models. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:2415-2424, 2017.

The Applications of Finite Element Analysis in Proximal Humeral Fractures

Proximal humeral fractures are common and most challenging, due to the complexity of the glenohumeral joint, especially in the geriatric population with impacted fractures, that the development of implants continues because currently the problems with their fixation are not solved. Pre-, intra-, and postoperative assessments are crucial in management of those patients. Finite element analysis, as one of the valuable tools, has been implemented as an effective and noninvasive method to analyze proximal humeral fractures, providing solid evidence for management of troublesome patients. However, no review article about the applications and effects of finite element analysis in assessing proximal humeral fractures has been reported yet. This review article summarized the applications, contribution, and clinical significance of finite element analysis in assessing proximal humeral fractures. Furthermore, the limitations of finite element analysis, the difficulties of more realistic simulation, and the validation and also the creation of validated FE models were discussed. We concluded that although some advancements in proximal humeral fractures researches have been made by using finite element analysis, utility of this powerful tool for routine clinical management and adequate simulation requires more state-of-the-art studies to provide evidence and bases.

Minimizing Pedicle Screw Pullout Risks: A Detailed Biomechanical Analysis of Screw Design and Placement

STUDY DESIGN

Detailed biomechanical analysis of the anchorage performance provided by different pedicle screw designs and placement strategies under pullout loading.

OBJECTIVE

To biomechanically characterize the specific effects of surgeon-specific pedicle screw design parameters on anchorage performance using a finite element model.

SUMMARY OF BACKGROUND DATA

Pedicle screw fixation is commonly used in the treatment of spinal pathologies. However, there is little consensus on the selection of an optimal screw type, size, and insertion trajectory depending on vertebra dimension and shape.

METHODS

Different screw diameters and lengths, threads, and insertion trajectories were computationally tested using a design of experiment approach. A detailed finite element model of an L3 vertebra was created including elastoplastic bone properties and contact interactions with the screws. Loads and boundary conditions were applied to the screws to simulate axial pullout tests. Force-displacement responses and internal stresses were analyzed to determine the specific effects of each parameter.

RESULTS

The design of experiment analysis revealed significant effects (P<0.01) for all tested principal parameters along with the interactions between diameter and trajectory. Screw diameter had the greatest impact on anchorage performance. The best insertion trajectory to resist pullout involved placing the screw threads closer to the pedicle walls using the straightforward insertion technique, which showed the importance of the cortical layer grip. The simulated cylindrical single-lead thread screws presented better biomechanical anchorage than the conical dual-lead thread screws in axial loading conditions.

CONCLUSIONS

The model made it possible to quantitatively measure the effects of both screw design characteristics and surgical choices, enabling to recommend strategies to improve single pedicle screw performance under axial loading.

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.

Influence of trabecular bone quality and implantation direction on press-fit mechanics

Achieving primary stability of uncemented press-fit prostheses in patients with poor quality bone can involve axial implantation forces large enough to cause bone fracture. Radial implantation eliminates intraoperative impaction forces and could prevent this damage. Platens of two commercial implant surfaces ("Beaded" and "Flaked") were implanted onto trabecular bone specimens of varying quality in a press-fit simulator. Samples were implanted with varying interference, either axially (shear) or radially (normal). Push-in and pull-out forces were measured to assess stability. Microstructural changes in the bone were determined from μCT analysis. For force-defined implantation analysis, push-in and pull-out forces both increased proportionally with increasing radial force, independent of implantation direction, bone quality or implant surface. For position-defined implantation analysis, pull-out forces were generally found to increase with interference and to be greater for radial than axial implantation direction, and to be lower for poor quality bone. Bone density increased locally at the tested interface due to implantation, in particular for the Beaded surface under axial implantation. If a safe radial stress can be determined for cortical bone in a particular patient, the associated implantation force, and pull-out force which represents primary stability, can be directly derived, regardless of implantation direction, bone quality or implant surface. Radial implantation delivers primary stability that is no worse than that for axial implantation and may eliminate potentially damaging impaction forces. Development of implant designs based on this principal might improve implant fixation. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:224-233, 2017.

Implicit modeling of screw threads for efficient finite element analysis of complex bone-implant systems

Finite element analysis is commonly used to assist in the development and evaluation of orthopedic devices. The physics of these models are simplified through approximations that enable more efficient simulations, without compromising the accuracy of the relative comparisons between implant designs or configurations. This study developed and evaluated a technique to approximate the behavior of a finely threaded screw using a smooth cylinder with the threads implicitly represented through interfacial contact conditions. This pseudo-threaded model was calibrated by comparing to simulations that explicitly modeled the thread geometry with frictional contact. A parametric analysis was performed with a single screw-in-bone system, five loading directions, and three Young׳s moduli that span the range of cancellous bone (200, 600, and 1,000MPa). Considering that screw cut-out from cancellous bone is a critical clinical issue in the osteoporotic proximal humerus, the pseudo-threaded model was compared with a bonded interface to examine three different screw configurations in a 3-part proximal humerus fracture across 10 patients. In the single screw-in-bone system, the pseudo-threaded model predicted the screw displacement of the explicitly threaded model with 1-5% difference and estimated the strain distributions and magnitudes more accurately than a bonded interface. Yet, the relative comparisons of implant stability across the three different screw configurations in the proximal humerus were not affected by the modeling choice for the bone-screw interface. Therefore, the bonded interface could serve as a more efficient methodology for making relative comparisons between implants that utilize the same thread profile.

Bone density and anisotropy affect periprosthetic cement and bone stresses after anatomical glenoid replacement: A micro finite element analysis

Glenoid loosening is still a main complication for shoulder arthroplasty. We hypothesize that cement and bone stresses potentially leading to fixation failure are related not only to glenohumeral conformity, fixation design or eccentric loading, but also to bone volume fraction, cortical thickness and degree of anisotropy in the glenoid. In this study, periprosthetic bone and cement stresses were computed with micro finite element models of the replaced glenoid depicting realistic bone microstructure. These models were used to quantify potential effects of bone microstructural parameters under loading conditions simulating different levels of glenohumeral conformity and eccentric loading simulating glenohumeral instability. Results show that peak cement stresses were achieved near the cement-bone interface in all loading schemes. Higher stresses within trabecular bone tissue and cement mantle were obtained within specimens of lower bone volume fraction and in regions of low anisotropy, increasing with decreasing glenohumeral conformity and reaching their maxima below the keeled design when the load is shifted superiorly. Our analyses confirm the combined influences of eccentric load shifts with reduced bone volume fraction and anisotropy on increasing periprosthetic stresses. They finally suggest that improving fixation of glenoid replacements must reduce internal cement and bone tissue stresses, in particular in glenoids of low bone density and heterogeneity.

Peri-implant stress correlates with bone and cement morphology: Micro-FE modeling of implanted cadaveric glenoids

Aseptic loosening of cemented joint replacements is a complex biological and mechanical process, and remains a clinical concern especially in patients with poor bone quality. Utilizing high resolution finite element analysis of a series of implanted cadaver glenoids, the objective of this study was to quantify relationships between construct morphology and resulting mechanical stresses in cement and trabeculae. Eight glenoid cadavers were implanted with a cemented central peg implant. Specimens were imaged by micro-CT, and subject-specific finite element models were developed. Bone volume fraction, glenoid width, implant-cortex distance, cement volume, cement-cortex contact, and cement-bone interface area were measured. Axial loading was applied to the implant of each model and stress distributions were characterized. Correlation analysis was completed across all specimens for pairs of morphological and mechanical variables. The amount of trabecular bone with high stress was strongly negatively correlated with both cement volume and contact between the cement and cortex (r = -0.85 and -0.84, p < 0.05). Bone with high stress was also correlated with both glenoid width and implant-cortex distance. Contact between the cement and underlying cortex may dramatically reduce trabecular bone stresses surrounding the cement, and this contact depends on bone shape, cement amount, and implant positioning.

The influence of bone density and anisotropy in finite element models of distal radius fracture osteosynthesis: Evaluations and comparison to experiments

Continuum-level finite element (FE) models can be used to analyze and improve osteosynthesis procedures for distal radius fractures (DRF) from a biomechanical point of view. However, previous models oversimplified the bone material and lacked thorough experimental validation. The goal of this study was to assess the influence of local bone density and anisotropy in FE models of DRF osteosynthesis for predictions of axial stiffness, implant plate stresses, and screw loads. Experiments and FE analysis were conducted in 25 fresh frozen cadaveric radii with DRFs treated by volar locking plate osteosynthesis. Specimen specific geometries were captured using clinical quantitative CT (QCT) scans of the prepared samples. Local bone material properties were computed based on high resolution CT (HR-pQCT) scans of the intact radii. The axial stiffness and individual screw loads were evaluated in FE models, with (1) orthotropic inhomogeneous (OrthoInhom), (2) isotropic inhomogeneous (IsoInhom), and (3) isotropic homogeneous (IsoHom) bone material and compared to the experimental axial stiffness and screw-plate interface failures. FE simulated and experimental axial stiffness correlated significantly (p<0.0001) for all three model types. The coefficient of determination was similar for OrthoInhom (R(2)=0.807) and IsoInhom (R(2)=0.816) models but considerably lower for IsoHom models (R(2)=0.500). The peak screw loads were in qualitative agreement with experimental screw-plate interface failure. Individual loads and implant plate stresses of IsoHom models differed significantly (p<0.05) from OrthoInhom and IsoInhom models. In conclusion, including local bone density in FE models of DRF osteosynthesis is essential whereas local bone anisotropy hardly effects the models׳ predictive abilities.

The relationship between porosity and specific surface in human cortical bone is subject specific

A characteristic relationship for bone between bone volume fraction (BV/TV) and specific surface (BS/TV) has previously been proposed based on 2D histological measurements. This relationship has been suggested to be bone intrinsic, i.e., to not depend on bone type, bone site and health state. In these studies, only limited data comes from cortical bone. The aim of this paper was to investigate the relationship between BV/TV and BS/TV in human cortical bone using high-resolution micro-CT imaging and the correlations with subject-specific biometric data such as height, weight, age and sex. Images from femoral cortical bone samples of the Melbourne Femur Collection were obtained using synchrotron radiation micro-CT (SPring8, Japan). Sixteen bone samples from thirteen individuals were analysed in order to find bone volume fraction values ranging from 0.20 to 1. Finally, morphological models of the tissue microstructure were developed to help explain the relationship between BV/TV and BS/TV. Our experimental findings indicate that the BV/TV vs BS/TV relationship is subject specific rather than intrinsic. Sex and pore density were statistically correlated with the individual curves. However no correlation was found with body height, weight or age. Experimental cortical data points deviate from interpolating curves previously proposed in the literature. However, these curves are largely based on data points from trabecular bone samples. This finding challenges the universality of the curve: highly porous cortical bone is significantly different to trabecular bone of the same porosity. Finally, our morphological models suggest that changes in BV/TV within the same sample can be explained by an increase in pore area rather than in pore density. This is consistent with the proposed mechanisms of age-related endocortical bone loss. In addition, these morphological models highlight that the relationship between BV/TV and BS/TV is not linear at high BV/TV as suggested in the literature but is closer to a square root function.

Finite element analysis of dental implant loading on atrophic and non-atrophic cancellous and cortical mandibular bone - a feasibility study

The first aim of this study was to assess displacements and micro-strain induced on different grades of atrophic cortical and trabecular mandibular bone by axially loaded dental implants using finite element analysis (FEA). The second aim was to assess the micro-strain induced by different implant geometries and the levels of bone-to-implant contact (BIC) on the surrounding bone. Six mandibular bone segments demonstrating different grades of mandibular bone atrophy and various bone volume fractions (from 0.149 to 0.471) were imaged using a micro-CT device. The acquired bone STL models and implant (Brånemark, Straumann, Ankylos) were merged into a three-dimensional finite elements structure. The mean displacement value for all implants was 3.1 ±1.2 µm. Displacements were lower in the group with a strong BIC. The results indicated that the maximum strain values of cortical and cancellous bone increased with lower bone density. Strain distribution is the first and foremost dependent on the shape of bone and architecture of cancellous bone. The geometry of the implant, thread patterns, grade of bone atrophy and BIC all affect the displacement and micro-strain on the mandible bone. Preoperative finite element analysis could offer improved predictability in the long-term outlook of dental implant restorations.

Bone Anchors—A Preliminary Finite Element Study of Some Factors Affecting Pullout

Bone anchors (or suture anchors) are used to provide attachment points for sutures to connect tissue such as tendons or ligaments to bone, and work by engaging a threaded portion—sometimes tapered—to the cancellous and/or cortical bone. Such repair is often needed after trauma, or as part of reconstructive surgery. This paper uses the finite element method to compare the pullout characteristics of one common type of bone anchor in different cancellous bone structures. Finite element models are created by using computed tomography (CT) scans of cancellous bone and building computer-aided design (CAD) models to define the cancellous bone geometry. Orthopedic surgeons will sometimes remove parts of the cortical shell and this paper also examines the mechanical effects of decortication. Furthermore, the importance of the connection between anchor and cortical layer is examined. One of the key outcomes from the model is that the coefficient of friction between bone and anchor determines potential mechanisms of pullout. The stiffness of anchors and the effect of the cortical layer are presented for different pullout angles to obtain the theoretical response. The results show the detailed modeling that includes the micro-architecture of the cancellous bone is necessary to capture the large variations that can exist.

Influence of object location in different FOVs on trabecular bone microstructure measurements of human mandible: a cone beam CT study

The aim of this study was to assess the influence of different object locations in different fields of view (FOVs) of two cone beam CT (CBCT) systems on trabecular bone microstructure measurements of a human mandible. A block of dry human mandible was scanned at five different locations (centre, left, right, anterior and posterior) using five different FOVs of two CBCT systems (NewTom™ 5G; QR Verona, Verona, Italy and Accuitomo 170; Morita, Kyoto, Japan). Image analysis software (CTAn software v. 1.1; SkyScan, Kontich, Belgium) was used to assess the trabecular bone microstructural parameters (thickness, Tb.Th; spacing, Tb.Sp; number, Tb.N; bone volume density, BV/TV). All measurements were taken twice by one trained observer. Tb.Th, Tb.Sp and Tb.N varied significantly across different FOVs in the NewTom 5G (p < 0.001) and the Accuitomo 170 (p < 0.001). For location, a significant difference was observed only when measuring BV/TV (p = 0.03) using the NewTom 5G. The trabecular bone microstructural measurements obtained from CBCT systems are influenced by the size of FOVs. Not all trabecular bone parameters measured using different CBCT systems are affected when varying the object location within the FOVs.

Patient-specific modelling of bone and bone-implant systems: the challenges

In the past three decades, finite element (FE) modelling has provided considerable understanding to the area of musculoskeletal biomechanics. However, most of this understanding has been generated using generic, standardised or idealised models. Patient-specific modelling (PSM) is almost never used for making clinical decisions. Imaging technologies have made it possible to create patient-specific geometries and FE meshes for modelling. While these have brought us closer to PSM, several challenges associated with the definition of material properties, loads, boundary conditions and interaction between components still need to be overcome. This study reviews the current status of PSM with respect to defining material behaviour and prescribing boundary conditions and interactions. With regard to the constitutive modelling of bone, it is seen that imaging is being increasingly used to define elastic properties (isotropic as well as anisotropic). However, the post-elastic and time-dependent behaviour, important for several modelling situations, is mostly obtained from in vitro experiments. Strain-based plasticity, not commonly available in FE codes, appears to have the potential of reducing an element of patient-specificity in modelling the yielding behaviour of bone. PSM of real boundary conditions that include muscles and ligaments continues to remain a challenge; many clinically relevant questions can be, however, answered without their inclusion. Simulation techniques to undertake PSM of interactions between bone and uncemented implants are available. Interference fit employed in both joint replacement fracture treatments induces considerable preload whose inclusion in models is important for the prediction of interface behaviour.