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.
Continuum finite element (FE) models of bones are commonly generated based on CT scans. Element material properties in such models are usually derived from bone density values using some empirical relationships. However, many different empirical relationships have been proposed. Most of these will provide isotropic material properties but relationships that can provide a full orthotropic elastic stiffness tensor have been proposed as well. Presently it is not clear which of these relationships best describes the material behavior of bone in continuum models, nor is it clear to what extent anisotropic models can improve upon isotropic models. The best way to determine the accuracy of such relationships for continuum analyses would be by quantifying the accuracy of the calculated stress/strain distribution, but this requires an accurate reference distribution that does not depend on such empirical relationships. In the present study, we propose a novel approach to generate such a reference stress distribution. With this approach, stress results obtained from a micro-FE model of a whole bone, that can represent the bone trabecular architecture in detail, are homogenized and the homogenized stresses are then used as a reference for stress results obtained from continuum models. The goal of the present study was to demonstrate this new approach and to provide examples of comparing continuum models with anisotropic versus isotropic material properties. Continuum models that implemented isotropic and orthotropic material definitions were generated for two proximal femurs for which micro-FE results were available as well, one representing a healthy and the other an osteoporotic femur. It was found that the continuum FE stress distributions calculated for the healthy femur compared well to the homogenized results of the micro-FE although slightly better for the orthotropic model (r=0.83) than for the isotropic model (r=0.79). For the osteoporotic bone also, the orthotropic model did better (r=0.83) than the isotropic model (r=0.77). We propose that this approach will enable a more relevant and accurate validation of different material models than experimental methods used so far.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.