The importance of modeling the human cerebral vasculature in blunt trauma

The Significance of Cross-Sectional Shape Accuracy and Non-Linear Elasticity on the Numerical Modelling of Cerebral Veins under Tensile Loading

Traumatic brain injury (TBI) is a serious global health issue, leading to serious disabilities. One type of TBI is acute subdural haematoma (ASDH), which occurs when a bridging vein ruptures. Many numerical models of these structures, mainly based on the finite element method, have been developed. However, most rely on linear elasticity (without validation) and others on simplifications at the geometrical level. An example of the latter is the assumption of a regular cylinder with a constant radius, or the geometry of the vein acquired from medical images. Unfortunately, these do not replicate the real conditions of a mechanical tensile test. In this work, the main goal is to evaluate the influence of the vein's geometry in its mechanical behaviour under tensile loading, simulating the real conditions of experimental tests. The second goal is to implement a hyperelastic model of the bridging veins where it would be possible to observe its non-linear elastic behaviour. The results of the developed finite element models were compared to experimental data available in the literature and other models. It was possible to conclude that the geometry of the vein structure influences the tensile stress-strain curve, which means that flattened specimens should be modelled when validating constitutive models for bridging veins. Additionally, the implementation of hyperelastic material models has been verified, highlighting the potential application of the Marlow and reduced polynomial (of fourth and sixth orders) constitutive models.

An atlas of the heterogeneous viscoelastic brain with local power-law attenuation synthesised using Prony-series

This review addresses the acute need to acknowledge the mechanical heterogeneity of brain matter and to accurately calibrate its local viscoelastic material properties accordingly. Specifically, it is important to compile the existing and disparate literature on attenuation power-laws and dispersion to make progress in wave physics of brain matter, a field of research that has the potential to explain the mechanisms at play in diffuse axonal injury and mild traumatic brain injury in general. Currently, viscous effects in the brain are modelled using Prony-series, i.e., a sum of decaying exponentials at different relaxation times. Here we collect and synthesise the Prony-series coefficients appearing in the literature for twelve regions: brainstem, basal ganglia, cerebellum, corona radiata, corpus callosum, cortex, dentate gyrus, hippocampus, thalamus, grey matter, white matter, homogeneous brain, and for eight different mammals: pig, rat, human, mouse, cow, sheep, monkey and dog. Using this data, we compute the fractional-exponent attenuation power-laws for different tissues of the brain, the corresponding dispersion laws resulting from causality, and the averaged Prony-series coefficients. STATEMENT OF SIGNIFICANCE: Traumatic brain injuries are considered a silent epidemic and finite element methods (FEMs) are used in modelling brain deformation, requiring access to viscoelastic properties of brain. To the best of our knowledge, this work presents 1) the first multi-frequency viscoelastic atlas of the heterogeneous brain, 2) the first review focusing on viscoelastic modelling in both FEMs and experimental works, 3) the first attempt to conglomerate the disparate existing literature on the viscoelastic modelling of the brain and 4) the largest collection of viscoelastic parameters for the brain (212 different Prony-series spanning 12 different tissues and 8 different animal surrogates). Furthermore, this work presents the first brain atlas of attenuation power-laws essential for modelling shear waves in brain.

Age effects on the mechanical behavior of human cerebral bridging veins

BACKGROUND

It is well established that the probability of occurrence of acute subdural hematomas in traumatic situations increases with age, since the main cause of such hematomas is the mechanical failure of cerebral blood vessels known as bridging veins. This research aims to determine whether there is an effect of age on the mechanical properties of these cerebral vessels, because previous reported studies were conflicting.

METHODS

This study used mechanical tests blue of cerebral bridging veins from post-mortem human subjects. In particular, a series of in vitro tensile tests were performed on a balanced sample of bridging veins from different human subjects.

FINDINGS

The mechanical parameters measured from the tests were analyzed by means of regression analysis looking for age related effects. The results show that there is a significant effect on both the ultimate strength, maximum stress and strain that the specimens can withstand. The quantitative analysis shows reductions of nearly 50% in ultimate stress, and almost 35% in ultimate strain.

INTERPRETATION

Mechanical deterioration of the mechanical strength of cerebral blood vessels seems to be a major factor involved in the increase of frequency of acute subdural hematoma in elderly people in a wide range of life-threatening traumatic situations.

A biomechanical-based approach to scale blast-induced molecular changes in the brain

Animal studies provide valuable insights on how the interaction of blast waves with the head may injure the brain. However, there is no acceptable methodology to scale the findings from animals to humans. Here, we propose an experimental/computational approach to project observed blast-induced molecular changes in the rat brain to the human brain. Using a shock tube, we exposed rats to a range of blast overpressures (BOPs) and used a high-fidelity computational model of a rat head to correlate predicted biomechanical responses with measured changes in glial fibrillary acidic protein (GFAP) in rat brain tissues. Our analyses revealed correlates between model-predicted strain rate and measured GFAP changes in three brain regions. Using these correlates and a high-fidelity computational model of a human head, we determined the equivalent BOPs in rats and in humans that induced similar strain rates across the two species. We used the equivalent BOPs to project the measured GFAP changes in the rat brain to the human. Our results suggest that, relative to the rat, the human requires an exposure to a blast wave of a higher magnitude to elicit similar brain-tissue responses. Our proposed methodology could assist in the development of safety guidelines for blast exposure.

A Finite Element Model of Cerebral Vascular Injury for Predicting Microbleeds Location

Finite Element (FE) models of brain mechanics have improved our understanding of the brain response to rapid mechanical loads that produce traumatic brain injuries. However, these models have rarely incorporated vasculature, which limits their ability to predict the response of vessels to head impacts. To address this shortcoming, here we used high-resolution MRI scans to map the venous system anatomy at a submillimetre resolution. We then used this map to develop an FE model of veins and incorporated it in an anatomically detailed FE model of the brain. The model prediction of brain displacement at different locations was compared to controlled experiments on post-mortem human subject heads, yielding over 3,100 displacement curve comparisons, which showed fair to excellent correlation between them. We then used the model to predict the distribution of axial strains and strain rates in the veins of a rugby player who had small blood deposits in his white matter, known as microbleeds, after sustaining a head collision. We hypothesised that the distribution of axial strain and strain rate in veins can predict the pattern of microbleeds. We reconstructed the head collision using video footage and multi-body dynamics modelling and used the predicted head accelerations to load the FE model of vascular injury. The model predicted large axial strains in veins where microbleeds were detected. A region of interest analysis using white matter tracts showed that the tract group with microbleeds had 95th percentile peak axial strain and strain rate of 0.197 and 64.9 s⁻¹ respectively, which were significantly larger than those of the group of tracts without microbleeds (0.163 and 57.0 s⁻¹). This study does not derive a threshold for the onset of microbleeds as it investigated a single case, but it provides evidence for a link between strain and strain rate applied to veins during head impacts and structural damage and allows for future work to generate threshold values. Moreover, our results suggest that the FE model has the potential to be used to predict intracranial vascular injuries after TBI, providing a more objective tool for TBI assessment and improving protection against it.

A 3-D Finite-Element Minipig Model to Assess Brain Biomechanical Responses to Blast Exposure

Despite years of research, it is still unknown whether the interaction of explosion-induced blast waves with the head causes injury to the human brain. One way to fill this gap is to use animal models to establish “scaling laws” that project observed brain injuries in animals to humans. This requires laboratory experiments and high-fidelity mathematical models of the animal head to establish correlates between experimentally observed blast-induced brain injuries and model-predicted biomechanical responses. To this end, we performed laboratory experiments on Göttingen minipigs to develop and validate a three-dimensional (3-D) high-fidelity finite-element (FE) model of the minipig head. First, we performed laboratory experiments on Göttingen minipigs to obtain the geometry of the cerebral vasculature network and to characterize brain-tissue and vasculature material properties in response to high strain rates typical of blast exposures. Next, we used the detailed cerebral vasculature information and species-specific brain tissue and vasculature material properties to develop the 3-D high-fidelity FE model of the minipig head. Then, to validate the model predictions, we performed laboratory shock-tube experiments, where we exposed Göttingen minipigs to a blast overpressure of 210 kPa in a laboratory shock tube and compared brain pressures at two locations. We observed a good agreement between the model-predicted pressures and the experimental measurements, with differences in maximum pressure of less than 6%. Finally, to evaluate the influence of the cerebral vascular network on the biomechanical predictions, we performed simulations where we compared results of FE models with and without the vasculature. As expected, incorporation of the vasculature decreased brain strain but did not affect the predictions of brain pressure. However, we observed that inclusion of the cerebral vasculature in the model changed the strain distribution by as much as 100% in regions near the interface between the vasculature and the brain tissue, suggesting that the vasculature does not merely decrease the strain but causes drastic redistributions. This work will help establish correlates between observed brain injuries and predicted biomechanical responses in minipigs and facilitate the creation of scaling laws to infer potential injuries in the human brain due to exposure to blast waves.

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.

Mechanical Behavior of Blood Vessels: Elastic and Viscoelastic Contributions

The mechanical properties of the cerebral bridging veins (CBVs) were studied using advanced microtensile equipment. Detailed high-quality curves were obtained at different strain rates, showing a clearly nonlinear stress-strain response. In addition, the tissue of the CBVs exhibits stress relaxation and a preconditioning effect under cyclic loading, unequivocal indications of viscoelastic behavior. Interestingly, most previous literature that conducts uniaxial tensile tests had not found significant viscoelastic effects in CBVs, but the use of more sensitive tests allowed to observe the viscoelastic effects. For that reason, a careful mathematical analysis is presented, clarifying why in uniaxial tests with moderate strain rates, it is difficult to observe any viscoelastic effect. The analysis provides a theoretical explanation as to why many recent studies that investigated mechanical properties did not find a significant viscoelastic effect, even though in other circumstances, the CBV tissue would clearly exhibit viscoelastic behavior. Finally, this study provides reference values for the usual mechanical properties, as well as calculations of constitutive parameters for nonlinear elastic and viscoelastic models that would allow more accurate numerical simulation of CBVs in Finite Element-based computational models in future works.