Morphology–elasticity relationships using decreasing fabric information of human trabecular bone from three major anatomical locations

Super-resolution of clinical CT: Revealing microarchitecture in whole bone clinical CT image data

Osteoporotic fractures, prevalent in the elderly, pose a significant health and economic burden. Current methods for predicting fracture risk, primarily relying on bone mineral density, provide only modest accuracy. If better spatial resolution of trabecular bone in a clinical scan were available, a more complete assessment of fracture risk would be obtained using microarchitectural measures of bone (i.e. trabecular thickness, trabecular spacing, bone volume fraction, etc.). However, increased resolution comes at the cost of increased radiation or can only be applied at small volumes of distal skeletal locations. This study explores super-resolution (SR) technology to enhance clinical CT scans of proximal femurs and better reveal the trabecular microarchitecture of bone. Using a deep-learning-based (i.e. subset of artificial intelligence) SR approach, low-resolution clinical CT images were upscaled to higher resolution and compared to corresponding MicroCT-derived images. SR-derived 2-dimensional microarchitectural measurements, such as degree of anisotropy, bone volume fraction, trabecular spacing, and trabecular thickness were within 16 % error compared to MicroCT data, whereas connectivity density exhibited larger error (as high as 1094 %). SR-derived 3-dimensional microarchitectural metrics exhibited errors <18 %. This work showcases the potential of SR technology to enhance clinical bone imaging and holds promise for improving fracture risk assessments and osteoporosis detection. Further research, including larger datasets and refined techniques, can advance SR's clinical utility, enabling comprehensive microstructural assessment across whole bones, thereby improving fracture risk predictions and patient-specific treatment strategies.

QCT-based computational bone strength assessment updated with MRI-derived 'hidden' microporosity

Microdamage accumulated through sustained periods of cyclic loading or single overloading events contributes to bone fragility through a reduction in stiffness and strength. Monitoring microdamage in vivo remains unattainable by clinical imaging modalities. As such, there are no established computational methods for clinical fracture risk assessment that account for microdamage that exists in vivo at any specific timepoint. We propose a method that combines multiple clinical imaging modalities to identify an indicative surrogate, which we term 'hidden porosity', that incorporates pre-existing bone microdamage in vivo. To do so, we use the third metacarpal bone of the equine athlete as an exemplary model for fatigue induced microdamage, which coalesces in the subchondral bone. N = 10 metacarpals were scanned by clinical quantitative computed tomography (QCT) and magnetic resonance imaging (MRI). We used a patch-based similarity method to quantify the signal intensity of a fluid sensitive MRI sequence in bone regions where microdamage coalesces. The method generated MRI-derived pseudoCT images which were then used to determine a pre-existing damage (Dpex) variable to quantify the proposed surrogate and which we incorporate into a nonlinear constitutive model for bone tissue. The minimum, median, and maximum detected Dpex of 0.059, 0.209, and 0.353 reduced material stiffness by 5.9%, 20.9%, and 35.3% as well as yield stress by 5.9%, 20.3%, and 35.3%. Limb-specific voxel-based finite element meshes were equipped with the updated material model. Lateral and medial condyles of each metacarpal were loaded to simulate physiological joint loading during gallop. The degree of detected Dpex correlated with a relative reduction in both condylar stiffness (p = 0.001, R2 > 0.74) and strength (p < 0.001, R2 > 0.80). Our results illustrate the complementary value of looking beyond clinical CT, which neglects the inclusion of microdamage due to partial volume effects. As we use clinically available imaging techniques, our results may aid research beyond the equine model on fracture risk assessment in human diseases such as osteoarthritis, bone cancer, or osteoporosis.

Hip joint load prediction using inverse bone remodeling with homogenized FE models: Comparison to micro-FE and influence of material modeling strategy

BACKGROUND AND OBJECTIVE

Measuring physiological loading conditions in vivo can be challenging, as methods are invasive or pose a high modeling effort. However, the physiological loading of bones is also imprinted in the bone microstructure due to bone (re)modeling. This information can be retrieved by inverse bone remodeling (IBR). Recently, an IBR method based on micro-finite-element (µFE) modeling was translated to homogenized-FE (hFE) to decrease computational effort and tested on the distal radius. However, this bone has a relatively simple geometry and homogeneous microstructure. Therefore, the objective of this study was to assess the agreement of hFE-based IBR with µFE-based IBR to predict hip joint loading from the head of the femur; a bone with more complex loading as well as more heterogeneous microstructure.

METHODS

hFE-based IBR was applied to a set of 19 femoral heads using four different material mapping laws. One model with a single homogeneous material for both trabecular and cortical volume and three models with a separated cortex and either homogeneous, density-dependent inhomogeneous, or density and fabric-dependent orthotropic material. Three different evaluation regions (full bone, trabecular bone only, head region only) were defined, in which IBR was applied. µFE models were created for the same bones, and the agreement of the predicted hip joint loading history obtained from hFE and µFE models was evaluated. The loading history was discretized using four unit load cases.

RESULTS

The computational time for FE solving was decreased on average from 500 h to under 1 min (CPU time) when using hFE models instead of µFE models. Using more information in the material model in the hFE models led to a better prediction of hip joint loading history. Inhomogeneous and inhomogeneous orthotropic models gave the best agreement to µFE-based IBR (RMSE% <14%). The evaluation region only played a minor role.

CONCLUSIONS

hFE-based IBR was able to reconstruct the dominant joint loading of the femoral head in agreement with µFE-based IBR and required considerably lower computational effort. Results indicate that cortical and trabecular bone should be modeled separately and at least density-dependent inhomogeneous material properties should be used with hFE models of the femoral head to predict joint loading.

Accuracy of osseointegrated screw-bone construct stiffness and peri-implant loading predicted by homogenized FE models relative to micro-FE models

Computational predictions of stiffness and peri-implant loading of screw-bone constructs are highly relevant to investigate and improve bone fracture fixations. Homogenized finite element (hFE) models have been used for this purpose in the past, but their accuracy has been questioned given the numerous simplifications, such as neglecting screw threads and modelling the trabecular bone structure as a continuum. This study aimed to investigate the accuracy of hFE models of an osseointegrated screw-bone construct when compared to micro-FE models considering the simplified screw geometry and different trabecular bone material models. Micro-FE and hFE models were created from 15 cylindrical bone samples with a virtually inserted, osseointegrated screw (fully bonded interface). Micro-FE models were created including the screw with threads (=reference models) and without threads to quantify the error due to screw geometry simplification. In the hFE models, the screws were modelled without threads and four different trabecular bone material models were used, including orthotropic and isotropic material derived from homogenization with kinematic uniform boundary conditions (KUBC), as well as from periodicity-compatible mixed uniform boundary conditions (PMUBC). Three load cases were simulated (pullout, shear in two directions) and errors in the construct stiffness and the volume average strain energy density (SED) in the peri-implant region were evaluated relative to the micro-FE model with a threaded screw. The pooled error caused by only omitting screw threads was low (max: 8.0%) compared to the pooled error additionally including homogenized trabecular bone material (max: 92.2%). Stiffness was predicted most accurately using PMUBC-derived orthotropic material (error: -0.7 ± 8.0%) and least accurately using KUBC-derived isotropic material (error: +23.1 ± 24.4%). Peri-implant SED averages were generally well correlated (R² ≥ 0.76), but slightly over- or underestimated by the hFE models and SED distributions were qualitatively different between hFE and micro-FE models. This study suggests that osseointegrated screw-bone construct stiffness can be predicted accurately using hFE models when compared to micro-FE models and that volume average peri-implant SEDs are well correlated. However, the hFE models are highly sensitive to the choice of trabecular bone material properties. PMUBC-derived isotropic material properties represented the best trade-off between model accuracy and complexity in this study.

The impact of the parameters of the constitutive model on the distribution of strain in the femoral head

The rapid spread of the finite element method has caused that it has become, among other methods, the standard tool for pre-clinical estimates of bone properties. This paper presents an application of this method for the calculation and prediction of strain and stress fields in the femoral head. The aim of the work is to study the influence of the considered anisotropy and heterogeneity of the modeled bone on the mechanical fields during a typical gait cycle. Three material models were tested with different properties of porous bone carried out in literature: a homogeneous isotropic model, a heterogeneous isotropic model, and a heterogeneous anisotropic model. In three cases studied, the elastic properties of the bone were determined basing on the Zysset-Curnier approach. The tensor of elastic constants defining the local properties of porous bone is correlated with a local porosity and a second order fabric tensor describing the bone microstructure. In the calculations, a model of the femoral head generated from high-resolution tomographic scans was used. Experimental data were drawn from publicly available database "Osteoporotic Virtual Physiological Human Project." To realistically reflect the load on the femoral head, main muscles were considered, and their contraction forces were determined based on inverse kinematics. For this purpose, the results from OpenSim packet were used. The simulations demonstrated that differences between the results predicted by these material models are significant. Only the anisotropic model allowed for the plausible distribution of stresses along the main trabecular groups. The outcomes also showed that the precise evaluation of the mechanical fields is critical in the context of bone tissue remodeling under mechanical stimulations.

The predictive ability of a QCT-FE model of the proximal femoral stiffness under multiple load cases is strongly influenced by experimental uncertainties

Despite significant improvements in terms of the predictive ability of Quantitative Computed Tomography based Finite Element (QCT-FE) models in estimating femoral strength (fracture load and stiffness), no substantial clinical adoption of this method has taken place to date. Narrowing the wide variability of FE results by standardizing the methodology and validation protocols, as well as reducing the uncertainties in the FEA process have been proposed as routes towards improved reliability. The aim of this study was to: First, validate a QCT-FE model of proximal femoral stiffness in multiple stance load cases, and second, using a parametric approach, determine the influence of select experimental and modeling parameters on the predictive ability of our model. Ten fresh frozen human femoral samples were tested in neutral stance, 15° adducted and 15° abducted load cases. Voxel-based linear-elastic QCT-FE models of the samples were generated to predict the models' stiffness values in all load cases. The base FE models were validated against the experimental results using linear regression. Thirty six deviated models were created using the minimum and maximum values of experiment-based "plausible range" for 18 parameters in 4 categories of embedding, loading, material, and segmentation. The predictive ability of the models were compared in terms of the coefficient of determination (R²) of the linear regression between the measured and predicted stiffness values in all load cases. Our model was capable of capturing 90% of the variation in the experimental stiffness of the samples in neutral stance position (R² = 0.9, concordance correlation coefficient (CCC) = 0.93, percent root mean squared error (RMSE%) = 8.4%, slope and intercept not significantly different from unity and zero, respectively). Embedding and loading categories strongly affected the predictive ability of the models with an average percent difference in R² of 4.36% ± 2.77 and 2.96% ± 1.69 for the stance-neutral load case, respectively. The performance of the models were significantly different in adducted and abducted load cases with their R² dropping to 71% and 70%, respectively. Similarly, off-axes load cases were affected by the parameters differently compared to the neutral load case, with the loading parameter category imposing more than 10% difference on their R², larger than all other categories. We also showed that automatically selecting the best performing plausible value for each parameter and each sample would result in a perfectly linear correlation (R²> 0.99) between the "tuned" model's predicted stiffness and experimental results. Based on our results, high sensitivity of the model performance to experimental parameters requires extra diligence in modeling the embedding geometry and the loading angles since these sources of uncertainty could dwarf the effects of material modeling and image processing parameters. The results of this study could help in improving the robustness of the QCT-FE models of proximal femur by limiting the uncertainties in the experimental and modeling steps.

A Density-Dependent Target Stimulus for Inverse Bone (Re)modeling with Homogenized Finite Element Models

Inverse bone (re)modeling (IBR) can infer physiological loading conditions from the bone microstructure. IBR scales unit loads, imposed on finite element (FE) models of a bone, such that the trabecular microstructure is homogeneously loaded and the difference to a target stimulus is minimized. Micro-FE (µFE) analyses are typically used to model the microstructure, but computationally more efficient, homogenized FE (hFE) models, where the microstructure is replaced by an equivalent continuum, could be used instead. However, also the target stimulus has to be translated from the tissue to the continuum level. In this study, a new continuum-level target stimulus relating relative bone density and strain energy density is proposed. It was applied using different types of hFE models to predict the physiological loading of 21 distal radii sections, which was subsequently compared to µFE-based IBR. The hFE models were able to correctly identify the dominant load direction and showed a high correlation of the predicted forces, but mean magnitude errors ranged from − 14.7 to 26.6% even for the best models. While µFE-based IBR can still be regarded as a gold standard, hFE-based IBR enables faster predictions, the usage of more sophisticated boundary conditions, and the usage of clinical images.

Photobiomodulation Therapy Affects the Elastic Modulus, Cytoskeletal Rearrangement and Migration Capability of Human Osteosarcoma Cells

Photobiomodulation (PBM) therapy utilizes low-power lasers to modulate the viability of living human cells and leads to changes in proliferation, differentiation, adhesion and gene expression, even though the rearrangement of cytoskeleton was not previously studied. The present study aims to evaluate the photobiological effects on the elastic behavior of human osteosarcoma cells (MG-63) and their morphological changes. Fluorescence staining, confocal imaging and atomic force microscopy (AFM) topography were performed to study the effects of PBM therapy with the exposure of 532 nm-25mW, 650 nm-3mW, 650 nm-150mW and 780 nm-70mW beams following the 5-min continuous irradiation. The area of each beam was 3.14cm² with a source-surface distance of 20 cm. Besides the cell proliferation assessment, the migratory potential of MG-63 was determined with the wound healing technique. The results indicated an increase in stiffness and shape index of radiation-induced cells 24 h after exposure along with the obvious F-actins changes. But, cell stiffening was not observed 72 h after 532 nm laser irradiation. Also, a decrease in the migration rate was seen in all of the groups after 72 h of irradiation except cells treated with 532 nm wavelength. However, 532 nm laser beams increase the migratory potential 24 h after exposure. Within 72 h after irradiation, the cell proliferation was only affected by applying 532 nm and 650 nm-150mW laser beams. It was concluded that applying photobiomodulation with wavelengths of 650 nm (at both utilized powers) and 780 nm alters the migration capability and provides a quantitative description of cytoskeletal changes. Moreover, membrane stiffening can be considered as the biological marker of PBM treatments.

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.

Structural Design Method for Constructions: Simulation, Manufacturing and Experiment

The development of additive manufacturing technology leads to new concepts for design implants and prostheses. The necessity of such approaches is fueled by patient-oriented medicine. Such a concept involves a new way of understanding material and includes complex structural geometry, lattice constructions, and metamaterials. This leads to new design concepts. In the article, the structural design method is presented. The general approach is based on the separation of the micro- and macro-mechanical parameters. For this purpose, the investigated region as a complex of the basic cells was considered. Each basic cell can be described by a parameters vector. An initializing vector was introduced to control the changes in the parameters vector. Changing the parameters vector according to the stress-strain state and the initializing vector leads to changes in the basic cells and consequently to changes in the microarchitecture. A medium with a spheroidal pore was considered as a basic cell. Porosity and ellipticity were used for the parameters vector. The initializing vector was initialized and depended on maximum von Mises stress. A sample was designed according to the proposed method. Then, solid and structurally designed samples were produced by additive manufacturing technology. The samples were scanned by computer tomography and then tested by structural loads. The results and analyses were presented.

Bone mineral density modeling via random field: Normality, stationarity, sex and age dependence

BACKGROUND AND OBJECTIVE

Capturing the population variability of bone properties is of paramount importance to biomedical engineering. The aim of the present paper is to describe variability and correlations in bone mineral density with a spatial random field inferred from routine computed tomography data.

METHODS

Random fields were simulated by transforming pairwise uncorrelated Gaussian random variables into correlated variables through the spectral decomposition of an age-detrended correlation matrix. The validity of the random field model was demonstrated in the spatiotemporal analysis of bone mineral density. The similarity between the computed tomography samples and those generated via random fields was analyzed with the energy distance metric.

RESULTS

The random field of bone mineral density was found to be approximately Gaussian/slightly left-skewed/strongly right-skewed at various locations. However, average bone density could be simulated well with the proposed Gaussian random field for which the energy distance, i.e., a measure that quantifies discrepancies between two distribution functions, is convergent with respect to the number of correlation eigenpairs.

CONCLUSIONS

The proposed random field model allows the enhancement of computational biomechanical models with variability in bone mineral density, which could increase the usability of the model and provides a step forward in in-silico medicine.

Can DXA image-based deep learning model predict the anisotropic elastic behavior of trabecular bone?

3D image-based finite element (FE) and bone volume fraction (BV/TV)/fabric tensor modeling techniques are currently used to determine the apparent stiffness tensor of trabecular bone for assessing its anisotropic elastic behavior. Inspired by the recent success of deep learning (DL) techniques, we hypothesized that DL modeling techniques could be used to predict the apparent stiffness tensor of trabecular bone directly using dual-energy X-ray absorptiometry (DXA) images. To test the hypothesis, a convolutional neural network (CNN) model was trained and validated to predict the apparent stiffness tensor of trabecular bone cubes using their DXA images. Trabecular bone cubes obtained from human cadaver proximal femurs were used to obtain simulated DXA images as input, and the apparent stiffness tensor of the trabecular cubes determined by using micro-CT based FE simulations was used as output (ground truth) to train the DL model. The prediction accuracy of the DL model was evaluated by comparing it with the micro-CT based FE models, histomorphometric parameter based multiple linear regression models, and BV/TV/fabric tensor based multiple linear regression models. The results showed that DXA image-based DL model achieved high fidelity in predicting the apparent stiffness tensor of trabecular bone cubes (R² = 0.905-0.973), comparable to or better than the histomorphometric parameter based multiple linear regression and BV/TV/fabric tensor based multiple linear regression models, thus supporting the hypothesis of this study. The outcome of this study could be used to help develop DXA image-based DL techniques for clinical assessment of bone fracture risk.

Loss of longitudinal superiority marks the microarchitecture deterioration of osteoporotic cancellous bones

Osteoporosis (OP), a skeletal disease making bone mechanically deteriorate and easily fracture, is a global public health issue due to its high prevalence. It has been well recognized that besides bone loss, microarchitecture degradation plays a crucial role in the mechanical deterioration of OP bones, but the specific role of microarchitecture in OP has not been well clarified and quantified from mechanics perspective. Here, we successfully decoupled and identified the specific roles of microarchitecture, bone mass and tissue property in the failure properties of cancellous bones, through μCT-based digital modeling and finite element method simulations on bone samples from healthy and ovariectomy-induced osteoporotic mice. The results show that the microarchitecture of healthy bones exhibits longitudinal superiority in mechanical properties such as the effective stiffness, strength and toughness, which fits them well to bearing loads along their longitudinal direction. OP does not only reduce bone mass but also impair the microarchitecture topology. The former is mainly responsible for the mechanical degradation of bones in magnitude, wherever the latter accounts for the breakdown of their function-favorable anisotropy, the longitudinal superiority. Hence, we identified the microarchitecture-deterioration-induced directional mismatch between material and loading as a hazardous feature of OP and defined a longitudinal superiority index as measurement of the health status of bone microarchitecture. These findings provide useful insights and guidelines for OP diagnosis and treat assessment.

Towards optimization of volar plate fixations of distal radius fractures: Using finite element analyses to reduce the number of screws

BACKGROUND

Using fewer distal screws in volar plate fixation of distal radius fractures could reduce treatment costs and complications. However, there is currently no consensus on the ideal screw configuration, likely due to experimental limitations and its subject-specific nature. In this study, finite element analysis was used to investigate (1) if reducing the number of screws is biomechanically feasible and (2) if an optimal screw configuration is subject-specific.

METHODS

Validated subject-specific finite element models of 16 human radii with extra articular distal radius fractures and volar plate fixation with six distal screws were used as a baseline. 41 additional configurations with three to six distal screws were simulated for each subject. Axial stiffness and peri-implant strains around the distal screws were evaluated. Subject-specific optimum configurations were determined using a lower bound for the axial stiffness and minimizing peri-implant strains.

FINDINGS

Even using three distal screws led to only minor deterioration of the biomechanical properties in the best configuration (axial stiffness: -11.2%, peri-implant strains: -35.0%), but a considerable deterioration in the worst configuration (axial stiffness: -46.2%, peri-implant strains: +112.4%). The optimization showed that the ideal screw configuration is subject-specific and on average 1.9 screws could be saved based on the herein used optimization criterion.

INTERPRETATION

This study highlights that not only how many, but which screws are used in volar plate fixation of distal radius fractures is critical. Using a patient-specific selection of distal screws bears potential to save costs and reduce complications.

Empirical relationships between bone density and ultimate strength: A literature review

INTRODUCTION

Ultimate strength-density relationships for bone have been reported with widely varying results. Reliable bone strength predictions are crucial for many applications that aim to assess bone failure. Bone density and bone morphology have been proposed to explain most of the variance in measured bone strength. If this holds true, it could lead to the derivation of a single ultimate strength-density-morphology relationship for all anatomical sites.

METHODS

All relevant literature was reviewed. Ultimate strength-density relationships derived from mechanical testing of human bone tissue were included. The reported relationships were translated to ultimate strength-apparent density relationships and normalized with respect to strain rate. Results were grouped based on bone tissue type (cancellous or cortical), anatomical site, and loading mode (tension vs. compression). When possible, the relationships were compared to existing ultimate strength-density-morphology relationships.

RESULTS

Relationships that considered bone density and morphology covered the full spectrum of eight-fold inter-study difference in reported compressive ultimate strength-density relationships for trabecular bone. This was true for studies that tested specimens in different loading direction and tissue from different anatomical sites. Sparse data was found for ultimate strength-density relationships in tension and for cortical bone properties transverse to the main loading axis of the bone.

CONCLUSIONS

Ultimate strength-density-morphology relationships could explain measured strength across anatomical sites and loading directions. We recommend testing of bone specimens in other directions than along the main trabecular alignment and to include bone morphology in studies that investigate bone material properties. The lack of tensile strength data did not allow for drawing conclusions on ultimate strength-density-morphology relationships. Further studies are needed. Ideally, these studies would investigate both tensile and compressive strength-density relationships, including morphology, to close this gap and lead to more accurate evaluation of bone failure.

Comparison of different microCT-based morphology assessment tools using human trabecular bone

MicroCT-based morphological parameters are often used to quantify the structural properties of trabecular bone. Various software tools are available for calculating these parameters. Studies that examine the comparability of their results are rare. Four different software tools were used to analyse a set of 701 microCT images from human trabecular bone samples. Bone volume to total volume (BV/TV), bone surface (BS), trabecular thickness (Tb. Th.) and degree of anisotropy (DA) were evaluated. BV/TV shows very low difference (-0.18 ± 0.15%). The difference in BS could be reduced below 5% if artificial cut surfaces are not included. Tb. Th. and Tb. Sp. show differences of maximal -12% although the same theoretical background is used. DA is most critical with differences from 4.75 ± 3.70% (medtool vs. Scanco), over -38.61 ± 13.15% (BoneJ vs. Scanco), up to 80.52 ± 50.04% (medtool vs. BoneJ). Quantitative results should be considered with caution, especially when comparing different studies. Introducing standardization procedures and the disclosure of underlying algorithms and their respective implementations could improve this issue.

An explicit micro-FE approach to investigate the post-yield behaviour of trabecular bone under large deformations

Homogenised finite element (FE) analyses are able to predict osteoporosis-related bone fractures and become useful for clinical applications. The predictions of FE analyses depend on the apparent, heterogeneous, anisotropic, elastic, and yield material properties, which are typically determined by implicit micro-FE (μFE) analyses of trabecular bone. The objective of this study is to explore an explicit μFE approach to determine the apparent post-yield behaviour of trabecular bone, beyond the elastic and yield properties. The material behaviour of bone tissue was described by elasto-plasticity with a von Mises yield criterion closed by a planar cap for positive hydrostatic stresses to distinguish the post-yield behaviour in tension and compression. Two ultimate strains for tension and compression were calibrated to trigger element deletion and reproduce damage of trabecular bone. A convergence analysis was undertaken to assess the role of the mesh. Thirteen load cases using periodicity-compatible mixed uniform boundary conditions were applied to three human trabecular bone samples of increasing volume fractions. The effect of densification in large strains was explored. The convergence study revealed a strong dependence of the apparent ultimate stresses and strains on element size. An apparent quadric strength surface for trabecular bone was successfully fitted in a normalised stress space. The effect of densification was reproduced and correlated well with former experimental results. This study demonstrates the potential of the explicit FE formulation and the element deletion technique to reproduce damage in trabecular bone using μFE analyses. The proper account of the mesh sensitivity remains challenging for practical computing times.

The Multi-Axial Failure Response of Porcine Trabecular Skull Bone Estimated Using Microstructural Simulations

The development of a multi-axial failure criterion for trabecular skull bone has many clinical and biological implications. This failure criterion would allow for modeling of bone under daily loading scenarios that typically are multi-axial in nature. Some yield criteria have been developed to evaluate the failure of trabecular bone, but there is a little consensus among them. To help gain deeper understanding of multi-axial failure response of trabecular skull bone, we developed 30 microstructural finite element models of porous porcine skull bone and subjected them to multi-axial displacement loading simulations that spanned three-dimensional (3D) stress and strain space. High-resolution microcomputed tomography (microCT) scans of porcine trabecular bone were obtained and used to develop the meshes used for finite element simulations. In total, 376 unique multi-axial loading cases were simulated for each of the 30 microstructure models. Then, results from the total of 11,280 simulations (approximately 135,360 central processing unit-hours) were used to develop a mathematical expression, which describes the average three-dimensional yield surface in strain space. Our results indicate that the yield strain of porcine trabecular bone under multi-axial loading is nearly isotropic and despite a spread of yielding points between the 30 different microstructures, no significant relationship between the yield strain and bone volume fraction is observed. The proposed yield equation has simple format and it can be implemented into a macroscopic model for the prediction of failure of whole bones.

Mapping anisotropy improves QCT-based finite element estimation of hip strength in pooled stance and side-fall load configurations

Hip fractures are one of the most severe consequences of osteoporosis. Compared to the clinical standard of DXA-based aBMD at the femoral neck, QCT-based FEA delivers a better surrogate of femoral strength and gains acceptance for the calculation of hip fracture risk when a CT reconstruction is available. Isotropic, homogenised voxel-based, finite element (hvFE) models are widely used to estimate femoral strength in cross-sectional and longitudinal clinical studies. However, fabric anisotropy is a classical feature of the architecture of the proximal femur and the second determinant of the homogenised mechanical properties of trabecular bone. Due to the limited resolution, fabric anisotropy cannot be derived from clinical CT reconstructions. Alternatively, fabric anisotropy can be extracted from HR-pQCT images of cadaveric femora. In this study, fabric anisotropy from HR-pQCT images was mapped onto QCT-based hvFE models of 71 human proximal femora for which both HR-pQCT and QCT images were available. Stiffness and ultimate load computed from anisotropic hvFE models were compared with previous biomechanical tests in both stance and side-fall configurations. The influence of using the femur-specific versus a mean fabric distribution on the hvFE predictions was assessed. Femur-specific and mean fabric enhance the prediction of experimental ultimate force for the pooled, i.e. stance and side-fall, (isotropic: r²=0.81, femur-specific fabric: r²=0.88, mean fabric: r²=0.86,p<0.001) but not for the individual configurations. Fabric anisotropy significantly improves bone strength prediction for the pooled configurations, and mapped fabric provides a comparable prediction to true fabric. The mapping of fabric anisotropy is therefore expected to help generate more accurate QCT-based hvFE models of the proximal femur for personalised or multiple load configurations.

Direct estimation of human trabecular bone stiffness using cone beam computed tomography

OBJECTIVES

The aim of this study was to evaluate the possibility of estimating the biomechanical properties of trabecular bone through finite element simulations by using dental cone beam computed tomography data.

STUDY DESIGN

Fourteen human radius specimens were scanned in 3 cone beam computed tomography devices: 3-D Accuitomo 80 (J. Morita MFG., Kyoto, Japan), NewTom 5 G (QR Verona, Verona, Italy), and Verity (Planmed, Helsinki, Finland). The imaging data were segmented by using 2 different methods. Stiffness (Young modulus), shear moduli, and the size and shape of the stiffness tensor were studied. Corresponding evaluations by using micro-CT were regarded as the reference standard.

RESULTS

The 3-D Accuitomo 80 (J. Morita MFG., Kyoto, Japan) showed good performance in estimating stiffness and shear moduli but was sensitive to the choice of segmentation method. NewTom 5 G (QR Verona, Verona, Italy) and Verity (Planmed, Helsinki, Finland) yielded good correlations, but they were not as strong as Accuitomo 80 (J. Morita MFG., Kyoto, Japan). The cone beam computed tomography devices overestimated both stiffness and shear compared with the micro-CT estimations.

CONCLUSIONS

Finite element-based calculations of biomechanics from cone beam computed tomography data are feasible, with strong correlations for the Accuitomo 80 scanner (J. Morita MFG., Kyoto, Japan) combined with an appropriate segmentation method. Such measurements might be useful for predicting implant survival by in vivo estimations of bone properties.

Influence of the shape of the micro-finite element model on the mechanical properties calculated from micro-finite element analysis

Assessing the biomechanical properties of trabecular bone is of major biological and clinical significance for the research of bone diseases, fractures and their treatments. Micro-finite element (µFE) models are becoming increasingly popular for investigating the biomechanical properties of trabecular bone. The shapes of µFE models typically include cube and cylinder. Whether there are differences between cubic and cylindrical µFE models has not yet been studied. In the present study, cubic and cylindrical µFE models of human vertebral trabecular bone were constructed. A 1% strain was prescribed to the model along the superior-inferior direction. E values were calculated from these models, and paired t-tests were performed to determine whether these were any differences between E values obtained from cubic and cylindrical models. The results demonstrated that there were no statistically significant differences in the E values between cubic and cylindrical models, and there were no significant differences in Von Mises stress distributions between the two models. These findings indicated that, to construct µFE models of vertebral trabecular bone, cubic or cylindrical models were both feasible. Choosing between the cubic or cylindrical µFE model is dependent upon the specific study design.

Accounting for spatial variation of trabecular anisotropy with subject-specific finite element modeling moderately improves predictions of local subchondral bone stiffness at the proximal tibia

INTRODUCTION

Previously, a finite element (FE) model of the proximal tibia was developed and validated against experimentally measured local subchondral stiffness. This model indicated modest predictions of stiffness (R²=0.77, normalized root mean squared error (RMSE%)=16.6%). Trabecular bone though was modeled with isotropic material properties despite its orthotropic anisotropy. The objective of this study was to identify the anisotropic FE modeling approach which best predicted (with largest explained variance and least amount of error) local subchondral bone stiffness at the proximal tibia.

METHODS

Local stiffness was measured at the subchondral surface of 13 medial/lateral tibial compartments using in situ macro indentation testing. An FE model of each specimen was generated assuming uniform anisotropy with 14 different combinations of cortical- and tibial-specific density-modulus relationships taken from the literature. Two FE models of each specimen were also generated which accounted for the spatial variation of trabecular bone anisotropy directly from clinical CT images using grey-level structure tensor and Cowin's fabric-elasticity equations. Stiffness was calculated using FE and compared to measured stiffness in terms of R² and RMSE%.

RESULTS

The uniform anisotropic FE model explained 53-74% of the measured stiffness variance, with RMSE% ranging from 12.4 to 245.3%. The models which accounted for spatial variation of trabecular bone anisotropy predicted 76-79% of the variance in stiffness with RMSE% being 11.2-11.5%.

CONCLUSIONS

Of the 16 evaluated finite element models in this study, the combination of Synder and Schneider (for cortical bone) and Cowin's fabric-elasticity equations (for trabecular bone) best predicted local subchondral bone stiffness.

Mimetization of the elastic properties of cancellous bone via a parameterized cellular material

Bone tissue mechanical properties and trabecular microarchitecture are the main factors that determine the biomechanical properties of cancellous bone. Artificial cancellous microstructures, typically described by a reduced number of geometrical parameters, can be designed to obtain a mechanical behavior mimicking that of natural bone. In this work, we assess the ability of the parameterized microstructure introduced by Kowalczyk (Comput Methods Biomech Biomed Eng 9:135-147, 2006. doi: 10.1080/10255840600751473 ) to mimic the elastic response of cancellous bone. Artificial microstructures are compared with actual bone samples in terms of elasticity matrices and their symmetry classes. The capability of the parameterized microstructure to combine the dominant isotropic, hexagonal, tetragonal and orthorhombic symmetry classes in the proportions present in the cancellous bone is shown. Based on this finding, two optimization approaches are devised to find the geometrical parameters of the artificial microstructure that better mimics the elastic response of a target natural bone specimen: a Sequential Quadratic Programming algorithm that minimizes the norm of the difference between the elasticity matrices, and a Pattern Search algorithm that minimizes the difference between the symmetry class decompositions. The pattern search approach is found to produce the best results. The performance of the method is demonstrated via analyses for 146 bone samples.

μCT-based trabecular anisotropy can be reproducibly computed from HR-pQCT scans using the triangulated bone surface

The trabecular structure can be assessed at the wrist or tibia via high-resolution peripheral quantitative computed tomography (HR-pQCT). Yet on this modality, the performance of the existing methods, evaluating trabecular anisotropy is usually overlooked, especially in terms of reproducibility. We thus proposed to compare the TRI routine used by SCANCO Medical AG (Brüttisellen, Switzerland), the classical mean intercept length (MIL), and the grey-level structure tensor (GST) to the mean surface length (MSL), a new method for evaluating a second-order fabric tensor based on the triangulation of the bone surface. The distal radius of 24 fresh-frozen human forearms was scanned three times via HR-pQCT protocols (61μm, 82μm nominal voxel size), dissected, and imaged via micro computed tomography (μCT) at 16μm nominal voxel size. After registering the scans, we compared for each resolution the fabric tensors, determined by the mentioned techniques for 182 trabecular regions of interest. We then evaluated the reproducibility of the fabric information measured by HR-pQCT via precision errors. On μCT, TRI and GST were respectively the best and worst surrogates for MIL_μCT (MIL computed on μCT) in terms of eigenvalues and main direction of anisotropy. On HR-pQCT, however, MSL provided the best approximation of MIL_μCT. Surprisingly, surface-based approaches (TRI, MSL) also proved to be more precise than both MIL and GST. Our findings confirm that MSL can reproducibly estimate MIL_μCT, the current gold standard. MSL thus enables the direct mapping of the fabric-dependent material properties required in homogenised HR-pQCT-based finite element models.

Finite Element-Based Mechanical Assessment of Bone Quality on the Basis of In Vivo Images

Beyond bone mineral density (BMD), bone quality designates the mechanical integrity of bone tissue. In vivo images based on X-ray attenuation, such as CT reconstructions, provide size, shape, and local BMD distribution and may be exploited as input for finite element analysis (FEA) to assess bone fragility. Further key input parameters of FEA are the material properties of bone tissue. This review discusses the main determinants of bone mechanical properties and emphasizes the added value, as well as the important assumptions underlying finite element analysis. Bone tissue is a sophisticated, multiscale composite material that undergoes remodeling but exhibits a rather narrow band of tissue mineralization. Mechanically, bone tissue behaves elastically under physiologic loads and yields by cracking beyond critical strain levels. Through adequate cell-orchestrated modeling, trabecular bone tunes its mechanical properties by volume fraction and fabric. With proper calibration, these mechanical properties may be incorporated in quantitative CT-based finite element analysis that has been validated extensively with ex vivo experiments and has been applied increasingly in clinical trials to assess treatment efficacy against osteoporosis.

Relationship between sample volumes and modulus of human vertebral trabecular bone in micro-finite element analysis

Micro-finite element (μFE) models have been widely used to assess the biomechanical properties of trabecular bone. How to choose a proper sample volume of trabecular bone, which could predict the real bone biomechanical properties and reduce the calculation time, was an interesting problem. Therefore, the purpose of this study was to investigate the relationship between different sample volumes and apparent elastic modulus (E) calculated from μFE model. 5 Human lumbar vertebral bodies (L1-L5) were scanned by micro-CT. Cubic concentric samples of different lengths were constructed as the experimental groups and the largest possible volumes of interest (VOI) were constructed as the control group. A direct voxel-to-element approach was used to generate μFE models and steel layers were added to the superior and inferior surface to mimic axial compression tests. A 1% axial strain was prescribed to the top surface of the model to obtain the E values. ANOVA tests were performed to compare the E values from the different VOIs against that of the control group. Nonlinear function curve fitting was performed to study the relationship between volumes and E values. The larger cubic VOI included more nodes and elements, and more CPU times were needed for calculations. E values showed a descending tendency as the length of cubic VOI decreased. When the volume of VOI was smaller than (7.34mm(3)), E values were significantly different from the control group. The fit function showed that E values approached an asymptotic values with increasing length of VOI. Our study demonstrated that apparent elastic modulus calculated from μFE models were affected by the sample volumes. There was a descending tendency of E values as the length of cubic VOI decreased. Sample volume which was not smaller than (7.34mm(3)) was efficient enough and timesaving for the calculation of E.

A patient-specific model of total knee arthroplasty to estimate patellar strain: A case study

BACKGROUND

Inappropriate patellar cut during total knee arthroplasty can lead to patellar complications due to increased bone strain. In this study, we evaluated patellar bone strain of a patient who had a deeper patellar cut than the recommended.

METHODS

A patient-specific model based on patient preoperative data was created. The model was decoupled into two levels: knee and patella. The knee model predicted kinematics and forces on the patella during squat movement. The patella model used these values to predict bone strain after total knee arthroplasty. Mechanical properties of the patellar bone were identified with micro-finite element modeling testing of cadaveric samples. The model was validated with a robotic knee simulator and postoperative X-rays. For this patient, we compared the deeper patellar cut depth to the recommended one, and evaluated patellar bone volume with octahedral shear strain above 1%.

FINDINGS

Model predictions were consistent with experimental measurements of the robotic knee simulator and postoperative X-rays. Compared to the recommended cut, the deeper cut increased the critical strain bone volume, but by less than 3% of total patellar volume.

INTERPRETATION

We thus conclude that the predicted increase in patellar strain should be within an acceptable range, since this patient had no complaints 8 months after surgery. This validated patient-specific model will later be used to address other questions on groups of patients, to eventually improve surgical planning and outcome of total knee arthroplasty.

Comparison of mixed and kinematic uniform boundary conditions in homogenized elasticity of femoral trabecular bone using microfinite element analyses

Mechanical properties of human trabecular bone play an important role in age-related bone fragility and implant stability. Microfinite element (lFE) analysis allows computing the apparent elastic properties of trabecular bone for use in homogenized FE (hFE) analysis,but the results depend unfortunately on the type of applied boundary conditions(BCs). In this study, 167 human femoral trabecular cubic regions with a side length of 5.3mm were extracted from three proximal femora and analyzed using lFE analysis to compare systematically their stiffness with kinematic uniform BCs (KUBCs) and periodicity-compatible mixed uniform BCs (PMUBCs). The obtained elastic constants were then used in the volume fraction and fabric-based orthotropic Zysset–Curnier model to identify their respective model parameters. As expected, PMUBCs lead to more compliant apparent elastic properties than KUBCs, especially in shear. The differences in stiffness decreased with bone volume fraction and mean intercept length (MIL). Unlike KUBCs, PMUBCs were sensitive to heterogeneity of the biopsies. The Zysset–Curnier model fitted the apparent elastic constants successfully in both cases with adjusted coefficients of determination (r2adj) of 0.986 for KUBCs and 0.975 for PMUBCs. The proper use of these BCs for hFE analysis of whole bones will need to be investigated in future work.

Elastic Anisotropy of Trabecular Bone in the Elderly Human Vertebra

Knowledge of the nature of the elastic symmetry of trabecular bone is fundamental to the study of bone adaptation and failure. Previous studies have classified human vertebral trabecular bone as orthotropic or transversely isotropic but have typically obtained samples from only selected regions of the centrum. In this study, the elastic symmetry of human vertebral trabecular bone was characterized using microfinite element (μFE) analyses performed on 1019 cubic regions of side length equal to 5 mm, obtained via thorough sampling of the centrums of 18 human L1 vertebrae (age = 81.17 ± 7.7 yr; eight males and ten females). An optimization procedure was used to find the closest orthotropic representation of the resulting stiffness tensor for each cube. The orthotropic elastic constants and orientation of the principal elastic axes were then recorded for each cube and were compared to the constants predicted from Cowin's fabric-based constitutive model (Cowin, 1985, "The Relationship Between the Elasticity Tensor and the Fabric Tensor," Mech. Mater., 4(2), pp. 137-147.) and the orientation of the principal axes of the fabric tensor, respectively. Deviations from orthotropy were quantified by the "orthotropic error" (van Rietbergen et al., 1996, "Direct Mechanics Assessment of Elastic Symmetries and Properties of Trabecular Bone Architecture," J. Biomech., 29(12), pp. 1653-1657), and deviations from transverse isotropy were determined by statistical comparison of the secondary and tertiary elastic moduli. The orthotropic error was greater than 50% for nearly half of the cubes, and the secondary and tertiary moduli differed from one another (p < 0.0001). Both the orthotropic error and the difference between secondary and tertiary moduli decreased with increasing bone volume fraction (BV/TV; p ≤ 0.007). Considering only the cubes with an orthotropic error less than 50%, only moderate correlations were observed between the fabric-based and the μFE-computed elastic moduli (R2 ≥ 0.337; p < 0.0001). These results indicate that when using a criterion of 5 mm for a representative volume element (RVE), transverse isotropy or orthotropy cannot be assumed for elderly human vertebral trabecular bone. Particularly at low values of BV/TV, this criterion does not ensure applicability of theories of continuous media. In light of the very sparse and inhomogeneous microstructure found in the specimens analyzed in this study, further work is needed to establish guidelines for selecting a RVE within the aged vertebral centrum.

The influence of bone density and anisotropy in finite element models of distal radius fracture osteosynthesis: Evaluations and comparison to experiments

Continuum-level finite element (FE) models can be used to analyze and improve osteosynthesis procedures for distal radius fractures (DRF) from a biomechanical point of view. However, previous models oversimplified the bone material and lacked thorough experimental validation. The goal of this study was to assess the influence of local bone density and anisotropy in FE models of DRF osteosynthesis for predictions of axial stiffness, implant plate stresses, and screw loads. Experiments and FE analysis were conducted in 25 fresh frozen cadaveric radii with DRFs treated by volar locking plate osteosynthesis. Specimen specific geometries were captured using clinical quantitative CT (QCT) scans of the prepared samples. Local bone material properties were computed based on high resolution CT (HR-pQCT) scans of the intact radii. The axial stiffness and individual screw loads were evaluated in FE models, with (1) orthotropic inhomogeneous (OrthoInhom), (2) isotropic inhomogeneous (IsoInhom), and (3) isotropic homogeneous (IsoHom) bone material and compared to the experimental axial stiffness and screw-plate interface failures. FE simulated and experimental axial stiffness correlated significantly (p<0.0001) for all three model types. The coefficient of determination was similar for OrthoInhom (R(2)=0.807) and IsoInhom (R(2)=0.816) models but considerably lower for IsoHom models (R(2)=0.500). The peak screw loads were in qualitative agreement with experimental screw-plate interface failure. Individual loads and implant plate stresses of IsoHom models differed significantly (p<0.05) from OrthoInhom and IsoInhom models. In conclusion, including local bone density in FE models of DRF osteosynthesis is essential whereas local bone anisotropy hardly effects the models׳ predictive abilities.

Clinical Use of Quantitative Computed Tomography-Based Finite Element Analysis of the Hip and Spine in the Management of Osteoporosis in Adults: the 2015 ISCD Official Positions-Part II

The International Society for Clinical Densitometry (ISCD) has developed new official positions for the clinical use of quantitative computed tomography (QCT)-based finite element analysis of the spine and hip. The ISCD task force for QCT reviewed the evidence for clinical applications and presented a report with recommendations at the 2015 ISCD Position Development Conference. Here we discuss the agreed upon ISCD official positions with supporting medical evidence, rationale, controversy, and suggestions for further study. Parts I and III address the clinical use of QCT of the hip, and the clinical feasibility of existing techniques for opportunistic screening of osteoporosis using CT scans obtained for other diagnosis such as colonography was addressed.

Bone volume fraction and fabric anisotropy are better determinants of trabecular bone stiffness than other morphological variables

As our population ages, more individuals suffer from osteoporosis. This disease leads to impaired trabecular architecture and increased fracture risk. It is essential to understand how morphological and mechanical properties of the cancellous bone are related. Morphology-elasticity relationships based on bone volume fraction (BV/TV) and fabric anisotropy explain up to 98% of the variation in elastic properties. Yet, other morphological variables such as individual trabeculae segmentation (ITS) and trabecular bone score (TBS) could improve the stiffness predictions. A total of 743 micro-computed tomography (μCT) reconstructions of cubic trabecular bone samples extracted from femur, radius, vertebrae, and iliac crest were analyzed. Their morphology was assessed via 25 variables and their stiffness tensor (CFE) was computed from six independent load cases using micro finite element (μFE) analyses. Variance inflation factors were calculated to evaluate collinearity between morphological variables and decide upon their inclusion in morphology-elasticity relationships. The statistically admissible morphological variables were included in a multiple linear regression model of the dependent variable CFE. The contribution of each independent variable was evaluated (ANOVA). Our results show that BV/TV is the best determinant of CFE(r(2) adj  = 0.889), especially in combination with fabric anisotropy (r(2) adj  = 0.968). Including the other independent predictors hardly affected the amount of variance explained by the model (r(2) adj  = 0.975). Across all anatomical sites, BV/TV explained 87% of the variance of the bone elastic properties. Fabric anisotropy further described 10% of the bone stiffness, but the improvement in variance explanation by adding other independent factors was marginal (<1%). These findings confirm that BV/TV and fabric anisotropy are the best determinants of trabecular bone stiffness and show, against common belief, that other morphological variables do not bring any further contribution. These overall conclusions remain to be confirmed for specific bone diseases and postelastic properties.

A novel approach to estimate trabecular bone anisotropy using a database approach

Continuum finite element (FE) models of bones have become a standard pre-clinical tool to estimate bone strength. These models are usually based on clinical CT scans and material properties assigned are chosen as isotropic based only on the density distribution. It has been shown, however, that trabecular bone elastic behavior is best described as orthotropic. Unfortunately, the use of orthotropic models in FE analysis derived from CT scans is hampered by the fact that the measurement of a trabecular orientation (fabric) is not possible from clinical CT images due to the low resolution of such images. In this study, we explore the concept of using a database (DB) of high-resolution bone models to derive the fabric information that is missing in clinical images. The goal of this study was to investigate if models with fabric derived from a relatively small database can already produce more accurate results than isotropic models. A DB of 33 human proximal femurs was generated from micro-CT scans with a nominal isotropic resolution of 82 µm. Continuum FE models were generated from the images using a pre-defined mesh template in combination with an iso-anatomic mesh morphing tool. Each element within the mesh template is at a specific anatomical location. For each element within the cancellous bone, a spherical region around the element centroid with a radius of 2mm was defined. Bone volume fraction and the mean-intercept-length fabric tensor were analyzed for that region. Ten femurs were used as test cases. For each test femur, four different models were generated: (1) an orthotropic model based on micro-CT fabric measurements (gold standard), (2) an orthotropic model based on the fabric derived from the best-matched database model, (3) an isotropic-I model in which the fabric tensor was set to the identity tensor, and (4) a second isotropic-II model with its total bone stiffness fitted to the gold standard. An elastic-plastic damage model was used to simulate failure and post failure behavior during a fall to the side. The results show that all models produce a similar stress distribution. However, compared to the gold standard, both isotropic-I and II models underestimated the stress/damage distributions significantly. We found no significant difference between DB-derived and gold standard models. Compared to the gold standard, the isotropic-I models further underestimated whole bone stiffness by 26.3% and ultimate load by 14.5%, while these differences for the DB-derived orthotropic models were only 4.9% and 3.1% respectively. The results indicate that the concept of using a DB to estimate patient-specific anisotropic material properties can considerably improve the results. We expect that this approach can lead to more accurate results in particular for cases where bone anisotropy plays an important role, such as in osteoporotic patients and around implants.

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.