Reliable simulations of the human proximal femur by high-order finite element analysis validated by experimental observations

Background: The mechanical response of patient-specific bone to various load conditions is of major clinical importance in orthopedics. Herein we enhance the methods presented in Yosibash et al. [2007. A CT-based high-order finite element analysis of the human proximal femur compared to in-vitro experiments. ASME Journal of Biomechanical Engineering 129(3), 297–309.] for the reliable simulations of the human proximal femur by high-order finite elements (FEs) and validate the simulations by experimental observations.

Method of approach: A fresh-frozen human femur was scanned by quantitative computed tomography (QCT) and thereafter loaded (in vitro experiments) by a quasi-static force of up to 1250 N. QCT scans were manipulated to generate a high-order FE bone model with distinct cortical and trabecular regions having inhomogeneous isotropic elastic properties with Young's modulus represented by continuous spatial functions. Sensitivity analyses were performed to quantify parameters that mostly influence the mechanical response. FE results were compared to displacements and strains measured in the experiments.

Results: Young moduli correlated to QCT Hounsfield Units by relations in Keyak and Falkinstein [2003. Comparison of in situ and in vitro CT scan-based finite element model predictions of proximal femoral fracture load. Medical Engineering and Physics 25, 781–787.] were found to provide predictions that match the experimental results closely. Excellent agreement was found for both the displacements and strains. The presented study demonstrates that reliable and validated high-order patient-specific FE simulations of human femurs based on QCT data are achievable for clinical computer-aided decision making.     

A CT-based high-order finite element analysis of the human proximal femur compared to in-vitro experiments

The prediction of patient-specific proximal femur mechanical response to various load conditions is of major clinical importance in orthopaedics. This paper presents a novel, empirically validated high-order finite element method (FEM) for simulating the bone response to loads. A model of the bone geometry was constructed from a quantitative computerized tomography (QCT) scan using smooth surfaces for both the cortical and trabecular regions. Inhomogeneous isotropic elastic properties were assigned to the finite element model using distinct continuous spatial fields for each region. The Young's modulus was represented as a continuous function computed by a least mean squares method. p-FEMs were used to bound the simulation numerical error and to quantify the modeling assumptions. We validated the FE results with in-vitro experiments on a fresh-frozen femur loaded by a quasi-static force of up to 1500 N at four different angles. We measured the vertical displacement and strains at various locations and investigated the sensitivity of the simulation. Good agreement was found for the displacements, and a fair agreement found in the measured strain in some of the locations. The presented study is a first step toward a reliable p-FEM simulation of human femurs based on QCT data for clinical computer aided decision making.     

Prediction of strength and strain of the proximal femur by a CT-based finite element method

Hip fractures are the most serious complication of osteoporosis and have been recognized as a major public health problem. In elderly persons, hip fractures occur as a result of increased fragility of the proximal femur due to osteoporosis. It is essential to precisely quantify the strength of the proximal femur in order to estimate the fracture risk and plan preventive interventions. CT-based finite element analysis could possibly achieve precise assessment of the strength of the proximal femur. The purpose of this study was to create a simulation model that could accurately predict the strength and surface strains of the proximal femur using a CT-based finite element method and to verify the accuracy of our model by load testing using fresh frozen cadaver specimens. Eleven right femora were collected. The axial CT scans of the proximal femora were obtained with a calibration phantom, from which the 3D finite element models were constructed. Materially nonlinear finite element analyses were performed. The yield and fracture loads were calculated, while the sites where elements failed and the distributions of the principal strains were determined. The strain gauges were attached to the proximal femoral surfaces. A quasi-static compression test of each femur was conducted. The yield loads, fracture loads and principal strains of the prediction significantly correlated with those measured (r=0.941, 0.979, 0.963). Finite element analysis showed that the solid elements and shell elements in undergoing compressive failure were at the same subcapital region as the experimental fracture site.     

Prediction of the mechanical response of the femur with uncertain elastic properties

A mandatory requirement for any reliable prediction of the mechanical response of bones, based on quantitative computer tomography, is an accurate relationship between material properties (usually Young's modulus E) and bone density ρ. Many such E-ρ relationships are available based on different experiments on femur specimens with a large spread due to uncertainties. The first goal of this study is to pool and analyze the relevant available experimental data and develop a stochasticE-ρ relationship. This analysis highlights that there is no experimental data available to cover the entire density range of the human femur and that some "popular" E-ρ relationships are based on data that contains extreme scatter, while others are based on a very limited amount of information. The second goal is to use the newly developed stochastic E-ρ relationship in high-order finite element analyses (FEAs) for the computation of strains and displacements in two human proximal femurs, mimicking in vitro experiments. When compared with the experimental observations, the FEA predictions using the median of the stochastic E-ρ relationship follow the underlying distribution of the stochastic E-ρ relationship. Thus, most deviations of the FEA predictions from experimental observations can possibly be explained by uncertain elastic properties of the femur.     

Use of a statistical model of the whole femur in a large scale,multi-model study of femoral neck fracture risk

Interpatient variability is often overlooked in orthopaedic computational studies due to the substantial challenges involved in sourcing and generating large numbers of bone models. A statistical model of the whole femur incorporating both geometric and material property variation was developed as a potential solution to this problem. The statistical model was constructed using principal component analysis, applied to 21 individual computer tomography scans. To test the ability of the statistical model to generate realistic, unique, finite element (FE) femur models it was used as a source of 1000 femurs to drive a study on femoral neck fracture risk. The study simulated the impact of an oblique fall to the side, a scenario known to account for a large proportion of hip fractures in the elderly and have a lower fracture load than alternative loading approaches. FE model generation, application of subject specific loading and boundary conditions, FE processing and post processing of the solutions were completed automatically. The generated models were within the bounds of the training data used to create the statistical model with a high mesh quality, able to be used directly by the FE solver without remeshing. The results indicated that 28 of the 1000 femurs were at highest risk of fracture. Closer analysis revealed the percentage of cortical bone in the proximal femur to be a crucial differentiator between the failed and non-failed groups. The likely fracture location was indicated to be intertrochantic. Comparison to previous computational, clinical and experimental work revealed support for these findings.     

To what extent can linear finite element models of human femora predict failure under stance and fall loading configurations?

Proximal femur strength estimates from computed tomography (CT)-based finite element (FE) models are finding clinical application. Published models reached a high in-vitro accuracy, yet many of them rely on nonlinear methodologies or internal best-fitting of parameters. The aim of the present study is to verify to what extent a linear FE modelling procedure, fully based on independently determined parameters, can predict the failure characteristics of the proximal femur in stance and sideways fall loading configurations.     

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.     

Patient-specific finite-element analyses of the proximal femur with orthotropic material properties validated by experiments

Patient-specific high order finite-element (FE) models of human femurs based on quantitative computer tomography (QCT) with inhomogeneous orthotropic and isotropic material properties are addressed. The point-wise orthotropic properties are determined by a micromechanics (MM) based approach in conjunction with experimental observations at the osteon level, and two methods for determining the material trajectories are proposed (along organs outer surface, or along principal strains). QCT scans on four fresh-frozen human femurs were performed and high-order FE models were generated with either inhomogeneous MM-based orthotropic or empirically determined isotropic properties. In vitro experiments were conducted on the femurs by applying a simple stance position load on their head, recording strains on femurs' surface and head's displacements. After verifying the FE linear elastic analyses that mimic the experimental setting for numerical accuracy, we compared the FE results to the experimental observations to identify the influence of material properties on models' predictions. The strains and displacements computed by FE models having MM-based inhomogeneous orthotropic properties match the FE-results having empirically based isotropic properties well, and both are in close agreement with the experimental results. When only the strains in the femoral neck are being compared a more pronounced difference is noticed between the isotropic and orthotropic FE result. These results lay the foundation for applying more realistic inhomogeneous orthotropic material properties in FEA of femurs.     

Sensitivity of periprosthetic stress-shielding to load and the bone density-modulus relationship in subject-specific finite element models

Subject-specific finite element (FE) computer models of the proximal femur in hip replacement could potentially predict stress-shielding and subsequent bone loss in individual patients. Before such predictions can be made, it is important first to determine if between subject differences in stress-shielding are sensitive to poorly defined parameters such as the load and the bone material properties. In this study we investigate if subject-specific FE models provide consistent stress-shielding patterns in the bone, independent of the choice of the loading conditions and the bone density-modulus relationship used in the computer model. FE models of two right canine femurs with and without implants were constructed based on contiguous computed tomography (CT) scans so that subject-specific estimates of stress-shielding could be calculated. Four different loading conditions and two bone density-modulus relationships were tested. Stress-shielding was defined as the decrease of strain energy per gram bone mass in the femur with the implant in place relative to the intact femur.The analyses showed that for the four loading conditions and two bone density-modulus relationships the difference in stress-shielding between the two subjects was essentially constant (1% variation) when the same loading condition and density-modulus relationship was used for both subjects. The severity of stress-shielding within a subject was sensitive to these input parameters, varying up to 20% in specific regions with a change in loading conditions and up to 10% for a change in the assumed density-modulus relationship. We conclude that although the choice of input parameters can substantially affect stress-shielding in an individual, this choice had virtually no effect on the relative differences in femoral periprosthetic stress-shielding between individuals. Thus, while care should be taken in the interpretation of the absolute value of stress-shielding calculated with these type of models, subject-specific FE models may be useful for explaining the variation in bone adaptation responsiveness between different subjects in experimental or clinical studies.     

The effect of boundary condition on the biomechanics of a human pelvic joint under an axial compressive load: a three-dimensional finite element model

The finite element (FE) model of the pelvic joint is helpful for clinical diagnosis and treatment of pelvic injuries. However, the effect of an FE model boundary condition on the biomechanical behavior of a pelvic joint has not been well studied. The objective of this study was to study the effect of boundary condition on the pelvic biomechanics predictions. A 3D FE model of a pelvis using subject-specific estimates of intact bone structures, main ligaments and bone material anisotropy by computed tomography (CT) gray value was developed and validated by bone surface strains obtained from rosette strain gauges in an in vitro pelvic experiment. Then three FE pelvic models were constructed to analyze the effect of boundary condition, corresponding to an intact pelvic joint, a pelvic joint without sacroiliac ligaments and a pelvic joint without proximal femurs, respectively. Vertical load was applied to the same pelvis with a fixed prosthetic femoral stem and the same load was simulated in the FE model. A strong correlation coefficient (R(2)=0.9657) was calculated, which indicated a strong correlation between the FE analysis and experimental results. The effect of boundary condition changes on the biomechanical response depended on the anatomical location and structure of the pelvic joint. It was found that acetabulum fixed in all directions with the femur removed can increase the stress distribution on the acetabular inner plate (approximately double the original values) and decrease that on the superior of pubis (from 7 MPa to 0.6 MPa). Taking sacrum and ilium as a whole, instead of sacroiliac and iliolumber ligaments, can influence the stress distribution on ilium and pubis bone vastly. These findings suggest pelvic biomechanics is very dependent on the boundary condition in the FE model.     

Incorporating uncertainty in mechanical properties for finite element-based evaluation of bone mechanics

Finite element (FE) models of bone, developed from computed tomography (CT) scan data, are used to evaluate stresses and strains, load transfer and fixation of implants, and potential for fracture. The experimentally derived relationships used to transform CT scan data in Hounsfield unit to modulus and strength contain substantial scatter. The scatter in these relationships has potential to impact the results and conclusions of bone studies. The objectives of this study were to develop a computationally efficient probabilistic FE-based platform capable of incorporating uncertainty in bone property relationships, and to apply the model to a representative analysis; variability in stresses and fracture risk was predicted in five proximal femurs under stance loading conditions. Based on published variability in strength and modulus relationships derived in the proximal femur, the probabilistic analysis predicted the distributions of stress and risk. For the five femurs analyzed, the 1 and 99 percentile bounds varied by an average of 17.3 MPa for stress and by 0.28 for risk. In each femur, the predicted variability in risk was greater than 50% of the mean risk calculated, with obvious implications for clinical assessment. Results using the advanced mean value (AMV) method required only seven analysis trials (1h) and differed by less than 2% when compared to a 1000-trial Monte-Carlo simulation (400 h). The probabilistic modeling platform developed has broad applicability to bone studies and can be similarly implemented to investigate other loading conditions, structures, sources of uncertainty, or output measures of interest.     

The use of sparse CT datasets for auto-generating accurate FE models of the femur and pelvis

The finite element (FE) method when coupled with computed tomography (CT) is a powerful tool in orthopaedic biomechanics. However, substantial data is required for patient-specific modelling. Here we present a new method for generating a FE model with a minimum amount of patient data. Our method uses high order cubic Hermite basis functions for mesh generation and least-square fits the mesh to the dataset. We have tested our method on seven patient data sets obtained from CT assisted osteodensitometry of the proximal femur. Using only 12 CT slices we generated smooth and accurate meshes of the proximal femur with a geometric root mean square (RMS) error of less than 1 mm and peak errors less than 8 mm. To model the complex geometry of the pelvis we developed a hybrid method which supplements sparse patient data with data from the visible human data set. We tested this method on three patient data sets, generating FE meshes of the pelvis using only 10 CT slices with an overall RMS error less than 3 mm. Although we have peak errors about 12 mm in these meshes, they occur relatively far from the region of interest (the acetabulum) and will have minimal effects on the performance of the model. Considering that linear meshes usually require about 70-100 pelvic CT slices (in axial mode) to generate FE models, our method has brought a significant data reduction to the automatic mesh generation step. The method, that is fully automated except for a semi-automatic bone/tissue boundary extraction part, will bring the benefits of FE methods to the clinical environment with much reduced radiation risks and data requirement.     

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.     

Strain and micromotion in intact and resurfaced composite femurs: Experimental and numerical investigations

Understanding the load transfer within a resurfaced femur is necessary to determine the influence of mechanical factors on potential failure mechanisms such as early femoral neck fractures and stress shielding. In this study, an attempt has been made to measure the stem-bone micromotion and implant cup-bone relative displacements (along medial-lateral and anterior-posterior direction), in addition to surface strains at different locations and orientations on the proximal femur and to compare these measurements with those predicted by equivalent FE models. The loading and the support conditions of the experiment were closely replicated in the FE models. A new experimental set-up has been developed, with specially designed fixtures and load application mechanism, which can effectively impose bending and deflection of the tested femurs, almost in any direction. High correlation coefficient (0.92–0.95), low standard error of the estimate (170–379 με) and low percentage error in regression slope (12.8–17.5%), suggested good agreement between the numerical and measured strains. The effect of strain shielding was observed in two (out of eight) strain gauges located on the posterior side. A pronounced strain increase occurred in strain gauges located on the anterior head and neck regions after implantation. Experimentally measured stem-bone micromotion and implant cup-bone relative displacements (0–13.7 μm) were small and similar in trends predicted by the FE models (0–25 μm). Despite quantitative deviations in the measured and numerical results, it appears that the FE model can be used as a valid predictor of the actual strain and stem-bone micromotion.     

A novel methodology for generating 3D finite element models of the hip from 2D radiographs

Finite element (FE) modelling has been proposed as a tool for estimating fracture risk and patient-specific FE models are commonly based on computed tomography (CT). Here, we present a novel method to automatically create personalised 3D models from standard 2D hip radiographs. A set of geometrical parameters of the femur were determined from seven a–p hip radiographs and compared to the 3D femoral shape obtained from CT as training material; the error in reconstructing the 3D model from the 2D radiographs was assessed. Using the geometry parameters as the input, the 3D shape of another 21 femora was built and meshed, separating a cortical and trabecular compartment. The material properties were derived from the homogeneity index assessed by texture analysis of the radiographs, with focus on the principal tensile and compressive trabecular systems. The ability of these FE models to predict failure load as determined by experimental biomechanical testing was evaluated and compared to the predictive ability of DXA. The average reconstruction error of the 3D models was 1.77 mm (±1.17 mm), with the error being smallest in the femoral head and neck, and greatest in the trochanter. The correlation of the FE predicted failure load with the experimental failure load was r²=64% for the reconstruction FE model, which was significantly better (p<0.05) than that for DXA (r²=24%). This novel method for automatically constructing a patient-specific 3D finite element model from standard 2D radiographs shows encouraging results in estimating patient-specific failure loads.     

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).     

Computer-integrated finite element modeling of human middle ear

 The objective of this study was to produce an improved finite element (FE) model of the human middle ear and to compare the model with human data. We began with a systematic and accurate geometric modeling technique for reconstructing the middle ear from serial sections of a freshly frozen temporal bone. A geometric model of a human middle ear was constructed in a computer-aided design (CAD) environment with particular attention to geometry and microanatomy. Using the geometric model, a working FE model of the human middle ear was created using previously published material properties of middle ear components. This working FE model was finalized by a cross-calibration technique, comparing its predicted stapes footplate displacements with laser Doppler interferometry measurements from fresh temporal bones. The final FE model was shown to be reasonable in predicting the ossicular mechanics of the human middle ear. Received: 18 February 2002 / Accepted: 6 June 2002 The preparation of temporal bone histological sections of Robert K. Dyer, Jr., MD is gratefully recognized. The Whitaker Foundation supported this work (Research Grant RG-98-0305).     

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.     

Computer aided stress analysis of long bones utilizing computed tomography

A computer aided analysis method has been developed which utilizes computed tomography (CT) and a finite element (FE) computer program to determine the stress-displacement pattern in a long bone section. The CT data file provides the geometry, the apparent density and the elastic properties for the three-dimensional FE model. A developed pre-processor generates the FE model of a human diaphyseal tibia section which is then analyzed by the SAP IV finite element program. The results obtained are sorted and displayed by a developed post-processor and compared with stresses and deformations from the literature. The model generation method was verified by applying it to a model of simple geometry and boundary conditions, then comparing the results with the analytical solution of the same problem. The convergence behavior of nodal displacements was tested as a function of mesh refinement. This method provides an automatic, versatile, non-invasive and accurate tool of long bone modeling for finite element stress analysis.     

Finite element analysis of the impact of bone nanostructure on its piezoelectric response

Biomechanics and Modeling in Mechanobiology - The piezoelectric response of bone at the submicron scale is analyzed under mechanical loadings using the finite element (FE) method. A new algorithm...