skip to main content


Title: Incorporating tissue anisotropy and heterogeneity in finite element models of trabecular bone altered predicted local stress distributions
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
Award ID(s):
1643164
NSF-PAR ID:
10049929
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Biomechanics and Modeling in Mechanobiology
ISSN:
1617-7959
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Abstract

    Meniscal tears are a common, painful, and debilitating knee injury with limited treatment options. Computational models that predict meniscal tears may help advance injury prevention and repair, but first these models must be validated using experimental data. Here we simulated meniscal tears with finite element analysis using continuum damage mechanics (CDM) in a transversely isotropic hyperelastic material. Finite element models were built to recreate the coupon geometry and loading conditions of forty uniaxial tensile experiments of human meniscus that were pulled to failure either parallel or perpendicular to the preferred fiber orientation. Two damage criteria were evaluated for all experiments: von Mises stress and maximum normal Lagrange strain. After we successfully fit all models to experimental force–displacement curves (grip-to-grip), we compared model predicted strains in the tear region at ultimate tensile strength to the strains measured experimentally with digital image correlation (DIC). In general, the damage models underpredicted the strains measured in the tear region, but models using von Mises stress damage criterion had better overall predictions and more accurately simulated experimental tear patterns. For the first time, this study has used DIC to expose strengths and weaknesses of using CDM to model failure behavior in soft fibrous tissue.

     
    more » « less
  3. Abstract

    A novel finite element method (FEM) is developed to study mechanical response of axons embedded in extra cellular matrix (ECM) when subjected to harmonic uniaxial stretch under purely non-affine kinematic boundary conditions. The proposed modeling approach combines hyper-elastic (such as Ogden model) and time/frequency domain viscoelastic constitutive models to evaluate the effect of parametrically varying oligodendrocyte-axon tethering under harmonic stretch at 50Hz. A hybrid hyper-viscoelastic material (HVE) model enabled the analysis of repeated uniaxial load on stress propagation and damage accumulation in white matter.

    In the proposed FEM, oligodendrocyte connections to axons are depicted via a spring-dashpot model. This tethering technique facilitates contact definition at various locations, parameterizes connection points and varies stiffness of connection hubs. Results from a home-grown FE submodel configuration of a single oligodendrocyte tethered to axons at various locations are presented. Root mean square deviation (RMSD) are computed between stress-strain plots to depict trends in mechanical response. Steady-state dynamic (SSD) simulations show stress relaxation in axons. Gradual axonal softening under repetitive loads is illustrated employing Prony series - HVE models. Representative von-Mises stress plots indicate that undulated axons experience bending stresses along their tortuous path, suggesting greater susceptibility to damage accumulation and fatigue failure due to repeated strains.

     
    more » « less
  4. This publication contains a finite element model for the analysis of bone core under consideration of bone tissue heterogeneity and tissue anisotropy. The model for bone tissue heterogeneity and anisotropy follows: Hammond, M.A., Wallace, J.M., Allen, M.R. and Siegmund, T., 2018. Incorporating tissue anisotropy and heterogeneity in finite element models of trabecular bone altered predicted local stress distributions. Biomechanics and Modeling in Mechanobiology, 17(2), pp.605-614. In this publication the finite element model, material set assignment and local orientations are provided. This dataset contains an inp file in the syntax of Abaqus/Standard software v2017. 
    more » « less
  5. This publication contains a finite element model for the analysis of single bone trabeculae under consideration of bone tissue heterogeneity and tissue anisotropy. The model for bone tissue heterogeneity and anisotropy follows: Hammond, M.A., Wallace, J.M., Allen, M.R. and Siegmund, T., 2018. Incorporating tissue anisotropy and heterogeneity in finite element models of trabecular bone altered predicted local stress distributions. Biomechanics and Modeling in Mechanobiology, 17(2), pp.605-614. In this publication the finite element model, material set assignment and local orientations are provided. This dataset contains an inp file in the syntax of Abaqus/Standard software v2017. 
    more » « less