An official website of the United States government
In image-based finite element analysis of bone, partial volume effects (PVEs) arise from image blur at tissue boundaries and as a byproduct of geometric reconstruction and meshing during model creation. In this study, we developed and validated a material assignment approach to mitigate partial volume effects. Our validation data consisted of physical torsion testing of intact tibiae from N = 20 Swiss alpine sheep. We created finite element models from micro-CT scans of these tibiae using three popular element types (10-node tetrahedral, 8-node hexahedral, and 20-node hexahedral). Without partial volume management, the models over-predicted the torsional rigidity compared to physical biomechanical tests. To address this problem, we implemented a dual-zone material model to treat elements that overlap low-density surface voxels as soft tissue rather than bone. After in situ inverse optimization, the dual-zone material model produced strong correlations and high absolute agreement between the virtual and physical tests. This suggests that with appropriate partial volume management, virtual mechanical testing can be a reliable surrogate for physical biomechanical testing. For maximum flexibility in partial volume management regardless of element type, we recommend the use of the following dual-zone material model for ovine tibiae: soft-tissue cutoff density of 665 mgHA/cm3 with a soft tissue modulus of 50 MPa (below cutoff) and a density-modulus conversion slope of 10,225 MPa-cm3/mgHA for bone (above cutoff).
more »
« less
Biomechanical duality of fracture healing captured using virtual mechanical testing and validated in ovine bones
Abstract Bone fractures commonly repair by forming a bridging structure called callus, which begins as soft tissue and gradually ossifies to restore rigidity to the bone. Virtual mechanical testing is a promising technique for image-based assessment of structural bone healing in both preclinical and clinical settings, but its accuracy depends on the validity of the material model used to assign tissue mechanical properties. The goal of this study was to develop a constitutive model for callus that captures the heterogeneity and biomechanical duality of the callus, which contains both soft tissue and woven bone. To achieve this, a large-scale optimization analysis was performed on 2363 variations of 3D finite element models derived from computed tomography (CT) scans of 33 osteotomized sheep under normal and delayed healing conditions. A piecewise material model was identified that produced high absolute agreement between virtual and physical tests by differentiating between soft and hard callus based on radiodensity. The results showed that the structural integrity of a healing long bone is conferred by an internal architecture of mineralized hard callus that is supported by interstitial soft tissue. These findings suggest that with appropriate material modeling, virtual mechanical testing is a reliable surrogate for physical biomechanical testing.
more »
« less
- Award ID(s):
- 1943287
- PAR ID:
- 10363048
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Scientific Reports
- Volume:
- 12
- Issue:
- 1
- ISSN:
- 2045-2322
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Bone healing has been traditionally described as a four-phase process: inflammatory response, soft callus formation, hard callus development, and remodeling. The remodeling phase has been largely neglected in most numerical mechanoregulation models of fracture repair in favor of capturing early healing using a pre-defined callus domain. However, in vivo evidence suggests that remodeling occurs concurrently with repair and causes changes in cortical bone adjacent to callus that are typically neglected in numerical models of bone healing. The objective of this study was to use image processing techniques to quantify this early-stage remodeling in ovine osteotomies. To accomplish this, we developed a numerical method for radiodensity profilometry with optimization-based curve fitting to mathematically model the bone density gradients in the radial direction across the cortical wall and callus. After assessing data from 26 sheep, we defined a dimensionless density fitting function that revealed significant remodeling occurring in the cortical wall adjacent to callus during early healing, a 23% average reduction in density compared to intact. This fitting function is robust for modeling radial density gradients in both intact bone and fracture repair scenarios and can capture a wide variety of the healing responses. The fitting function can also be scaled easily for comparison to numerical model predictions and may be useful for validating future mechanoregulatory models of coupled fracture repair and remodeling.more » « less
-
In large animal studies, the mechanical reintegration of the bone fragments is measured using postmortem physical testing, but these assessments can only be performed once, after sacrifice. Image‐based virtual mechanical testing is an attractive alternative because it could be used to monitor healing longitudinally. However, the procedures and software required to perform finite element analysis (FEA) on subject‐specific models for virtual mechanical testing can be time consuming and costly. Accordingly, the goal of this study was to determine whether a simpler image‐based geometric measure—the torsion constant, sometimes known as polar moment of inertia—can be reliably used as a surrogate measure of bone healing in large animals. To achieve this, postmortem biomechanical testing and microCT scans were analyzed for a total of 33 operated and 20 intact ovine tibiae. An image‐processing procedure to compute the attenuation‐weighted torsion constant from the microCT scans was developed in MATLAB and this code has been made freely available. Linear regression analysis was performed between the postmortem biomechanical data, the results of virtual mechanical testing using FEA, and the torsion constants measured from the scans. The results showed that virtual mechanical testing is the most reliable surrogate measure of postmortem torsional rigidity, having strong correlations and high absolute agreement. However, when FEA is not practical, the torsion constant is a viable alternative surrogate measure that is moderately correlated with postmortem torsional rigidity and can be readily calculated.more » « less
-
Trabecular bone is composed of organized mineralized collagen fibrils, which results in heterogeneous and anisotropic mechanical properties at the tissue level. Recently, biomechanical models computing stresses and strains in trabecular bone have indicated a significant effect of tissue heterogeneity on predicted stresses and strains. How-ever, the effect of the tissue-level mechanical anisotropy on the trabecular bone biomechanical response is unknown. Here, a computational method was established to automatically impose physiologically relevant orientation inherent in trabecular bone tissue on a trabecular bone microscale finite element model. Spatially varying tissue-level anisotropic elastic properties were then applied according to the bone mineral density and the local tissue orientation. The model was used to test the hypothesis that anisotropy in both homogeneous and heterogeneous models alters the predicted distribution of stress invariants. Linear elastic finite element computations were performed on a 3 mm cube model isolated from a microcomputed tomography scan of human trabecular bone from the distal femur. Hydrostatic stress and von Mises equivalent stress were recorded at every element, and the distributions of these values were analyzed. Anisotropy reduced the range of hydrostatic stress in both tension and compression more strongly than the associated increase in von Mises equivalent stress. The effect of anisotropy was independent of the spatial redistribution high compressive stresses due to tissue elastic heterogeneity. Tissue anisotropy and heterogeneity are likely important mechanisms to protect bone from failure and should be included for stress analyses in trabecular bone.more » « less
-
At the nanoscale bone is composed of aligned mineralized collagen fibrils organized into packets along the surface of trabecular bone creating an anisotropic tissue microstructure. Newer packets at the trabecular surfaces are usually less mineralized than older bone in the interior of the trabeculae, which along with irregular mineral deposition within packets, forms a heterogeneous material across the span of a trabeculae. However, finite element (FE) models of bone typically use homogenous isotropic material properties, because it is challenging to build anisotropy and heterogeneity into a model in a way that is applicable to the complex geometries of trabecular bone. Both the material anisotropy and heterogeneity may influence the stress state of trabecular bone, and it is important to understand the implications of such differences for determining bone biomechanical failure. It was hypothesized that taking into consideration both the tissue anisotropy and heterogeneity of bone’s biomechanical properties would alter the expected failure locations by reducing tensile stress on near surface elements of an FE model of canine trabecular bone. The objective of this study was to test this hypothesis and to develop a method to apply anisotropic and heterogeneous material properties to a model automatically from micro-computed tomography (μCT) data.more » « less
