skip to main content

Title: Finite Element Models of Linear Microcracks in Trabecular Bone with Simulated Bisphosphonate and Raloxifene Treatment
At the nanoscale bone is composed of aligned heterogeneously mineralized collagen fibrils. While raloxifene (Ral) and bisphosphonate (BP) treatment preserve bone mass, they also affect bone quality through changes in collagen hydration and mineral density/heterogeneity, respectively. It was hypothesized that the effects of pharmacological treatment on the tissue would alter linear microcracking in finite element (FE) models of trabeculae reflecting control (Ctrl), Ral and BP. A FE mesh of a single canine vertebral body trabecula was generated from a micro- CT scan using ScanIP. A custom MATLAB code imposed tissue property heterogeneity and a collagen fibril orientation parallel to the trabecular surface. Ctrl was heterogeneous (based on vBMD) in both modulus and strength, and BP was homogenous (+25% of Ctrl mean modulus and strength). Ctrl and BP models had identical microcracking toughness. Ral had increased microcracking toughness (+25%) and the same modulus and strength heterogeneity as Ctrl. Transverse deflections were applied to simulate bending of the trabeculae, microcrack formation and propagation was simulated with the imposed orientation using the extended FE method in Abaqus/Standard, and the energy dissipated by the microcrack was assessed.
; ; ;
Award ID(s):
Publication Date:
Journal Name:
Journal of bone and mineral research
Page Range or eLocation-ID:
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.more »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.« less
  2. Abstract Background The pia arachnoid complex (PAC) is a cerebrospinal fluid-filled tissue conglomerate that surrounds the brain and spinal cord. Pia mater adheres directly to the surface of the brain while the arachnoid mater adheres to the deep surface of the dura mater. Collagen fibers, known as subarachnoid trabeculae (SAT) fibers, and microvascular structure lie intermediately to the pia and arachnoid meninges. Due to its structural role, alterations to the biomechanical properties of the PAC may change surface stress loading in traumatic brain injury (TBI) caused by sub-concussive hits. The aim of this study was to quantify the mechanical andmore »morphological properties of ovine PAC. Methods Ovine brain samples (n = 10) were removed from the skull and tissue was harvested within 30 min post-mortem. To access the PAC, ovine skulls were split medially from the occipital region down the nasal bone on the superior and inferior aspects of the skull. A template was used to remove arachnoid samples from the left and right sides of the frontal and occipital regions of the brain. 10 ex-vivo samples were tested with uniaxial tension at 2 mm s −1 , average strain rate of 0.59 s −1 , until failure at < 5 h post extraction. The force and displacement data were acquired at 100 Hz. PAC tissue collagen fiber microstructure was characterized using second-harmonic generation (SHG) imaging on a subset of n = 4 stained tissue samples. To differentiate transverse blood vessels from SAT by visualization of cell nuclei and endothelial cells, samples were stained with DAPI and anti-von Willebrand Factor, respectively. The Mooney-Rivlin model for average stress–strain curve fit was used to model PAC material properties. Results The elastic modulus, ultimate stress, and ultimate strain were found to be 7.7 ± 3.0, 2.7 ± 0.76 MPa, and 0.60 ± 0.13, respectively. No statistical significance was found across brain dissection locations in terms of biomechanical properties. SHG images were post-processed to obtain average SAT fiber intersection density, concentration, porosity, tortuosity, segment length, orientation, radial counts, and diameter as 0.23, 26.14, 73.86%, 1.07 ± 0.28, 17.33 ± 15.25 µm, 84.66 ± 49.18°, 8.15%, 3.46 ± 1.62 µm, respectively. Conclusion For the sizes, strain, and strain rates tested, our results suggest that ovine PAC mechanical behavior is isotropic, and that the Mooney-Rivlin model is an appropriate curve-fitting constitutive equation for obtaining material parameters of PAC tissues.« less
  3. 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 densitymore »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.« less
  4. 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.
  5. Abstract

    Corneal stiffness plays a critical role in shaping the cornea with respect to intraocular pressure and physical interventions. However, it remains difficult to measure the mechanical properties noninvasively. Here, we report the first measurement of shear modulus in human corneas in vivo using optical coherence elastography (OCE) based on surface elastic waves. In a pilot study of 12 healthy subjects aged between 25 and 67, the Rayleigh-wave speed was 7.86 ± 0.75 m/s, corresponding to a shear modulus of 72 ± 14 kPa. Our data reveal two unexpected trends: no correlation was found between the wave speed and IOP between 13–18 mmHg, and shear modulus decreasesmore »with age (− 0.32 ± 0.17 m/s per decade). We propose that shear stiffness is governed by the interfibrillar matrix, whereas tensile strength is dominated by collagen fibrils. Rayleigh-wave OCE may prove useful for clinical diagnosis, refractive surgeries, and treatment monitoring.

    « less