Development and validation of an atlas-based finite element brain model

The Impact of Drop Test Conditions on Brain Strain Location and Severity: A Novel Approach Using a Deep Learning Model

In contact sports such as rugby, players are at risk of sustaining traumatic brain injuries (TBI) due to high-intensity head impacts that generate high linear and rotational accelerations of the head. Previous studies have established a clear link between high-intensity head impacts and brain strains that result in concussions. This study presents a novel approach to investigating the effect of a range of laboratory controlled drop test parameters on regional peak and mean maximum principal strain (MPS) predictions within the brain using a trained convolutional neural network (CNN). The CNN is publicly available at https://github.com/Jilab-biomechanics/CNN-brain-strains . The results of this study corroborate previous findings that impacts to the side of the head result in significantly higher regional MPS than forehead impacts. Forehead impacts tend to result in the lowest region-averaged MPS values for impacts where the surface angle was at 0° and 45°, while side impacts tend to result in higher regional peak and mean MPS. The absence of a neck in drop tests resulted in lower regional peak and mean MPS values. The results indicated that the relationship between drop test parameters and resulting regional peak and mean MPS predictions is complex. The study's findings offer valuable insights into how deep learning models can be used to provide more detailed insights into how drop test conditions impact regional MPS. The novel approach used in this paper to predict brain strains can be applied in the development of better methods to reduce the brain strain resulting from head accelerations such as protective sports headgear.

A comparison of brain retraction mechanisms using finite element analysis and the effects of regionally heterogeneous material properties

Finite element (FE) simulations of the brain undergoing neurosurgical procedures present us with the great opportunity to better investigate, understand, and optimize surgical techniques and equipment. FE models provide access to data such as the stress levels within the brain that would otherwise be inaccessible with the current medical technology. Brain retraction is often a dangerous but necessary part of neurosurgery, and current research focuses on minimizing trauma during the procedure. In this work, we present a simulation-based comparison of different types of retraction mechanisms. We focus on traditional spatulas and tubular retractors. Our results show that tubular retractors result in lower average predicted stresses, especially in the subcortical structures and corpus callosum. Additionally, we show that changing the location of retraction can greatly affect the predicted stress results. As the model predictions highly depend on the material model and parameters used for simulations, we also investigate the importance of using region-specific hyperelastic and viscoelastic material parameters when modelling a three-dimensional human brain during retraction. Our investigations demonstrate how FE simulations in neurosurgical techniques can provide insight to surgeons and medical device manufacturers. They emphasize how further work into this direction could greatly improve the management and prevention of injury during surgery. Additionally, we show the importance of modelling the human brain with region-dependent parameters in order to provide useful predictions for neurosurgical procedures.

Validation of a computational biomechanical mouse brain model for rotational head acceleration

Recent mouse brain injury experiments examine diffuse axonal injury resulting from accelerative head rotations. Evaluating brain deformation during these events would provide valuable information on tissue level thresholds for brain injury, but there are many challenges to imaging the brain's mechanical response during dynamic loading events, such as a blunt head impact. To address this shortcoming, we present an experimentally validated computational biomechanics model of the mouse brain that predicts tissue deformation, given the motion of the mouse head during laboratory experiments. First, we developed a finite element model of the mouse brain that computes tissue strains, given the same head rotations as previously conducted in situ hemicephalic mouse brain experiments. Second, we calibrated the model using a single brain segment, and then validated the model based on the spatial and temporal strain responses of other regions. The result is a computational tool that will provide researchers with the ability to predict brain tissue strains that occur during mouse laboratory experiments, and to link the experiments to the resulting neuropathology, such as diffuse axonal injury.

Brain Material Properties and Integration of Arachnoid Complex for Biofidelic Impact Response for Human Head Finite Element Model

Finite element head models offer great potential to study brain-related injuries; however, at present may be limited by geometric and material property simplifications required for continuum-level human body models. Specifically, the mechanical properties of the brain tissues are often represented with simplified linear viscoelastic models, or the material properties have been optimized to specific impact cases. In addition, anatomical structures such as the arachnoid complex have been omitted or implemented in a simple lumped manner. Recent material test data for four brain regions at three strain rates in three modes of loading (tension, compression, and shear) was used to fit material parameters for a hyper-viscoelastic constitutive model. The material model was implemented in a contemporary detailed head finite element model. A detailed representation of the arachnoid trabeculae was implemented with mechanical properties based on experimental data. The enhanced head model was assessed by re-creating 11 ex vivo head impact scenarios and comparing the simulation results with experimental data. The hyper-viscoelastic model faithfully captured mechanical properties of the brain tissue in three modes of loading and multiple strain rates. The enhanced head model showed a high level of biofidelity in all re-created impacts in part due to the improved brain-skull interface associated with implementation of the arachnoid trabeculae. The enhanced head model provides an improved predictive capability with material properties based on tissue level data and is positioned to investigate head injury and tissue damage in the future.

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.

A Morphologically Individualized Deep Learning Brain Injury Model

The brain injury modeling community has recommended improving model subject specificity and simulation efficiency. Here, we extend an instantaneous (< 1 sec) convolutional neural network (CNN) brain model based on the anisotropic Worcester Head Injury Model (WHIM) V1.0 to account for strain differences due to individual morphological variations. Linear scaling factors relative to the generic WHIM along the three anatomical axes are used as additional CNN inputs. To generate training samples, the WHIM is randomly scaled to pair with augmented head impacts randomly generated from real-world data for simulation. An estimation of voxelized peak maximum principal strain of the whole-brain is said to be successful when the linear regression slope and Pearson's correlation coefficient relative to directly simulated do not deviate from 1.0 (when identical) by more than 0.1. Despite a modest training dataset (N = 1363 vs. ∼5.7 k previously), the individualized CNN achieves a success rate of 86.2% in cross-validation for scaled model responses, and 92.1% for independent generic model testing for impacts considered as complete capture of kinematic events. Using 11 scaled subject-specific models (with scaling factors determined from pre-established regression models based on head dimensions and sex and age information, and notably, without neuroimages), the morphologically individualized CNN remains accurate for impacts that also yield successful estimations for the generic WHIM. The individualized CNN instantly estimates subject-specific and spatially detailed peak strains of the entire brain and thus, supersedes others that report a scalar peak strain value incapable of informing the location of occurrence. This tool could be especially useful for youths and females due to their anticipated greater morphological differences relative to the generic model, even without the need for individual neuroimages. It has potential for a wide range of applications for injury mitigation purposes and the design of head protective gears. The voxelized strains also allow for convenient data sharing and promote collaboration among research groups.

A sensor-enabled cloud-based computing platform for computational brain biomechanics

BACKGROUND AND OBJECTIVES

Driven by the risk of repetitive head trauma, sensors have been integrated into mouthguards to measure head impacts in contact sports and military activities. These wearable devices, referred to as "instrumented" or "smart" mouthguards are being actively developed by various research groups and organizations. These instrumented mouthguards provide an opportunity to further study and understand the brain biomechanics due to impact. In this study, we present a brain modeling service that can use information from these sensors to predict brain injury metrics in an automated fashion.

METHODS

We have built a brain modeling platform using several of Amazon's Web Services (AWS) to enable cloud computing and scalability. We use a custom-built cloud-based finite element modeling code to compute the physics-based nonlinear response of the intracranial brain tissue and provide a frontend web application and an application programming interface for groups working on head impact sensor technology to include simulated injury predictions into their research pipeline.

RESULTS

The platform results have been validated against experimental data available in literature for brain-skull relative displacements, brain strains and intracranial pressure. The parallel processing capability of the platform has also been tested and verified. We also studied the accuracy of the custom head surfaces generated by Avatar 3D.

CONCLUSION

We present a validated cloud-based computational brain modeling platform that uses sensor data as input for numerical brain models and outputs a quantitative description of brain tissue strains and injury metrics. The platform is expected to generate transparent, reproducible, and traceable brain computing results.

In silico approach towards neuro-occlusal rehabilitation for the early correction of asymmetrical development in a unilateral crossbite patient

Neuro-occlusal rehabilitation (N.O.R.) is a discipline of the stomatognathic medicine that defends early treatments of functional malocclusions, such as unilateral crossbite, for the correction of craniofacial development, avoiding surgical procedures later in life. Nevertheless, N.O.R.'s advances have not been proved analytically yet due to the difficulties of evaluate the mechanical response after the treatment. This study aims to evaluate computationally the effect of N.O.R.'s treatments during childhood. Therefore, bilateral chewing and maximum intercuspation occlusion were modelled through a detailed finite element model of a paediatric craniofacial complex, before and after different selective grinding-alternatives. This model was subjected to the muscular forces derived from a musculoskeletal model and was validated by the occlusal contacts recorded experimentally. This approach yielded errors below 2% and reproduced successfully the occlusal, muscular, functional and mechanical imbalance before the therapies. Treatment strategies balanced the occlusal plane and reduced the periodontal overpressure (>4.7 kPa) and the mandibular over deformation (>0.002 ε) on the crossed side. Based on the principles of the mechanostat theory of bone remodelling and the pressure-tension theory of tooth movement, these findings could also demonstrate how N.O.R.'s treatments correct the malocclusion and the asymmetrical development of the craniofacial complex. Besides, N.O.R.'s treatments slightly modified the stress state and functions of the temporomandibular joints, facilitating the chewing by the unaccustomed side. These findings provide important biomechanical insights into the use of N.O.R.'s treatments for the correction of unilateral crossbite, but also encourage the application of computing methods in biomedical research and clinical practise.

Spinal Cord Boundary Conditions Affect Brain Tissue Strains in Impact Simulations

Brain and spinal cord injuries have devastating consequences on quality of life but are challenging to assess experimentally due to the traumatic nature of such injuries. Finite element human body models (HBM) have been developed to investigate injury but are limited by a lack of biofidelic spinal cord implementation. In many HBM, brain models terminate with a fixed boundary condition at the brain stem. The goals of this study were to implement a comprehensive representation of the spinal cord into a contemporary head and neck HBM, and quantify the effect of the spinal cord on brain deformation during simulated impacts. Spinal cord tissue geometries were developed, based on 3D medical imaging and literature data, meshed, and implemented into the GHBMC 50th percentile male model. The model was evaluated in frontal, lateral, rear, and oblique impact conditions, and the resulting maximum principal strains in the brain tissue were compared, with and without the spinal cord. A new cumulative strain curve metric was proposed to quantify brain strain distribution. Presence of the spinal cord increased brain tissue strains in all simulated cases, owing to a more compliant boundary condition, highlighting the importance of the spinal cord to assess brain response during impact.

A Review of Cyclist Head Injury, Impact Characteristics and the Implications for Helmet Assessment Methods

Head injuries are common for cyclists involved in collisions. Such collision scenarios result in a range of injuries, with different head impact speeds, angles, locations, or surfaces. A clear understanding of these collision characteristics is vital to design high fidelity test methods for evaluating the performance of helmets. We review literature detailing real-world cyclist collision scenarios and report on these key characteristics. Our review shows that helmeted cyclists have a considerable reduction in skull fracture and focal brain pathologies compared to non-helmeted cyclists, as well as a reduction in all brain pathologies. The considerable reduction in focal head pathologies is likely to be due to helmet standards mandating thresholds of linear acceleration. The less considerable reduction in diffuse brain injuries is likely to be due to the lack of monitoring head rotation in test methods. We performed a novel meta-analysis of the location of 1809 head impacts from ten studies. Most studies showed that the side and front regions are frequently impacted, with one large, contemporary study highlighting a high proportion of occipital impacts. Helmets frequently had impact locations low down near the rim line. The face is not well protected by most conventional bicycle helmets. Several papers determine head impact speed and angle from in-depth reconstructions and computer simulations. They report head impact speeds from 5 to 16 m/s, with a concentration around 5 to 8 m/s and higher speeds when there was another vehicle involved in the collision. Reported angles range from 10° to 80° to the normal, and are concentrated around 30°-50°. Our review also shows that in nearly 80% of the cases, the head impact is reported to be against a flat surface. This review highlights current gaps in data, and calls for more research and data to better inform improvements in testing methods of standards and rating schemes and raise helmet safety.

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.

Biomechanics of Traumatic Head and Neck Injuries on Women: A State-of-the-Art Review and Future Directions

The biomechanics of traumatic injuries of the human body as a consequence of road crashes, falling, contact sports, and military environments have been studied for decades. In particular, traumatic brain injury (TBI), the so-called "silent epidemic", is the traumatic insult responsible for the greatest percentage of death and disability, justifying the relevance of this research topic. Despite its great importance, only recently have research groups started to seriously consider the sex differences regarding the morphology and physiology of women, which differs from men and may result in a specific outcome for a given traumatic event. This work aims to provide a summary of the contributions given in this field so far, from clinical reports to numerical models, covering not only the direct injuries from inertial loading scenarios but also the role sex plays in the conditions that precede an accident, and post-traumatic events, with an emphasis on neuroendocrine dysfunctions and chronic traumatic encephalopathy. A review on finite element head models and finite element neck models for the study of specific traumatic events is also performed, discussing whether sex was a factor in validating them. Based on the information collected, improvement perspectives and future directions are discussed.

Evaluation of Deep Brain Stimulation (DBS) Lead Biomechanical Interaction with Brain Tissue

The current study aims to examine the effect of material properties on implanted leads used for deep brain stimulation (DBS) using finite element (FE) analysis to investigate brain deformation around an implanted DBS lead in response to daily head accelerations. FE analysis was used to characterize the relative motion of the DBS lead in a suite of fifteen cases sampled from a previously derived kinematic envelope representative of everyday activities describing translational and rotational pulse shape, magnitude, and duration. Load curves were applied to the atlas-based brain model (ABM) with a scaled Haversine acceleration pulse in four directions of rotation: + X, - Y, + Y, and + Z. In addition to the fifteen sampled cases, six experimental cases taken from a previous literature review were also simulated for comparison. The current investigation found that there was very little difference in brain response for the DBS leads with two different material properties. In general, the brain and DBS lead experienced the greatest deformation during rotation about the Z axis for similar load cases. In conclusion, this study showed that there was no significant difference in implanted DBS lead deformation based on lead material properties.

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.

Development and validation of subject-specific 3D human head models based on a nonlinear visco-hyperelastic constitutive framework

Computational head models are promising tools for understanding and predicting traumatic brain injuries. Most available head models are developed using inputs (i.e. head geometry, material properties and boundary conditions) from experiments on cadavers or animals and employ hereditary integral-based constitutive models that assume linear viscoelasticity in part of the rate-sensitive material response. This leads to high uncertainty and poor accuracy in capturing the nonlinear brain tissue response. To resolve these issues, a framework for the development of subject-specific three-dimensional head models is proposed, in which all inputs are derived in vivo from the same living human subject: head geometry via magnetic resonance imaging (MRI), brain tissue properties via magnetic resonance elastography (MRE), and full-field strain-response of the brain under rapid head rotation via tagged MRI. A nonlinear, viscous dissipation-based visco-hyperelastic constitutive model is employed to capture brain tissue response. Head models are validated using quantitative metrics that compare spatial strain distribution, temporal strain evolution, and the magnitude of strain maxima, with the corresponding experimental observations from tagged MRI. Results show that our head models accurately capture the strain-response of the brain. Further, employment of the nonlinear visco-hyperelastic constitutive framework provides improvements in the prediction of peak strains and temporal strain evolution over hereditary integral-based models.

Concussion Prone Scenarios: A Multi-Dimensional Exploration in Impact Directions, Brain Morphology, and Network Architectures Using Computational Models

While individual susceptibility to traumatic brain injury (TBI) has been speculated, past work does not provide an analysis considering how physical features of an individual's brain (e.g., brain size, shape), impact direction, and brain network features can holistically contribute to the risk of suffering a TBI from an impact. This work investigated each of these features simultaneously using computational modeling and analyses of simulated functional connectivity. Unlike the past studies that assess the severity of TBI based on the quantification of brain tissue damage (e.g., principal strain), we approached the brain as a complex network in which neuronal oscillations orchestrate to produce normal brain function (estimated by functional connectivity) and, to this end, both the anatomical damage location and its topological characteristics within the brain network contribute to the severity of brain function disruption and injury. To represent the variations in the population, we analyzed a publicly available database of brain imaging data and selected five distinct network architectures, seven different brain sizes, and three uniaxial head rotational conditions to study the consequences of 74 virtual impact scenarios. Results show impact direction produces the most significant change in connections across brain areas (structural connectome) and the functional coupling of activity across these brain areas (functional connectivity). Axial rotations were more injurious than those with sagittal and coronal rotations when the head kinematics were the same for each condition. When the impact direction was held constant, brain network architecture showed a significantly different vulnerability across axial and sagittal, but not coronal rotations. As expected, brain size significantly affected the expected change in structural and functional connectivity after impact. Together, these results provided groupings of predicted vulnerability to impact-a subgroup of male brain architectures exposed to axial impacts were most vulnerable, while a subgroup of female brain architectures was the most tolerant to the sagittal impacts studied. These findings lay essential groundwork for subject-specific analyses of concussion and provide invaluable guidance for designing personalized protection equipment.

American Football Helmet Effectiveness Against a Strain-Based Concussion Mechanism

Brain strain is increasingly being used in helmet design and safety performance evaluation as it is generally considered as the primary mechanism of concussion. In this study, we investigate whether different helmet designs can meaningfully alter brain strains using two commonly used metrics, peak maximum principal strain (MPS) of the whole brain and cumulative strain damage measure (CSDM). A convolutional neural network (CNN) that instantly produces detailed brain strains is first tested for accuracy for helmeted head impacts. Based on N = 144 impacts in 12 impact conditions from three random and representative helmet models, we conclude that the CNN is sufficiently accurate for helmet testing applications, for elementwise MPS (success rate of 98.6%), whole-brain peak MPS and CSDM (coefficient of determination of 0.977 and 0.980, with root mean squared error of 0.015 and 0.029, respectively). We then apply the technique to 23 football helmet models (N = 1104 impacts) to reproduce elementwise MPS. Assuming a concussion would occur when peak MPS or CSDM exceeds a threshold, we sweep their thresholds across the value ranges to evaluate the number of predicted hypothetical concussions that different helmets sustain across the impact conditions. Relative to the 12 impact conditions tested, we find that the "best" and "worst" helmets differ by an average of 22.5% in terms of predicted concussions, ranging from 0 to 42% (the latter achieved at the threshold value of 0.28 for peak MPS and 0.4 for CSDM, respectively). Such a large variation among helmets in strain-based concussion predictions demonstrate that helmet designs can still be optimized in a clinically meaningful way. The robustness and accuracy of the CNN tool also suggest its potential for routine use for helmet design and safety performance evaluation in the future. The CNN is freely available online at https://github.com/Jilab-biomechanics/CNN-brain-strains .

Real-time dynamic simulation for highly accurate spatiotemporal brain deformation from impact

Real-time dynamic simulation remains a significant challenge for spatiotemporal data of high dimension and resolution. In this study, we establish a transformer neural network (TNN) originally developed for natural language processing and a separate convolutional neural network (CNN) to estimate five-dimensional (5D) spatiotemporal brain-skull relative displacement resulting from impact (isotropic spatial resolution of 4 mm with temporal resolution of 1 ms). Sequential training is applied to train (N = 5184 samples) the two neural networks for estimating the complete 5D displacement across a temporal duration of 60 ms. We find that TNN slightly but consistently outperforms CNN in accuracy for both displacement and the resulting voxel-wise four-dimensional (4D) maximum principal strain (e.g., root mean squared error (RMSE) of ~1.0% vs. ~1.6%, with coefficient of determination, R ² >0.99 vs. >0.98, respectively, and normalized RMSE (NRMSE) at peak displacement of 2%-3%, based on an independent testing dataset; N = 314). Their accuracies are similar for a range of real-world impacts drawn from various published sources (dummy, helmet, football, soccer, and car crash; average RMSE/NRMSE of ~0.3 mm/~4%-5% and average R ² of ~0.98 at peak displacement). Sequential training is effective for allowing instantaneous estimation of 5D displacement with high accuracy, although TNN poses a heavier computational burden in training. This work enables efficient characterization of the intrinsically dynamic brain strain in impact critical for downstream multiscale axonal injury model simulation. This is also the first application of TNN in biomechanics, which offers important insight into how real-time dynamic simulations can be achieved across diverse engineering fields.

The Presence of the Temporal Horn Exacerbates the Vulnerability of Hippocampus During Head Impacts

Hippocampal injury is common in traumatic brain injury (TBI) patients, but the underlying pathogenesis remains elusive. In this study, we hypothesize that the presence of the adjacent fluid-containing temporal horn exacerbates the biomechanical vulnerability of the hippocampus. Two finite element models of the human head were used to investigate this hypothesis, one with and one without the temporal horn, and both including a detailed hippocampal subfield delineation. A fluid-structure interaction coupling approach was used to simulate the brain-ventricle interface, in which the intraventricular cerebrospinal fluid was represented by an arbitrary Lagrangian-Eulerian multi-material formation to account for its fluid behavior. By comparing the response of these two models under identical loadings, the model that included the temporal horn predicted increased magnitudes of strain and strain rate in the hippocampus with respect to its counterpart without the temporal horn. This specifically affected cornu ammonis (CA) 1 (CA1), CA2/3, hippocampal tail, subiculum, and the adjacent amygdala and ventral diencephalon. These computational results suggest that the presence of the temporal horn exacerbate the vulnerability of the hippocampus, highlighting the mechanobiological dependency of the hippocampus on the temporal horn.

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.

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.

Header biomechanics in youth and collegiate female soccer

Concerns about the effects of intentional heading in soccer have led to regulatory restrictions on headers for youth players. However, there is limited data describing how header exposure varies across age levels, and few studies have attempted to compare head impact exposure across different levels of play with the same sensor. Additionally, little is known about the biomechanical response of the brain to header impacts. The objective of this study was to evaluate head kinematics and the resulting tissue-level brain strain associated with intentional headers among youth and collegiate female soccer players. Six youth and 13 collegiate participants were instrumented with custom mouthpiece-based sensors measuring six-degree-of-freedom head kinematics of headers during practices and games. Kinematics of film-verified headers were used to drive impact simulations with a detailed brain finite element model to estimate tissue-level strain. Linear and rotational head kinematics and strain metrics, specifically 95th percentile maximum principal strain (ε_1,95) and the area under the cumulative strain damage measure curve (VSM₁), were compared across levels of play (i.e., youth vs. collegiate) while adjusting for session type and ball delivery method. A total of 483 headers (n = 227 youth, n = 256 collegiate) were analyzed. Level of play was significantly associated with linear acceleration, rotational acceleration, rotational velocity, ε_1,95, and VSM₁. Headers performed by collegiate players had significantly greater mean head kinematics and strain metrics compared to those performed by youth players (all p < .001). Targeted interventions aiming to reduce head impact magnitude in soccer should consider factors associated with the level of play.

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.

Mechanical characterisation of the human dura mater, falx cerebri and superior sagittal sinus

The cranial meninges have been shown to play a pivotal role in traumatic brain injury mechanopathology. However, while the mechanical response of the brain and its many subregions have been studied extensively, the meninges have conventionally been overlooked. This paper presents the first comparative mechanical analysis of human dura mater, falx cerebri and superior sagittal sinus tissues. Biaxial tensile analysis identified that these tissues are mechanically heterogeneous, in contrast to the assumption that the tissues are mechanically homogeneous which is typically employed in FE model design. A thickness of 0.91 ± 0.05 (standard error) mm for the falx cerebri was also identified. This data can aid in improving the biofidelity of the influential falx structure in FE models. Additionally, the use of a collagen hybridizing peptide on the superior sagittal sinus suggests this structure is particularly susceptible to the effects of circumferential stretch, which may have important implications for clinical treatment of dural venous sinus pathologies. Collectively, this research progresses understanding of meningeal mechanical and structural characteristics and may aid in elucidating the behaviour of these tissues in healthy and diseased conditions. STATEMENT OF SIGNIFICANCE: This study presents the first evaluation of human falx cerebri and superior sagittal sinus mechanical, geometrical and structural properties, along with a comparison to cranial dura mater. To mechanically characterise the tissues, biaxial tensile testing is conducted on the tissues. This analysis identifies, for the first time, mechanical stiffness differences between these tissues. Additionally, geometrical analysis identifies that there are thickness differences between the tissues. The evaluation of human meningeal tissues allows for direct implementation of the novel data to finite element head injury models to enable improved biofidelity of these influential structures in traumatic brain injury simulations. This work also identifies that the superior sagittal sinus may be easily damaged during clinical angioplasty procedures, which may inform the treatment of dural sinus pathologies.

MR Imaging of Human Brain Mechanics In Vivo: New Measurements to Facilitate the Development of Computational Models of Brain Injury

Computational models of the brain and its biomechanical response to skull accelerations are important tools for understanding and predicting traumatic brain injuries (TBIs). However, most models have been developed using experimental data collected on animal models and cadaveric specimens, both of which differ from the living human brain. Here we describe efforts to noninvasively measure the biomechanical response of the human brain with MRI-at non-injurious strain levels-and generate data that can be used to develop, calibrate, and evaluate computational brain biomechanics models. Specifically, this paper reports on a project supported by the National Institute of Neurological Disorders and Stroke to comprehensively image brain anatomy and geometry, mechanical properties, and brain deformations that arise from impulsive and harmonic skull loadings. The outcome of this work will be a publicly available dataset ( http://www.nitrc.org/projects/bbir ) that includes measurements on both males and females across an age range from adolescence to older adulthood. This article describes the rationale and approach for this study, the data available, and how these data may be used to develop new computational models and augment existing approaches; it will serve as a reference to researchers interested in using these data.

Regional White Matter Diffusion Changes Associated with the Cumulative Tensile Strain and Strain Rate in Nonconcussed Youth Football Players

The purpose of this study is to assess the relationship between regional white matter diffusion imaging changes and finite element strain measures in nonconcussed youth football players. Pre- and post-season diffusion-weighted imaging was performed in 102 youth football subject-seasons, in which no concussions were diagnosed. The diffusion data were normalized to the IXI template. Percent change in fractional anisotropy (%ΔFA) images were generated. Using data from the head impact telemetry system, the cumulative maximum principal strain one times strain rate (CMPS1 × SR), a measure of the cumulative tensile brain strain and strain rate for one season, was calculated for each subject. Two linear regression analyses were performed to identify significant positive or inverse relationships between CMPS1 × SR and %ΔFA within the international consortium for brain mapping white matter mask. Age, body mass index, days between pre- and post-season imaging, previous brain injury, attention disorder diagnosis, and imaging protocol were included as covariates. False discovery rate correction was used with corrected alphas of 0.025 and voxel thresholds of zero. Controlling for all covariates, a significant, positive linear relationship between %ΔFA and CMPS1 × SR was identified in the bilateral cingulum, fornix, internal capsule, external capsule, corpus callosum, corona radiata, corticospinal tract, cerebral and middle cerebellar peduncle, superior longitudinal fasciculus, and right superior fronto-occipital fasciculus. Post hoc analyses further demonstrated significant %ΔFA differences between high-strain football subjects and noncollision control athletes, no significant %ΔFA differences between low-strain subjects and noncollision control athletes, and that CMPS1 × SR significantly explained more %ΔFA variance than number of head impacts alone.

Potential Mechanisms of Acute Standing Balance Deficits After Concussions and Subconcussive Head Impacts: A Review

Standing balance deficits are prevalent after concussions and have also been reported after subconcussive head impacts. However, the mechanisms underlying such deficits are not fully understood. The objective of this review is to consolidate evidence linking head impact biomechanics to standing balance deficits. Mechanical energy transferred to the head during impacts may deform neural and sensory components involved in the control of standing balance. From our review of acute balance-related changes, concussions frequently resulted in increased magnitude but reduced complexity of postural sway, while subconcussive studies showed inconsistent outcomes. Although vestibular and visual symptoms are common, potential injury to these sensors and their neural pathways are often neglected in biomechanics analyses. While current evidence implies a link between tissue deformations in deep brain regions including the brainstem and common post-concussion balance-related deficits, this link has not been adequately investigated. Key limitations in current studies include inadequate balance sampling duration, varying test time points, and lack of head impact biomechanics measurements. Future investigations should also employ targeted quantitative methods to probe the sensorimotor and neural components underlying balance control. A deeper understanding of the specific injury mechanisms will inform diagnosis and management of balance deficits after concussions and subconcussive head impact exposure.

Toward subject-specific evaluation: methods of evaluating finite element brain models using experimental high-rate rotational brain motion

Computational models of the brain have become the gold standard in biomechanics to understand, predict, and mitigate traumatic brain injuries. Many models have been created and evaluated with limited experimental data and without accounting for subject-specific morphometry of the specimens in the dataset. Recent advancements in the measurement of brain motion using sonomicrometry allow for a comprehensive evaluation of brain model biofidelity using a high-rate, rotational brain motion dataset. In this study, four methods were used to determine the best technique to compare nodal displacement to experimental brain motion, including a new morphing method to match subject-specific inner skull geometry. Three finite element brain models were evaluated in this study: the isotropic GHBMC and SIMon models, as well as an anisotropic model with explicitly embedded axons (UVA-EAM). Using a weighted cross-correlation score (between 0 and 1), the anisotropic model yielded the highest average scores across specimens and loading conditions ranging from 0.53 to 0.63, followed by the isotropic GHBMC with average scores ranging from 0.46 to 0.58, and then the SIMon model with average scores ranging from 0.36 to 0.51. The choice of comparison method did not significantly affect the cross-correlation score, and differences of global strain up to 0.1 were found for the morphed geometry relative to baseline models. The morphed or scaled geometry is recommended when evaluating computational brain models to capture the subject-specific skull geometry of the experimental specimens.

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.

An anatomically detailed and personalizable head injury model: Significance of brain and white matter tract morphological variability on strain

Finite element head (FE) models are important numerical tools to study head injuries and develop protection systems. The generation of anatomically accurate and subject-specific head models with conforming hexahedral meshes remains a significant challenge. The focus of this study is to present two developmental works: first, an anatomically detailed FE head model with conforming hexahedral meshes that has smooth interfaces between the brain and the cerebrospinal fluid, embedded with white matter (WM) fiber tracts; second, a morphing approach for subject-specific head model generation via a new hierarchical image registration pipeline integrating Demons and Dramms deformable registration algorithms. The performance of the head model is evaluated by comparing model predictions with experimental data of brain-skull relative motion, brain strain, and intracranial pressure. To demonstrate the applicability of the head model and the pipeline, six subject-specific head models of largely varying intracranial volume and shape are generated, incorporated with subject-specific WM fiber tracts. DICE similarity coefficients for cranial, brain mask, local brain regions, and lateral ventricles are calculated to evaluate personalization accuracy, demonstrating the efficiency of the pipeline in generating detailed subject-specific head models achieving satisfactory element quality without further mesh repairing. The six head models are then subjected to the same concussive loading to study the sensitivity of brain strain to inter-subject variability of the brain and WM fiber morphology. The simulation results show significant differences in maximum principal strain and axonal strain in local brain regions (one-way ANOVA test, p < 0.001), as well as their locations also vary among the subjects, demonstrating the need to further investigate the significance of subject-specific models. The techniques developed in this study may contribute to better evaluation of individual brain injury and the development of individualized head protection systems in the future. This study also contains general aspects the research community may find useful: on the use of experimental brain strain close to or at injury level for head model validation; the hierarchical image registration pipeline can be used to morph other head models, such as smoothed-voxel models.

Brain Strain: Computational Model-Based Metrics for Head Impact Exposure and Injury Correlation

Athletes participating in contact sports are exposed to repetitive subconcussive head impacts that may have long-term neurological consequences. To better understand these impacts and their effects, head impacts are often measured during football to characterize head impact exposure and estimate injury risk. Despite widespread use of kinematic-based metrics, it remains unclear whether any single metric derived from head kinematics is well-correlated with measurable changes in the brain. This shortcoming has motivated the increasing use of finite element (FE)-based metrics, which quantify local brain deformations. Additionally, quantifying cumulative exposure is of increased interest to examine the relationship to brain changes over time. The current study uses the atlas-based brain model (ABM) to predict the strain response to impacts sustained by 116 youth football athletes and proposes 36 new, or derivative, cumulative strain-based metrics that quantify the combined burden of head impacts over the course of a season. The strain-based metrics developed and evaluated for FE modeling and presented in the current study present potential for improved analytics over existing kinematically-based and cumulative metrics. Additionally, the findings highlight the importance of accounting for directional dependence and expand the techniques to explore spatial distribution of the strain response throughout the brain.

Mechanical Properties of the Cranial Meninges: A Systematic Review

The meninges are membranous tissues that are pivotal in maintaining homeostasis of the central nervous system. Despite the importance of the cranial meninges in nervous system physiology and in head injury mechanics, our knowledge of the tissues' mechanical behavior and structural composition is limited. This systematic review analyzes the existing literature on the mechanical properties of the meningeal tissues. Publications were identified from a search of Scopus, Academic Search Complete, and Web of Science and screened for eligibility according to Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines. The review details the wide range of testing techniques employed to date and the significant variability in the observed experimental findings. Our findings identify many gaps in the current literature that can serve as a guide for future work for meningeal mechanics investigators. The review identifies no peer-reviewed mechanical data on the falx and tentorium tissues, both of which have been identified as key structures in influencing brain injury mechanics. A dearth of mechanical data for the pia-arachnoid complex also was identified (no experimental mechanics studies on the human pia-arachnoid complex were identified), which is desirable for biofidelic modeling of human head injuries. Finally, this review provides recommendations on how experiments can be conducted to allow for standardization of test methodologies, enabling simplified comparisons and conclusions on meningeal mechanics.

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.

Integrating material properties from magnetic resonance elastography into subject-specific computational models for the human brain

Advances in brain imaging and computational methods have facilitated the creation of subject-specific computational brain models that aid researchers in investigating brain trauma using simulated impacts. The emergence of magnetic resonance elastography (MRE) as a non-invasive mechanical neuroimaging tool has enabled in vivo estimation of material properties at low-strain, harmonic loading. An open question in the field has been how this data can be integrated into computational models. The goals of this study were to use a novel MRI dataset acquired in human volunteers to generate models with subject-specific anatomy and material properties, and then to compare simulated brain deformations to subject-specific brain deformation data under non-injurious loading. Models of five subjects were simulated with linear viscoelastic (LVE) material properties estimated directly from MRE data. Model predictions were compared to experimental brain deformation acquired in the same subjects using tagged MRI. Outcomes from the models matched the spatial distribution and magnitude of the measured peak strain components as well as the 95^th percentile in-plane peak strains within 0.005 mm/mm and maximum principal strain within 0.012 mm/mm. Sensitivity to material heterogeneity was also investigated. Simulated brain deformations from a model with homogenous brain properties and a model with brain properties discretized with up to ten regions were very similar (a mean absolute difference less than 0.0015 mm/mm in peak strains). Incorporating material properties directly from MRE into a biofidelic subject-specific model is an important step toward future investigations of higher-order model features and simulations under more severe loading conditions.

Variations in User Implementation of the CORA Rating Metric

The CORA rating metric is frequently used in the field of injury biomechanics to compare the similarity of response time histories. However, subjectivity exists within the CORA metric in the form of user-customizable parameters that give the metric the flexibility to be used for a variety of applications. How these parameters are customized is not always reported in the literature, and it is unknown how these customizations affect the CORA scores. Therefore, the purpose of this study was to evaluate how variations in the CORA parameters affect the resulting similarity scores. A literature review was conducted to determine how the CORA parameters are commonly customized within the literature. Then, CORA scores for two datasets were calculated using the most common parameter customizations and the default parameters. Differences between the CORA scores using customized and default parameters were statistically significant for all customizations. Furthermore, most customizations produced score increases relative to the default settings. The use of standard deviation corridors and exclusion of the corridor component were found to produce the largest score differences. The observed differences demonstrated the need for researchers to exercise transparency when using customized parameters in CORA analyses.

A Review of Validation Methods for the Intracranial Response of FEHM to Blunt Impacts

The following is a review of the processes currently employed when validating the intracranial response of Finite Element Head Models (FEHM) against blunt impacts. The authors aim to collate existing validation tools, their applications and findings on their effectiveness to aid researchers in the validation of future FEHM and potential efforts in improving procedures. In this vain, publications providing experimental data on the intracranial pressure, relative brain displacement and brain strain responses to impacts in human subjects are surveyed and key data are summarised. This includes cases that have previously been used in FEHM validation and alternatives with similar potential uses. The processes employed to replicate impact conditions and the resulting head motion are reviewed, as are the analytical techniques used to judge the validity of the models. Finally, publications exploring the validation process and factors affecting it are critically discussed. Reviewing FEHM validation in this way highlights the lack of a single best practice, or an obvious solution to create one using the tools currently available. There is clear scope to improve the validation process of FEHM, and the data available to achieve this. By collecting information from existing publications, it is hoped this review can help guide such developments and provide a point of reference for researchers looking to validate or investigate FEHM in the future, enabling them to make informed choices about the simulation of impacts, how they are generated numerically and the factors considered during output assessment, whilst being aware of potential limitations in the process.

An anatomically detailed and personalizable head injury model: Significance of brain and white matter tract morphological variability on strain

Finite element head (FE) models are important numerical tools to study head injuries and develop protection systems. The generation of anatomically accurate and subject-specific head models with conforming hexahedral meshes remains a significant challenge. The focus of this study is to present two developmental works: first, an anatomically detailed FE head model with conforming hexahedral meshes that has smooth interfaces between the brain and the cerebrospinal fluid, embedded with white matter (WM) fiber tracts; second, a morphing approach for subject-specific head model generation via a new hierarchical image registration pipeline integrating Demons and Dramms deformable registration algorithms. The performance of the head model is evaluated by comparing model predictions with experimental data of brain-skull relative motion, brain strain, and intracranial pressure. To demonstrate the applicability of the head model and the pipeline, six subject-specific head models of largely varying intracranial volume and shape are generated, incorporated with subject-specific WM fiber tracts. DICE similarity coefficients for cranial, brain mask, local brain regions, and lateral ventricles are calculated to evaluate personalization accuracy, demonstrating the efficiency of the pipeline in generating detailed subject-specific head models achieving satisfactory element quality without further mesh repairing. The six head models are then subjected to the same concussive loading to study the sensitivity of brain strain to inter-subject variability of the brain and WM fiber morphology. The simulation results show significant differences in maximum principal strain and axonal strain in local brain regions (one-way ANOVA test, p < 0.001), as well as their locations also vary among the subjects, demonstrating the need to further investigate the significance of subject-specific models. The techniques developed in this study may contribute to better evaluation of individual brain injury and the development of individualized head protection systems in the future. This study also contains general aspects the research community may find useful: on the use of experimental brain strain close to or at injury level for head model validation; the hierarchical image registration pipeline can be used to morph other head models, such as smoothed-voxel models.

An Image Registration-Based Morphing Technique for Generating Subject-Specific Brain Finite Element Models

Finite element (FE) models of the brain are crucial for investigating the mechanisms of traumatic brain injury (TBI). However, FE brain models are often limited to a single neuroanatomy because the manual development of subject-specific models is time consuming. The objective of this study was to develop a pipeline to automatically generate subject-specific FE brain models using previously developed nonlinear image registration techniques, preserving both external and internal neuroanatomical characteristics. To verify the morphing-induced mesh distortions did not influence the brain deformation response, strain distributions predicted using the morphed model were compared to those from manually created voxel models of the same subject. Morphed and voxel models were generated for 44 subjects ranging in age, and simulated using head kinematics from a football concussion case. For each subject, brain strain distributions predicted by each model type were consistent, and differences in strain prediction was less than 4% between model type. This automated technique, taking approximately 2 h to generate a subject-specific model, will facilitate interdisciplinary research between the biomechanics and neuroimaging fields and could enable future use of biomechanical models in the clinical setting as a tool for improving diagnosis.

A machine learning approach for magnetic resonance image-based mouse brain modeling and fast computation in controlled cortical impact

Computational modeling of the brain is crucial for the study of traumatic brain injury. An anatomically accurate model with refined details could provide the most accurate computational results. However, computational models with fine mesh details could take prolonged computation time that impedes the clinical translation of the models. Therefore, a way to construct a model with low computational cost while maintaining a computational accuracy comparable with that of the high-fidelity model is desired. In this study, we constructed magnetic resonance (MR) image-based finite element (FE) models of a mouse brain for simulations of controlled cortical impact. The anatomical details were kept by mapping each image voxel to a corresponding FE mesh element. We constructed a super-resolution neural network that could produce computational results of a refined FE model with a mesh size of 70 μm from a coarse FE model with a mesh size of 280 μm. The peak signal-to-noise ratio of the reconstructed results was 33.26 dB, while the computational speed was increased by 50-fold. This proof-of-concept study showed that using machine learning techniques, MR image-based computational modeling could be applied and evaluated in a timely fashion. This paved ways for fast FE modeling and computation based on MR images. Results also support the potential clinical applications of MR image-based computational modeling of the human brain in a variety of scenarios such as brain impact and intervention.Graphical abstract MR image-based FE models with different mesh sizes were generated for CCI. The training and testing data sets were computed with 5 different impact locations and 3 different impact velocities. High-resolution strain maps were estimated using a SR neural network with greatly reduced computational cost.

Simulation of the Strain Amplification in Sulci Due to Blunt Impact to the Head

Traumatic brain injury (TBI) has become a concern in sports, automobile accidents and combat operations. A better understanding of the mechanics leading to a TBI is required to cope with both the short-term life-threatening effects and long-term effects of TBIs, such as the development chronic traumatic encephalopathy (CTE). Kornguth et al. (1) proposed that an inflammatory and autoimmune process initiated by a water hammer effect at the bases of the sulci of the brain is a mechanism of TBI leading to CTE. A major objective of this study is to investigate whether the water hammer effect is present due to blunt impacts through the use of computational models. Frontal blunt impacts were simulated with 2D finite element models developed to capture the biofidelic geometry of a human head. The models utilized the Arbitrary Lagrangian Eulerian (ALE) method to model a layer of cerebrospinal fluid (CSF) as a deforming fluid allowing for CSF to move in and out of sulci. During the simulated impacts, CSF was not observed to be driven into the sulci during the transient response. However, elevated shear strain levels near the base of the sulci were exhibited. Further, increased shear strain was present when differentiation between white and gray matter was taken into account. Both of the results support clinical observations of (1).

Biomechanics of Periventricular Injury

Periventricular injury is frequently noted as one aspect of severe traumatic brain injury (TBI) and the presence of the ventricles has been hypothesized to be a primary pathogenesis associated with the prevalence of periventricular injury in patients with TBI. Although substantial endeavors have been made to elucidate the potential mechanism, a thorough explanation for this hypothesis appears lacking. In this study, a three-dimensional (3D) finite element (FE) model of the human head with an accurate representation of the cerebral ventricles is developed accounting for the fluid properties of the intraventricular cerebrospinal fluid (CSF) as well as its interaction with the brain. An additional model is developed by replacing the intraventricular CSF with a substitute with brain material. Both models are subjected to rotational accelerations with magnitudes suspected to induce severe diffuse axonal injury. The results reveal that the presence of the ventricles leads to increased strain in the periventricular region, providing a plausible explanation for the vulnerability of the periventricular region. In addition, the strain-exacerbation effect associated with the presence of the ventricles is also noted in the paraventricular region, although less pronounced than that in the periventricular region. The current study advances the understanding of the periventricular injury mechanism as well as the detrimental effects that the ventricles exert on the periventricular and paraventricular brain tissue.

Displacement- and Strain-Based Discrimination of Head Injury Models across a Wide Range of Blunt Conditions

Successful validation of a head injury model is critical to ensure its biofidelity. However, there is an ongoing debate on what experimental data are suitable for model validation. Here, we report that CORrelation and Analysis (CORA) scores based on the commonly adopted relative brain-skull displacements or recent marker-based strains from cadaveric head impacts may not be effective in discriminating model-simulated whole-brain strains across a wide range of blunt conditions. We used three versions of the Worcester Head Injury Model (WHIM; isotropic and anisotropic WHIM V1.0, and anisotropic WHIM V1.5) to simulate 19 experiments, including eight high-rate cadaveric impacts, seven mid-rate cadaveric pure rotations simulating impacts in contact sports, and four in vivo head rotation/extension tests. All WHIMs achieved similar average CORA scores based on cadaveric displacement (~ 0.70; 0.47-0.88) and strain (V1.0: 0.86; 0.73-0.97 vs. V1.5: 0.78; 0.62-0.96), using the recommended settings. However, WHIM V1.5 produced ~ 1.17-2.69 times strain of the two V1.0 variants with substantial differences in strain distribution as well (Pearson correlation of ~ 0.57-0.92) when comparing their whole-brain strains across the range of blunt conditions. Importantly, their strain magnitude differences were similar to that in cadaveric marker-based strain (~ 1.32-3.79 times). This suggests that cadaveric strains are capable of discriminating head injury models for their simulated whole-brain strains (e.g., by using CORA magnitude sub-rating alone or peak strain magnitude ratio), although the aggregated CORA may not. This study may provide fresh insight into head injury model validation and the harmonization of simulation results from diverse head injury models. It may also facilitate future experimental designs to improve model validation.

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

Cerebral vasculature is several orders of magnitude stiffer than the brain tissue. However, only a handful of studies have investigated its potential stiffening effect on dynamic brain strains; yet, they report contradictory findings. Here, we reanalyze the cerebrovascular stiffening effect by incorporating vasculature derived from the latest neuroimaging atlases into a re-meshed Worcester Head Injury Model using an embedded element method. Regional brain strains with and without vasculature were simulated using a reconstructed, predominantly sagittal head impact. Using the two previously adopted linear or non-linear vessel material models, we reproduced the earlier conflicting results (~40% vs. ~1-6% in regional strain reductions). Nevertheless, with refitted non-linear material models chosen to represent the average dynamic tension behaviors of arteries and veins, respectively, inclusion of vasculature reduced regional brain strains by ~13-36% relative to the baselines without vasculature. Compared to the whole brain baseline response, inclusion of vasculature led to an element-wise linear regression slope of 0.8 and a Pearson correlation coefficient of 0.8. The vascular stiffening effect appears mild for the whole brain but more significant locally, which should not be ignored in head injury models. Nevertheless, more work is necessary to investigate the cerebrovascular mechanical behaviors and loading environment to allow for more biofidelic modeling of the brain in the future.

A network-based response feature matrix as a brain injury metric

Conventional brain injury metrics are scalars that treat the whole head/brain as a single unit but do not characterize the distribution of brain responses. Here, we establish a network-based "response feature matrix" to characterize the magnitude and distribution of impact-induced brain strains. The network nodes and edges encode injury risks to the gray matter regions and their white matter interconnections, respectively. The utility of the metric is illustrated in injury prediction using three independent, real-world datasets: two reconstructed impact datasets from the National Football League (NFL) and Virginia Tech, respectively, and measured concussive and non-injury impacts from Stanford University. Injury predictions with leave-one-out cross-validation are conducted using the two reconstructed datasets separately, and then by combining all datasets into one. Using support vector machine, the network-based injury predictor consistently outperforms four baseline scalar metrics including peak maximum principal strain of the whole brain (MPS), peak linear/rotational acceleration, and peak rotational velocity across all five selected performance measures (e.g., maximized accuracy of 0.887 vs. 0.774 and 0.849 for MPS and rotational acceleration with corresponding positive predictive values of 0.938, 0.772, and 0.800, respectively, using the reconstructed NFL dataset). With sufficient training data, real-world injury prediction is similar to leave-one-out in-sample evaluation, suggesting the potential advantage of the network-based injury metric over conventional scalar metrics. The network-based response feature matrix significantly extends scalar metrics by sampling the brain strains more completely, which may serve as a useful framework potentially allowing for other applications such as characterizing injury patterns or facilitating targeted multi-scale modeling in the future.

Development and Multi-Scale Validation of a Finite Element Football Helmet Model

Head injury is a growing concern within contact sports, including American football. Computational tools such as finite element (FE) models provide an avenue for researchers to study, and potentially optimize safety tools, such as helmets. The goal of this study was to develop an accurate representative helmet model that could be used in further study of head injury to mitigate the toll of concussions in contact sports. An FE model of a Schutt Air XP Pro football helmet was developed through three major steps: geometry development, material characterization, and model validation. The fully assembled helmet model was fit onto a Hybrid III dummy head–neck model and National Operating Committee on Standards for Athletic Equipment (NOCSAE) head model and validated through a series of 67 representative impacts similar to those experienced by a football player. The kinematic and kinetic response of the model was compared to the response of the physical experiments, which included force, head linear acceleration, head angular velocity, and carriage acceleration. The outputs between the model and the physical tests were quantitatively evaluated using CORelation and Analysis (CORA), amounting to an overall averaged score of 0.76. The model described in this study has been extensively validated and can function as a building block for innovation in player safety.