skip to main content


Title: Brain White Matter Model of Orthotropic Viscoelastic Properties in Frequency Domain
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
Award ID(s):
1763005
NSF-PAR ID:
10190337
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
ASME INternational Mechanical Engineering Congress and Exposition
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Material properties of brain white matter (BWM) show high anisotropy due to the complicated internal three-dimensional microstructure and variant interaction between heterogeneous brain-tissue (axon, myelin, and glia). From our previous study, finite element methods were used to merge micro-scale Representative Volume Elements (RVE) with orthotropic frequency domain viscoelasticity to an integral macro-scale BWM. Quantification of the micro-scale RVE with anisotropic frequency domain viscoelasticity is the core challenge in this study.

    The RVE behavior is expressed by a viscoelastic constitutive material model, in which the frequency-related viscoelastic properties are imparted as storage modulus and loss modulus for the composite comprised of axonal fibers and extracellular glia. Using finite elements to build RVEs with anisotropic frequency domain viscoelastic material properties is computationally very consuming and resource-draining. Additionally, it is very challenging to build every single RVE using finite elements since the architecture of each RVE is arbitrary in an infinite data set. The architecture information encoded in the voxelized location is employed as input data and is consequently incorporated into a deep 3D convolution neural network (CNN) model that cross-references the RVEs’ material properties (output data). The output data (RVEs’ material properties) is calculated in parallel using an in-house developed finite element method, which models RVE samples of axon-myelin-glia composites. This novel combination of the CNN-RVE method achieved a dramatic reduction in the computation time compared with directly using finite element methods currently present in the literature. 

    more » « less
  2. 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
  3. Abstract

    Numerical simulations using non-linear hyper-elastic material models to describe interactions between brain white matter (axons and extra cellular matrix (ECM)) have enabled high-fidelity characterization of stress-strain response. In this paper, a novel finite element model (FEM) has been developed to study mechanical response of axons embedded in ECM when subjected to tensile loads under purely non-affine kinematic boundary conditions. FEM leveraging Ogden hyper-elastic material model is deployed to understand impact of parametrically varying oligodendrocyte-axon tethering and analyze influence of aging material characteristics on stress propagation. In proposed FEM, oligodendrocyte connections to axons are represented via spring-dashpot model, such tethering technique facilitates contact definition at various locations, parameterize connection points and vary stiffness of connection hubs. Two FE submodels are discussed: 1) multiple oligodendrocytes arbitrarily tethered to the nearest axons, and 2) single oligodendrocyte tethered to all axons at various locations. Root mean square deviation (RMSD) were computed between stress-strain plots to depict trends in mechanical response. Axonal stiffness was found to rise with increasing tethering, indicating role of oligodendrocytes in stress redistribution. Finally, stress state results for aging axon material, with varying stiffnesses and number of connections in FEM ensemble have also been discussed to demonstrate gradual softening of tissues.

     
    more » « less
  4. null (Ed.)
    Abstract

    Composites can be tailored to specific applications by adjusting process variables. These variables include those related to composition, such as volume fraction of the constituents and those associated with processing methods, methods that can affect composite topology. In the case of particle matrix composites, orientation of the inclusions affects the resulting composite properties, particularly so in instances where the particles can be oriented and arranged into structures. In this work, we study the effects of coupled electric and magnetic field processing with externally applied fields on those structures, and consequently on the resulting material properties that arise. The ability to vary these processing conditions with the goal of generating microstructures that yield target material properties adds an additional level of control to the design of composite material properties. Moreover, while analytical models allow for the prediction of resulting composite properties from constituents and composite topology, these models do not build upward from process variables to make these predictions.

    This work couples simulation of the formation of microscale architectures, which result from coupled electric and magnetic field processing of particulate filled polymer matrix composites, with finite element analysis of those structures to provide a direct and explicit linkages between process, structure, and properties. This work demonstrates the utility of these method as a tool for determining composite properties from constituent and processing parameters. Initial particle dynamics simulation incorporating electromagnetic responses between particles and between the particles and the applied fields, including dielectrophoresis, are used to stochastically generate representative volume elements for a given set of process variables. Next, these RVEs are analyzed as periodic structures using FEA yielding bulk material properties. The results are shown to converge for simulation size and discretization, validating the RVE as an appropriate representation of the composite volume. Calculated material properties are compared to traditional effective medium theory models. Simulations allow for mapping of composite properties with respect to not only composition, but also fundamentally from processing simulations that yield varying particle configurations, a step not present in traditional or more modern effective medium theories such as the Halpin Tsai or double-inclusion theories.

     
    more » « less
  5. 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