Embedded Finite Elements for Modeling Axonal Injury

Dynamic strain fields of the mouse brain during rotation

Mouse models are used to better understand brain injury mechanisms in humans, yet there is a limited understanding of biomechanical relevance, beginning with how the murine brain deforms when the head undergoes rapid rotation from blunt impact. This problem makes it difficult to translate some aspects of diffuse axonal injury from mouse to human. To address this gap, we present the two-dimensional strain field of the mouse brain undergoing dynamic rotation in the sagittal plane. Using a high-speed camera with digital image correlation measurements of the exposed mid-sagittal brain surface, we found that pure rotations (no direct impact to the skull) of 100-200 rad/s are capable of producing complex strain fields that evolve over time with respect to rotational acceleration and deceleration. At the highest rotational velocity tested, the largest tensile strains (≥ 21% elongation) in selected regions of the mouse brain approach strain thresholds previously associated with axonal injury in prior work. These findings provide a benchmark to validate the mechanical response in biomechanical computational models predicting diffuse axonal injury, but much work remains in correlating tissue deformation patterns from computational models with underlying neuropathology.

Histology-informed multiscale modeling of human brain white matter

In this study, we propose a novel micromechanical model for the brain white matter, which is described as a heterogeneous material with a complex network of axon fibers embedded in a soft ground matrix. We developed this model in the framework of RVE-based multiscale theories in combination with the finite element method and the embedded element technique for embedding the fibers. Microstructural features such as axon diameter, orientation and tortuosity are incorporated into the model through distributions derived from histological data. The constitutive law of both the fibers and the matrix is described by isotropic one-term Ogden functions. The hyperelastic response of the tissue is derived by homogenizing the microscopic stress fields with multiscale boundary conditions to ensure kinematic compatibility. The macroscale homogenized stress is employed in an inverse parameter identification procedure to determine the hyperelastic constants of axons and ground matrix, based on experiments on human corpus callosum. Our results demonstrate the fundamental effect of axon tortuosity on the mechanical behavior of the brain's white matter. By combining histological information with the multiscale theory, the proposed framework can substantially contribute to the understanding of mechanotransduction phenomena, shed light on the biomechanics of a healthy brain, and potentially provide insights into neurodegenerative processes.

Approximating subject-specific brain injury models via scaling based on head-brain morphological relationships

Most human head/brain models represent a generic adult male head/brain. They may suffer in accuracy when investigating traumatic brain injury (TBI) on a subject-specific basis. Subject-specific models can be developed from neuroimages; however, neuroimages are not typically available in practice. In this study, we establish simple and elegant regression models between brain outer surface morphology and head dimensions measured from neuroimages along with age and sex information (N = 191; 141 males and 50 females with age ranging 14-25 years). The regression models are then used to approximate subject-specific brain models by scaling a generic counterpart, without using neuroimages. Model geometrical accuracy is assessed using adjusted [Formula: see text] and absolute percentage error (e.g., 0.720 and 3.09 ± 2.38%, respectively, for brain volume when incorporating tragion-to-top). For a subset of 11 subjects (from smallest to largest in brain volume), impact-induced brain strains are compared with those from "morphed models" derived from neuroimage-based mesh warping. We find that regional peak strains from the scaled subject-specific models are comparable to those of the morphed counterparts but could be considerably different from those of the generic model (e.g., linear regression slope of 1.01-1.03 for gray and white matter regions versus 1.16-1.19, or up to ~ 20% overestimation for the smallest brain studied). These results highlight the importance of incorporating brain morphological variations in impact simulation and demonstrate the feasibility of approximating subject-specific brain models without neuroimages using age, sex, and easily measurable head dimensions. The scaled models may improve subject specificity for future TBI investigations.

Cerebral Vasculature Influences Blast-Induced Biomechanical Responses of Human Brain Tissue

Multiple finite-element (FE) models to predict the biomechanical responses in the human brain resulting from the interaction with blast waves have established the importance of including the brain-surface convolutions, the major cerebral veins, and using non-linear brain-tissue properties to improve model accuracy. We hypothesize that inclusion of a more detailed network of cerebral veins and arteries can further enhance the model-predicted biomechanical responses and help identify correlates of blast-induced brain injury. To more comprehensively capture the biomechanical responses of human brain tissues to blast-wave exposure, we coupled a three-dimensional (3-D) detailed-vasculature human-head FE model, previously validated for blunt impact, with a 3-D shock-tube FE model. Using the coupled model, we computed the biomechanical responses of a human head facing an incoming blast wave for blast overpressures (BOPs) equivalent to 68, 83, and 104 kPa. We validated our FE model, which includes the detailed network of cerebral veins and arteries, the gyri and the sulci, and hyper-viscoelastic brain-tissue properties, by comparing the model-predicted intracranial pressure (ICP) values with previously collected data from shock-tube experiments performed on cadaver heads. In addition, to quantify the influence of including a more comprehensive network of brain vessels, we compared the biomechanical responses of our detailed-vasculature model with those of a reduced-vasculature model and a no-vasculature model for the same blast-loading conditions. For the three BOPs, the predicted ICP values matched well with the experimental results in the frontal lobe, with peak-pressure differences of 4-11% and phase-shift differences of 9-13%. As expected, incorporating the detailed cerebral vasculature did not influence the ICP, however, it redistributed the peak brain-tissue strains by as much as 30% and yielded peak strain differences of up to 7%. When compared to existing reduced-vasculature FE models that only include the major cerebral veins, our high-fidelity model redistributed the brain-tissue strains in most of the brain, highlighting the importance of including a detailed cerebral vessel network in human-head FE models to more comprehensively account for the biomechanical responses induced by blast exposure.

Toward a Comprehensive Delineation of White Matter Tract-Related Deformation

Finite element (FE) models of the human head are valuable instruments to explore the mechanobiological pathway from external loading, localized brain response, and resultant injury risks. The injury predictability of these models depends on the use of effective criteria as injury predictors. The FE-derived normal deformation along white matter (WM) fiber tracts (i.e., tract-oriented strain) recently has been suggested as an appropriate predictor for axonal injury. However, the tract-oriented strain only represents a partial depiction of the WM fiber tract deformation. A comprehensive delineation of tract-related deformation may improve the injury predictability of the FE head model by delivering new tract-related criteria as injury predictors. Thus, the present study performed a theoretical strain analysis to comprehensively characterize the WM fiber tract deformation by relating the strain tensor of the WM element to its embedded fiber tract. Three new tract-related strains with exact analytical solutions were proposed, measuring the normal deformation perpendicular to the fiber tracts (i.e., tract-perpendicular strain), and shear deformation along and perpendicular to the fiber tracts (i.e., axial-shear strain and lateral-shear strain, respectively). The injury predictability of these three newly proposed strain peaks along with the previously used tract-oriented strain peak and maximum principal strain (MPS) were evaluated by simulating 151 impacts with known outcome (concussion or non-concussion). The results preliminarily showed that four tract-related strain peaks exhibited superior performance than MPS in discriminating concussion and non-concussion cases. This study presents a comprehensive quantification of WM tract-related deformation and advocates the use of orientation-dependent strains as criteria for injury prediction, which may ultimately contribute to an advanced mechanobiological understanding and enhanced computational predictability of brain injury.

Investigation of the Compressive Viscoelastic Properties of Brain Tissue Under Time and Frequency Dependent Loading Conditions

The mechanical characterization of brain tissue has been generally analyzed in the frequency and time domain. It is crucial to understand the mechanics of the brain under realistic, dynamic conditions and convert it to enable mathematical modelling in a time domain. In this study, the compressive viscoelastic properties of brain tissue were investigated under time and frequency domains with the same physical conditions and the theory of viscoelasticity was applied to estimate the prediction of viscoelastic response in the time domain based on frequency-dependent mechanical moduli through Finite Element models. Storage and loss modulus were obtained from white and grey matter, of bovine brains, using dynamic mechanical analysis and time domain material functions were derived based on a Prony series representation. The material models were evaluated using brain testing data from stress relaxation and hysteresis in the time dependent analysis. The Finite Element models were able to represent the trend of viscoelastic characterization of brain tissue under both testing domains. The outcomes of this study contribute to a better understanding of brain tissue mechanical behaviour and demonstrate the feasibility of deriving time-domain viscoelastic parameters from frequency-dependent compressive data for biological tissue, as validated by comparing experimental tests with computational simulations.

Multiscale modelling of cerebrovascular injury reveals the role of vascular anatomy and parenchymal shear stresses

Neurovascular injury is often observed in traumatic brain injury (TBI). However, the relationship between mechanical forces and vascular injury is still unclear. A key question is whether the complex anatomy of vasculature plays a role in increasing forces in cerebral vessels and producing damage. We developed a high-fidelity multiscale finite element model of the rat brain featuring a detailed definition of the angioarchitecture. Controlled cortical impacts were performed experimentally and in-silico. The model was able to predict the pattern of blood-brain barrier damage. We found strong correlation between the area of fibrinogen extravasation and the brain area where axial strain in vessels exceeds 0.14. Our results showed that adjacent vessels can sustain profoundly different axial stresses depending on their alignment with the principal direction of stress in parenchyma, with a better alignment leading to larger stresses in vessels. We also found a strong correlation between axial stress in vessels and the shearing component of the stress wave in parenchyma. Our multiscale computational approach explains the unrecognised role of the vascular anatomy and shear stresses in producing distinct distribution of large forces in vasculature. This new understanding can contribute to improving TBI diagnosis and prevention.

Effect of axonal fiber architecture on mechanical heterogeneity of the white matter-a statistical micromechanical model

A diffusion tensor imaging (DTI) -based statistical micromechanical model was developed to study the effect of axonal fiber architecture on the inter- and intra-regional mechanical heterogeneity of the white matter. Three characteristic regions within the white matter, i.e., corpus callosum, brain stem, and corona radiata, were studied considering the previous observations of locations of diffuse axonal injury. The embedded element technique was used to create a fiber-reinforced model, where the fiber was characterized by a Holzapfel hyperelastic material model with variable dispersion of axonal orientations. A relationship between the fractional anisotropy and the dispersion parameter of the hyperelastic model was used to introduce the statistical DTI data into the representative volume element. The FA-informed statistical micromechanical models of three characteristic regions of white matter were developed by deriving the corresponding probabilistic measures of FA variations. Comparison of the model predictions and experimental data indicated a good agreement, suggesting that the model could reasonably capture the inter-regional heterogeneity of white matter. Moreover, the standard deviations of experimental results correlated well with the model predictions, suggesting that the model could capture the intra-regional mechanical heterogeneity for different regions of white matter.

Calibration of a Heterogeneous Brain Model Using a Subject-Specific Inverse Finite Element Approach

Central to the investigation of the biomechanics of traumatic brain injury (TBI) and the assessment of injury risk from head impact are finite element (FE) models of the human brain. However, many existing FE human brain models have been developed with simplified representations of the parenchyma, which may limit their applicability as an injury prediction tool. Recent advances in neuroimaging techniques and brain biomechanics provide new and necessary experimental data that can improve the biofidelity of FE brain models. In this study, the CAB-20MSym template model was developed, calibrated, and extensively verified. To implement material heterogeneity, a magnetic resonance elastography (MRE) template image was leveraged to define the relative stiffness gradient of the brain model. A multi-stage inverse FE (iFE) approach was used to calibrate the material parameters that defined the underlying non-linear deviatoric response by minimizing the error between model-predicted brain displacements and experimental displacement data. This process involved calibrating the infinitesimal shear modulus of the material using low-severity, low-deformation impact cases and the material non-linearity using high-severity, high-deformation cases from a dataset of in situ brain displacements obtained from cadaveric specimens. To minimize the geometric discrepancy between the FE models used in the iFE calibration and the cadaveric specimens from which the experimental data were obtained, subject-specific models of these cadaveric brain specimens were developed and used in the calibration process. Finally, the calibrated material parameters were extensively verified using independent brain displacement data from 33 rotational head impacts, spanning multiple loading directions (sagittal, coronal, axial), magnitudes (20–40 rad/s), durations (30–60 ms), and severity. Overall, the heterogeneous CAB-20MSym template model demonstrated good biofidelity with a mean overall CORA score of 0.63 ± 0.06 when compared to in situ brain displacement data. Strains predicted by the calibrated model under non-injurious rotational impacts in human volunteers (N = 6) also demonstrated similar biofidelity compared to in vivo measurements obtained from tagged magnetic resonance imaging studies. In addition to serving as an anatomically accurate model for further investigations of TBI biomechanics, the MRE-based framework for implementing material heterogeneity could serve as a foundation for incorporating subject-specific material properties in future models.

Instantaneous Whole-Brain Strain Estimation in Dynamic Head Impact

Head injury models are notoriously time consuming and resource demanding in simulations, which prevents routine application. Here, we extend a convolutional neural network (CNN) to instantly estimate element-wise distribution of peak maximum principal strain (MPS) of the entire brain (>36 k speedup accomplished on a low-end computing platform). To achieve this, head impact rotational velocity and acceleration temporal profiles are combined into two-dimensional images to serve as CNN input for training and prediction of MPS. Compared with the directly simulated counterparts, the CNN-estimated responses (magnitude and distribution) are sufficiently accurate for 92.1% of the cases via 10-fold cross-validation using impacts drawn from the real world (n = 5661; range of peak rotational velocity in augmented data extended to 2-40 rad/sec). The success rate further improves to 97.1% for "in-range" impacts (n = 4298). When using the same CNN architecture to train (n = 3064) and test on an independent, reconstructed National Football League (NFL) impact dataset (n = 53; 20 concussions and 33 non-injuries), 51 out of 53, or 96.2% of the cases, are sufficiently accurate. The estimated responses also achieve virtually identical concussion prediction performances relative to the directly simulated counterparts, and they often outperform peak MPS of the whole brain (e.g., accuracy of 0.83 vs. 0.77 via leave-one-out cross-validation). These findings support the use of CNN for accurate and efficient estimation of spatially detailed brain strains across the vast majority of head impacts in contact sports. Our technique may hold the potential to transform traumatic brain injury (TBI) research and the design and testing standards of head protective gears by facilitating the transition from acceleration-based approximation to strain-based design and analysis. This would have broad implications in the TBI biomechanics field to accelerate new scientific discoveries. The pre-trained CNN is freely available online at https://github.com/Jilab-biomechanics/CNN-brain-strains.

Rapid Estimation of Entire Brain Strain Using Deep Learning Models

OBJECTIVE

Many recent studies have suggested that brain deformation resulting from a head impact is linked to the corresponding clinical outcome, such as mild traumatic brain injury (mTBI). Even though several finite element (FE) head models have been developed and validated to calculate brain deformation based on impact kinematics, the clinical application of these FE head models is limited due to the time-consuming nature of FE simulations. This work aims to accelerate the process of brain deformation calculation and thus improve the potential for clinical applications.

METHODS

We propose a deep learning head model with a five-layer deep neural network and feature engineering, and trained and tested the model on 2511 total head impacts from a combination of head model simulations and on-field college football and mixed martial arts impacts.

RESULTS

The proposed deep learning head model can calculate the maximum principal strain (Green Lagrange) for every element in the entire brain in less than 0.001s with an average root mean squared error of 0.022, and with a standard deviation of 0.001 over twenty repeats with random data partition and model initialization.

CONCLUSION

Trained and tested using the dataset of 2511 head impacts, this model can be applied to various sports in the calculation of brain strain with accuracy, and its applicability can even further be extended by incorporating data from other types of head impacts.

SIGNIFICANCE

In addition to the potential clinical application in real-time brain deformation monitoring, this model will help researchers estimate the brain strain from a large number of head impacts more efficiently than using FE models.

White Matter Tract-Oriented Deformation Is Dependent on Real-Time Axonal Fiber Orientation

Traumatic axonal injury (TAI) is a critical public health issue with its pathogenesis remaining largely elusive. Finite element (FE) head models are promising tools to bridge the gap between mechanical insult, localized brain response, and resultant injury. In particular, the FE-derived deformation along the direction of white matter (WM) tracts (i.e., tract-oriented strain) has been shown to be an appropriate predictor for TAI. The evolution of fiber orientation in time during the impact and its potential influence on the tract-oriented strain remains unknown, however. To address this question, the present study leveraged an embedded element approach to track real-time fiber orientation during impacts. A new scheme to calculate the tract-oriented strain was proposed by projecting the strain tensors from pre-computed simulations along the temporal fiber direction instead of its static counterpart directly obtained from diffuse tensor imaging. The results revealed that incorporating the real-time fiber orientation not only altered the direction but also amplified the magnitude of the tract-oriented strain, resulting in a generally more extended distribution and a larger volume ratio of WM exposed to high deformation along fiber tracts. These effects were exacerbated with the impact severities characterized by the acceleration magnitudes. Results of this study provide insights into how best to incorporate fiber orientation in head injury models and derive the WM tract-oriented deformation from computational simulations, which is important for furthering our understanding of the underlying mechanisms of TAI.

The importance of modeling the human cerebral vasculature in blunt trauma

BACKGROUND

Multiple studies describing human head finite element (FE) models have established the importance of including the major cerebral vasculature to improve the accuracy of the model predictions. However, a more detailed network of cerebral vasculature, including the major veins and arteries as well as their branch vessels, can further enhance the model-predicted biomechanical responses and help identify correlates to observed blunt-induced brain injury.

METHODS

We used an anatomically accurate three-dimensional geometry of a 50th percentile U.S. male head that included the skin, eyes, sinuses, spine, skull, brain, meninges, and a detailed network of cerebral vasculature to develop a high-fidelity model. We performed blunt trauma simulations and determined the intracranial pressure (ICP), the relative displacement (RD), the von Mises stress, and the maximum principal strain. We validated our detailed-vasculature model by comparing the model-predicted ICP and RD values with experimental measurements. To quantify the influence of including a more comprehensive network of brain vessels, we compared the biomechanical responses of our detailed-vasculature model with those of a reduced-vasculature model and a no-vasculature model.

RESULTS

For an inclined frontal impact, the predicted ICP matched well with the experimental results in the fossa, frontal, parietal, and occipital lobes, with peak-pressure differences ranging from 2.4% to 9.4%. For a normal frontal impact, the predicted ICP matched the experimental results in the frontal lobe and lateral ventricle, with peak-pressure discrepancies equivalent to 1.9% and 22.3%, respectively. For an offset parietal impact, the model-predicted RD matched well with the experimental measurements, with peak RD differences of 27% and 24% in the right and left cerebral hemispheres, respectively. Incorporating the detailed cerebral vasculature did not influence the ICP but redistributed the brain-tissue stresses and strains by as much as 30%. In addition, our detailed-vasculature model predicted strain reductions by as much as 28% when compared to current reduced-vasculature FE models that only include the major cerebral vessels.

CONCLUSIONS

Our study highlights the importance of including a detailed representation of the cerebral vasculature in FE models to more accurately estimate the biomechanical responses of the human brain to blunt impact.

Multi-Scale White Matter Tract Embedded Brain Finite Element Model Predicts the Location of Traumatic Diffuse Axonal Injury

Finite element models (FEMs) are used increasingly in the traumatic brain injury (TBI) field to provide an estimation of tissue responses and predict the probability of sustaining TBI after a biomechanical event. However, FEM studies have mainly focused on predicting the absence/presence of TBI rather than estimating the location of injury. In this study, we created a multi-scale FEM of the pig brain with embedded axonal tracts to estimate the sites of acute (≤6 h) traumatic axonal injury (TAI) after rapid head rotation. We examined three finite element (FE)-derived metrics related to the axonal bundle deformation and three FE-derived metrics based on brain tissue deformation for prediction of acute TAI location. Rapid head rotations were performed in pigs, and TAI neuropathological maps were generated and colocalized to the FEM. The distributions of the FEM-derived brain/axonal deformations spatially correlate with the locations of acute TAI. For each of the six metric candidates, we examined a matrix of different injury thresholds (th_x) and distance to actual TAI sites (d_s) to maximize the average of two optimization criteria. Three axonal deformation-related TAI candidates predicted the sites of acute TAI within 2.5 mm, but no brain tissue metric succeeded. The optimal range of thresholds for maximum axonal strain, maximum axonal strain rate, and maximum product of axonal strain and strain rate were 0.08-0.14, 40-90, and 2.0-7.5 s^-1, respectively. The upper bounds of these thresholds resulted in higher true-positive prediction rate. In summary, this study confirmed the hypothesis that the large axonal-bundle deformations occur on/close to the areas that sustained TAI.

Embedded axonal fiber tracts improve finite element model predictions of traumatic brain injury

With the growing rate of traumatic brain injury (TBI), there is an increasing interest in validated tools to predict and prevent brain injuries. Finite element models (FEM) are valuable tools to estimate tissue responses, predict probability of TBI, and guide the development of safety equipment. In this study, we developed and validated an anisotropic pig brain multi-scale FEM by explicitly embedding the axonal tract structures and utilized the model to simulate experimental TBI in piglets undergoing dynamic head rotations. Binary logistic regression, survival analysis with Weibull distribution, and receiver operating characteristic curve analysis, coupled with repeated k-fold cross-validation technique, were used to examine 12 FEM-derived metrics related to axonal/brain tissue strain and strain rate for predicting the presence or absence of traumatic axonal injury (TAI). All 12 metrics performed well in predicting of TAI with prediction accuracy rate of 73-90%. The axonal-based metrics outperformed their rival brain tissue-based metrics in predicting TAI. The best predictors of TAI were maximum axonal strain times strain rate (MASxSR) and its corresponding optimal fraction-based metric (AF-MASxSR7.5) that represents the fraction of axonal fibers exceeding MASxSR of 7.5 s-1. The thresholds compare favorably with tissue tolerances found in in-vitro/in-vivo measurements in the literature. In addition, the damaged volume fractions (DVF) predicted using the axonal-based metrics, especially MASxSR (DVF = 0.05-4.5%), were closer to the actual DVF obtained from histopathology (AIV = 0.02-1.65%) in comparison with the DVF predicted using the brain-related metrics (DVF = 0.11-41.2%). The methods and the results from this study can be used to improve model prediction of TBI in humans.

Biomechanics of Blast TBI With Time-Resolved Consecutive Primary, Secondary, and Tertiary Loads

Blast-induced traumatic brain injury (bTBI) has become a signature casualty of recent military operations. In spite of significant clinical and preclinical TBI research, current understanding of injury mechanisms and short- and long-term outcomes is limited. Mathematical models of bTBI biomechanics may help in better understanding of injury mechanisms and in the development of improved neuroprotective strategies. Until present, bTBI has been analyzed as a single event of a blast pressure wave propagating through the brain. In many bTBI events, the loads on the body and the head are spatially and temporarily distributed, involving the primary intracranial pressure wave, followed by the head rotation and then by head impact on the ground. In such cases, the brain microstructures may experience time/space distributed (consecutive) damage and recovery events. The paper presents a novel multiscale simulation framework that couples the body/brain scale biomechanics with micro-scale mechanobiology to study the effects of micro-damage to neuro-axonal structures. Our results show that the micro-mechanical responses of neuro-axonal structures occur sequentially in time with “damage” and “relaxation” periods in different parts of the brain. A new integrated computational framework is described coupling the brain-scale biomechanics with micro-mechanical damage to axonal and synaptic structures.

Material and Structural Modeling Aspects of Brain Tissue Deformation under Dynamic Loads

The aim of this work was to assess the numerous approaches to structural and material modeling of brain tissue under dynamic loading conditions. The current technological improvements in material modeling have led to various approaches described in the literature. However, the methods used for the determination of the brain’s characteristics have not always been stated or clearly defined and material data are even more scattered. Thus, the research described in this paper explicitly underlines directions for the development of numerical brain models. An important element of this research is the development of a numerical model of the brain based on medical imaging methods. This approach allowed the authors to assess the changes in the mechanical and geometrical parameters of brain tissue caused by the impact of mechanical loads. The developed model was verified through comparison with experimental studies on post-mortem human subjects described in the literature, as well as through numerical tests. Based on the current research, the authors identified important aspects of the modeling of brain tissue that influence the assessment of the actual biomechanical response of the brain for dynamic analyses.