Accurate characterization of the mechanical properties of the human brain at both microscopic and macroscopic length scales is a critical requirement for modeling of traumatic brain injury and brain folding. To date, most experimental studies that employ classical tension/compression/shear tests report the mechanical properties of the brain averaged over both the gray and white matter within the macroscopic regions of interest. As a result, there is a missing correlation between the independent mechanical properties of the microscopic constituent elements and the composite bulk macroscopic mechanical properties of the tissue. This microstructural computational study aims to inversely predict the hyperelastic mechanical properties of the axonal fibers and their surrounding extracellular matrix (ECM) from the bulk tissue's mechanical properties. We develop a representative volume element (RVE) model of the bulk tissue consisting of axonal fibers and ECM with the embedded element technique. A multiobjective optimization technique is implemented to calibrate the model and establish the independent mechanical properties of axonal fibers and ECM based on seven previously reported experimental mechanical tests for bulk white matter tissue from the corpus callosum. The result of the study shows that the discrepancy between the reported values for the elastic behavior of white matter in literature stems from the anisotropy of the tissue at the microscale. The shear modulus of the axonal fiber is seven times larger than the ECM, with axonal fibers that also show greater nonlinearity, contrary to the common assumption that both components exhibit identical nonlinear characteristics.
Statement of significance
The reported mechanical properties of white matter microstructure used in traumatic brain injury or brain mechanics studies vary widely, in some cases by up to two orders of magnitude. Currently, the material parameters of the white matter microstructure are identified by a single loading mode or ultimately two modes of the bulk tissue. The presented material models only define the response of the bulk and homogenized white matter at a macroscopic scale and cannot explicitly capture the connection between the material properties of microstructure and bulk structure. To fill this knowledge gap, our study characterizes the hyperelastic material properties of axonal fibers and ECM using microscale computational modeling and multiobjective optimization. The hyperelastic material properties for axonal fibers and ECM presented in this study are more accurate than previously proposed because they have been optimized using seven or six loading modes of the bulk tissue, which were previously limited to only two of the seven possible loading modes. As such, the predicted values with high accuracy could be used in various computational modeling studies. The systematic characterization of the material properties of the human brain tissue at both macro- and microscales will lead to more accurate computational predictions, which will enable a better understanding of injury criteria, and has a positive impact on the improved development of smart protection systems, and more accurate prediction of brain development and disease progression.
more »
« less
A theoretical framework for predicting the heterogeneous stiffness map of brain white matter tissue
Finding the stiffness map of biological tissues is of great importance in evaluating their healthy or pathological conditions. However, due to the heterogeneity and anisotropy of biological fibrous tissues, this task presents challenges and significant uncertainty when characterized only by single-mode loading experiments. In this study, we propose a new theoretical framework to map the stiffness landscape of fibrous tissues, specifically focusing on brain white matter tissue. Initially, a finite element (FE) model of the fibrous tissue was subjected to six loading cases, and their corresponding stress–strain curves were characterized. By employing multiobjective optimization, the material constants of an equivalent anisotropic material model were inversely extracted to best fit all six loading modes simultaneously. Subsequently, large-scale FE simulations were conducted, incorporating various fiber volume fractions and orientations, to train a convolutional neural network capable of predicting the equivalent anisotropic material properties solely based on the fibrous architecture of any given tissue. The proposed method, leveraging brain fiber tractography, was applied to a localized volume of white matter, demonstrating its effectiveness in precisely mapping the anisotropic behavior of fibrous tissue. In the long-term, the proposed method may find applications in traumatic brain injury, brain folding studies, and neurodegenerative diseases, where accurately capturing the material behavior of the tissue is crucial for simulations and experiments.
more » « less- PAR ID:
- 10552921
- Publisher / Repository:
- IOP Publishing
- Date Published:
- Journal Name:
- Physical Biology
- Volume:
- 21
- Issue:
- 6
- ISSN:
- 1478-3967
- Format(s):
- Medium: X Size: Article No. 066004
- Size(s):
- Article No. 066004
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Carbon Fiber Reinforced Plastics (CFRPs) are widely used due to their high stiffness to weight ratios. A common process manufacturers use to increase the strength to weight ratio is debulking. Debulking is the process of compacting a dry fibrous reinforcement prior to resin infusion. This process is meant to decrease the average inter-fiber distance, effectively increasing the fiber volume fraction of the sample. While this process is widely understood macroscopically its effects on fibrous microstructures have not yet been well characterized. The aim of this work is to compare the microstructures of three CFRP laminates, varying only the debulking step in the manufacturing process. High resolution serial sections of all three laminates were taken for analysis. Using these scans, the fiber positions were reconstructed. Statistical descriptors such as local fiber and void volume fractions, fiber orientation, and void distribution and morphology were then generated for each sample. Fiber clusters present within the material were identified and analyzed for each level of debulking applied. Using these descriptors, the effects of debulking on the morphology and organization of the composite microstructure was evaluated.more » « less
-
Wrinkling is a ubiquitous surface phenomenon in many biological tissues and is believed to play an important role in arterial health. As arteries are highly nonlinear, anisotropic, multilayered composite systems, it is necessary to investigate wrinkling incorporating these material characteristics. Several studies have examined surface wrinkling mechanisms with nonlinear isotropic material relationships. Nevertheless, wrinkling associated with anisotropic constitutive models such as Ogden–Gasser–Holzapfel (OGH), which is suitable for soft biological tissues, and in particular arteries, still requires investigation. Here, the effects of OGH parameters such as fibers’ orientation, stiffness, and dispersion on the onset of wrinkling, wrinkle wavelength and amplitude are elucidated through analysis of a bilayer system composed of a thin, stiff neo-Hookean membrane and a soft OGH substrate subjected to compression. Critical contractile strain at which wrinkles occur is predicted using both finite element analysis and analytical linear perturbation approach. Results suggest that besides stiffness mismatch, anisotropic features associated with fiber stiffness and distribution might be used in natural layered systems to adjust wrinkling and subsequent folding behaviors. Further analysis of a bilayer system with fibers in the (x–y) plane subjected to compression in the x direction shows a complex dependence of wrinkling strain and wavelength on fiber angle, stiffness, and dispersion. This behavior is captured by an approximation utilizing the linearized anisotropic properties derived from OGH model. Such understanding of wrinkling in this artery wall-like system will help identify the role of wrinkling mechanisms in biological artery in addition to the design of its synthetic counterparts.more » « less
-
more » « less
Abstract Tissues and engineered biomaterials exhibit exquisite local variation in stiffness that defines their function. Conventional elastography quantifies stiffness in soft (e.g. brain, liver) tissue, but robust quantification in stiff (e.g. musculoskeletal) tissues is challenging due to dissipation of high frequency shear waves. We describe new development of
finite deformation elastography that utilizes magnetic resonance imaging of low frequency, physiological-level (large magnitude) displacements, coupled to an iterative topology optimization routine to investigate stiffness heterogeneity, including spatial gradients and inclusions. We reconstruct 2D and 3D stiffness distributions in bilayer agarose hydrogels and silicon materials that exhibit heterogeneous displacement/strain responses. We map stiffness in porcine and sheep articular cartilage deep within the bony articular joint spacein situ for the first time. Elevated cartilage stiffness localized to the superficial zone is further related to collagen fiber compaction and loss of water content during cyclic loading, as assessed by independentT 2measurements. We additionally describe technical challenges needed to achievein vivo elastography measurements. Our results introduce new functional imaging biomarkers, which can be assessed nondestructively, with clinical potential to diagnose and track progression of disease in early stages, including osteoarthritis or tissue degeneration. -
more » « less
In tissues and organs, the extracellular matrix (ECM) helps maintain inter- and intracellular architectures that sustain the structure–function relationships defining physiological homeostasis. Combining fiber scaffolds and cells to form engineered tissues is a means of replicating these relationships. Engineered tissues' fiber scaffolds are designed to mimic the topology and chemical composition of the ECM network. Here, we asked how cells found in the heart compare in their propensity to align their cytoskeleton and self-organize in response to topological cues in fibrous scaffolds. We studied cardiomyocytes, valvular interstitial cells, and vascular endothelial cells as they adapted their inter- and intracellular architectures to the extracellular space. We used focused rotary jet spinning to manufacture aligned fibrous scaffolds to mimic the length scale and three-dimensional (3D) nature of the native ECM in the muscular, valvular, and vascular tissues of the heart. The representative cardiovascular cell types were seeded onto fiber scaffolds and infiltrated the fibrous network. We measured different cell types' propensity for cytoskeletal alignment in response to fiber scaffolds with differing levels of anisotropy. The results indicated that valvular interstitial cells on moderately anisotropic substrates have a higher propensity for cytoskeletal alignment than cardiomyocytes and vascular endothelial cells. However, all cell types displayed similar levels of alignment on more extreme (isotropic and highly anisotropic) fiber scaffold organizations. These data suggest that in the hierarchy of signals that dictate the spatiotemporal organization of a tissue, geometric cues within the ECM and cellular networks may homogenize behaviors across cell populations and demographics.
