The Effect of Material Heterogeneity, Element Type, and Down-Sampling on Trabecular Stiffness in Micro Finite Element Models

A statistical analysis of human talar shape and bone density distribution

The shape of the talus, its internal structure, and its mechanical properties are important in determining talar behavior during loading, which may be significant for the design of surgical tools and implants. Although recent studies using statistical shape modeling have described quantitative talar external shape variation, no similar quantitative study exists to describe the density distribution of internal talar structure. The goal of this study is to quantify statistical variation in talar shape and density to benefit the design of talar implants. To this end, weight-bearing computed tomography (CT) scans of the ankle were collected in neutral, bilateral standing posture, and three-dimensional models were generated for each talus. Local density derived from the Hounsfield unit of each CT voxel was extracted. A weighted spherical harmonic analysis was performed to quantify the talar external shape. One hundred and seventy-nine volumes of interest were placed in the same relative position within each talus to quantify the talar density. Additionally, a finite element analysis (FEA) was conducted on a talus with both heterogeneous and homogeneous material properties to observe the effect of these properties on the stress and strain response. Significant differences were found in the talar density in sex and age, as well as in the stress and strain response between homogeneous and heterogeneous FEA. These differences show the importance of considering heterogeneity when examining the load response of tarsal bones.

Experimental DVC validation of heterogeneous micro finite element models applied to subchondral trabecular bone of the humeral head

Subchondral trabecular bone (STB) undergoes adaptive changes during osteoarthritic (OA) disease progression. These changes alter both the mineralization patterns and structure of bone and may contribute to variations in the mechanical properties. Similarly, when images are downsampled - as is often performed in micro finite element model (microFEM) generation - the morphological and mineralization patterns may further alter the mechanical properties due to partial volume effects. MicroFEMs accounting for material heterogeneity can account for these tissue variations, but no studies have validated these with robust full-field testing methods. As such, this study compared homogeneous and heterogeneous microFEMs to experimentally loaded trabecular bone cores from the humeral head combined with digital volume correlation (DVC). These microFEMs were used to compare apparent mechanical properties between normal and OA STB. Morphological and mineralization patterns between groups were also compared. There were no significant differences in tissue or bone mineral density between groups. The only significant differences in morphometric parameters were in trabecular thickness between groups. There were no significant differences in linear regression parameters between normal and OA STB apparent mechanical properties estimated using heterogeneous microFEMs with an element-wise bilinear elastic-plastic constitutive model. Clinical significance: Validated heterogeneous microFEMs applied to STB of the humeral head have the potential to significantly improve our understanding of mechanical variations in the bone that occur during OA progression.

Fracture Risk Evaluation of Bone Metastases: A Burning Issue

Major progress has been achieved to treat cancer patients and survival has improved considerably, even for stage-IV bone metastatic patients. Locomotive health has become a crucial issue for patient autonomy and quality of life. The centerpiece of the reflection lies in the fracture risk evaluation of bone metastasis to guide physician decision regarding physical activity, antiresorptive agent prescription, and local intervention by radiotherapy, surgery, and interventional radiology. A key mandatory step, since bone metastases may be asymptomatic and disseminated throughout the skeleton, is to identify the bone metastasis location by cartography, especially within weight-bearing bones. For every location, the fracture risk evaluation relies on qualitative approaches using imagery and scores such as Mirels and spinal instability neoplastic score (SINS). This approach, however, has important limitations and there is a need to develop new tools for bone metastatic and myeloma fracture risk evaluation. Personalized numerical simulation qCT-based imaging constitutes one of these emerging tools to assess bone tumoral strength and estimate the femoral and vertebral fracture risk. The next generation of numerical simulation and artificial intelligence will take into account multiple loadings to integrate movement and obtain conditions even closer to real-life, in order to guide patient rehabilitation and activity within a personalized-medicine approach.

Mechanical regulation of bone formation and resorption around implants in a mouse model of osteopenic bone

Although mechanical stimulation is considered a promising approach to accelerate implant integration, our understanding of load-driven bone formation and resorption around implants is still limited. This lack of knowledge may delay the development of effective loading protocols to prevent implant loosening, especially in osteoporosis. In healthy bone, formation and resorption are mechanoregulated processes. In the intricate context of peri-implant bone regeneration, it is not clear whether bone (re)modelling can still be load-driven. Here, we investigated the mechanical control of peri-implant bone (re)modelling with a well-controlled mechanobiological experiment. We applied cyclic mechanical loading after implant insertion in tail vertebrae of oestrogen depleted mice and we monitored peri-implant bone response by in vivo micro-CT. Experimental data were combined with micro-finite element simulations to estimate local tissue strains in (re)modelling locations. We demonstrated that a substantial increase in bone mass around the implant could be obtained by loading the entire bone. This augmentation could be attributed to a large reduction in bone resorption rather than to an increase in bone formation. We also showed that following implantation, mechanical regulation of bone (re)modelling was transiently lost. Our findings should help to clarify the role of mechanical stimulation on the maintenance of peri-implant bone mass.

MRI-based assessment of proximal femur strength compared to mechanical testing

Half of the women who sustain a hip fracture would not qualify for osteoporosis treatment based on current DXA-estimated bone mineral density criteria. Therefore, a better approach is needed to determine if an individual is at risk of hip fracture from a fall. The objective of this study was to determine the association between radiation-free MRI-derived bone strength and strain simulations compared to results from direct mechanical testing of cadaveric femora. Imaging was conducted on a 3-Tesla MRI scanner using two sequences: one balanced steady-state free precession sequence with 300 μm isotropic voxel size and one spoiled gradient echo with anisotropic voxel size of 234 × 234 × 1500 μm. Femora were dissected free of soft-tissue and 4350-ohm strain-gauges were securely applied to surfaces at the femoral shaft, inferior neck, greater trochanter, and superior neck. Cadavers were mechanically tested with a hydraulic universal test frame to simulate loading in a sideways fall orientation. Sideways fall forces were simulated on MRI-based finite element meshes and bone stiffness, failure force, and force for plastic deformation were computed. Simulated bone strength metrics from the 300 μm isotropic sequence showed strong agreement with experimentally obtained values of bone strength, with stiffness (r = 0.88, p = 0.0002), plastic deformation point (r = 0.89, p < 0.0001), and failure force (r = 0.92, p < 0.0001). The anisotropic sequence showed similar trends for stiffness, plastic deformation point, and failure force (r = 0.68, 0.70, 0.84; p = 0.02, 0.01, 0.0006, respectively). Surface strain-gauge measurements showed moderate to strong agreement with simulated magnitude strain values at the greater trochanter, superior neck, and inferior neck (r = -0.97, -0.86, 0.80; p ≤0.0001, 0.003, 0.03, respectively). The findings from this study support the use of MRI-based FE analysis of the hip to reliably predict the mechanical competence of the human femur in clinical settings.

Morphological and Apparent-Level Stiffness Variations Between Normal and Osteoarthritic Bone in the Humeral Head

Osteoarthritis (OA) is characterized by morphological changes that alter bone structure and mechanical properties. This study compared bone morphometric parameters and apparent modulus between humeral heads excised from end-stage OA patients undergoing total shoulder arthroplasty (n = 28) and non-pathologic normal cadavers (n = 28). Morphometric parameters were determined in central cores, with regional variations compared in four medial to lateral regions. Linear regression compared apparent modulus, morphometric parameters, and age. Micro finite element models estimated trabecular apparent modulus and derived density-modulus relationships. Significant differences were found for bone volume fraction (p < 0.001) and trabecular thickness (p < 0.001) in the most medial regions. No significant differences occurred between morphometric parameters and apparent modulus or age, except in slope between groups for apparent modulus versus trabecular number (p = 0.021), and in intercept for trabecular thickness versus age (p = 0.040). Significant differences occurred in both slope and intercept between density-modulus regression fits for each group (p ≤ 0.001). The normal group showed high correlations in the power-fit (r²  = 0.87), with a lower correlation (r²  = 0.61) and a more linear relationship, in the OA group. This study suggests that alterations in structure and apparent modulus persist mainly in subchondral regions of end-stage OA bone. As such, if pathologic regions are removed during joint replacement, computational models that utilize modeling parameters from non-pathologic normal bone may be applied to end-stage OA bone. An improved understanding of humeral trabecular bone variations has potential to improve the surgical management of end-stage OA patients. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 38:503-509, 2020.