Development of a crushable foam model for human trabecular bone

The effect of bone plasticity models on simulations of primary fixation in total knee arthroplasty

Predictions of primary fixation in total knee arthroplasty (TKA) can aid in implant design, optimizing long-term survival. Finite element (FE) simulations are commonly used to assess micromotions at the bone-implant interface during daily activities, requiring accurate computational models. A key factor is the material model used to simulate bone properties. This study evaluated two plastic material models-Isotropic Crushable Foam (ICF) and softening Von-Mises (sVM)-for predicting primary fixation in femoral TKA components. Mechanical tests on human femoral trabecular bone samples under cyclic loading were replicated using FE simulations with ICF and sVM models. These models were then applied to FE simulations of three femoral TKA reconstructions, representing patients aged 57, 73, and 90 years. The ICF model replicated the gradual plastic deformation observed in experiments, unlike the sVM model, and more closely matched the permanent deformation of bone samples. In primary fixation simulations, micromotions at the bone-implant interface averaged 27 µm with ICF and 17 µm with sVM. While both predictions fall within acceptable ranges, the ICF model, as a pressure-dependent material model, more accurately replicates experimental energy dissipation and plastic deformation patterns, closely mirroring human bone's plastic behavior. This makes it better suited for simulating implant-bone interface micromotions in primary TKA fixation.

Development of a continuum-based, meshless, finite element modeling approach for representation of trabecular bone indentation

Implant subsidence into the underlying trabecular bone is a common problem in orthopaedic surgeries; however, the ability to pre-operatively predict implant subsidence remains limited. Current state-of-the-art computational models for predicting subsidence have issues addressing this clinical problem, often resulting from the size and complexity of existing subject-specific, image-based finite element (FE) models. The current study aimed to develop a simplified approach to FE modeling of subject-specific trabecular bone indentation resulting from implant penetration. Confined indentation experiments of human trabecular bone with flat- and sharp-tip indenters were simulated using FE analysis. A generalized continuum-level approach using a meshless smoothed particle hydrodynamics (SPH) approach and an isotropic crushable foam (CF) material model was developed for the trabecular bone specimens. Five FE models were generated with CF material parameters calibrated to cadaveric specimens spanning a range of bone mineral densities (BMD). Additionally, an alternative model configuration was developed that included consideration of bone marrow, with bone and marrow material parameters assigned to elements randomly according to bone volume (BV%) measurements of experimental specimens, owing to the non-uniform nature of trabecular bone tissue microstructure. Statistical analysis found significant correlation between the shapes of the numerical and experimental force-displacement curves. FE models accurately captured the bone densification patterns observed experimentally. Inclusion of marrow elements offered improved response prediction of the flat-tip indenter tests. Ultimately, the developed approach demonstrates the ability of a generalizable continuum-level SPH approach to capture bone variability using clinical bone imaging metrics without needing detailed image-based geometries, a significant step towards simplified subject-specific modeling of implant subsidence.

Comparative finite element analysis of posterior short segment fixation constructs with or without intermediate screws in the fractured vertebrae for the treatment of type a thoracolumbar fracture

Six-screw short-segment posterior fixation for thoracolumbar fractures, which involves intermediate screws at the fractured vertebrae has been proposed to reduce the rates of kyphosis recurrence and implant failure. Yet, little is known about the mechanisms and biomechanical responses by which intermediate screws at the fracture vertebrae enhance fixation strength. The objective of this study was to investigate the biomechanical properties that are associated with the augmentation of intermediate screws in relation to the severity of type A thoracolumbar fracture using finite element analysis. Short-segment stabilization models with or without augmentation screws at fractured vertebrae were established based on finite element model of moderate compressive fractures, severe compressive fractures and burst fractures. The spinal stiffness, stresses at the implanted hardware, and axial displacement of the bony defect were measured and compared under mechanical loading conditions. All six-screw stabilization showed a decreased range of motion in extension, lateral bending, and axial rotation compared to the traditional four-screw fixation models. Burst thoracolumbar fracture benefited more from augmentation of intermediate screws at the fracture vertebrae. The stress of the rod in six-screw models increased while decreased that of pedicle screws. Our results suggested that patients with more unstable fractures might achieve greater benefits from augmentation of intermediate screws at the fracture vertebrae. Augmentation of intermediate screws at the fracture vertebrae is recommended for patients with higher wedge-shaped or burst fractures to reduce the risk of hardware failure and postoperative re-collapse of injured vertebrae.

The effect of periprosthetic bone loss on the failure risk of tibial total knee arthroplasty

The effect of long-term periprosthetic bone loss on the process of aseptic loosening of tibial total knee arthroplasty (TKA) is subject to debate. Contradicting studies can be found in literature, reporting either bone resorption or bone formation before failure of the tibial tray. The aim of the current study was to investigate the effects of bone resorption on failure of tibial TKA, by simulating clinical postoperative bone density changes in finite element analysis (FEA) models and FEA models were created of two tibiae representing cases with good and poor initial bone quality which were subjected to a walking configuration and subsequently to a traumatic stumbling load. Bone failure was simulated using a crushable foam model incorporating progressive yielding. Repetitive loading under a level walking load did not result in failure of the periprosthetic bone in neither the good nor poor bone quality tibia at the baseline bone densities. When applying a stumble load, a collapse of the tibial reconstruction was noticed in the poor bone quality model. Incorporating postoperative bone loss led to a significant increase of the failure risk, particularly for the poor bone quality model in which subsidence of the tibial component was substantial. Our results suggest bone loss can lead to an increased risk of a collapse of the tibial component, particularly in case of poor bone quality at the time of surgery. The study also examined the probability of medial or lateral subsidence of the implant and aimed to improve clinical implications. The FEA model simulated plastic deformation of the bone and implant subsidence, with further validation required via mechanical experiments.

The Porosity Design and Deformation Behavior Analysis of Additively Manufactured Bone Scaffolds through Finite Element Modelling and Mechanical Property Investigations

Additively manufactured synthetic bone scaffolds have emerged as promising candidates for the replacement and regeneration of damaged and diseased bones. By employing optimal pore architecture, including pore morphology, sizes, and porosities, 3D-printed scaffolds can closely mimic the mechanical properties of natural bone and withstand external loads. This study aims to investigate the deformation pattern exhibited by polymeric bone scaffolds fabricated using the PolyJet (PJ) 3D printing technique. Cubic and hexagonal closed-packed uniform scaffolds with porosities of 30%, 50%, and 70% are utilized in finite element (FE) models. The crushable foam plasticity model is employed to analyze the scaffolds' mechanical response under quasi-static compression. Experimental validation of the FE results demonstrates a favorable agreement, with an average percentage error of 12.27% ± 7.1%. Moreover, the yield strength and elastic modulus of the scaffolds are evaluated and compared, revealing notable differences between cubic and hexagonal closed-packed designs. The 30%, 50%, and 70% porous cubic pore-shaped bone scaffolds exhibit significantly higher yield strengths of 46.89%, 58.29%, and 66.09%, respectively, compared to the hexagonal closed-packed bone scaffolds at percentage strains of 5%, 6%, and 7%. Similarly, the elastic modulus of the 30%, 50%, and 70% porous cubic pore-shaped bone scaffolds is 42.68%, 59.70%, and 58.18% higher, respectively, than the hexagonal closed-packed bone scaffolds at the same percentage strain levels. Furthermore, it is observed in comparison with our previous study the μSLA-printed bone scaffolds demonstrate 1.5 times higher elastic moduli and yield strengths compared to the PJ-printed bone scaffolds.

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.

Application of vertebral body compression osteotomy in pedicle subtraction osteotomy on ankylosing spondylitis kyphosis: Finite element analysis and retrospective study

BACKGROUND

Ankylosing spondylitis (AS) is a chronic inflammatory rheumatic disease, with pathological characteristics of bone erosion, inflammation of attachment point, and bone ankylosis. Due to the ossified intervertebral disc and ligament, pedicle subtraction osteotomy (PSO) is one of the mainstream surgeries of AS-related thoracolumbar kyphosis, but the large amount of blood loss and high risk of instrumental instability limit its clinical application. The purpose of our study is to propose a new transpedicular vertebral body compression osteotomy (VBCO) in PSO to reduce blood loss and improve stability.

METHODS

A retrospective analysis was performed on patients with AS-related thoracolumbar kyphosis who underwent one-level PSO in our hospital from February 2009 to May 2019. A total of 31 patients were included in this study; 6 received VBCO and 25 received eggshell vertebral body osteotomy. We collected demographic data containing gender and age at diagnosis. Surgical data contained operation time, estimated blood loss (EBL), and complications. Radiographic data contained pre-operative and follow-up sagittal parameters including chin brow-vertical angle (CBVA), global kyphosis (GK), thoracic kyphosis (TK), and lumbar lordosis (LL). A typical case with L2-PSO was used to establish a finite element model. The mechanical characteristics of the internal fixation device, vertebral body, and osteotomy plane of the two osteotomy models were analyzed under different working conditions.

RESULTS

The VBCO could provide comparable restoring of CBVA, GK, TK, and LL in the eggshell osteotomy procedure (all p > 0.05). The VBCO significantly reduced EBL compared to those with eggshell osteotomy [800.0 ml (500.0-1,439.5 ml) vs. 1,455.5 ml (1,410.5-1,497.8 ml), p = 0.033]. Compared with the eggshell osteotomy, VBCO showed better mechanical property. For the intra-pedicular screw fixation, the VBCO group had a more average distributed and lower stress condition on both nails and connecting rod. VBCO had a flattened osteotomy plane than the pitted osteotomy plane of the eggshell group, showing a lower and more average distributed maximum stress and displacement of osteotomy plane.

CONCLUSION

In our study, we introduced VBCO as an improved method in PSO, with advantages in reducing blood loss and providing greater stability. Further investigation should focus on clinical research and biomechanical analysis for the application of VBCO.