Explicit finite element analysis can predict the mechanical response of conical implant press-fit in homogenized trabecular bone

3-D finite element model of the impaction of a press-fitted femoral stem under various biomechanical environments

BACKGROUND

Uncemented femoral stem insertion into the bone is achieved by applying successive impacts on an inserter tool called "ancillary". Impact analysis has shown to be a promising technique to monitor the implant insertion and to improve its primary stability.

METHOD

This study aims to provide a better understanding of the dynamic phenomena occurring between the hammer, the ancillary, the implant and the bone during femoral stem insertion, to validate the use of impact analyses for implant insertion monitoring. A dynamic 3-D finite element model of the femoral stem insertion via an impaction protocol is proposed. The influence of the trabecular bone Young's modulus (Et), the interference fit (IF), the friction coefficient at the bone-implant interface (μ) and the impact velocity (v0) on the implant insertion and on the impact force signal is evaluated.

RESULTS

For all configurations, a decrease of the time difference between the two first peaks of the impact force signal is observed throughout the femoral stem insertion, up to a threshold value of 0.23 ms. The number of impacts required to reach this value depends on Et, v0 and IF and varies between 3 and 8 for the set of parameters considered herein. The bone-implant contact ratio reached after ten impacts varies between 60% and 98%, increases as a function of v0 and decreases as a function of IF, μ and Et.

CONCLUSION

This study confirms the potential of an impact analyses-based method to monitor implant insertion and to retrieve bone-implant contact properties.

A convolutional neural network-based method for the generation of super-resolution 3D models from clinical CT images

BACKGROUND AND OBJECTIVE

The accurate evaluation of bone mechanical properties is essential for predicting fracture risk based on clinical computed tomography (CT) images. However, blurring and noise in clinical CT images can compromise the accuracy of these predictions, leading to incorrect diagnoses. Although previous studies have explored enhancing trabecular bone CT images to super-resolution (SR), none of these studies have examined the possibility of using clinical CT images from different instruments, typically of lower resolution, as a basis for analysis. Additionally, previous studies rely on 2D SR images, which may not be sufficient for accurate mechanical property evaluation, due to the complex nature of the 3D trabecular bone structures. The objective of this study was to address these limitations.

METHODS

A workflow was developed that utilizes convolutional neural networks to generate SR 3D models across different clinical CT instruments. The morphological and finite-element-derived mechanical properties of these SR models were compared with ground truth models obtained from micro-CT scans.

RESULTS

A significant improvement in analysis accuracy was demonstrated, where the new SR models increased the accuracy by up to 700 % compared with the low-resolution data, i.e. clinical CT images. Additionally, we found that the mixture of different CT image datasets may improve the SR model performance.

CONCLUSIONS

SR images, generated by convolutional neural networks, outperformed clinical CT images in the determination of morphological and mechanical properties. The developed workflow could be implemented for fracture risk prediction, potentially leading to improved diagnoses and subsequent clinical decision making.

Hydroxyapatite particle shape affects screw attachment in cancellous bone when augmented with hydroxyapatite-containing hydrogels

Screw-bone construct failures are a true challenge in orthopaedic implant fixation, particularly in poor quality bone. Whilst augmentation with bone cement can improve the primary stability of screws, there are cements, e.g. PMMA, that may impede blood flow and nutrients and hamper bone remodelling. In this study, soft, non-setting biomaterials based on Hyalectin gels and hydroxyapatite (HA) particles with different morphological parameters were evaluated as potential augmentation materials, using a lapine ex vivo bone model. The pull-out force, stiffness, and work to fracture were considered in evaluating screw attachment. The pull-out force of constructs reinforced with Hyalectin containing irregularly shaped nano-HA and spherically shaped micro-HA particles were found to be significantly higher than the control group (no augmentation material). The pull-out stiffness increased for the micro-HA particles and the work to fracture increased for the irregular nano-HA particles. However, there were no significant augmentation effect found for the spherical shaped nano-HA particles. In conclusion, injectable Hyalectin gel loaded with hydroxyapatite particles was found to have a potentially positive effect on the primary stability of screws in trabecular bone, depending on the HA particle shape and size.

Finite element analysis of the angle range in trans-inferior alveolar nerve implantation at the mandibular second molar

BACKGROUND

Trans- inferior alveolar nerve (IAN) implantation technique was wildly used while the potential appropriate angle range in which the residual alveolar bone can bear the stress without absorption are currently unclear. This study aimed to evaluate the stress distribution pattern of the interface between bone and implant by finite element analysis (FEA) to determine the appropriate range of the implant tilt angle.

METHODS

Cone beam computed tomography (CBCT) images of 120 patients with missing mandibular second molars and vertical bone height < 9 mm in the edentulous area were selected. The distances from the mandibular nerve canal to the buccal cortex, the lingual cortex and the alveolar ridge crest were measured by using a combination of software. The angular ranges of the buccal-lingual inclination of simulated trans-IAN implants were measured and three-dimensional finite element models were constructed in the mandibular second molar area according to the differences of the inclination angles. A vertical load (200N) was then applied to analyze the biomechanical conditions of the implant-bone interface during median occlusion.

RESULTS

The distance at the second molar from the nerve canal to the buccal cortex, lingual cortex and alveolar crest were 6.861 ± 1.194 mm, 2.843 ± 0.933 mm and 7.944 ± 0.77 mm. Trans-IAN implantation was feasible in 73.33% of patients. The minimum angle and maximum angles of the buccal-lingual inclination of the simulated implant were 19.135 ± 6.721° and 39.282 ± 6.581°. When a vertical static load of 200N was applied, the tensile stress in cortical bone gradually increased with the increase of the implant tilt angle. When the inclination angle reached 30°, the tensile stress (105.9 MPa) exceeded the yield strength (104 MPa) of cortical bone. Compared with the conventional implants, the stress peak value of the vertical ultra-short implant in cortical bone was greater than the stress peak value of the conventional implants at 10°(79.81 MPa) and 20°(82.83 MPa) and was smaller than the stress of the implant at 30°(105.9 MPa) and 40°(107.8 MPa). Therefore, when the bone mass allows, conventional-length implants should be selected whenever possible, and an operative range of the trans-IAN implantation in the mandibular second molar could be retained with an inclination angle of < 30°.

CONCLUSIONS

The mandibular nerve canal at the mandibular second molar was obviously biased to the lingual side, which ensured sufficient bone mass at the buccal side. In most patients with severe mandibular atrophy, it was possible to maintain a safe distance from the nerve canal with conventional-length implants via the trans-IAN implantation technique.

The application of an isotropic crushable foam model to predict the femoral fracture risk

For biomechanical simulations of orthopaedic interventions, it is imperative to implement a material model that can realistically reproduce the nonlinear behavior of the bone structure. However, a proper material model that adequately combines the trabecular and cortical bone response is not yet widely identified. The current paper aims to investigate the possibility of using an isotropic crushable foam (ICF) model dependent on local bone mineral density (BMD) for simulating the femoral fracture risk. The elastoplastic properties of fifty-nine human femoral trabecular cadaveric bone samples were determined and combined with existing cortical bone properties to characterize two forms of the ICF model, a continuous and discontinuous model. Subsequently, the appropriateness of this combined material model was evaluated by simulating femoral fracture experiments, and a comparison with earlier published results of a softening Von-Mises (sVM) material model was made. The obtained mechanical properties of the trabecular bone specimens were comparable to previous findings. Furthermore, the ultimate failure load predicted by the simulations of femoral fractures was on average 79% and 90% for the continuous and discontinuous forms of the ICF model and 82% of the experimental value for the sVM material model. Also, the fracture locations predicted by ICF models were comparable to the experiments. In conclusion, a nonlinear material model dependent on BMD was characterized for human femoral bone. Our findings indicate that the ICF model could predict the femoral bone strength and reproduce the variable fracture locations in the experiments.

Detection of periprosthetic fractures around the femoral stem by resonance frequency analysis: An in vitro study

Periprosthetic femoral bone fractures are frequent complications of Total Hip Arthroplasty (THA) and may occur during the insertion of uncemented Femoral Stems (FS), due to the nature of the press-fit fixation. Such fracture may lead to the surgical failure of the THA and require a revision surgery, which may have dramatic consequences. Therefore, an early detection of intra-operative fractures is important to avoid worsening the fracture and/or to enable a peroperative treatment. The aim of this in vitro study is to determine the sensitivity of a method based on resonance frequency analysis of the bone-stem-ancillary system for periprosthetic fractures detection. A periprosthetic fracture was artificially created close to the lesser-trochanter of 10 femoral bone mimicking phantoms. The bone-stem-ancillary resonance frequencies in the range (2-12) kHz were measured on an ancillary instrumented with piezoelectric sensors, which was fixed to the femoral stem. The measurements were repeated for different fracture lengths from 4 to 55 mm. The results show a decrease of the resonance frequencies due to the fracture occurrence and propagation. The frequency shift reached up to 170 Hz. The minimum fracture length that can be detected varies from 3.1±1.7 mm to 5.9±1.9 mm according to the mode and to the specimen. A significantly higher sensitivity (p = 0.011) was obtained for a resonance frequency around 10.6 kHz, corresponding to a mode vibrating in a plane perpendicular to the fracture. This study opens new paths toward the development of non-invasive vibration-based methods for intra-operative periprosthetic fractures detection.

Influence of the biomechanical environment on the femoral stem insertion and vibrational behavior: a 3-D finite element study

The long-term success of cementless surgery strongly depends on the implant primary stability. The femoral stem initial fixation relies on multiple geometrical and material factors, but their influence on the biomechanical phenomena occurring during the implant insertion is still poorly understood, as they are difficult to quantify in vivo. The aim of the present study is to evaluate the relationship between the resonance frequencies of the bone-implant-ancillary system and the stability of the femoral stem under various biomechanical environments. The interference fit IF, the trabecular bone Young's modulus [Formula: see text] and the bone-implant contact friction coefficient [Formula: see text] are varied to investigate their influence on the implant insertion phenomena and on the system vibration behavior. The results exhibit for all the configurations, a nonlinear increase in the bone-implant contact throughout femoral stem insertion, until the proximal contact is reached. While the pull-out force increases with [Formula: see text], IF and [Formula: see text], the bone-implant contact ratio decreases, which shows that a compromise on the set of parameters could be found in order to achieve the largest bone-implant contact while maintaining sufficient pull-out force. The modal analysis on the range [2-7] kHz shows that the resonance frequencies of the bone-implant-ancillary system increase with the bone-implant contact ratio and the trabecular bone Young's modulus, with a sensitivity that varies over the modes. Both the pull-out forces and the vibration behavior are consistent with previous experimental studies. This study demonstrates the potential of using vibration methods to guide the surgeons for optimizing implant stability in various patients and surgical configurations.

Development of robust finite element models to investigate the stability of osteochondral grafts within porcine femoral condyles

Osteoarthritis (OA) is the most prevalent chronic rheumatic disease worldwide with knee OA having an estimated lifetime risk of approximately 14%. Autologous osteochondral grafting has demonstrated positive outcomes in some patients, however, understanding of the biomechanical function and how treatments can be optimised remains limited. Increased short-term stability of the grafts allows cartilage surfaces to remain congruent prior to graft integration. In this study methods for generating specimen specific finite element (FE) models of osteochondral grafts were developed, using parallel experimental data for calibration and validation. Experimental testing of the force required to displace osteochondral grafts by 2 mm was conducted on three porcine knees, each with four grafts. Specimen specific FE models of the hosts and grafts were created from registered μCT scans captured from each knee (pre- and post-test). Material properties were based on the μCT background with a conversion between μCT voxel brightness and Young's modulus. This conversion was based on the results of the separate testing of eight porcine condyles and optimization of specimen specific FE models. The comparison between the experimental and computational push-in forces gave a strong agreement with a concordance correlation coefficient (CCC) = 0.75, validating the modelling approach. The modelling process showed that homogenous material properties based on whole bone BV/TV calculations are insufficient for accurate modelling and that an intricate description of the density distribution is required. The robust methodology can provide a method of testing different treatment options and can be used to investigate graft stability in full tibiofemoral joints.

Predicting the trabecular bone apparent stiffness tensor with spherical convolutional neural networks

The apparent stiffness tensor is relevant for characterizing trabecular bone quality. Previous studies have used morphology-stiffness relationships for estimating the apparent stiffness tensor. In this paper, we propose to train spherical convolutional neural networks (SphCNNs) to estimate this tensor. Information of the edges, trabecular thickness, and spacing are summarized in functions on the unitary sphere used as inputs for the SphCNNs. The concomitant dimensionality reduction makes it possible to train neural networks on relatively small datasets. The predicted tensors were compared to the stiffness tensors computed by using the micro-finite element method (μFE), which was considered as the gold standard, and models based on fourth-order fabric tensors. Combining edges and trabecular thickness yields significant improvements in the accuracy compared to the methods based on fourth-order fabric tensors. From the results, SphCNNs are promising for replacing the more expensive μFE stiffness estimations.

•

Characteristic stiffness tensors are derived from trabecular bone micro-CT samples.

•

Previous approximation methods fall short on heterogeneous data-sets.

•

The gradient, trabecular thickness and spacing are mapped to a spherical domain.

•

Spherical convolutional neural networks are used for the prediction.

•

The prediction error is significantly reduced compared to the state-of-the-art.

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.

Determinants of the primary stability of cementless acetabular cup implants: A 3D finite element study

Primary stability of cementless implants is crucial for the surgical success and long-term stability. However, primary stability is difficult to quantify in vivo and the biomechanical phenomena occurring during the press-fit insertion of an acetabular cup (AC) implant are still poorly understood. The aim of this study is to investigate the influence of the cortical and trabecular bone Young's moduli E_c and E_t, the interference fit IF and the sliding friction coefficient of the bone-implant interface μ on the primary stability of an AC implant. For each parameter combination, the insertion of the AC implant into the hip cavity and consequent pull-out are simulated with a 3D finite element model of a human hemi-pelvis. The primary stability is assessed by determining the polar gap and the maximum pull-out force. The polar gap increases along with all considered parameters. The pull-out force shows a continuous increase with E_c and E_t and a non-linear variation as a function of μ and IF is obtained. For μ > 0.6 and IF > 1.4 mm the primary stability decreases, and a combination of smaller μ and IF lead to a better fixation. Based on the patient's bone stiffness, optimal combinations of μ and IF can be identified. The results are in good qualitative agreement with previous studies and provide a better understanding of the determinants of the AC implant primary stability. They suggest a guideline for the optimal choice of implant surface roughness and IF based on the patient's bone quality.

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.