Quantitative Computed Tomography Protocols Affect Material Mapping and Quantitative Computed Tomography-Based Finite-Element Analysis Predicted Stiffness

Design and validation of finite element models for the assessment of post-fixation distal femur fracture motion

Fracture site motion is thought to play an important role in the healing of complex fractures of the distal femur via mechanotransduction. Measuring this motion in vivo is challenging, and this has led researchers to turn to finite element modeling approaches to gain insights into the mechanical environment at the fracture site. Developing a systematic understanding of the effect of different model choices for distal femur fractures may allow more accurate prediction of fracture site motion from these types of simulations. In this study, we aim to assess the effect of four different modeling choices and parameters. We looked at the effect of using bone specific density distributions vs generic values, employing landmark-based geometry generation, varying fracture alignment within clinically relevant ranges, and determining whether direct apposition of the fracture to the plate was achieved. For validation, five cadaveric femurs had fractures created and repaired with plated constructs, and these were then loaded and fracture site motion was directly measured. We found that using landmark based bone geometry and patient-specific bone density distributions had a minimal effect on the overall model predictions. Changing the alignment, particularly into varus and procurvatum could have a large (>50%) effect on predicted shear motion, as could direct apposition of the bone to the plate. These findings demonstrate that modeling choices can play an important role in simulating distal femur fracture mechanics, and it is particularly critical that patient customized models attempt to accurately represent alignment of the bone fragments and lateral plate apposition.

Image-based finite element model stiffness and vBMD by single and dual energy CT reconstruction kernel

Single-energy quantitative computed tomography (SEQCT) provides volumetric bone mineral density (vBMD) measures for bone analysis and input to image-based finite element models (FEMs). Dual-energy CT (DECT) improves vBMD by accounting for voxel-specific material variations utilizing scans at multiple x-ray energies. vBMD is also altered by reconstruction kernel that cannot be accounted for using calibration phantoms. This study compared vBMD and FEM stiffness derived from SEQCT and DECT images reconstructed with two common kernels. SEQCT and DECT images of cadaveric shoulders (n = 10) were collected using standard (STD) and boneplus (BONE) kernels. Hounsfield Units were converted to vBMD using specimen-specific calibrations. DECT STD and BONE images were generated using an established material decomposition method with 40 and 90 keV simulated monochromatic images. A proximal humerus bone section below the anatomic neck was used for vBMD analysis and FEM generation. FEMs were loaded to 1% apparent strain for stiffness measurements. Between STD and BONE kernel images, average vBMD differed 0.9 mgK2HPO4/cc and 4.1 mg K2HPO4/cc, in SEQCT and DECT images, respectively. Significant differences occurred in DECT images (p = 0.001). BONE reconstructed images produced higher vBMD measures across both SEQCT and DECT images. The difference between STD and BONE in both SEQCT- and DECT-based FEMs persisted, with larger estimated stiffness in BONE models. For six of the models DECT-based had higher stiffness than SEQCT-based models using the same kernel, although these models differed between STD and BONE kernels. Differences in stiffness between STD and BONE derived models were similar across image types (DECT: 17.5 kN/mm; SEQCT: 19.0 kN/mm). Stiffness values were significantly different within SECT kernels and between SEQCT BONE and DECT STD models. This study shows important differences in vBMD and FEM stiffness that occur due to CT-based imaging parameters alone. These results indicate that consistent imaging parameters should be used for vBMD analysis and FEM input to avoid systematic measurement errors.

Evaluation of bone-related mechanical properties in female patients with long-term remission of Cushing's syndrome using quantitative computed tomography-based finite element analysis

BACKGROUND

Hypercortisolism in Cushing's syndrome (CS) is associated with bone loss, skeletal fragility, and altered bone quality. No studies evaluated bone geometric and strain-stress values in CS patients after remission thus far.

PATIENTS AND METHODS

Thirty-two women with CS in remission (mean age [±SD] 51 ± 11; body mass index [BMI], 27 ± 4 kg/m2; mean time of remission, 120 ± 90 months) and 32 age-, BMI-, and gonadal status-matched female controls. Quantitative computed tomography (QCT) was used to assess volumetric bone mineral density (vBMD) and buckling ratio, cross-sectional area, and average cortical thickness at the level of the proximal femur. Finite element (FE) models were generated from QCT to calculate strain and stress values (maximum principal strain [MPE], maximum strain energy density [SED], maximum Von Mises [VM], and maximum principal stress [MPS]). Areal BMD (aBMD) and trabecular bone score (TBS) were assessed by dual-energy X-ray absorptiometry (2D DXA).

RESULTS

Trabecular vBMD at total hip and trochanter were lower in CS as compared with controls (P < .05). Average cortical thickness was lower, and buckling ratio was greater in CS vs controls (P < .01). All strain and stress values were higher in CS patients vs controls (P < .05). 2D DXA-derived measures were similar between patients and controls (P > .05). Prior hypercortisolism predicted both VM (β .30, P = .014) and MPS (β .30, P = .015), after adjusting for age, BMI, menopause, delay to diagnosis, and duration of remission.

CONCLUSIONS

Women with prior hypercortisolism have reduced trabecular vBMD and impaired bone geometrical and mechanical properties, which may contribute to an elevated fracture risk despite long-term remission.

Biomechanical MRI detects reduced bone strength in subjects with vertebral fractures

Vertebral fracture is one of the most serious consequences of osteoporosis. Estimation of vertebral strength from magnetic resonance imaging (MRI) scans may provide a new approach for the prediction of vertebral fractures. To that end, we sought to establish a biomechanical MRI (BMRI) method to compute vertebral strength and test its ability to distinguish fracture from non-fracture subjects. This case-control study included 30 subjects without vertebral fractures and 15 subjects with vertebral fractures. All subjects underwent MRI with a mDIXON-Quant sequence and quantitative computed tomography (QCT), from which proton fat fraction-based bone marrow adipose tissue (BMAT) content and volumetric bone mineral density (vBMD) were measured, respectively. Nonlinear finite element analysis was applied to MRI and QCT scans of L2 vertebrae to compute vertebral strength (BMRI- and BCT-strength). The differences in BMAT content, vBMD, BMRI-strength and BCT-strength between the two groups were examined by t-tests. Receiver operating characteristic (ROC) analysis was performed to assess the ability of each measured parameter to distinguish fracture from non-fracture subjects. Results showed that the fracture group had 23 % lower BMRI-strength (P < .001) and 19 % higher BMAT content (P < .001) than the non-fracture group, whereas no significant difference in vBMD was detected between the two groups. A poor correlation was found between vBMD and BMRI-strength (R² = 0.33). Compared to vBMD and BMAT content, BMRI- and BCT-strength had the larger area under the curve (0.82 and 0.84, respectively) and provided better sensitivity and specificity in separating fracture from non-fracture subjects. In conclusion, BMRI is capable of detecting reduced bone strength in patients with vertebral fracture, and may serve as a new approach for risk assessment of vertebral fracture.

Influence of image reconstruction kernel on computed tomography-based finite element analysis in the clinical opportunistic screening of osteoporosis-A preliminary result

PURPOSE

This study aimed to evaluate the difference in vertebral mechanical properties estimated by finite element analysis (FEA) with different computed tomography (CT) reconstruction kernels and evaluate their accuracy in the screening and classification of osteoporosis.

METHODS

There were 31 patients enrolled retrospectively from the quantitative CT database of our hospital, uniformly covering the range from osteoporosis to normal. All subjects' CT raw data were reconstructed both with a smooth standard convolution kernel (B40f) and a sharpening bone convolution kernel (B70f), and FEA was performed on L1 of each subject based on two reconstructed images to obtain vertebral estimated strength and stiffness. The trabecular volumetric bone mineral density (vBMD) of the same vertebral body was also measured. FEA measurements between two kernels and their accuracy for osteoporosis screening were compared.

RESULTS

The vertebral stiffness and strength measured in FEA-B40f were significantly lower compared with those of FEA-B70f (12.0%, p = 0.000 and 10.7%, p = 0.000, respectively). The correlation coefficient between FEA-B70F and vBMD was slightly higher than that of FEA-B40F in both vertebral strength and stiffness (strength: r ²-B40f = 0.21, p = 0.009 vs. r ²-B70f = 0.27, p = 0.003; stiffness: r ²-B40f = 0.37, p = 0.002 vs. r ²-B70f = 0.45, p=0.000). The receiver operator characteristic curve showed little difference in the classification of osteoporosis between FEA-B40f and FEA-B70f.

CONCLUSION

Two kernels both seemed to be applicable to the opportunistic screening of osteoporosis by CT-FEA despite variance in FE-estimated bone strength and bone stiffness. A protocol for CT acquisition and FEA is still required to guarantee the reproducibility of clinical use.

A Patient-Specific Fracture Risk Assessment Tool for Femoral Bone Metastases: Using the Bone Strength (BOS) Score in Clinical Practice

Patients with femoral metastases are at risk of fracturing bones. It is important to prevent fractures in order to maintain mobility and quality of life. The BOne Strength (BOS) score is based on a computed tomography (CT)-based patient-specific finite element (FE) computer model that objectively calculates bone strength. In this pilot study, the added clinical value of the BOS score towards treatment-related decision making was assessed. In December 2019, the BOS score was implemented in four radiotherapy centers. The BOS scores and fracture risks of individual patients were calculated and returned to the physician to assist in treatment decisions. The physicians filled out a questionnaire, which was qualitatively analyzed. A follow-up to identify fractures and/or death was performed after six months. Until June 2021, 42 BOS scores were delivered (20 high, 9 moderate, and 13 low fracture risk). In 48%, the BOS score led to an adaptation of treatment plans. Physicians indicated that the BOS score provided objective insight into fracture risk, was reassuring for physicians and patients, and improved multidisciplinary discussions and shared decision making. In conclusion, the BOS score is an objective tool to assess fracture risk in femoral bone metastases and aids physicians and patients in making a more informed decision regarding the most appropriate treatment.

Biomechanical analysis of vertebral wedge deformity in elderly women with quantitative CT-based finite element analysis

BACKGROUND

To explore the vertebral deformity angle (VD angle) of 1^st lumbar vertebral body (L1) in elderly women, investigate the influence of VD on vertebral stiffness (VS) by biomechanical analysis using quantitative computed tomography-based finite element analysis (QCT-FEA).

METHODS

Two hundred seventy eight participants were recruited, and underwent QCT scan. Measured VD angles of L1, and constructed QCT-FEA models of L1 with the minimum (0.59°), median (5.79°) and maximum (11.15°) VD angles, respectively. Loads in two directions were applied on the upper edge of L1 with a force of 700 N, and vertebral stiffness (VS) was defined as the ratio of 700 N and displacement at the superior reference point: (1) perpendicular to the upper edge of L1 (defined as VS-U); (2) perpendicular to the lower edge of L1(defined as VS-L).

RESULTS

Age was very weak positively correlated with VD angle, moderate negatively correlated with vBMD, and moderate negatively correlated with VS (P < 0.05). VS-U was significantly different among three VD angles, so was VS-L (P < 0.001). VS-U was higher than VS-L in 5.79° and 11.15° VD angles (P < 0.05), however no difference in 0.59° VD angles (P > 0.10).

CONCLUSIONS

VD angle of L1 was slightly increased with age and not correlated with vBMD, and VS was moderate negatively correlated with age, showing that the vertebral body was more likely to fracture with aging. VS-U and VS-L were gradually decreased with the increase of VD angle, and VS-L was lower than VS-U with the increase of VD angle, which showed that vertebral body was more prone to fracture when the load was perpendicular to the lower edge of the vertebral body as the VD angle increasing.

Sensitivity of the stress field of the proximal femur predicted by CT-based FE analysis to modeling uncertainties

Proximal femur anatomy and bone mineral density vary widely among individuals, precluding the use of one predefined finite element (FE) model to determine the stress field for all proximal femurs. This variability poses a challenge in current prosthetic hip design approach. Given the numerous options for generating computed tomography (CT)-based FE models, selecting the best methods for defining the mechanical behavior of the proximal femur is difficult. In this study, a combination of computational and experimental approaches was used to explore the susceptibility of the predicted stress field of the proximal femur to different combinations of density-elasticity relationships, element type, element size, and calibration error. Our results suggest that FE models with first-order voxelized elements generated by the Keyak and Falkinstein density-elasticity relationship or quadratic tetrahedral elements generated by the Morgan density-elasticity relationship lead to accurate estimations of the mechanical behavior of human femurs. Other combinations of element size, element type, and mathematical relationships produce less accurate results, especially in the cortical bone of the femoral neck and calcar region. The voxelized model was more susceptible to variation of element size and density-elasticity relationships than FE models with quadratic tetrahedral elements. Regardless of element type, the stress fields predicted by the Keyak and Falkinstein and the Morgan relationships were the most robust to calibration error when deriving material density from CT-generated Hounsfield data. These results provide insight into the implementation of a robust platform for designing patient-specific implants capable of maintaining or modifying the stress in bones.

Material properties of human vertebral trabecular bone under compression can be predicted based on quantitative computed tomography

Background

The prediction of the stability of bones is becoming increasingly important. Especially osteoporotic vertebral body fractures are a growing problem and an increasing burden on the health system. Therefore, the aim of this study was to provide the best possible description of the relationship between the material properties of human vertebral trabecular bone measured under the most physiological conditions possible and the bone mineral density (BMD) determined by clinical quantitative computed tomography (QCT).

Methods

Forty eight cylindric cancellous bone samples with a diameter of 7.2 mm obtained from 13 human fresh-frozen lumbar vertebrae from 5 donors (3 men, 2 women) have been used for this study. After the specimens were temporarily reinserted into the vertebral body, the QCT was performed. For mechanical testing, the samples were embedded in a load-free manner using polymethylmetacrylate (PMMA). The surrounding test chamber was filled with phosphate buffered saline (PBS) and heated to 37 °C during the test. After 10 preconditioning load cycles, destructive testing was performed under axial compression. After determining the fracture site, BMD has been evaluated in this region only. Regression analyses have been performed.

Results

Fracture site had an average length of 2.4 (±1.4) mm and a position of 43.9 (±10.9) percent of the measurement length from the cranial end. No fracture reached the embedding. The average BMD at the fracture site was 80.2 (±28.7 | min. 14.5 | max. 137.8) mgCaHA/ml.

In summary the results of the regression analyses showed for all three parameters a very good quality of fit by a power regression.

Conclusion

The results of this study show that QCT-based bone density measurements have a good predictive power for the material properties of the vertebral cancellous bone measured under near to physiological conditions. The mechanical bone properties of vertebral cancellous bone could be modelled with high accuracy in the investigated bone density range.

Experimental testing and biomechanical CT analysis of Chinese cadaveric vertebrae with different modeling approaches

Osteoporosis is characterized by reduced bone strength predisposing to an increased risk of fracture. Biomechanical computed tomography (BCT), predicting bone strength via CT-based finite element analysis (FEA), is now clinically available in the USA for diagnosing osteoporosis or assessing fracture risk. However, it has not been previously validated using a cohort of only Chinese subjects. Additionally, the effect of various modeling approaches on BCT outcomes remains elusive. To address these issues, we performed DXA and QCT scanning, compression testing, and BCT analyses on thirteen vertebrae derived from Chinese donors. Three BCT models were created (voxBCT and tetBCT: voxel-based and tetrahedral element-based FE models generated by a commercial software; matBCT: tetrahedral element-based FE model generated by a custom MATLAB program). BCT-computed outcomes were compared with experimental measures or between different BCT models. Results showed that, DXA-measured areal bone mineral density (aBMD) showed weak correlations with experimentally-measured vertebral stiffness (R² = 0.28) and strength (R² = 0.34). Compared to DXA-aBMD, BCT-computed stiffness provided improved correlations with experimentally-measured stiffness (voxBCT: R² = 0.82; tetBCT: R² = 0.77; matBCT: R² = 0.76) and strength (voxBCT: R² = 0.55; tetBCT: R² = 0.57; matBCT: R² = 0.53); BCT-computed mechanical parameters (stiffness, stress and strain) of the three different models were highly correlated with each other, with coefficient of determination (R²) values of 0.89-0.98. These results, based on a cohort of Chinese vertebral cadavers, suggest that BCT is superior over aBMD to consistently predict vertebral mechanical characteristics, regardless of the modeling approaches of choice.

Evaluation of patient tissue selection methods for deriving equivalent density calibration for femoral bone quantitative CT analyses

Osteoporosis affects an increasing number of people every year and patient specific finite element analysis of the femur has been proposed to identify patients that could benefit from preventative treatment. The aim of this study was to demonstrate, verify, and validate an objective process for selecting tissues for use as the basis of phantomless calibration to enable patient specific finite element analysis derived hip fracture risk prediction. Retrospective reanalysis of patient computed tomography (CT) scans has the potential to yield insights into more accurate prediction of osteoporotic fracture. Bone mineral density (BMD) specific calibration scans are not typically captured during routine clinical practice. Tissue-based BMD calibration can therefore empower the retrospective study of patient CT scans captured during routine clinical practice. Together the method for selecting tissues as the basis for phantomless calibration coupled with the post-processing steps for deriving a calibration equation using the selected tissues provide an estimation of quantitative equivalent density results derived using calibration phantoms. Patient tissues from a retrospective cohort of 211 patients were evaluated. The best phantomless calibration resulted in a femoral strength (FS) [N] bias of 0.069 ± 0.07% over FS derived from inline calibration and a BMD [kg/cm³] bias of 0.038 ± 0.037% over BMD derived from inline calibration. The phantomless calibration slope for the best method presented was within the range of patient specific calibration curves available for comparison and demonstrated a small bias of 0.028 ± 0.054 HU/(mg/cm³), assuming the Mindways Model 3 BMD inline calibration phantom as the gold standard. The presented method of estimating a calibration equation from tissues showed promise for CT-based femoral fracture analyses of retrospective cohorts without readily available calibration data.

Subject-specific rib finite element models with material data derived from coupon tests under bending loading

Rib fractures are common thoracic injuries in motor vehicle crashes. Several human finite element (FE) human models have been created to numerically assess thoracic injury risks. However, the accurate prediction of rib biomechanical response has shown to be challenging due to human variation and modeling approaches. The main objective of this study was to better understand the role of modeling approaches on the biomechanical response of human ribs in anterior-posterior bending. Since the development of subject specific rib models is a time-consuming process, the second objective of this study was to develop an accurate morphing approach to quickly generate high quality subject specific rib meshes. The exterior geometries and cortical-trabecular boundaries of five human 6th-level ribs were extracted from CT-images. One rib mesh was developed in a parametric fashion and the other four ribs were developed with an in-house morphing algorithm. The morphing algorithm automatically defined landmarks on both the periosteal and endosteal boundaries of the cortical layer, which were used to morph the template nodes to target geometries. Three different cortical bone material models were defined based on the stress-strain data obtained from subject-specific tensile coupon tests for each rib. Full rib anterior-posterior bending tests were simulated based on data recorded in testing. The results showed similar trends to test data with some sensitivity relative to the material modeling approach. Additionally, the FE models were substantially more resistant to failure, highlighting the need for better techniques to model rib fracture. Overall, the results of this work can be used to improve the biofidelity of human rib finite element models.

Opportunistic Computed Tomography and Spine Surgery: A Narrative Review

STUDY DESIGN

Narrative review.

OBJECTIVE

This article seeks to provide a narrative review regarding the ability of opportunistic information available from computed tomography (CT) scans to guide decisions in spine surgery related to patient bone quality.

METHODS

A review of the literature (limited to human and English language) was performed via PubMed and Google Scholar using the search terms; "osteoporosis" AND "opportunistic" AND "computed tomography" AND "spine surgery." The titles and then abstracts of all identified citations were reviewed for inclusion by 2 of the authors (MS, BAF). Relevant articles were then studied in full text.

RESULTS

A review of the literature found 25 articles that were selected for inclusion in this narrative review. These articles were broadly divided into 4 subcategories: (1) opportunistic CT (oCT) and osteoporosis detection, (2) oCT data and the quality of screw fixation, (3) utilization of Hounsfield units to assess clinical and/or radiographic outcomes following spine fusion, and (4) virtual stress testing in spine surgery.

CONCLUSION

The literature on oCT, as well as associated virtual stress-testing techniques, demonstrate the potential to enhance spine surgery outcomes by preoperatively identifying at-risk patients in need of bone health optimization and informing best techniques for performing spinal fusion surgery on patients with diminished bone quality. While our narrative summary of the limited literature to date suggests a promising future for oCT data, significant additional research and/or radiographic workflow standardization is needed to validate these methods for clinical use.

Density-Dependent Material and Failure Criteria Equations Highly Affect the Accuracy and Precision of QCT/FEA-Based Predictions of Osteoporotic Vertebral Fracture Properties

About 700,000 vertebral fractures occur in the US as a result of bone loss. Quantitative computed tomography (QCT)-based finite element analysis (FEA) is a promising tool for fracture risk prediction that is becoming attractive in the clinical setting. The goals of this study were (1) to perform individual and pooled specimen optimization using inverse QCT/FEA modeling to obtain ash density-elastic modulus equations incorporating the whole vertebral body and accounting for all variables used during FE modeling, and (2) to determine the effect of material equations and failure criteria on the accuracy and precision of mechanical properties. Fifty-four (54) human vertebrae were used to optimize material equations based on experimental outcomes and, together with a previously proposed material equation, were implemented in our models using three different failure criteria to obtain fracture loads. Our robust QCT/FEA approach predicted 78% of the failure loads. Material equations resulted in poor accuracy in the predicted stiffness, yet yielded good precision and, more importantly, strong correlations with fracture loads. Both material and fracture criterion equations are equally important in estimating accurate and precise QCT/FEA predictions. Results suggest that both elastic modulus and fracture criterion equations should be validated against experimental outcomes to better explain the response of the tissue under various conditions.

Fracture risk assessment and clinical decision making for patients with metastatic bone disease

Metastatic breast, prostate, lung, and other cancers often affect bone, causing pain, increasing fracture risk, and decreasing function. Management of metastatic bone disease (MBD) is clinically challenging when there is potential but uncertain risk of pathological fracture. Management of MBD has become a major focus within orthopedic oncology with respect to fracture and impending fracture care. If impending skeletal-related events (SREs), particularly pathologic fracture, could be predicted, increasing evidence suggests that prophylactic surgical treatment improves patient outcomes. However, current fracture risk assessment and radiographic metrics do not have high accuracy and have not been combined with relevant patient survival tools. This review first explores the prevalence, incidence, and morbidity of MBD and associated SREs for different cancer types. Strengths and limitations of current fracture risk scoring systems for spinal stability and long bone fracture are highlighted. More recent computed tomography (CT)-based structural rigidity analysis (CTRA) and finite element (FE) analysis methods offer advantages of increased specificity (true negative rate), but are limited in availability. Other fracture prediction approaches including parametric response mapping and positron emission tomography/computed tomography measures show early promise. Substantial new information to inform clinical decision-making includes measures of survival, clinical benefits, and economic analysis of prophylactic treatment compared to after-fracture stabilization. Areas of future research include use of big data and machine learning to predict SREs, greater access and refinement of CTRA/FE approaches, combination of clinical survival prediction tools with radiographically based fracture risk assessment, and net benefit analysis for fracture risk assessment and prophylactic treatment.

Deciphering the "Art" in Modeling and Simulation of the Knee Joint: Overall Strategy

Recent explorations of knee biomechanics have benefited from computational modeling, specifically leveraging advancements in finite element analysis and rigid body dynamics of joint and tissue mechanics. A large number of models have emerged with different levels of fidelity in anatomical and mechanical representation. Adapted modeling and simulation processes vary widely, based on justifiable choices in relation to anticipated use of the model. However, there are situations where modelers' decisions seem to be subjective, arbitrary, and difficult to rationalize. Regardless of the basis, these decisions form the "art" of modeling, which impact the conclusions of simulation-based studies on knee function. These decisions may also hinder the reproducibility of models and simulations, impeding their broader use in areas such as clinical decision making and personalized medicine. This document summarizes an ongoing project that aims to capture the modeling and simulation workflow in its entirety-operation procedures, deviations, models, by-products of modeling, simulation results, and comparative evaluations of case studies and applications. The ultimate goal of the project is to delineate the art of a cohort of knee modeling teams through a publicly accessible, transparent approach and begin to unravel the complex array of factors that may lead to a lack of reproducibility. This manuscript outlines our approach along with progress made so far. Potential implications on reproducibility, on science, engineering, and training of modeling and simulation, on modeling standards, and on regulatory affairs are also noted.

Effect of metastatic lesion size and location on the load-bearing capacity of vertebrae using an optimized ash density-modulus equation

About 1.8 million new cancer cases are estimated in the US in 2019 from which 50-85% might metastasize to the thoracic and lumbar spines. Subject-specific quantitative computed tomography-based finite element analysis (QCT/FEA) is a promising used tool to predict vertebral fracture properties. The aims of this study were twofold: First, to develop an optimized equation for the elastic modulus accounting for all input parameters in FE modeling of fracture properties. Second, to assess the effect of lesion size and location on the predicted fracture loads. An inverse QCT/FEA method was implemented to determine optimal coefficients for the modulus equation as a function of ash density. Lesions of 16 and 20 mm were then virtually located at the center, off-centered, anterior, and posterior regions of the vertebrae. A total of 6426 QCT/FEA models were run to optimize the coefficients and evaluate the effect of lesions on fracture properties. QCT/FEA predicted stiffness showed high correlations (50%) with the experimentally measured values. Compared to a 16 mm lesion size, a 20 mm lesion had a reduction in failure load of 55%, 57%, 52%, and 44% at the center, off-centered, anterior cortex, and pedicle, respectively (p < 0.001). Lesions affecting mostly trabecular bone showed the largest reduction in predicted failure loads (about 55%), and females presented weaker outcomes than males. An optimal elastic modulus equation resulted in accurate vertebral stiffness predictions. A deterioration of the trabecular bone due to the presence of a lesion highly affected the predicted fracture loads, and this reduction was significantly higher in females compared to males.

Dual-energy estimates of volumetric bone mineral densities in the lumbar spine using quantitative computed tomography better correlate with fracture properties when compared to single-energy BMD outcomes

It is estimated that over 200 million people worldwide are affected by osteoporosis. Vertebral fracture risk prediction using dual energy x-ray absorptiometry (DXA) is confounded by limitations of the technology, such as 2D measurements of bone mineral density (BMD), inability to measure bone distribution and heterogeneity, and potential overestimations of BMD due to degenerative diseases. To overcome these shortcomings, single energy (SE) quantitative computed tomography (QCT) imaging estimates of Hounsfield units (HU) and volumetric BMD have been implemented as alternative methodologies for assessing fracture risk. However, marrow fat within the vertebrae can highly affect the vBMD and fracture properties estimations. To address this issue, 54 vertebrae were dissected from nine cadaveric spines and scanned using SE-QCT (120kVp) and dual energy (DE)-QCT (80/140 kVp), with the latter accounting for marrow fat within the vertebrae. The vertebrae were then scanned using DXA and subjected to mechanical testing to obtain fracture properties. aBMD outcomes from DXA showed a better correlation with DE-QCT vBMD versus SE outcomes [DE: aBMD vs. vBMD (R²: 0.61); SE: aBMD vs. vBMD (R²: 0.27)]. SE-QCT underestimated vertebral vBMD by -56% (p<0.0001) when compared to DE-QCT. vBMD estimates from SE-QCT could predict 45% and 37% of the vertebral failure loads and stiffness, respectively, compared to 67% and 46% from DE-QCT. DE-QCT vBMD outcomes highly correlated with fracture properties of vertebrae as compared to SE-QCT metrics. As DE scanning has the ability to correct for the effects of bone marrow fat, estimated vBMD from SE-QCT were significantly underestimated compared to DE-QCT. Dual energy CT scanning has the potential to more accurately predict vertebral failure and aid the clinician in the evaluation of appropriate interventions. Future studies should consider implementing DE-QCT in their fracture assessment.

The Influence of Reconstruction Kernel on Bone Mineral and Strength Estimates Using Quantitative Computed Tomography and Finite Element Analysis

Quantitative computed tomography has been posed as an alternative imaging modality to investigate osteoporosis. We examined the influence of computed tomography convolution back-projection reconstruction kernels on the analysis of bone quantity and estimated mechanical properties in the proximal femur. Eighteen computed tomography scans of the proximal femur were reconstructed using both a standard smoothing reconstruction kernel and a bone-sharpening reconstruction kernel. Following phantom-based density calibration, we calculated typical bone quantity outcomes of integral volumetric bone mineral density, bone volume, and bone mineral content. Additionally, we performed finite element analysis in a standard sideways fall on the hip loading configuration. Significant differences for all outcome measures, except integral bone volume, were observed between the 2 reconstruction kernels. Volumetric bone mineral density measured using images reconstructed by the standard kernel was significantly lower (6.7%, p < 0.001) when compared with images reconstructed using the bone-sharpening kernel. Furthermore, the whole-bone stiffness and the failure load measured in images reconstructed by the standard kernel were significantly lower (16.5%, p < 0.001, and 18.2%, p < 0.001, respectively) when compared with the image reconstructed by the bone-sharpening kernel. These data suggest that for future quantitative computed tomography studies, a standardized reconstruction kernel will maximize reproducibility, independent of the use of a quantitative calibration phantom.

Effect of different CT scanners and settings on femoral failure loads calculated by finite element models

In a multi-center patient study, using different CT scanners, CT-based finite element (FE) models are utilized to calculate failure loads of femora with metastases. Previous studies showed that using different CT scanners can result in different outcomes. This study aims to quantify the effects of (i) different CT scanners; (ii) different CT protocols with variations in slice thickness, field of view (FOV), and reconstruction kernel; and (iii) air between calibration phantom and patient, on Hounsfield Units (HU), bone mineral density (BMD), and FE failure load. Six cadaveric femora were scanned on four CT scanners. Scans were made with multiple CT protocols and with or without an air gap between the body model and calibration phantom. HU and calibrated BMD were determined in cortical and trabecular regions of interest. Non-linear isotropic FE models were constructed to calculate failure load. Mean differences between CT scanners varied up to 7% in cortical HU, 6% in trabecular HU, 6% in cortical BMD, 12% in trabecular BMD, and 17% in failure load. Changes in slice thickness and FOV had little effect (≤4%), while reconstruction kernels had a larger effect on HU (16%), BMD (17%), and failure load (9%). Air between the body model and calibration phantom slightly decreased the HU, BMD, and failure loads (≤8%). In conclusion, this study showed that quantitative analysis of CT images acquired with different CT scanners, and particularly reconstruction kernels, can induce relatively large differences in HU, BMD, and failure loads. Additionally, if possible, air artifacts should be avoided. © 2018 Orthopaedic Research Society. © 2018 The Authors. Journal of Orthopaedic Research® Published by Wiley Periodicals, Inc. on behalf of the Orthopaedic Research Society. J Orthop Res.

Effects of Scan Resolutions and Element Sizes on Bovine Vertebral Mechanical Parameters from Quantitative Computed Tomography-Based Finite Element Analysis

Quantitative computed tomography-based finite element analysis (QCT/FEA) has been developed to predict vertebral strength. However, QCT/FEA models may be different with scan resolutions and element sizes. The aim of this study was to explore the effects of scan resolutions and element sizes on QCT/FEA outcomes. Nine bovine vertebral bodies were scanned using the clinical CT scanner and reconstructed from datasets with the two-slice thickness, that is, 0.6 mm (PA resolution) and 1 mm (PB resolution). There were significantly linear correlations between the predicted and measured principal strains (R² > 0.7, P < 0.0001), and the predicted vertebral strength and stiffness were modestly correlated with the experimental values (R² > 0.6, P < 0.05). Two different resolutions and six different element sizes were combined in pairs, and finite element (FE) models of bovine vertebral cancellous bones in the 12 cases were obtained. It showed that the mechanical parameters of FE models with the PB resolution were similar to those with the PA resolution. The computational accuracy of FE models with the element sizes of 0.41 × 0.41 × 0.6 mm³ and 0.41 × 0.41 × 1 mm³ was higher by comparing the apparent elastic modulus and yield strength. Therefore, scan resolution and element size should be chosen optimally to improve the accuracy of QCT/FEA.