Incorporation of vasculature in a head injury model lowers local mechanical strains in dynamic impact

Embedded finite element modeling of the mechanics of brain axonal fiber tracts under head impact conditions

Investigating and understanding the biomechanical kinematics and kinetics of human brain axonal fibers during head impact process is crucial to study the mechanisms of Traumatic Axonal Injury (TAI). Such a study may require the explicit incorporation of brain fiber tracts into the host brain in order to distinguish the mechanical states of axonal fibers and brain tissue. Herein we extend our previously developed human head model by using an embedded element method to include fiber tracts reconstructed from diffusion tensor images in a host brain with the purpose of numerically tracking the deformation state of axonal fiber tracts during a head impact simulation. The updated model is validated by comparing its prediction of intracranial pressures with experimental data, followed by a thorough study of the effects of element types used for fiber tracts and the stiffness ratios of fiber to host brain. The validated model is also used to predict and visualize the damaged region of fiber tracts during the head impact process based on different injury criteria. The model is promising in tracking the state of fiber tracts and can add more objective functions such as axonal fiber deformation if used in the future design optimization of head protective equipment such as a football helmet.

Anatomical Features and Material Properties of Human Surrogate Head Models Affect Spatial and Temporal Brain Motion under Blunt Impact

Traumatic brain injury (TBI) is a biomechanical problem where the initiating event is dynamic loading (blunt, inertial, blast) to the head. To understand the relationship between the mechanical parameters of the injury and the deformation patterns in the brain, we have previously developed a surrogate head (SH) model capable of measuring spatial and temporal deformation in a surrogate brain under blunt impact. The objective of this work was to examine how material properties and anatomical features affect the motion of the brain and the development of injurious deformations. The SH head model was modified to study six variables independently under blunt impact: surrogate brain stiffness, surrogate skull stiffness, inclusion of cerebrospinal fluid (CSF), head/skull size, inclusion of vasculature, and neck stiffness. Each experimental SH was either crown or frontally impacted at 1.3 m/s (3 mph) using a drop tower system. Surrogate brain material, the Hybrid III neck stiffness, and skull stiffness were measured and compared to published properties. Results show that the most significant variables affecting changes in brain deformation are skull stiffness, inclusion of CSF and surrogate brain stiffness. Interestingly, neck stiffness and SH size significantly affected the strain rate only suggesting these parameters are less important in blunt trauma. While the inclusion of vasculature locally created strain concentrations at the interface of the artery and brain, overall deformation was reduced.

Vascular microforamina and endocranial surface: Normal variation and distribution in adult humans: Vascular biology

The term craniovascular traits refers to the imprints left by arteries and veins on the skull bones. These features can be used in biological anthropology and archaeology to investigate the morphology of the vascular network in extinct species and past populations. Generally, the term refers to macrovascular features of the endocranial cavity, like those associated with the middle meningeal artery, venous sinuses, emissary foramina, and diploic channels. However, small vascular passages (here called microforamina or microchannels) have been occasionally described on the endocranial surface. The larger ones (generally with a diameter between 0.5 and 2.0 mm) can be detected through medical scanners on osteological collections. In this study, we describe and quantify the number and distribution of these microforamina in adult humans (N = 45) and, preliminarily, in a small sample of children (N = 7). Adults display more microchannels than juvenile skulls. A higher frequency in females is also observed, although this result is not statistically significant and might be associated with allometric cranial variations. The distribution of the microforamina is particularly concentrated on the top of the vault, in particular along the sagittal, metopic, and coronal sutures, matching the course of major venous sinuses and parasagittal bridging veins. Nonetheless, the density is lower in the region posterior to bregma. Beyond oxygenation, these vessels are likely involved in endocranial thermal regulation, infection, inflammation, and immune responses, and their distribution and prevalence can hence be of interest in human biology, evolutionary anthropology, and medicine.

A computational pipeline towards large-scale and multiscale modeling of traumatic axonal injury

Contemporary biomechanical modeling of traumatic brain injury (TBI) focuses on either the global brain as an organ or a representative tiny section of a single axon. In addition, while it is common for a global brain model to employ real-world impacts as input, axonal injury models have largely been limited to inputs of either tension or compression with assumed peak strain and strain rate. These major gaps between global and microscale modeling preclude a systematic and mechanistic investigation of how tissue strain from impact leads to downstream axonal damage throughout the white matter. In this study, a unique subject-specific multimodality dataset from a male ice-hockey player sustaining a diagnosed concussion is used to establish an efficient and scalable computational pipeline. It is then employed to derive voxelized brain deformation, maximum principal strains and white matter fiber strains, and finally, to produce diverse fiber strain profiles of various shapes in temporal history necessary for the development and application of a deep learning axonal injury model in the future. The pipeline employs a structured, voxelized representation of brain deformation with adjustable spatial resolution independent of model mesh resolution. The method can be easily extended to other head impacts or individuals. The framework established in this work is critical for enabling large-scale (i.e., across the entire white matter region, head impacts, and individuals) and multiscale (i.e., from organ to cell length scales) modeling for the investigation of traumatic axonal injury (TAI) triggering mechanisms. Ultimately, these efforts could enhance the assessment of concussion risks and design of protective headgear. Therefore, this work contributes to improved strategies for concussion detection, mitigation, and prevention.

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.

Use of Brain Biomechanical Models for Monitoring Impact Exposure in Contact Sports

Head acceleration measurement sensors are now widely deployed in the field to monitor head kinematic exposure in contact sports. The wealth of impact kinematics data provides valuable, yet challenging, opportunities to study the biomechanical basis of mild traumatic brain injury (mTBI) and subconcussive kinematic exposure. Head impact kinematics are translated into brain mechanical responses through physics-based computational simulations using validated brain models to study the mechanisms of injury. First, this article reviews representative legacy and contemporary brain biomechanical models primarily used for blunt impact simulation. Then, it summarizes perspectives regarding the development and validation of these models, and discusses how simulation results can be interpreted to facilitate injury risk assessment and head acceleration exposure monitoring in the context of contact sports. Recommendations and consensus statements are presented on the use of validated brain models in conjunction with kinematic sensor data to understand the biomechanics of mTBI and subconcussion. Mainly, there is general consensus that validated brain models have strong potential to improve injury prediction and interpretation of subconcussive kinematic exposure over global head kinematics alone. Nevertheless, a major roadblock to this capability is the lack of sufficient data encompassing different sports, sex, age and other factors. The authors recommend further integration of sensor data and simulations with modern data science techniques to generate large datasets of exposures and predicted brain responses along with associated clinical findings. These efforts are anticipated to help better understand the biomechanical basis of mTBI and improve the effectiveness in monitoring kinematic exposure in contact sports for risk and injury mitigation purposes.

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 Strain-Rate Dependent Constitutive Model for Göttingen Minipig Cerebral Arteries

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 the present 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.

Cerebral vascular strains in dynamic head impact using an upgraded model with brain material property heterogeneity

Cerebral vascular injury (CVI) is a frequent consequence of traumatic brain injury but has often been neglected. Substantial experimental work exists on vascular material properties and failure/subfailure thresholds. However, little is known about vascular in vivo loading conditions in dynamic head impact, which is necessary to investigate the risk, severity, and extent of CVI. In this study, we resort to the Worcester Head Injury Model (WHIM) V2.1 for investigation. The model embeds the cerebral vasculature network and is further upgraded to incorporate brain material property heterogeneity based on magnetic resonance elastography. The brain material property is calibrated to match with the previously validated anisotropic V1.0 version in terms of whole-brain strains against six experimental datasets of a wide range of blunt impact conditions. The upgraded WHIM is finally used to simulate five representative real-world head impacts drawn from contact sports and automotive crashes. We find that peak strains in veins are considerably higher than those in arteries and that peak circumferential strains are also higher than peak axial strains. For a typical concussive head impact, cerebral vascular axial strains reach the lowest reported yield strain of ∼7-8%. For severe automotive impacts, axial strains could reach ∼20%, which is on the order of the lowest reported ultimate failure strain of ∼24%. These results suggest in vivo mechanical loading conditions of the cerebral vasculature (excluding bridging veins not assessed here) due to rapid head rotation are at the lower end of failure/subfailure thresholds established from ex vivo experiments. This study provides some first insight into the risk, severity, and extent of CVI in real-world head impacts.

The Effect of 3D Whole, Major, and Small Vasculature On Mouse Brain Strain Under Both Diffuse and Focal Brain Injury Loading

Blood vessels are much stiffer than brain parenchyma and their effects in finite element (FE) brain models need to be investigated. Despite the publication of some comprehensive three-dimensional (3D) brain vasculature models, no mechanical model exists for the mouse brain vasculature. Moreover, how the vasculature affects the mechanical behavior of brain tissue remains controversial. Therefore, we developed FE mouse brain models with detailed 3D vasculature to investigate the effect of the vasculature on brain strains under both diffuse (closed-head impact) and focal injury (controlled cortical impact (CCI)) loading, two commonly laboratory models of traumatic brain injury. The effect of the vasculature was examined by comparing maximum principal strain in mouse brain FE models with and without the vasculature. On average, modeling comprehensive vasculature under diffuse injury loading reduced average brain strain predictions by 32% with non-linear elastic properties. Nearly three-fourths of the 32% strain reduction was attributable to the effects of the major branches of the vasculature. Meanwhile, during focal open-skull CCI injury loading, the contribution of the vasculature was limited, producing a less than 5% reduction in all cases. Overall, the vasculature, especially the major branches, increased the load-bearing capacity of the brain FE model and thus reduced brain strain predictions.

Displacement voxelization to resolve mesh-image mismatch: Application in deriving dense white matter fiber strains

BACKGROUND AND OBJECTIVE

It is common to combine biomechanical modeling and medical images for multimodal analyses. However, mesh-image mismatch may occur that prevents direct information exchange. To eliminate mesh-image mismatch, we develop a simple but elegant displacement voxelization technique based on image voxel corner nodes to achieve voxel-wise strain. We then apply the technique to derive dense white matter fiber strains along whole-brain tractography (∼35 k fiber tracts consisting of ∼3.3 million sampling points) resulting from head impact.

METHODS

Displacements at image voxel corner nodes are first obtained from model simulation via scattered interpolation. Each voxel is then scaled linearly to form a unit hexahedral element. This allows convenient and efficient voxel-wise strain tensor calculation and displacement interpolation at arbitrary fiber sampling points via shape functions. Fiber strains from displacement interpolation are then compared with those from the commonly used strain tensor projection using either voxel- or element-wise strain tensors.

RESULTS

Based on a synthetic displacement field, fiber strains interpolated from voxelized displacement are considerably more accurate than those from strain tensor projection relative to the prescribed ground-truth (determinant of coefficient (R²) of 1.00 and root mean squared error (RMSE) of 0.01 vs. 0.87 and 0.10, respectively). For a set of real-world reconstructed head impacts (N = 53), the strain tensor projection method performs similarly poorly (R² of 0.80-0.90 and RMSE of 0.03-0.07), with overestimation strongly correlated with strain magnitude (Pearson correlation coefficient >0.9). Up to ∼15% of the fiber strains are overestimated by more than the lower bound of a conservative injury threshold of 0.09. The percentage increases to ∼37% when halving the threshold. Voxel interpolation is also significantly more efficient (15 s vs. 40 s for element strain tensor projection, without parallelization).

CONCLUSIONS

Voxelized displacement interpolation is considerably more accurate and efficient in deriving dense white matter fiber strains than strain tensor projection. The latter generally overestimates with overestimation magnitude strongly correlating with fiber strain magnitude. Displacement voxelization is an effective technique to eliminate mesh-image mismatch and generates a convenient image representation of tissue deformation. This technique can be generalized to broadly facilitate a diverse range of image-related biomechanical problems for multimodal analyses. The convenient image format may also promote and facilitate biomechanical data sharing in the future.

Injury Metrics for Assessing the Risk of Acute Subdural Hematoma in Traumatic Events

Worldwide, the ocurrence of acute subdural hematomas (ASDHs) in road traffic crashes is a major public health problem. ASDHs are usually produced by loss of structural integrity of one of the cerebral bridging veins (CBVs) linking the parasagittal sinus to the brain. Therefore, to assess the risk of ASDH it is important to know the mechanical conditions to which the CBVs are subjected during a potentially traumatic event (such as a traffic accident or a fall from height). Recently, new studies on CBVs have been published allowing much more accurate prediction of the likelihood of mechanical failure of CBVs. These new data can be used to propose new damage metrics, which make more accurate predictions about the probability of occurrence of ASDH in road crashes. This would allow a better assessement of the effects of passive safety countermeasures and, consequently, to improve vehicle restraint systems. Currently, some widely used damage metrics are based on partially obsolete data and measurements of the mechanical behavior of CBVs that have not been confirmed by subsequent studies. This paper proposes a revision of some existing metrics and constructs a new metric based on more accurate recent data on the mechanical failure of human CBVs.

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.

Instantaneous Brain Strain Estimation for Automotive Head Impacts via Deep Learning

Efficient brain strain estimation is critical for routine application of a head injury model. Lately, a convolutional neural network (CNN) has been successfully developed to estimate spatially detailed brain strains instantly and accurately in contact sports. Here, we extend its application to automotive head impacts, where impact profiles are typically more complex with longer durations. Head impact kinematics (N=458) from two public databases were used to generate augmented impacts (N=2694). They were simulated using the anisotropic Worcester Head Injury Model (WHIM) V1.0, which provided baseline elementwise peak maximum principal strain (MPS). For each augmented impact, rotational velocity (v_rot) and the corresponding rotational acceleration (a_rot) profiles were concatenated as static images to serve as CNN input. Three training strategies were evaluated: 1) "baseline", using random initial weights; 2) "transfer learning", using weight transfer from a previous CNN model trained on head impacts drawn from contact sports; and 3) "combined training", combining previous training data from contact sports (N=5661) for training. The combined training achieved the best performances. For peak MPS, the CNN achieved a coefficient of determination (R²) of 0.932 and root mean squared error (RMSE) of 0.031 for the real-world testing dataset. It also achieved a success rate of 60.5% and 94.8% for elementwise MPS, where the linear regression slope, k, and correlation coefficient, r, between estimated and simulated MPS did not deviate from 1.0 (when identical) by more than 0.1 and 0.2, respectively. Cumulative strain damage measure (CSDM) from the CNN estimation was also highly accurate compared to those from direct simulation across a range of thresholds (R² of 0.899-0.943 with RMSE of 0.054-0.069). Finally, the CNN achieved an average k and r of 0.98±0.12 and 0.90±0.07, respectively, for six reconstructed car crash impacts drawn from two other sources independent of the training dataset. Importantly, the CNN is able to efficiently estimate elementwise MPS with sufficient accuracy while conventional kinematic injury metrics cannot. Therefore, the CNN has the potential to supersede current kinematic injury metrics that can only approximate a global peak MPS or CSDM. The CNN technique developed here may offer enhanced utility in the design and development of head protective countermeasures, including in the automotive industry. This is the first study aimed at instantly estimating spatially detailed brain strains for automotive head impacts, which employs >8.8 thousand impact simulations generated from ~1.5 years of nonstop computations on a high-performance computing platform.

A Machine Learning Approach to Investigate the Uncertainty of Tissue-Level Injury Metrics for Cerebral Contusion

Controlled cortical impact (CCI) on porcine brain is often utilized to investigate the pathophysiology and functional outcome of focal traumatic brain injury (TBI), such as cerebral contusion (CC). Using a finite element (FE) model of the porcine brain, the localized brain strain and strain rate resulting from CCI can be computed and compared to the experimentally assessed cortical lesion. This way, tissue-level injury metrics and corresponding thresholds specific for CC can be established. However, the variability and uncertainty associated with the CCI experimental parameters contribute to the uncertainty of the provoked cortical lesion and, in turn, of the predicted injury metrics. Uncertainty quantification via probabilistic methods (Monte Carlo simulation, MCS) requires a large number of FE simulations, which results in a time-consuming process. Following the recent success of machine learning (ML) in TBI biomechanical modeling, we developed an artificial neural network as surrogate of the FE porcine brain model to predict the brain strain and the strain rate in a computationally efficient way. We assessed the effect of several experimental and modeling parameters on four FE-derived CC injury metrics (maximum principal strain, maximum principal strain rate, product of maximum principal strain and strain rate, and maximum shear strain). Next, we compared the in silico brain mechanical response with cortical damage data from in vivo CCI experiments on pig brains to evaluate the predictive performance of the CC injury metrics. Our ML surrogate was capable of rapidly predicting the outcome of the FE porcine brain undergoing CCI. The now computationally efficient MCS showed that depth and velocity of indentation were the most influential parameters for the strain and the strain rate-based injury metrics, respectively. The sensitivity analysis and comparison with the cortical damage experimental data indicate a better performance of maximum principal strain and maximum shear strain as tissue-level injury metrics for CC. These results provide guidelines to optimize the design of CCI tests and bring new insights to the understanding of the mechanical response of brain tissue to focal traumatic brain injury. Our findings also highlight the potential of using ML for computationally efficient TBI biomechanics investigations.

Displacement Error Propagation From Embedded Markers to Brain Strain

Head injury model validation has evolved from against pressure to relative brain-skull displacement, and more recently, against marker-based strain. However, there are concerns on strain data quality. In this study, we parametrically investigate how displacement random errors and synchronization errors propagate into strain. Embedded markers from four representative configurations are used to form unique and nonoverlapping tetrahedrons, triangles, and linear elements. Marker displacements are then separately subjected to up to ±10% random displacement errors and up to ±2 ms synchronization errors. Based on 100 random trials in each perturbation test, we find that smaller strain errors relative to the baseline peak strains are significantly associated with larger element sizes (volume, area, or length; p < 0.05). When displacement errors are capped at the two extreme levels, the earlier "column" and "cluster" configurations provide few usable elements with relative strain error under an empirical threshold of 20%, while about 30-80% of elements in recent "repeatable" and "uniform" configurations are considered otherwise usable. Overall, denser markers are desired to provide exhaustive pairwise linear elements with a range of sizes to balance the need for larger elements to minimize strain error but smaller elements to increase the spatial resolution in strain sampling. Their signed strains also provide unique and unambiguous information on tissue tension and compression. This study may provide useful insights into the scrutinization of existing experimental data for head injury model strain validation and to inform how best to design new experiments in the future.

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.

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.

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.

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.

Validation and Comparison of Instrumented Mouthguards for Measuring Head Kinematics and Assessing Brain Deformation in Football Impacts

Because of the rigid coupling between the upper dentition and the skull, instrumented mouthguards have been shown to be a viable way of measuring head impact kinematics for assisting in understanding the underlying biomechanics of concussions. This has led various companies and institutions to further develop instrumented mouthguards. However, their use as a research tool for understanding concussive impacts makes quantification of their accuracy critical, especially given the conflicting results from various recent studies. Here we present a study that uses a pneumatic impactor to deliver impacts characteristic to football to a Hybrid III headform, in order to validate and compare five of the most commonly used instrumented mouthguards. We found that all tested mouthguards gave accurate measurements for the peak angular acceleration, the peak angular velocity, brain injury criteria values (mean average errors < 13, 8, 13%, respectively), and the mouthguards with long enough sampling time windows are suitable for a convolutional neural network-based brain model to calculate the brain strain (mean average errors < 9%). Finally, we found that the accuracy of the measurement varies with the impact locations yet is not sensitive to the impact velocity for the most part.