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 Nonlinear Reduced-Order Model of the Corpus Callosum Under Planar Coronal Excitation
Abstract Traumatic brain injury (TBI) is often associated with microstructural tissue damage in the brain, which results from its complex biomechanical behavior. Recent studies have shown that the deep white matter (WM) region of the human brain is susceptible to being damaged due to strain localization in that region. Motivated by these studies, in this paper, we propose a geometrically nonlinear dynamical reduced order model (ROM) to model and study the dynamics of the deep WM region of the human brain under coronal excitation. In this model, the brain hemispheres were modeled as lumped masses connected via viscoelastic links, resembling the geometry of the corpus callosum (CC). Employing system identification techniques, we determined the unknown parameters of the ROM, and ensured the accuracy of the ROM by comparing its response against the response of an advanced finite element (FE) model. Next, utilizing modal analysis techniques, we determined the energy distribution among the governing modes of vibration of the ROM and concluded that the demonstrated nonlinear behavior of the FE model might be predominantly due to the special geometry of the brain deep WM region. Furthermore, we observed that, for sufficiently high input energies, high frequency harmonics at approximately 45 Hz, were generated in the response of the CC, which, in turn, are associated with high-frequency oscillations of the CC. Such harmonics might potentially lead to strain localization in the CC. This work is a step toward understanding the brain dynamics during traumatic injury.
more »
« less
- NSF-PAR ID:
- 10168967
- Date Published:
- Journal Name:
- Journal of Biomechanical Engineering
- Volume:
- 142
- Issue:
- 9
- ISSN:
- 0148-0731
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)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 and 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.more » « less
-
Abstract Computational simulations of traumatic brain injury (TBI) are commonly used to advance understanding of the injury–pathology relationship, tissue damage thresholds, and design of protective equipment such as helmets. Both human and animal TBI models have developed substantially over recent decades, partially due to the inclusion of more detailed brain geometry and representation of tissues like cerebral blood vessels. Explicit incorporation of vessels dramatically affects local strain and enables researchers to investigate TBI-induced damage to the vasculature. While some studies have indicated that cerebral arteries are rate-dependent, no published experimentally based, rate-sensitive constitutive models of cerebral arteries exist. In this work, we characterize the mechanical properties of axially failed porcine arteries, both quasi-statically (0.01 s−1) and at high rate (>100 s−1), and propose a rate-sensitive model to fit the data. We find that the quasi-static and high-rate stress–stretch curves become significantly different (p < 0.05) above a stretch of 1.23. We additionally find a significant change in both failure stretch and stress as a result of strain rate. The stress–stretch curve is then modeled as a Holzapfel–Gasser–Ogden material, with a Prony series added to capture the effects of viscoelasticity. Ultimately, this paper demonstrates that rate dependence should be considered in the material properties of cerebral arteries undergoing high strain-rate deformations and provides a ready-to-use model for finite element implementation.more » « less
-
Finite element analysis is used to study brain axonal injury and develop Brain White Matter (BWM) models while accounting for both the strain magnitude and the strain rate. These models are becoming more sophisticated and complicated due to the complex nature of the BMW composite structure with different material properties for each constituent phase. State-of-the-art studies focus on employing techniques that combine information about the local axonal directionality in different areas of the brain with diagnostic tools such as Diffusion-Weighted Magnetic Resonance Imaging (Diffusion-MRI). The diffusion-MRI data offers localization and orientation information of axonal tracks which are analyzed in finite element models to simulate virtual loading scenarios. Here, a BMW biphasic material model comprised of axons and neuroglia is considered. The model’s architectural anisotropy represented by a multitude of axonal orientations, that depend on specific brain regions, adds to its complexity. During this effort, we develop a finite element method to merge micro-scale Representative Volume Elements (RVEs) with orthotropic frequency domain viscoelasticity to an integrated macro-scale BWM finite element model, which incorporates local axonal orientation. Previous studies of this group focused on building RVEs that combined different volume fractions of axons and neuroglia and simulating their anisotropic viscoelastic properties. Via the proposed model, we can assign material properties and local architecture on each element based on the information from the orientation of the axonal traces. Consecutively, a BWM finite element model is derived with fully defined both material properties and material orientation. The frequency-domain dynamic response of the BMW model is analyzed to simulate larger scale diagnostic modalities such as MRI and MRE.more » « less
-
Traumatic brain injury (TBI) is a leading cause of long-term neurological disability in the world and the strongest environmental risk factor for the development of dementia. Even mild TBI (resulting from concussive injuries) is associated with a greater than twofold increase in the risk of dementia onset. Little is known about the cellular mechanisms responsible for the progression of long-lasting cognitive deficits. The integrated stress response (ISR), a phylogenetically conserved pathway involved in the cellular response to stress, is activated after TBI, and inhibition of the ISR—even weeks after injury—can reverse behavioral and cognitive deficits. However, the cellular mechanisms by which ISR inhibition restores cognition are unknown. Here, we used longitudinal two-photon imaging in vivo after concussive injury in mice to study dendritic spine dynamics in the parietal cortex, a brain region involved in working memory. Concussive injury profoundly altered spine dynamics measured up to a month after injury. Strikingly, brief pharmacological treatment with the drug-like small-molecule ISR inhibitor ISRIB entirely reversed structural changes measured in the parietal cortex and the associated working memory deficits. Thus, both neural and cognitive consequences of concussive injury are mediated in part by activation of the ISR and can be corrected by its inhibition. These findings suggest that targeting ISR activation could serve as a promising approach to the clinical treatment of chronic cognitive deficits after TBI.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1115/1.4046503
