Finite element analysis of the hip and spine based on quantitative computed tomography

Biomechanical perspectives on image-based hip fracture risk assessment: advances and challenges

Hip fractures pose a significant health challenge, particularly in aging populations, leading to substantial morbidity and economic burden. Most hip fractures result from a combination of osteoporosis and falls. Accurate assessment of hip fracture risk is essential for identifying high-risk individuals and implementing effective preventive strategies. Current clinical tools, such as the Fracture Risk Assessment Tool (FRAX), primarily rely on statistical models of clinical risk factors derived from large population studies. However, these tools often lack specificity in capturing the individual biomechanical factors that directly influence fracture susceptibility. Consequently, image-based biomechanical approaches, primarily leveraging dual-energy X-ray absorptiometry (DXA) and quantitative computed tomography (QCT), have garnered attention for their potential to provide a more precise evaluation of bone strength and the impact forces involved in falls, thereby enhancing risk prediction accuracy. Biomechanical approaches rely on two fundamental components: assessing bone strength and predicting fall-induced impact forces. While significant advancements have been made in image-based finite element (FE) modeling for bone strength analysis and dynamic simulations of fall-induced impact forces, substantial challenges remain. In this review, we examine recent progress in these areas and highlight the key challenges that must be addressed to advance the field and improve fracture risk prediction.

Effect of septal extension graft on nasal tip support: A finite element analysis

BACKGROUND

Septal extension graft (SEG) is an effective method to control the projection, rotation, and shape of the nasal tip. However, the structural mechanics of SEG have not yet been adequately determined.

OBJECTIVES

The purpose of this study was to examine the effect of SEG parameters on nasal tip support using finite element analysis.

METHODS

A multicomponent nasal model was constructed from a computed tomographic scan. A control model without graft and a total of 15 models with different SEGs were created, regarding the direction, length, width, and piece of SEG. The nasal tip compression was simulated to analyze the von Mises stress, reaction force, and strain energy of the tip structure.

RESULTS

The SEG increased the max stress, reaction force, and strain energy of the nasal tip compared to the normal control. The SEG perpendicular to the nasal dorsum resulted in the highest maximum stress, reaction force, and strain energy for the same size of SEG. With the length increasing from 15 × 8 × 1 mm to 25 × 8 × 1 mm, the reaction force remained relatively stable, but the stress on the graft reduced significantly. Adding the width and pieces of the SEG increased the reaction force and strain energy of the tip.

CONCLUSION

The placement of SEG can strengthen the nasal tip support. The direction, length, width, and piece of SEG have an impact on the mechanics.

LEVEL OF EVIDENCE

Diagnostic, III.

Intradiscal cement leakage (ICL) increases the stress on adjacent vertebrae after kyphoplasty for osteoporotic vertebra compression fracture (OVCF): a finite-element study

This study aimed to explore the biomechanical effects on adjacent vertebra of thoracolumbar Osteoporotic Vertebra Compression Fracture (OVCF) after Percutaneous Kyphoplasty (PKP) with intraoperative intradiscal cement leakage (ICL) by applying a Finite-Element Analysis. We collected pre- and post-operative computer tomography (CT) images of a 71-year-old female patient with single T12 OVCF, who underwent an intraoperative cement leakage into the T12-L1 disc. Three-dimensional finite element models of thoracolumbar spine (T10-L2) were built with the support of Materialise Interactive Medical Image Control System (MIMICS) and ABAQUS software. The stress on adjacent vertebrae and endplates under the uniform compressive pressure (0.3 MPa) and during different loading moments were analyzed. The three-dimensional finite element models reveal an asymmetrical distribution of von Mises stresses on the adjacent endplate unaffected by the surgical intervention. The maximum von Mises stress on adjacent vertebral bodies increased during different loading conditions, especially for lateral bending and rotation loading conditions, whereas the maximum von Mises stress on distal non-treated vertebrae decreased on anteflexion and backward extension loading conditions. Post-operative adjacent vertebra compression fractures after PKP with intraoperative intradiscal cement leakage (ICL) may be closely related to the biomechanical changes of adjacent vertebrae of thoracolumbar OVCF, and it may increase the risk of postoperative fracture.

Imaging techniques to study diabetic bone disease

PURPOSE OF REVIEW

This review article presents the most recent research on bone fragility in individuals with diabetes from a medical imaging perspective.

RECENT FINDINGS

The widespread availability of dual-energy X-ray absorptiometry (DXA) and trabecular bone score (TBS) software has led to recent assessments of bone fragility with this texture parameter in several studies of type 2 diabetes mellitus (T2D), but in few of type 1 diabetes mellitus (T1D). Although most studies show a trend of reduced TBS values in T2D independent of areal bone mineral density (aBMD) of the lumbar spine, some studies also show the limitations of TBS in both T2D and T1D. Given the limitations of DXA to assess bone strength and investigate the etiology of bone fragility in diabetes, more investigators are incorporating three-dimensional (3D) medical imaging techniques in their studies. Recent use of 3D medical imaging to assess bone fragility in the setting of diabetes has been mostly limited to a few cross-sectional studies predominantly incorporating high-resolution peripheral quantitative computed tomography (HR-pQCT). Although HR-pQCT studies indicate higher tibial cortical porosity in subjects with T2D, results are inconsistent in T1D due to differences in study designs, sample sizes, and subject characteristics, among other factors. With respect to central CT, recent studies support a previous finding in the literature indicating femoral neck geometrical impairments in subjects with T2D and provide encouraging results for the incorporation of finite element analysis (FEA) to assess bone strength in studies of T2D. In the recent literature, there are no studies assessing bone fragility in T1D with QCT, and only two studies used pQCT reporting tibial and radial impairments in young women and children with T1D, respectively. Magnetic resonance imaging (MRI) has not been recently used in diabetic studies of bone fragility.

SUMMARY

As bone fragility in diabetes is not explained by DXA-derived aBMD and given the limitations of cross-sectional studies, it is imperative to use 3D imaging techniques for longitudinal assessments of the density, quality, and microenvironment of bone to improve our understanding of the effects of diabetes on bone and reduce the risk of fracture in this large and vulnerable population of subjects with diabetes.

Patient-Specific Finite Element Modeling of the Whole Lumbar Spine Using Clinical Routine Multi-Detector Computed Tomography (MDCT) Data-A Pilot Study

(1) Background: To study the feasibility of developing finite element (FE) models of the whole lumbar spine using clinical routine multi-detector computed tomography (MDCT) scans to predict failure load (FL) and range of motion (ROM) parameters. (2) Methods: MDCT scans of 12 subjects (6 healthy controls (HC), mean age ± standard deviation (SD): 62.16 ± 10.24 years, and 6 osteoporotic patients (OP), mean age ± SD: 65.83 ± 11.19 years) were included in the current study. Comprehensive FE models of the lumbar spine (5 vertebrae + 4 intervertebral discs (IVDs) + ligaments) were generated (L1−L5) and simulated. The coefficients of correlation (ρ) were calculated to investigate the relationship between FE-based FL and ROM parameters and bone mineral density (BMD) values of L1−L3 derived from MDCT (BMDQCT-L1-3). Finally, Mann−Whitney U tests were performed to analyze differences in FL and ROM parameters between HC and OP cohorts. (3) Results: Mean FE-based FL value of the HC cohort was significantly higher than that of the OP cohort (1471.50 ± 275.69 N (HC) vs. 763.33 ± 166.70 N (OP), p < 0.01). A strong correlation of 0.8 (p < 0.01) was observed between FE-based FL and BMDQCT-L1-L3 values. However, no significant differences were observed between ROM parameters of HC and OP cohorts (p = 0.69 for flexion; p = 0.69 for extension; p = 0.47 for lateral bending; p = 0.13 for twisting). In addition, no statistically significant correlations were observed between ROM parameters and BMDQCT- L1-3. (4) Conclusions: Clinical routine MDCT data can be used for patient-specific FE modeling of the whole lumbar spine. ROM parameters do not seem to be significantly altered between HC and OP. In contrast, FE-derived FL may help identify patients with increased osteoporotic fracture risk in the future.

Testosterone therapy and bone quality in men with diabetes and hypogonadism: Study design and protocol

Context

Type 2 diabetes mellitus (T2D) is often accompanied by male hypogonadism and both conditions are associated with increased risk for fractures. Testosterone (T) has been shown to improve the bone health of hypogonadal men but has not been tested in patients who also have T2D in addition to low T. To date, there is no treatment that is specifically recommended for bone disease among patients with T2D. This study will evaluate the effect of T therapy on the bone health of male veterans with low T who also have T2D.

Methods

This is a randomized double-blind placebo-controlled trial of 166 male veterans 35–65 years old, with T2D and hypogonadism, randomized to either T gel 1.62% or placebo for 12 months. We will evaluate the effect of T therapy on the following primary outcomes:1) changes in bone strength as measured by microfinite elements analysis (μFEA) using high-resolution peripheral quantitative computer tomography, 2) changes in bone turnover markers, and 3) changes in circulating osteoblast progenitors (COP) and osteoclast precursors cells.

Discussion

We anticipate that T therapy will result in improvement in bone strength owing to improvement in bone remodeling through an increase in osteoblastic differentiation and proliferation in patients with hypogonadism and T2D.

MDCT-Based Finite Element Analyses: Are Measurements at the Lumbar Spine Associated with the Biomechanical Strength of Functional Spinal Units of Incidental Osteoporotic Fractures along the Thoracolumbar Spine?

Assessment of osteoporosis-associated fracture risk during clinical routine is based on the evaluation of clinical risk factors and T-scores, as derived from measurements of areal bone mineral density (aBMD). However, these parameters are limited in their ability to identify patients at high fracture risk. Finite element models (FEMs) have shown to improve bone strength prediction beyond aBMD. This study aims to investigate whether FEM measurements at the lumbar spine can predict the biomechanical strength of functional spinal units (FSUs) with incidental osteoporotic vertebral fractures (VFs) along the thoracolumbar spine. Multi-detector computed tomography (MDCT) data of 11 patients (5 females and 6 males, median age: 67 years) who underwent MDCT twice (median interval between baseline and follow-up MDCT: 18 months) and sustained an incidental osteoporotic VF between baseline and follow-up scanning were used. Based on baseline MDCT data, two FSUs consisting of vertebral bodies and intervertebral discs (IVDs) were modeled: one standardly capturing L1-IVD-L2-IVD-L3 (FSU_L1-L3) and one modeling the incidentally fractured vertebral body at the center of the FSU (FSU_F). Furthermore, volumetric BMD (vBMD) derived from MDCT, FEM-based displacement, and FEM-based load of the single vertebrae L1 to L3 were determined. Statistically significant correlations (adjusted for a BMD ratio of fracture/L1-L3 segments) were revealed between the FSU_F and mean load of L1 to L3 (r = 0.814, p = 0.004) and the mean vBMD of L1 to L3 (r = 0.745, p = 0.013), whereas there was no statistically significant association between the FSU_F and FSU_L1-L3 or between FSU_F and the mean displacement of L1 to L3 (p > 0.05). In conclusion, FEM measurements of single vertebrae at the lumbar spine may be able to predict the biomechanical strength of incidentally fractured vertebral segments along the thoracolumbar spine, while FSUs seem to predict only segment-specific fracture risk.

Establishment and Simulation of 3D Geometric Models of Mini-Pig and Sheep Knee Joints Using Finite Element Analysis

Background

Our objective was to establish and compare three-dimensional models of knee joints of mini-pigs and sheep, the 2 most commonly used animal models of osteoarthritis.

Material/Methods

Three-dimensional geometric models of knee joints were used to assess their biomechanical properties by analysis of the three-dimensional finite element stress load for flexion at 30° and 60°.

Results

Analysis of multiple tissues indicated that the sheep knee had greater stress peaks than the mini-pig knee at 30° flexion (range: 12.5 to 30.4 Mpa for sheep vs. 11.1 to 20.2 Mpa for mini-pig) and at 60° flexion (range: 17.9 to 43.5 Mpa for sheep vs. 15.9 to 28.9 Mpa for mini-pig). In addition, there was uneven distribution of stress loads in the surrounding ligaments during flexion.

Conclusions

Our three-dimensional finite element analysis indicated that the mini-pig knee joint had stress values and changes of cartilage, meniscus, and peripheral ligaments that were similar to those of the human knee.

Case-specific non-linear finite element models to predict failure behavior in two functional spinal units

Current finite element (FE) models predicting failure behavior comprise single vertebrae, thereby neglecting the role of the posterior elements and intervertebral discs. Therefore, this study aimed to develop a more clinically relevant, case-specific non-linear FE model of two functional spinal units able to predict failure behavior in terms of (i) the vertebra predicted to fail; (ii) deformation of the specimens; (iii) stiffness; and (iv) load to failure. For this purpose, we also studied the effect of different bone density-mechanical properties relationships (material models) on the prediction of failure behavior. Twelve two functional spinal units (T6-T8, T9-T11, T12-L2, and L3-L5) with and without artificial metastases were destructively tested in axial compression. These experiments were simulated using CT-based case-specific non-linear FE models. Bone mechanical properties were assigned using four commonly used material models. In 10 of the 11 specimens our FE model was able to correctly indicate which vertebrae failed during the experiments. However, predictions of the three-dimensional deformations of the specimens were less promising. Whereas stiffness of the whole construct could be strongly predicted (R²  = 0.637-0.688, p < 0.01), we obtained weak correlations between FE predicted and experimentally determined load to failure, as defined by the total reaction force exhibiting a drop in force (R²  = 0.219-0.247, p > 0.05). Additionally, we found that the correlation between predicted and experimental fracture loads did not strongly depend on the material model implemented, but the stiffness predictions did. In conclusion, this work showed that, in its current state, our FE models may be used to identify the weakest vertebra, but that substantial improvements are required in order to quantify in vivo failure loads. © 2018 The Authors. Journal of Orthopaedic Research® Published by Wiley Periodical, Inc. on behalf of Orthopaedic Research Society. J Orthop Res 36:3208-3218, 2018.

Automation of a DXA-based finite element tool for clinical assessment of hip fracture risk

Finite element analysis of medical images is a promising tool for assessing hip fracture risk. Although a number of finite element models have been developed for this purpose, none of them have been routinely used in clinic. The main reason is that the computer programs that implement the finite element models have not been completely automated, and heavy training is required before clinicians can effectively use them. By using information embedded in clinical dual energy X-ray absorptiometry (DXA), we completely automated a DXA-based finite element (FE) model that we previously developed for predicting hip fracture risk. The automated FE tool can be run as a standalone computer program with the subject's raw hip DXA image as input. The automated FE tool had greatly improved short-term precision compared with the semi-automated version. To validate the automated FE tool, a clinical cohort consisting of 100 prior hip fracture cases and 300 matched controls was obtained from a local community clinical center. Both the automated FE tool and femoral bone mineral density (BMD) were applied to discriminate the fracture cases from the controls. Femoral BMD is the gold standard reference recommended by the World Health Organization for screening osteoporosis and for assessing hip fracture risk. The accuracy was measured by the area under ROC curve (AUC) and odds ratio (OR). Compared with femoral BMD (AUC = 0.71, OR = 2.07), the automated FE tool had a considerably improved accuracy (AUC = 0.78, OR = 2.61 at the trochanter). This work made a large step toward applying our DXA-based FE model as a routine clinical tool for the assessment of hip fracture risk. Furthermore, the automated computer program can be embedded into a web-site as an internet application.

Peripheral Quantitative Computed Tomography (pQCT) Measures Contribute to the Understanding of Bone Fragility in Older Patients With Low-trauma Fracture

Dual-energy X-ray absorptiometry (DXA) as currently used has limitations in identifying patients with osteoporosis and predicting occurrence of fracture. We aimed to express peripheral quantitative computed tomography (pQCT) variables of patients with low-trauma fracture as T-scores by using T-score scales obtained from healthy young women, and to evaluate the potential clinical utility of pQCT for the assessment of bone fragility. Fracture patients were recruited from a fracture liaison service at the Royal Melbourne Hospital. Reference pQCT data were obtained from studies on women's health conducted by our group. A study visit was arranged with fracture patients, during which DXA and pQCT were applied to measure their bone strength. A total of 59 fracture patients were recruited, and reference data were obtained from 78 healthy young females. All DXA variables and most pQCT variables were significantly different between healthy young females and fracture patients (p < 0.05), except polar stress-strain index (p = 0.34) and cortical bone density (p = 0.19). Fracture patients were divided into osteoporosis and non-osteoporosis groups according to their DXA T-scores. Significant differences were observed in most pQCT variables (p < 0.05), except trabecular area and cortical density (p > 0.9 and p = 0.5, respectively). By applying pQCT T-scores, 11 (27%) of patients who were classified as having low or medium risk of osteoporosis on DXA T-scores alone were reclassified as high risk. Results of logistic regression suggested trabecular bone density as an independent predictor of osteoporosis status. More patients can be identified with osteoporosis by applying pQCT T-score variables in older people with low-trauma fracture. Peripheral QCT T-scores contribute to the understanding of bone fragility in this population.

A biomechanical sorting of clinical risk factors affecting osteoporotic hip fracture

Osteoporotic fracture has been found associated with many clinical risk factors, and the associations have been explored dominantly by evidence-based and case-control approaches. The major challenges emerging from the studies are the large number of the risk factors, the difficulty in quantification, the incomplete list, and the interdependence of the risk factors. A biomechanical sorting of the risk factors may shed lights on resolving the above issues. Based on the definition of load-strength ratio (LSR), we first identified the four biomechanical variables determining fracture risk, i.e., the risk of fall, impact force, bone quality, and bone geometry. Then, we explored the links between the FRAX clinical risk factors and the biomechanical variables by looking for evidences in the literature. To accurately assess fracture risk, none of the four biomechanical variables can be ignored and their values must be subject-specific. A clinical risk factor contributes to osteoporotic fracture by affecting one or more of the biomechanical variables. A biomechanical variable represents the integral effect from all the clinical risk factors linked to the variable. The clinical risk factors in FRAX mostly stand for bone quality. The other three biomechanical variables are not adequately represented by the clinical risk factors. From the biomechanical viewpoint, most clinical risk factors are interdependent to each other as they affect the same biomechanical variable(s). As biomechanical variables must be expressed in numbers before their use in calculating LSR, the numerical value of a biomechanical variable can be used as a gauge of the linked clinical risk factors to measure their integral effect on fracture risk, which may be more efficient than to study each individual risk factor.

Bone Imaging and Fracture Risk after Spinal Cord Injury

Spinal cord injury (SCI) is characterized by marked bone loss and an increased risk of fracture with high complication rate. Recent research based on advanced imaging analysis, including quantitative computed tomography (QCT) and patient-specific finite element (FE) modeling, has provided new and important insights into the magnitude and temporal pattern of bone loss, as well as the associated changes to bone structure and strength, following SCI. This work has illustrated the importance of early therapeutic treatment to prevent bone loss after SCI and may someday serve as the basis for a clinical fracture risk assessment tool for the SCI population. This review provides an update on the epidemiology of fracture after SCI and discusses new findings and significant developments related to bone loss and fracture risk assessment in the SCI population based on QCT analysis and patient-specific FE modeling.

Finite element analysis for prediction of bone strength

Finite element (FE) analysis has been applied for the past 40 years to simulate the mechanical behavior of bone. Although several validation studies have been performed on specific anatomical sites and load cases, this study aims to review the predictability of human bone strength at the three major osteoporotic fracture sites quantified in recently completed in vitro studies at our former institute. Specifically, the performance of FE analysis based on clinical computer tomography (QCT) is compared with the ones of the current densitometric standards, bone mineral content, bone mineral density (BMD) and areal BMD (aBMD). Clinical fractures were produced in monotonic axial compression of the distal radii, vertebral sections and in side loading of the proximal femora. QCT-based FE models of the three bones were developed to simulate as closely as possible the boundary conditions of each experiment. For all sites, the FE methodology exhibited the lowest errors and the highest correlations in predicting the experimental bone strength. Likely due to the improved CT image resolution, the quality of the FE prediction in the peripheral skeleton using high-resolution peripheral CT was superior to that in the axial skeleton with whole-body QCT. Because of its projective and scalar nature, the performance of aBMD in predicting bone strength depended on loading mode and was significantly inferior to FE in axial compression of radial or vertebral sections but not significantly inferior to FE in side loading of the femur. Considering the cumulated evidence from the published validation studies, it is concluded that FE models provide the most reliable surrogates of bone strength at any of the three fracture sites.