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

Efficient Generation of Pretraining Samples for Developing a Deep Learning Brain Injury Model via Transfer Learning

The large amount of training samples required to develop a deep learning brain injury model demands enormous computational resources. Here, we study how a transformer neural network (TNN) of high accuracy can be used to efficiently generate pretraining samples for a convolutional neural network (CNN) brain injury model to reduce computational cost. The samples use synthetic impacts emulating real-world events or augmented impacts generated from limited measured impacts. First, we verify that the TNN remains highly accurate for the two impact types (N = 100 each;R 2 of 0.948-0.967 with root mean squared error, RMSE, ~ 0.01, for voxelized peak strains). The TNN-estimated samples (1000-5000 for each data type) are then used to pretrain a CNN, which is further finetuned using directly simulated training samples (250-5000). An independent measured impact dataset considered of complete capture of impact event is used to assess estimation accuracy (N = 191). We find that pretraining can significantly improve CNN accuracy via transfer learning compared to a baseline CNN without pretraining. It is most effective when the finetuning dataset is relatively small (e.g., 2000-4000 pretraining synthetic or augmented samples improves success rate from 0.72 to 0.81 with 500 finetuning samples). When finetuning samples reach 3000 or more, no obvious improvement occurs from pretraining. These results support using the TNN to rapidly generate pretraining samples to facilitate a more efficient training strategy for future deep learning brain models, by limiting the number of costly direct simulations from an alternative baseline model. This study could contribute to a wider adoption of deep learning brain injury models for large-scale predictive modeling and ultimately, enhancing safety protocols and protective equipment.

Head Impact Kinematics and Brain Tissue Strains in High School Lacrosse

Male lacrosse and female lacrosse have differences in history, rules, and equipment. There is current debate regarding the need for enhanced protective headwear in female lacrosse like that worn by male lacrosse players. To inform this discussion, 17 high school lacrosse players (6 female and 11 male) wore the Stanford Instrumented Mouthguard during 26 competitive games over the 2021 season. Time-windowing and video review were used to remove false-positive recordings and verify head acceleration events (HAEs). The HAE rate in high school female lacrosse (0.21 per athlete exposure and 0.24 per player hour) was approximately 35% lower than the HAE rate in high school male lacrosse (0.33 per athlete exposure and 0.36 per player hour). Previously collected kinematics data from the 2019 high school male and female lacrosse season were combined with the newly collected 2021 kinematics data, which were used to drive a finite element head model and simulate 42 HAEs. Peak linear acceleration (PLA), peak angular velocity (PAV), and 95th percentile maximum principal strain (MPS95) of brain tissue were compared between HAEs in high school female and male lacrosse. Median values for peak kinematics and MPS95 of HAEs in high school female lacrosse (PLA, 22.3 g; PAV, 10.4 rad/s; MPS95, 0.05) were lower than for high school male lacrosse (PLA, 24.2 g; PAV, 15.4 rad/s; MPS95, 0.07), but the differences were not statistically significant. Quantifying a lower HAE rate in high school female lacrosse compared to high school male lacrosse, but similar HAE magnitudes, provides insight into the debate regarding helmets in female lacrosse. However, due to the small sample size, additional video-verified data from instrumented mouthguards are required.

An Objective Injury Threshold for the Maximum Principal Strain Criterion for Brain Tissue in the Finite Element Head Model and Its Application

Although the finite element head model (FEHM) has been widely utilized to analyze injury locations and patterns in traumatic brain injury, significant controversy persists regarding the selection of a mechanical injury variable and its corresponding threshold. This paper aims to determine an objective injury threshold for maximum principal strain (MPS) through a novel data-driven method, and to validate and apply it. We extract the peak responses from all elements across 100 head impact simulations to form a dataset, and then determine the objective injury threshold by analyzing the relationship between the combined injury degree and the threshold according to the stationary value principle. Using an occipital impact case from a clinical report as an example, we evaluate the accuracy of the injury prediction based on the new threshold. The results show that the injury area predicted by finite element analysis closely matches the main injury area observed in CT images, without the issue of over- or underestimating the injury due to an unreasonable threshold. Furthermore, by applying this threshold to the finite element analysis of designed occipital impacts, we observe, for the first time, supra-tentorium cerebelli injury, which is related to visual memory impairment. This discovery may indicate the biomechanical mechanism of visual memory impairment after occipital impacts reported in clinical cases.

The White Matter Fiber Tract Deforms Most in the Perpendicular Direction During In Vivo Volunteer Impacts

White matter (WM) tract-related strains are increasingly used to quantify brain mechanical responses, but their dynamics in live human brains during in vivo impact conditions remain largely unknown. Existing research primarily looked into the normal strain along the WM fiber tracts (i.e., tract-oriented normal strain), but it is rarely the case that the fiber tract only endures tract-oriented normal strain during impacts. In this study, we aim to extend the in vivo measurement of WM fiber deformation by quantifying the normal strain perpendicular to the fiber tract (i.e., tract-perpendicular normal strain) and the shear strain along and perpendicular to the fiber tract (i.e., tract-oriented shear strain and tract-perpendicular shear strain, respectively). To achieve this, we combine the three-dimensional strain tensor from the tagged magnetic resonance imaging with the diffusion tensor imaging (DTI) from an open-access dataset, including 44 volunteer impacts under two head loading modes, i.e., neck rotations (N = 30) and neck extensions (N = 14). The strain tensor is rotated to the coordinate system with one axis aligned with DTI-revealed fiber orientation, and then four tract-related strain measures are calculated. The results show that tract-perpendicular normal strain peaks are the largest among the four strain types (p < 0.05, Friedman's test). The distribution of tract-related strains is affected by the head loading mode, of which laterally symmetric patterns with respect to the midsagittal plane are noted under neck extensions, but not under neck rotations. Our study presents a comprehensive in vivo strain quantification toward a multifaceted understanding of WM dynamics. We find that the WM fiber tract deforms most in the perpendicular direction, illuminating new fundamentals of brain mechanics. The reported strain images can be used to evaluate the fidelity of computational head models, especially those intended to predict fiber deformation under noninjurious conditions.

Uncovering the hidden effects of repetitive subconcussive head impact exposure: A mega-analytic approach characterizing seasonal brain microstructural changes in contact and collision sports athletes

Repetitive subconcussive head impacts (RSHI) are believed to induce sub-clinical brain injuries, potentially resulting in cumulative, long-term brain alterations. This study explores patterns of longitudinal brain white matter changes across sports with RSHI-exposure. A systematic literature search identified 22 datasets with longitudinal diffusion magnetic resonance imaging data. Four datasets were centrally pooled to perform uniform quality control and data preprocessing. A total of 131 non-concussed active athletes (American football, rugby, ice hockey; mean age: 20.06 ± 2.06 years) with baseline and post-season data were included. Nonparametric permutation inference (one-sample t tests, one-sided) was applied to analyze the difference maps of multiple diffusion parameters. The analyses revealed widespread lateralized patterns of sports-season-related increases and decreases in mean diffusivity (MD), radial diffusivity (RD), and axial diffusivity (AD) across spatially distinct white matter regions. Increases were shown across one MD-cluster (3195 voxels; mean change: 2.34%), one AD-cluster (5740 voxels; mean change: 1.75%), and three RD-clusters (817 total voxels; mean change: 3.11 to 4.70%). Decreases were shown across two MD-clusters (1637 total voxels; mean change: -1.43 to -1.48%), two RD-clusters (1240 total voxels; mean change: -1.92 to -1.93%), and one AD-cluster (724 voxels; mean change: -1.28%). The resulting pattern implies the presence of strain-induced injuries in central and brainstem regions, with comparatively milder physical exercise-induced effects across frontal and superior regions of the left hemisphere, which need further investigation. This article highlights key considerations that need to be addressed in future work to enhance our understanding of the nature of observed white matter changes, improve the comparability of findings across studies, and promote data pooling initiatives to allow more detailed investigations (e.g., exploring sex- and sport-specific effects).

Full-field, frequency-domain comparison of simulated and measured human brain deformation

We propose a robust framework for quantitatively comparing model-predicted and experimentally measured strain fields in the human brain during harmonic skull motion. Traumatic brain injuries (TBIs) are typically caused by skull impact or acceleration, but how skull motion leads to brain deformation and consequent neural injury remains unclear and comparison of model predictions to experimental data remains limited. Magnetic resonance elastography (MRE) provides high-resolution, full-field measurements of dynamic brain deformation induced by harmonic skull motion. In the proposed framework, full-field strain measurements from human brain MRE in vivo are compared to simulated strain fields from models with similar harmonic loading. To enable comparison, the model geometry and subject anatomy, and subsequently, the predicted and measured strain fields are nonlinearly registered to the same standard brain atlas. Strain field correlations (\(\:{C}_{v}\)), both global (over the brain volume) and local (over smaller sub-volumes), are then computed from the inner product of the complex-valued strain tensors from model and experiment at each voxel. To demonstrate our approach, we compare strain fields from MRE in six human subjects to predictions from two previously developed models. Notably, global \(\:{C}_{v}\) values are higher when comparing strain fields from different subjects (\(\:{C}_{v}\)~0.6-0.7) than when comparing strain fields from either of the two models to strain fields in any subject. The proposed framework provides a quantitative method to assess similarity (and to identify discrepancies) between model predictions and experimental measurements of brain deformation, and thus can aid in the development and evaluation of improved models of brain biomechanics.

American Football On-Field Head Impact Kinematics: Influence of Acceleration Signal Characteristics on Peak Maximal Principal Strain

Recorded head kinematics from head-impact measurement devices (HIMd) are pivotal for evaluating brain stress and strain through head finite element models (hFEM). The variability in kinematic recording windows across HIMd presents challenges as they yield inconsistent hFEM responses. Despite establishing an ideal recording window for maximum principal strain (MPS) in brain tissue, uncertainties persist about the impact characteristics influencing vulnerability when this window is shortened. This study aimed to scrutinize factors within impact kinematics affecting the reliability of different recording windows on whole-brain peak MPS using a validated hFEM. Utilizing 53 on-field head impacts recorded via an instrumented mouthguard during a Canadian varsity football game, 10 recording windows were investigated with varying pre- and post-impact-trigger durations. Tukey pair-wise comparisons revealed no statistically significant differences in MPS responses for the different recording windows. However, specific impacts showed marked variability up to 40%. It was found, through correlation analyses, that impacts with lower peak linear acceleration exhibited greater response variability across different pre-trigger durations. Signal shape, analyzed through spectral analysis, influenced the time required for MPS development, resulting in specific impacts requiring a prolonged post-trigger duration. This study adds to the existing consensus on standardizing HIMd acquisition time windows and sheds light on impact characteristics leading to peak MPS variation across different head impact kinematic recording windows. Considering impact characteristics in research assessments is crucial, as certain impacts, affected by recording duration, may lead to significant errors in peak MPS responses during cumulative longitudinal exposure assessments.

Cavitation in blunt impact traumatic brain injury

Traumatic brain injury (TBI) poses a major public health challenge. No proven therapies for the condition exist so protective equipment that prevents or mitigates these injuries plays a critical role in minimizing the societal burden of this condition. Our ability to optimize protective equipment depends on our capacity to relate the mechanics of head impact events to morbidity and mortality. This capacity, in turn, depends on correctly identifying the mechanisms of injury. For several decades, a controversial theory of TBI biomechanics has attributed important classes of injury to cavitation inside the cranial vault during blunt impact. This theory explains counter-intuitive clinical observations, including the coup-contre-coup pattern of injury. However, it is also difficult to validate experimentally in living subjects. Also, blunt impact TBI is a broad term that covers a range of different head impact events, some of which may be better described by cavitation theory than others. This review surveys what has been learned about cavitation through mathematical modeling, physical modeling, and experimentation with living tissues and places it in context with competing theories of blunt injury biomechanics and recent research activity in the field in an attempt to understand what the theory has to offer the next generation of innovators in TBI biomechanics.

When to Pull the Trigger: Conceptual Considerations for Approximating Head Acceleration Events Using Instrumented Mouthguards

Head acceleration events (HAEs) are acceleration responses of the head following external short-duration collisions. The potential risk of brain injury from a single high-magnitude HAE or repeated occurrences makes them a significant concern in sport. Instrumented mouthguards (iMGs) can approximate HAEs. The distinction between sensor acceleration events, the iMG datum for approximating HAEs and HAEs themselves, which have been defined as the in vivo event, is made to highlight limitations of approximating HAEs using iMGs. This article explores the technical limitations of iMGs that constrain the approximation of HAEs and discusses important conceptual considerations for stakeholders interpreting iMG data. The approximation of HAEs by sensor acceleration events is constrained by false positives and false negatives. False positives occur when a sensor acceleration event is recorded despite no (in vivo) HAE occurring, while false negatives occur when a sensor acceleration event is not recorded after an (in vivo) HAE has occurred. Various mechanisms contribute to false positives and false negatives. Video verification and post-processing algorithms offer effective means for eradicating most false positives, but mitigation for false negatives is less comprehensive. Consequently, current iMG research is likely to underestimate HAE exposures, especially at lower magnitudes. Future research should aim to mitigate false negatives, while current iMG datasets should be interpreted with consideration for false negatives when inferring athlete HAE exposure.

Brain Deformation Estimation With Transfer Learning for Head Impact Datasets Across Impact Types

OBJECTIVE

The machine-learning head model (MLHM) to accelerate the calculation of brain strain and strain rate, which are the predictors for traumatic brain injury (TBI), but the model accuracy was found to decrease sharply when the training/test datasets were from different head impacts types (i.e., car crash, college football), which limits the applicability of MLHMs to different types of head impacts and sports. Particularly, small sizes of target dataset for specific impact types with tens of impacts may not be enough to train an accurate impact-type-specific MLHM.

METHODS

To overcome this, we propose data fusion and transfer learning to develop a series of MLHMs to predict the maximum principal strain (MPS) and maximum principal strain rate (MPSR).

RESULTS

The strategies were tested on American football (338), mixed martial arts (457), reconstructed car crash (48) and reconstructed American football (36) and we found that the MLHMs developed with transfer learning are significantly more accurate in estimating MPS and MPSR than other models, with a mean absolute error (MAE) smaller than 0.03 in predicting MPS and smaller than [Formula: see text] in predicting MPSR on all target impact datasets. High performance in concussion detection was observed based on the MPS and MPSR estimated by the transfer-learning-based models.

CONCLUSION

The MLHMs can be applied to various head impact types for rapidly and accurately calculating brain strain and strain rate.

SIGNIFICANCE

This study enables developing MLHMs for the head impact type with limited availability of data, and will accelerate the applications of MLHMs.

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.

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.

Using dropout based active learning and surrogate models in the inverse viscoelastic parameter identification of human brain tissue

Inverse mechanical parameter identification enables the characterization of ultrasoft materials, for which it is difficult to achieve homogeneous deformation states. However, this usually involves high computational costs that are mainly determined by the complexity of the forward model. While simulation methods like finite element models can capture nearly arbitrary geometries and implement involved constitutive equations, they are also computationally expensive. Machine learning models, such as neural networks, can help mitigate this problem when they are used as surrogate models replacing the complex high fidelity models. Thereby, they serve as a reduced order model after an initial training phase, where they learn the relation of in- and outputs of the high fidelity model. The generation of the required training data is computationally expensive due to the necessary simulation runs. Here, active learning techniques enable the selection of the "most rewarding" training points in terms of estimated gained accuracy for the trained model. In this work, we present a recurrent neural network that can well approximate the output of a viscoelastic finite element simulation while significantly speeding up the evaluation times. Additionally, we use Monte-Carlo dropout based active learning to identify highly informative training data. Finally, we showcase the potential of the developed pipeline by identifying viscoelastic material parameters for human brain tissue.

Inverse identification of region-specific hyperelastic material parameters for human brain tissue

The identification of material parameters accurately describing the region-dependent mechanical behavior of human brain tissue is crucial for computational models used to assist, e.g., the development of safety equipment like helmets or the planning and execution of brain surgery. While the division of the human brain into different anatomical regions is well established, knowledge about regions with distinct mechanical properties remains limited. Here, we establish an inverse parameter identification scheme using a hyperelastic Ogden model and experimental data from multi-modal testing of tissue from 19 anatomical human brain regions to identify mechanically distinct regions and provide the corresponding material parameters. We assign the 19 anatomical regions to nine governing regions based on similar parameters and microstructures. Statistical analyses confirm differences between the regions and indicate that at least the corpus callosum and the corona radiata should be assigned different material parameters in computational models of the human brain. We provide a total of four parameter sets based on the two initial Poisson's ratios of 0.45 and 0.49 as well as the pre- and unconditioned experimental responses, respectively. Our results highlight the close interrelation between the Poisson's ratio and the remaining model parameters. The identified parameters will contribute to more precise computational models enabling spatially resolved predictions of the stress and strain states in human brains under complex mechanical loading conditions.

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.

Sex Differences in Axonal Dynamic Responses Under Realistic Tension Using Finite Element Models

Existing axonal finite element models do not consider sex morphological differences or the fidelity in dynamic input. To facilitate a systematic investigation into the micromechanics of diffuse axonal injury, we develop a parameterized modeling approach for automatic and efficient generation of sex-specific axonal models according to specified geometrical parameters. Baseline female and male axonal models in the corpus callosum with random microtubule (MT) gap configurations are generated for model calibration and evaluation. They are then used to simulate a realistic tensile loading consisting of both a loading and a recovery phase (to return to an initial undeformed state) generated from dynamic corpus callosum fiber strain in a real-world head impact simulation. We find that MT gaps and the dynamic recovery phase are both critical to successfully reproduce MT undulation as observed experimentally, which has not been reported before. This strengthens confidence in model dynamic responses. A statistical approach is further employed to aggregate axonal responses from a large sample of random MT gap configurations for both female and male axonal models (n = 10,000 each). We find that peak strains in MTs and the Ranvier node and associated neurofilament failures in female axons are substantially higher than those in male axons because there are fewer MTs in the former and also because of the random nature of MT gap locations. Despite limitations in various model assumptions as a result of limited experimental data currently available, these findings highlight the need to systematically characterize MT gap configurations and to ensure a realistic model input for axonal dynamic simulations. Finally, this study may offer fresh and improved insight into the biomechanical basis of sex differences in brain injury, and sets the stage for more systematic investigations at the microscale in the future, both numerically and experimentally.

Magnetic Resonance Elastography-Based Technique to Assess the Biomechanics of the Skull-Brain Interface: Repeatability and Age-Sex Characteristics

Increasing concerns have been raised about the long-term negative effects of subconcussive repeated head impact (RHI). To elucidate RHI injury mechanisms, many efforts have studied how head impacts affect the skull-brain biomechanics and have found that mechanical interactions at the skull-brain interface dampen and isolate brain motions by decoupling the brain from the skull. Despite intense interest, in vivo quantification of the functional state of the skull-brain interface remains difficult. This study developed a magnetic resonance elastography (MRE) based technique to non-invasively assess skull-brain mechanical interactions (i.e., motion transmission and isolation function) under dynamic loading. The full MRE displacement data were separated into rigid body motion and wave motion. The rigid body motion was used to calculate the brain-to-skull rotational motion transmission ratio (Rtr) to quantify skull-brain motion transmissibility, and the wave motion was used to calculate the cortical normalized octahedral shear strain (NOSS) (calculated based on a partial derivative computing neural network) to evaluate the isolation capability of the skull-brain interface. Forty-seven healthy volunteers were recruited to investigate the effects of age/sex on Rtr and cortical NOSS, and 17 of 47 volunteers received multiple scans to test the repeatability of the proposed techniques under different strain conditions. The results showed that both Rtr and NOSS were robust to MRE driver variations and had good repeatability, with intraclass correlation coefficient (ICC) values between 0.68 and 0.97 (fair to excellent). No age or sex dependence were observed with Rtr, whereas a significant positive correlation between age and NOSS was found in the cerebrum, frontal, temporal, and parietal lobes (all p < 0.05), but not in the occipital lobe (p = 0.99). The greatest change in NOSS with age was found in the frontal lobe, one of the most frequent locations of traumatic brain injury (TBI). Except for the temporal lobe (p = 0.0087), there was no significant difference in NOSS between men and women. This work provides motivation for utilizing MRE as a non-invasive tool for quantifying the biomechanics of the skull-brain interface. It evaluated the age and sex dependence and may lead to a better understanding of the protective role and mechanisms of the skull-brain interface in RHI and TBI, as well as improve the accuracy of computational models in simulating the skull-brain interface.

An Instrumented Mouthguard for Real-Time Measurement of Head Kinematics under a Large Range of Sport Specific Accelerations

BACKGROUND

Head impacts in sports can produce brain injuries. The accurate quantification of head kinematics through instrumented mouthguards (iMG) can help identify underlying brain motion during injurious impacts. The aim of the current study is to assess the validity of an iMG across a large range of linear and rotational accelerations to allow for on-field head impact monitoring.

METHODS

Drop tests of an instrumented helmeted anthropometric testing device (ATD) were performed across a range of impact magnitudes and locations, with iMG measures collected concurrently. ATD and iMG kinematics were also fed forward to high-fidelity brain models to predict maximal principal strain.

RESULTS

The impacts produced a wide range of head kinematics (16-171 g, 1330-10,164 rad/s2 and 11.3-41.5 rad/s) and durations (6-18 ms), representing impacts in rugby and boxing. Comparison of the peak values across ATD and iMG indicated high levels of agreement, with a total concordance correlation coefficient of 0.97 for peak impact kinematics and 0.97 for predicted brain strain. We also found good agreement between iMG and ATD measured time-series kinematic data, with the highest normalized root mean squared error for rotational velocity (5.47 ± 2.61%) and the lowest for rotational acceleration (1.24 ± 0.86%). Our results confirm that the iMG can reliably measure laboratory-based head kinematics under a large range of accelerations and is suitable for future on-field validity assessments.

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.

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.