Prediction of femoral strength using 3D finite element models reconstructed from DXA images: validation against experiments

Computed tomography (CT)-based finite element (FE) models may improve the current osteoporosis diagnostics and prediction of fracture risk by providing an estimate for femoral strength. However, the need for a CT scan, as opposed to the conventional use of dual-energy X-ray absorptiometry (DXA) for osteoporosis diagnostics, is considered a major obstacle. The 3D shape and bone mineral density (BMD) distribution of a femur can be reconstructed using a statistical shape and appearance model (SSAM) and the DXA image of the femur. Then, the reconstructed shape and BMD could be used to build FE models to predict bone strength. Since high accuracy is needed in all steps of the analysis, this study aimed at evaluating the ability of a 3D FE model built from one 2D DXA image to predict the strains and fracture load of human femora. Three cadaver femora were retrieved, for which experimental measurements from ex vivo mechanical tests were available. FE models were built using the SSAM-based reconstructions: using only the SSAM-reconstructed shape, only the SSAM-reconstructed BMD distribution, and the full SSAM-based reconstruction (including both shape and BMD distribution). When compared with experimental data, the SSAM-based models predicted accurately principal strains (coefficient of determination >0.83, normalized root-mean-square error <16%) and femoral strength (standard error of the estimate 1215 N). These results were only slightly inferior to those obtained with CT-based FE models, but with the considerable advantage of the models being built from DXA images. In summary, the results support the feasibility of SSAM-based models as a practical tool to introduce FE-based bone strength estimation in the current fracture risk diagnostics.     

Comparison between mechanical stress and bone mineral density in the femur after total hip arthroplasty by using subject-specific finite element analyses

The mechanism underling bone mineral density (BMD) loss that occurs in the femur after total hip arthroplasty (THA) remains unknown. We compared the equivalent stress and strain energy density (SED) to BMD in the femur after THA using subject-specific finite element analyses. Twenty-four patients who had undergone primary cementless THA were analysed. BMD was measured using dual-energy X-ray absorptiometry (DEXA) at 1 week and 3, 6 and 12 months after THA. Seven regions of interest (ROIs) were defined in accordance with Gruen's system (ROIs 1–7). Computed tomography images of the femurs were acquired pre- and postoperatively, and the images were converted into three-dimensional finite element (FE) models. Equivalent stress and SED were analysed and compared with DEXA data. BMD was maintained 1 year after THA in ROIs 3, 4, 5 and 6, whereas BMD decreased in ROIs 1, 2 and 7. FE analysis revealed that equivalent stress in ROIs 3, 4, 5 and 6 was much higher than that in ROIs 1, 2 and 7. A significant correlation was observed between the rate of changes in BMD and equivalent stress. Reduction of equivalent stress may contribute to decrease in BMD in the femur after THA.     

The use of preserved tissue in finite element modelling of fresh ovine vertebral behaviour

The aim of this study was to investigate whether the predicted finite element (FE) stiffness of vertebral bone is altered when using images of preserved rather than fresh tissue to generate specimen-specific FE models. Fresh ovine vertebrae were used to represent embalmed (n = 3) and macerated dry-bone (n = 3) specimens and treated accordingly. Specimens were scanned pre- and post-treatment using micro-computed tomography. A constant threshold level derived from these images was used to calculate the respective bone volume fraction (BV/TV) from which the conversion factor validated for fresh tissue was used to determine material properties that were assigned to corresponding FE models. Results showed a definite change in the BV/TV between the fresh and the preserved bone. However, the changes in the predicted FE stiffness were not generally greater than the variations expected from assignment of loading and boundary conditions. In conclusion, images of preserved tissue can be used to generate FE models that are representative of fresh tissue with a tolerable level of error.     

Validation of a voxel-based FE method for prediction of the uniaxial apparent modulus of human trabecular bone using macroscopic mechanical tests and nanoindentation

Assessment of the mechanical properties of trabecular bone is of major biological and clinical importance for the investigation of bone diseases, fractures and their treatments. Finite element (FE) methods are getting increasingly popular for quantifying the elastic and failure properties of trabecular bone. In particular, voxel-based FE methods have been previously used to calculate the effective elastic properties of trabecular microstructures. However, in most studies, bone tissue moduli were assumed or back-calculated to match the apparent elastic moduli from experiments, which often lead to surprisingly low values when compared to nanoindentation results. In this study, voxel-based FE analysis of trabecular bone is combined with physical measures of volume fraction, micro-CT (microCT) reconstructions, uniaxial mechanical tests and specimen-specific nanoindentation tests for proper validation of the method. Cylindrical specimens of cancellous bone were extracted from human femurs and their volume fraction determined with Archimede's method. Uniaxial apparent modulus of the specimens was measured with an improved tension-compression testing protocol that minimizes boundary artefacts. Their microCT reconstructions were segmented to match the measured bone volume fraction and used to create full-size voxel models with 30-45 microm element size. For each specimen, linear isotropic elastic material properties were defined based on specific nanoindentation measurements of its embedded bone tissue. Linear FE analyses were finally performed to simulate the uniaxial mechanical tests. Additional parametric analyses were performed to evaluate the potential errors on the predicted apparent modulus arising from variations in segmentation threshold, tissue modulus, and the use of 125-mm(3) cubic sub-regions. The results demonstrate an excellent correspondence between experimental measures and FE predictions of uniaxial apparent modulus. In conclusion, the adopted voxel-based FE approach is found to be a robust method to predict the linear elastic properties of human cancellous bone, provided segmentation of the microCT reconstructions is carefully calibrated, tissue modulus is known a priori and the entire region of interest is included in the analysis.     

An accurate estimation of bone density improves the accuracy of subject-specific finite element models

An experimental-numerical study was performed to investigate the relationships between computed tomography (CT)-density and ash density, and between ash density and apparent density for bone tissue, to evaluate their influence on the accuracy of subject-specific FE models of human bones. Sixty cylindrical bone specimens were examined. CT-densities were computed from CT images while apparent and ash densities were measured experimentally. The CT/ash-density and ash/apparent-density relationships were calculated. Finite element models of eight human femurs were generated considering these relationships to assess their effect on strain prediction accuracy. CT and ash density were linearly correlated (R(2)=0.997) over the whole density range but not equivalent (intercep t <0, slope >1). A constant ash/apparent-density ratio (0.598+/-0.004) was found for cortical bone. A lower ratio, with a larger dispersion, was found for trabecular bone (0.459+/-0.100), but it became less dispersed, and equal to that of cortical tissue, when testing smaller trabecular specimens (0.598+/-0.036). This suggests that an experimental error occurred in apparent-density measurements for large trabecular specimens and a constant ratio can be assumed valid for the whole density range. Introducing the obtained relationships in the FE modelling procedure improved strain prediction accuracy (R(2)=0.95, RMSE=7%). The results suggest that: (i) a correction of the densitometric calibration should be used when evaluating bone ash-density from clinical CT scans, to avoid ash-density underestimation and overestimation for low- and high-density bone tissue, respectively; (ii) the ash/apparent-density ratio can be assumed constant in human femurs and (iii) the correction improves significantly the model accuracy and should be considered in subject-specific bone modelling.     

Postfailure modulus strongly affects microcracking and mechanical property change in human iliac cancellous bone: a study using a 2D nonlinear finite element method

A two-dimensional (2D) finite element (FE) method was used to estimate the ability of bone tissue to sustain damage as a function of postfailure modulus. Briefly, 2D nonlinear compact-tension FE models were created from quantitative back-scattered electron images taken of human iliac crest bone specimens. The effects of different postfailure moduli on predicted microcrack propagation were examined. The 2D FE models were used as surrogates for real bone tissues. The crack number was larger in models with higher postfailure modulus, while mean crack length and area were smaller in these models. The rate of stiffness reduction was greater in the models with lower postfailure modulus. Hence, the current results supported the hypothesis that hard tissue postfailure properties have strong effects on bone microdamage morphology and the rate of change in apparent mechanical 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.     

Femoral strength is better predicted by finite element models than QCT and DXA.

Clinicians and patients would benefit if accurate methods of predicting and monitoring bone strength in-vivo were available. A group of 51 human femurs (age range 21-93; 23 females, 28 males) were evaluated for bone density and geometry using quantitative computed tomography (QCT) and dual energy X-ray absorptiometry (DXA). Regional bone density and dimensions obtained from QCT and DXA were used to develop statistical models to predict femoral strength ex vivo. The QCT data also formed the basis of a three-dimensional finite element (FE) models to predict structural stiffness. The femurs were separated into two groups; a model training set (n = 25) was used to develop statistical models to predict ultimate load, and a test set (n = 26) was used to validate these models. The main goal of this study was to test the ability of DXA, QCT and FE techniques to predict fracture load non-invasively, in a simple load configuration which produces predominantly femoral neck fractures. The load configuration simulated the single stance phase portion of normal gait; in 87% of the specimens, clinical appearing sub-capital fractures were produced. The training/test study design provided a tool to validate that the predictive models were reliable when used on specimens with "unknown" strength characteristics. The FE method explained at least 20% more of the variance in strength than the DXA models. Planned refinements of the FE technique are expected to further improve these results. Three-dimensional FE models are a promising method for predicting fracture load, and may be useful in monitoring strength changes in vivo.     

Determination of orthotropic bone elastic constants using FEA and modal analysis

Finite element models have been widely employed in an effort to quantify the stress and strain distribution around implanted prostheses and to explore the influence of these distributions on their long-term stability. In order to provide meaningful predictions, such models must contain an appropriate reflection of mechanical properties. Detailed geometrical and density information is now readily available from CT scanning. However, despite the use of phantoms, a method of determining mechanical properties (or elastic constants) from bone density has yet to be made available in a usable form.In this study, a cadaveric bone was CT scanned and its natural frequencies were measured using modal analysis. Using the geometry obtained from the CT scan data, a finite element mesh was created with the distribution of density established by matching the mass of the FE bone model with the mass of the cadaveric bone. The maximum values of the orthotropic elastic constants were then established by matching the predictions from FE modal analyses to the experimental natural frequencies, giving a maximum error of 7.8% over 4 modes of vibration. Finally, the elastic constants of the bone derived from the analyses were compared with those measured using ultrasound techniques. This produced a difference of <1% for both the maximum density and axial Young's Modulus. This study has thereby produced an orthotropic finite element model of a human femur. More importantly, however, is the implication that it is possible to create a valid FE model by simply comparing the FE results with the measured resonant frequency of the CT scanned bone.     

A new approach for assigning bone material properties from CT images into finite element models

Generation of subject-specific finite element (FE) models from computed tomography (CT) datasets is of significance for application of the FE analysis to bone structures. A great challenge that remains is the automatic assignment of bone material properties from CT Hounsfield Units into finite element models. This paper proposes a new assignment approach, in which material properties are directly assigned to each integration point. Instead of modifying the dataset of FE models, the proposed approach divides the assignment procedure into two steps: generating the data file of the image intensity of a bone in a MATLAB program and reading the file into ABAQUS via user subroutines. Its accuracy has been validated by assigning the density of a bone phantom into a FE model. The proposed approach has been applied to the FE model of a sheep tibia and its applicability tested on a variety of element types. The proposed assignment approach is simple and illustrative. It can be easily modified to fit users’ situations.     

Development and validation of subject-specific finite element models for blunt trauma study

This study developed and validated finite element (FE) models of swine and human thoraxes and abdomens that had subject-specific anatomies and could accurately and efficiently predict body responses to blunt impacts. Anatomies of the rib cage, torso walls, thoracic, and abdominal organs were reconstructed from X-ray computed tomography (CT) images and extracted into geometries to build FE meshes. The rib cage was modeled as an inhomogeneous beam structure with geometry and bone material parameters determined directly from CT images. Meshes of soft components were generated by mapping structured mesh templates representative of organ topologies onto the geometries. The swine models were developed from and validated by 30 animal tests in which blunt insults were applied to swine subjects and CT images, chest wall motions, lung pressures, and pathological data were acquired. A comparison of the FE calculations of animal responses and experimental measurements showed a good agreement. The errors in calculated response time traces were within 10% for most tests. Calculated peak responses showed strong correlations with the experimental values. The stress concentration inside the ribs, lungs, and livers produced by FE simulations also compared favorably to the injury locations. A human FE model was developed from CT images from the Visible Human project and was scaled to simulate historical frontal and side post mortem human subject (PMHS) impact tests. The calculated chest deformation also showed a good agreement with the measurements. The models developed in this study can be of great value for studying blunt thoracic and abdominal trauma and for designing injury prevention techniques, equipments, and devices.     

Development of a parametric finite element model of the proximal femur using statistical shape and density modelling

Skeletal fractures associated with bone mass loss are a major clinical problem and economic burden, and lead to significant morbidity and mortality in the ageing population. Clinical image-based measures of bone mass show only moderate correlative strength with bone strength. However, engineering models derived from clinical image data predict bone strength with significantly greater accuracy. Currently, image-based finite element (FE) models are time consuming to construct and are non-parametric. The goal of this study was to develop a parametric proximal femur FE model based on a statistical shape and density model (SSDM) derived from clinical image data. A small number of independent SSDM parameters described the shape and bone density distribution of a set of cadaver femurs and captured the variability affecting proximal femur FE strength predictions. Finally, a three-dimensional FE model of an 'unknown' femur was reconstructed from the SSDM with an average spatial error of 0.016 mm and an average bone density error of 0.037 g/cm(3).     

Comparison of micro-level and continuum-level voxel models of the proximal femur

Continuum-level finite element (FE) models became standard computational tools for the evaluation of bone mechanical behavior from in vivo computed tomography scans. Such scans do not account for the anisotropy of the bone. Instead, local mechanical properties in the continuum-level FE models are assumed isotropic and are derived from bone density, using statistical relationships. Micro-FE models, on the other hand, incorporate the anisotropic structure in detail. This study aimed to quantify the effects of assumed isotropy, by comparing continuum-level voxel models of a healthy and a severely osteoporotic proximal femur with recently analyzed micro-FE models of the same bones. The micro-model element size was coarsened to generate continuum FE models with two different element sizes (0.64 and 3.04 mm) and two different density–modulus relationships found in the literature for wet and ash density. All FE models were subjected to the same boundary conditions that simulated a fall to the side, and the stress and strain distributions, model stiffness and yield load were compared. The results indicated that the stress and strain distributions could be reproduced well with the continuum models. The smallest differences between the continuum-level model and micro-level model predictions of the stiffness and yield load were obtained with the coarsest element size. Better results were obtained for both continuum-element sizes when isotropic moduli were based on ash density rather than wet density.     

Development and validation of a subject-specific finite element model of a human clavicle

This study aimed to develop and validate a finite element (FE) model of a human clavicle which can predict the structural response and bone fractures under both axial compression and anterior–posterior three-point bending loads. Quasi-static non-injurious axial compression and three-point bending tests were first conducted on a male clavicle followed by a dynamic three-point bending test to fracture. Then, two types of FE models of the clavicle were developed using bone material properties which were set to vary with the computed tomography image density of the bone. A volumetric solid FE model comprised solely of hexahedral elements was first developed. A solid-shell FE model was then created which modelled the trabecular bone as hexahedral elements and the cortical bone as quadrilateral shell elements. Finally, simulations were carried out using these models to evaluate the influence of variations in cortical thickness, mesh density, bone material properties and modelling approach on the biomechanical responses of the clavicle, compared with experimental data. The FE results indicate that the inclusion of density-based bone material properties can provide a more accurate reproduction of the force–displacement response and bone fracture timing than a model with uniform bone material properties. Inclusion of a variable cortical thickness distribution also slightly improves the ability of the model to predict the experimental response. The methods developed in this study will be useful for creating subject-specific FE models to better understand the biomechanics and injury mechanism of the clavicle.     

Mapping anisotropy of the proximal femur for enhanced image based finite element analysis

Finite element (FE) models of bone derived from quantitative computed tomography (QCT) rely on realistic material properties to accurately predict bone strength. QCT cannot resolve bone microarchitecture, therefore QCT-based FE models lack the anisotropy apparent within the underlying bone tissue. This study proposes a method for mapping femoral anisotropy using high-resolution peripheral quantitative computed tomography (HR-pQCT) scans of human cadaver specimens. Femur HR-pQCT images were sub-divided into numerous overlapping cubic sub-volumes and the local anisotropy was quantified using a ‘direct-mechanics’ method. The resulting directionality reflected all the major stress lines visible within the trabecular lattice, and provided a realistic estimate of the alignment of Harvesian systems within the cortical compartment. QCT-based FE models of the proximal femur were constructed with isotropic and anisotropic material properties, with directionality interpolated from the map of anisotropy. Models were loaded in a sideways fall configuration and the resulting whole bone stiffness was compared to experimental stiffness and ultimate strength. Anisotropic models were consistently less stiff, but no statistically significant differences in correlation were observed between material models against experimental data. The mean difference in whole bone stiffness between model types was approximately 26%, suggesting that anisotropy can still effect considerable change in the mechanics of proximal femur models. The under prediction of whole bone stiffness in anisotropic models suggests that the orthotropic elastic constants require further investigation. The ability to map mechanical anisotropy from high-resolution images and interpolate information into clinical-resolution models will allow testing of new anisotropic material mapping strategies.     

The influence of material property and morphological parameters on specimen-specific finite element models of porcine vertebral bodies

The use of finite element (FE) methods in spinal research is increasing, but there is only limited information available on the influence of different input parameters on the model predictions. The aim of this study was to investigate the role of these parameters in FE models of the vertebral body. Experimental tests were undertaken on porcine lumbar vertebral bodies and scans of the specimens were used to create specimen-specific FE models. Three models were created for each specimen with combinations of generic and specimen-specific parameters. Stiffness and strength predictions were also made directly from the specimen trabecular bone volume fraction (BVF) and cross-sectional area (CSA). The agreement between the experimental results and the FE models with generic morphology was poorer (concordance coefficients = 0.058, 0.125 for stiffness, strength) than those made from the BVF and CSA (concordance coefficients = 0.638, 0.609). The greatest levels of agreement were found with the morphologically specific models including element-specific material properties (concordance coefficients = 0.881, 0.752). This indicates that highly specific models, both in terms of morphology and bone quality, are necessary if the FE tool is to be used effectively for spinal research and clinical practice.     

Finite element and experimental cortex strains of the intact and implanted tibia

Finite Element (FE) models for the simulation of intact and implanted bone find their main purpose in accurately reproducing the associated mechanical behavior. FE models can be used for preclinical testing of joint replacement implants, where some biomechanical aspects are difficult, if not possible, to simulate and investigate in vitro. To predict mechanical failure or damage, the models should accurately predict stresses and strains. Commercially available synthetic femur models have been extensively used to validate finite element models, but despite the vast literature available on the characteristics of synthetic tibia, numerical and experimental validation of the intact and implant assemblies of tibia are very limited or lacking. In the current study, four FE models of synthetic tibia, intact and reconstructed, were compared against experimental bone strain data, and an overall agreement within 10% between experimental and FE strains was obtained. Finite element and experimental (strain gauge) models of intact and implanted synthetic tibia were validated based on the comparison of cortex bone strains. The study also includes the analysis carried out on standard tibial components with cemented and noncemented stems of the P.F.C Sigma Modular Knee System. The overall agreement within 10% previously established was achieved, indicating that FE models could be successfully validated. The obtained results include a statistical analysis where the root-mean-square-error values were always <10%. FE models can successfully reproduce bone strains under most relevant acting loads upon the condylar surface of the tibia. Moreover, FE models, once properly validated, can be used for preclinical testing of tibial knee replacement, including misalignment of the implants in the proximal tibia after surgery, simulation of long-term failure according to the damage accumulation failure scenario, and other related biomechanical aspects.     

Experimental and probabilistic analysis of distal femoral periprosthetic fracture: a comparison of locking plate and intramedullary nail fixation. Part B: probabilistic investigation

The following is Part B of a two-part study. Part A evaluated, biomechanically, intramedullary (IM) nails versus locking plates for fixation of an extra-articular, metaphyseal wedge fracture in synthetic osteoporotic bone. Part B of this study introduces deterministic finite element (FE) models of each construct type in synthetic osteoporotic bone and investigates the probability of periprosthetic fracture of the locking plate compared with the retrograde IM nail using Monte Carlo simulation. Deterministic FE models of the fractured femur implanted with IM nail and locking plate, respectively, were developed and validated using experimental data presented in Part A of this study. The models were validated by comparing the load-displacement curve of the experimental data with the load-displacement curve of the FE simulation with a root-mean square error of less than 3?mm. The validated FE models were then modified by defining the cortical and cancellous bone modulus of elasticity as uncertain variables that could be assumed to vary randomly. Monte Carlo simulation was used to evaluate the probability of fracture (POF) of each fixation. The POF represents the cumulative probability that the predicted shear stresses in the cortical bone will exceed the expected shear strength of the cortical bone. This investigation provides information regarding the significance of post-operative damage accumulation on the POF of the implanted bones when the two fixations are used. The probabilistic analysis found the locking plate fixation to have a higher POF than the IM nail fixation under the applied loading conditions (locking plate 21.8% versus IM nail 0.019%).     

Adaptive glenoid bone remodeling simulation

Glenoid prosthesis loosening is the most common cause for revision total shoulder arthroplasty. Stress-induced bone remodeling may compromise long-term prosthesis fixation and significantly contribute to loosening. Realistic, robust analysis of bone-prosthesis constructs need to look beyond initial post-implantation mechanics provided by static finite element (FE) simulation. Adaptive bone remodeling simulations based on Wolff's law are needed for evaluating long-term glenoid prostheses fixation. The purpose of this study was to take a first step towards this goal and create and validate two-dimensional FE simulations, using the intact glenoid, for computing subject-specific adaptive glenoid remodeling. Two-dimensional glenoid FE models were created from scapulae computed tomography images. Two distinct processes, “element” and “node” simulations, used the forward-Euler method to compute bone remodeling. Initial bone density was homogeneous. Center and offset load combinations were iteratively applied. To validate the simulations we performed location-specific statistical comparisons between predicted and actual bone density, load combinations, and “element” and “node” processes. Visually and quantitatively “element” simulations produced better results (p>0.22), and correlation coefficients ranged 0.51–0.69 (p<0.001). Having met this initial work's goals, we expect subject-specific FE glenoid bone remodeling simulations together with static FE stress analyses to be effective tools for designing and evaluating glenoid prostheses.     

Predicting surface strains at the human distal radius during an in vivo loading task — Finite element model validation and application

Bone strains resulting from physical activity are thought to be a primary driver of bone adaptation, but cannot be directly noninvasively measured. Because bone adapts nonuniformly, physical activity may make an important independent structural contribution to bone strength that is independent of bone mass and density. Our objective was to create and validate methods for subject-specific finite element (FE) model generation that would accurately predict the surface strains experienced by the distal radius during an in vivo loading task, and to apply these methods to a group of 23 women aged 23–35 to examine variations in strain, bone mass and density, and physical activity. Four cadaveric specimens were experimentally tested and specimen-specific FE models were developed to accurately predict periosteal surface strains (root mean square error=16.3%). In the living subjects, when 300 N load was simulated, mean strains were significantly inversely correlated with BMC (r=−0.893), BMD (r=−0.892) and physical activity level (r=−0.470). Although the group of subjects was relatively homogenous, BMD varied by two-fold (range: 0.19–0.40 g/cm³) and mean energy-equivalent strain varied by almost six-fold (range: 226.79–1328.41 με) with a simulated 300 N load. In summary, we have validated methods for estimating surface strains in the distal radius that occur while leaning onto the palm of the hand. In our subjects, strain varied widely across individuals, and was inversely related to bone parameters that can be measured using clinical CT, and inversely related to physical activity history.