Performance Evaluation of a Pre-computed Brain Response Atlas in Dummy Head Impacts

Peaks and Distributions of White Matter Tract-related Strains in Bicycle Helmeted Impacts: Implication for Helmet Ranking and Optimization

Traumatic brain injury (TBI) in cyclists is a growing public health problem, with helmets being the major protection gear. Finite element head models have been increasingly used to engineer safer helmets often by mitigating brain strain peaks. However, how different helmets alter the spatial distribution of brain strain remains largely unknown. Besides, existing research primarily used maximum principal strain (MPS) as the injury parameter, while white matter fiber tract-related strains, increasingly recognized as effective predictors for TBI, have rarely been used for helmet evaluation. To address these research gaps, we used an anatomically detailed head model with embedded fiber tracts to simulate fifty-one helmeted impacts, encompassing seventeen bicycle helmets under three impact locations. We assessed the helmet performance based on four tract-related strains characterizing the normal and shear strain oriented along and perpendicular to the fiber tract, as well as the prevalently used MPS. Our results showed that both the helmet model and impact location affected the strain peaks. Interestingly, we noted that different helmets did not alter strain distribution, except for one helmet under one specific impact location. Moreover, our analyses revealed that helmet ranking outcome based on strain peaks was affected by the choice of injury metrics (Kendall's Tau coefficient: 0.58-0.93). Significant correlations were noted between tract-related strains and angular motion-based injury metrics. This study provided new insights into computational brain biomechanics and highlighted the helmet ranking outcome was dependent on the choice of injury metrics. Our results also hinted that the performance of helmets could be augmented by mitigating the strain peak and optimizing the strain distribution with accounting the selective vulnerability of brain subregions and more research was needed to develop region-specific injury criteria.

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.

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.

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 2 >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 2 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.

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.

Toward a Comprehensive Delineation of White Matter Tract-Related Deformation

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

Effective Head Impact Kinematics to Preserve Brain Strain

Conventional kinematics-based brain injury metrics often approximate peak maximum principal strain (MPS) of the whole brain but ignore the anatomical location of occurrence. In this study, we develop effective impact kinematics consisting of peak rotational velocity and the associated rotational axis to preserve not only peak MPS but also spatially detailed MPS. A pre-computed brain response atlas (pcBRA) serves as a common reference. A training dataset (N = 3069) is used to develop a convolutional neural network (CNN) to automate impact simplification. When preserving peak MPS alone, the CNN-estimated effective peak rotational velocity achieves a coefficient of determination ([Formula: see text]) of ~ 0.96 relative to the directly identified counterpart, far outperforming nominal peak velocity from the resultant profiles ([Formula: see text] of ~ 0.34). Impacts from a subset of data (N = 1900) are also successfully matched with pcBRA idealized impacts based on elementwise MPS, where their regression slope and Pearson correlation coefficient do not deviate from 1.0 (when identical) by more than 0.1. The CNN-estimated effective peak rotation velocity and rotational axis are sufficiently accurate for ~ 73.5% of the impacts. This is not possible for the nominal peak velocity or any other conventional injury metric. The performance may be further improved by expanding the pcBRA to include deceleration and focusing on region-wise strains. This study establishes a new avenue to reduce an arbitrary head impact into an idealized but actual "impact mode" characterized by triplets of basic kinematic variables. They retain specific physical interpretations of head impact and may be an advancement over state-of-the-art kinematics-based scalar metrics for more effective impact comparison in the future.

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.

Ranking and Rating Bicycle Helmet Safety Performance in Oblique Impacts Using Eight Different Brain Injury Models

Bicycle helmets are shown to offer protection against head injuries. Rating methods and test standards are used to evaluate different helmet designs and safety performance. Both strain-based injury criteria obtained from finite element brain injury models and metrics derived from global kinematic responses can be used to evaluate helmet safety performance. Little is known about how different injury models or injury metrics would rank and rate different helmets. The objective of this study was to determine how eight brain models and eight metrics based on global kinematics rank and rate a large number of bicycle helmets (n=17) subjected to oblique impacts. The results showed that the ranking and rating are influenced by the choice of model and metric. Kendall’s tau varied between 0.50 and 0.95 when the ranking was based on maximum principal strain from brain models. One specific helmet was rated as 2-star when using one brain model but as 4-star by another model. This could cause confusion for consumers rather than inform them of the relative safety performance of a helmet. Therefore, we suggest that the biomechanics community should create a norm or recommendation for future ranking and rating methods.

Low-Rank Representation of Head Impact Kinematics: A Data-Driven Emulator

Head motion induced by impacts has been deemed as one of the most important measures in brain injury prediction, given that the vast majority of brain injury metrics use head kinematics as input. Recently, researchers have focused on using fast approaches, such as machine learning, to approximate brain deformation in real time for early brain injury diagnosis. However, training such models requires large number of kinematic measurements, and therefore data augmentation is required given the limited on-field measured data available. In this study we present a principal component analysis-based method that emulates an empirical low-rank substitution for head impact kinematics, while requiring low computational cost. In characterizing our existing data set of 537 head impacts, each consisting of 6 degrees of freedom measurements, we found that only a few modes, e.g., 15 in the case of angular velocity, is sufficient for accurate reconstruction of the entire data set. Furthermore, these modes are predominantly low frequency since over 70% of the angular velocity response can be captured by modes that have frequencies under 40 Hz. We compared our proposed method against existing impact parametrization methods and showed significantly better performance in injury prediction using a range of kinematic-based metrics—such as head injury criterion (HIC), rotational injury criterion (RIC), and brain injury metric (BrIC)—and brain tissue deformation-based metrics—such as brain angle metric (BAM), maximum principal strain (MPS), and axonal fiber strains (FS). In all cases, our approach reproduced injury metrics similar to the ground truth measurements with no significant difference, whereas the existing methods obtained significantly different (p < 0.01) values as well as substantial differences in injury classification sensitivity and specificity. This emulator will enable us to provide the necessary data augmentation to build a head impact kinematic data set of any size. The emulator and corresponding examples are available on our website¹.

A Nonlinear Reduced-Order Model of the Corpus Callosum Under Planar Coronal Excitation

Traumatic brain injury (TBI) is often associated with microstructural tissue damage in the brain, which results from its complex biomechanical behavior. Recent studies have shown that the deep white matter (WM) region of the human brain is susceptible to being damaged due to strain localization in that region. Motivated by these studies, in this paper, we propose a geometrically nonlinear dynamical reduced order model (ROM) to model and study the dynamics of the deep WM region of the human brain under coronal excitation. In this model, the brain hemispheres were modeled as lumped masses connected via viscoelastic links, resembling the geometry of the corpus callosum (CC). Employing system identification techniques, we determined the unknown parameters of the ROM, and ensured the accuracy of the ROM by comparing its response against the response of an advanced finite element (FE) model. Next, utilizing modal analysis techniques, we determined the energy distribution among the governing modes of vibration of the ROM and concluded that the demonstrated nonlinear behavior of the FE model might be predominantly due to the special geometry of the brain deep WM region. Furthermore, we observed that, for sufficiently high input energies, high frequency harmonics at approximately 45 Hz, were generated in the response of the CC, which, in turn, are associated with high-frequency oscillations of the CC. Such harmonics might potentially lead to strain localization in the CC. This work is a step toward understanding the brain dynamics during traumatic injury.

A Multibody Model for Predicting Spatial Distribution of Human Brain Deformation Following Impact Loading

With an increasing focus on long-term consequences of concussive brain injuries, there is a new emphasis on developing tools that can accurately predict the mechanical response of the brain to impact loading. Although finite element models (FEM) estimate the brain response under dynamic loading, these models are not capable of delivering rapid (∼seconds) estimates of the brain's mechanical response. In this study, we develop a multibody spring-mass-damper model that estimates the regional motion of the brain to rotational accelerations delivered either about one anatomic axis or across three orthogonal axes simultaneously. In total, we estimated the deformation across 120 locations within a 50th percentile human brain. We found the multibody model (MBM) correlated, but did not precisely predict, the computed finite element response (average relative error: 18.4 ± 13.1%). We used machine learning (ML) to combine the prediction from the MBM and the loading kinematics (peak rotational acceleration, peak rotational velocity) and significantly reduced the discrepancy between the MBM and FEM (average relative error: 9.8 ± 7.7%). Using an independent sports injury testing set, we found the hybrid ML model also correlated well with predictions from a FEM (average relative error: 16.4 ± 10.2%). Finally, we used this hybrid MBM-ML approach to predict strains appearing in different locations throughout the brain, with average relative error estimates ranging from 8.6% to 25.2% for complex, multi-axial acceleration loading. Together, these results show a rapid and reasonably accurate method for predicting the mechanical response of the brain for single and multiplanar inputs, and provide a new tool for quickly assessing the consequences of impact loading throughout the brain.

Describing headform pose and impact location for blunt impact testing

Reproduction of anthropomorphic test device (ATD) head impact test methods is a critical element needed to develop guidance and technologies that reduce the risk for brain injury in sport. However, there does not appear to be a consensus for reporting ATD pose and impact location for industry and researchers to follow. Thus, the purpose of this article is to explore the various methods used to report impact location and ATD head pose for sport-related head impact testing and provide recommendations for standardizing these descriptions. A database search and exclusion process identified 137 articles that met the review criteria. Only 4 of the 137 articles provided a description similar to the method we propose to describe ATD pose and impact location. We thus propose a method to unambiguously convey the impact location and pose of the ATD based on the sequence, quantifiable design, and articulation of ATD mount joints. This reporting method has been used to a limited extent in the literature, but we assert that adoption of this method will help to standardize the reporting of ATD headform pose and impact location as well as aid in the replication of impact test protocols across laboratories.

A 6-year survey of road traffic accidents in Southwest China: Emphasis on traumatic brain injury

BACKGROUND

The objective of this study is to provide an up-to-date overview of the patterns of injuries, especially traumatic brain injury (TBI) caused by RTAs and to discuss some of the public health consequences.

METHODS

A scientific team was established to collect road traffic accidents occurring between 2013 and 2018 in Chongqing, Southwest China. For each accident, the environment-, vehicle-, and person- variables were analyzed and determined. The overall injury distribution and TBI patterns of four types of road users (driver, passenger, motorcyclist and pedestrian) were compared. The environmental and time distribution of accidents with TBI were shown by bar and pie chart. The risks of severe brain injury whether motorcyclist wearing helmets or not were compared and the risk factors of severe TBI in pedestrian were determined by odds ratio analysis.

RESULTS

This study enrolled 2131 accidents with 2741 persons of all kind of traffic participants, 1149 of them suffered AIS1+ head injury and 1598(58%) died in 7 days. The most common cause of deaths is due to head injury with 714(85%) and 1266(79%) persons died within 2 hours. Among 423 persons suffered both skull fracture and intracranial injury, 102 (24.1%) have an intracranial injury but no skull fractures, while none of the skull fractures without intracranial injury was found. Besides, motorcyclists without a helmet were at higher risks for all the brain injury categories. The risk of pedestrian suffering severe TBI at an impact speed of more than 70 km/h is 100 times higher than that with an impact speed of less than 40 km/h.

CONCLUSION

It is urgently needed to develop a more reliable brain injury evaluation criterion for better protection of the road users. We believe that strengthening the emergency care to head injury at the scene is the most effective way to reduce traffic fatality.

Convolutional neural network for efficient estimation of regional brain strains

Head injury models are important tools to study concussion biomechanics but are impractical for real-world use because they are too slow. Here, we develop a convolutional neural network (CNN) to estimate regional brain strains instantly and accurately by conceptualizing head rotational velocity profiles as two-dimensional images for input. We use two impact datasets with augmentation to investigate the CNN prediction performances with a variety of training-testing configurations. Three strain measures are considered, including maximum principal strain (MPS) of the whole brain, MPS of the corpus callosum, and fiber strain of the corpus callosum. The CNN is further tested using an independent impact dataset (N = 314) measured in American football. Based on 2592 training samples, it achieves a testing R² of 0.916 and root mean squared error (RMSE) of 0.014 for MPS of the whole brain. Combining all impact-strain response data available (N = 3069), the CNN achieves an R² of 0.966 and RMSE of 0.013 in a 10-fold cross-validation. This technique may enable a clinical diagnostic capability to a sophisticated head injury model, such as facilitating head impact sensors in concussion detection via a mobile device. In addition, it may transform current acceleration-based injury studies into focusing on regional brain strains. The trained CNN is publicly available along with associated code and examples at https://github.com/Jilab-biomechanics/CNN-brain-strains. They will be updated as needed in the future.

A reanalysis of football impact reconstructions for head kinematics and finite element modeling

BACKGROUND

Head kinematics generated by laboratory reconstructions of professional football helmet impacts have been applied to computational models to study the biomechanics of concussion. Since the original publication of this data, techniques for evaluating accelerometer consistency and error correction have been developed. This study applies these techniques to the original reconstruction data and reanalyzes the results given the current state of concussion biomechanics.

METHODS

Consistency checks were applied to the sensor data collected in the head of each test dummy. Inconsistent data were corrected using analytical techniques, and head kinematics were recalculated from the corrected data. Reconstruction videos were reviewed to identify artefactual impacts during the reconstruction to establish the region of applicability for simulations. Corrected head kinematics were input into finite element brain models to investigate strain response to the corrected dataset.

FINDINGS

Multiple reconstruction cases had inconsistent sensor arrays caused by a problematic sensor; corrections to the arrays caused changes in calculated rotational head motion. These corrections increased median peak angular velocity for the concussion cases from 35.6 to 41.5 rad/s. Using the original kinematics resulted in an average error of 20% in maximum principal strain results for each case. Simulations of the reconstructions also demonstrated that simulation lengths less than 40 ms did not capture the entire brain strain response and under-predicted strain.

INTERPRETATION

This study corrects data that were used to determine concussion risk, and indicates altered head angular motion and brain strain response for many reconstructions. Conclusions based on the original data should be re-examined based on this new study.

A review of impact testing methods for headgear in sports: Considerations for improved prevention of head injury through research and standards

Standards for sports headgear were introduced as far back as the 1960s and many have remained substantially unchanged to present day. Since this time, headgear has virtually eliminated catastrophic head injuries such as skull fractures and changed the landscape of head injuries in sports. Mild traumatic brain injury (mTBI) is now a prevalent concern and the effectiveness of headgear in mitigating mTBI is inconclusive for most sports. Given that most current headgear standards are confined to attenuating linear head mechanics and recent brain injury studies have underscored the importance of angular mechanics in the genesis of mTBI, new or expanded standards are needed to foster headgear development and assess headgear performance that addresses all types of sport-related head and brain injuries. The aim of this review is to provide a basis for developing new sports headgear impact tests for standards by summarizing and critiquing: 1) impact testing procedures currently codified in published headgear standards for sports and 2) new or proposed headgear impact test procedures in published literature and/or relevant conferences. Research areas identified as needing further knowledge to support standards test development include defining sports-specific head impact conditions, establishing injury and age appropriate headgear assessment criteria, and the development of headgear specific head and neck surrogates for at-risk populations.

Mesh Convergence Behavior and the Effect of Element Integration of a Human Head Injury Model

Numerous head injury models exist that vary in mesh density by orders of magnitude. A careful study of the mesh convergence behavior is necessary, especially in terms of strain most relevant to brain injury. To this end, as well as to investigate the effect of element integration scheme on simulated strains, we re-meshed the Worcester Head Injury Model at five mesh densities (~ 7.2-1000 k high-quality hexahedral elements of the brain). Results from explicit dynamic simulations of three cadaveric impacts and an in vivo head rotation were compared. First, scalar metrics of the whole brain only considering magnitude were used, including peak maximum principal strain and population-based median strain. They were further extended to deep white matter regions and the entire brain elements, respectively, to form two "response vectors" to account for both magnitude and distribution. Using benchmark enhanced full-integration elements (C3D8I), a minimum of 202.8 k brain elements were necessary to converge for response vectors of the deep white matter regions. This model was further used to simulate with reduced integration (C3D8R). We found that hourglass energy higher than the common rule of thumb (e.g., up to 44.38% vs. < 10% of internal energy) was necessary to maintain comparable strain relative to C3D8I. Based on these results, it is recommended that a human head injury model should have a minimum number of 202.8 k elements, or an average element size of no larger than 1.8 mm, for the brain. C3D8R formulation with relax stiffness hourglass control using a high scaling factor is also recommended to achieve sufficient accuracy without substantial computational cost.

Concussion classification via deep learning using whole-brain white matter fiber strains

Developing an accurate and reliable injury predictor is central to the biomechanical studies of traumatic brain injury. State-of-the-art efforts continue to rely on empirical, scalar metrics based on kinematics or model-estimated tissue responses explicitly pre-defined in a specific brain region of interest. They could suffer from loss of information. A single training dataset has also been used to evaluate performance but without cross-validation. In this study, we developed a deep learning approach for concussion classification using implicit features of the entire voxel-wise white matter fiber strains. Using reconstructed American National Football League (NFL) injury cases, leave-one-out cross-validation was employed to objectively compare injury prediction performances against two baseline machine learning classifiers (support vector machine (SVM) and random forest (RF)) and four scalar metrics via univariate logistic regression (Brain Injury Criterion (BrIC), cumulative strain damage measure of the whole brain (CSDM-WB) and the corpus callosum (CSDM-CC), and peak fiber strain in the CC). Feature-based machine learning classifiers including deep learning, SVM, and RF consistently outperformed all scalar injury metrics across all performance categories (e.g., leave-one-out accuracy of 0.828-0.862 vs. 0.690-0.776, and .632+ error of 0.148-0.176 vs. 0.207-0.292). Further, deep learning achieved the best cross-validation accuracy, sensitivity, AUC, and .632+ error. These findings demonstrate the superior performances of deep learning in concussion prediction and suggest its promise for future applications in biomechanical investigations of traumatic brain injury.

Material properties of the brain in injury-relevant conditions - Experiments and computational modeling

Material properties of the brain have been extensively studied but remain poorly characterized. The vast variations in constitutive models and material constants are well documented. However, no study exists to translate the variations into disparities in impact-induced brain strains most relevant to brain injury. Here, we reviewed a subset of injury-relevant brain material properties either characterized in experiments or adopted in recent head injury models. To highlight how variations in measured brain material properties manifested in simulated brain strains, we selected six experiments that have provided a complete set of brain material model and constants to implement a common head injury model. Responses resulting from two extreme events representing a high-rate cadaveric head impact and a low-rate in vivo head rotation, respectively, varied substantially. We hypothesized, and further confirmed, that the time-varying shear moduli at the appropriate time scales (e.g., ~5 ms and ~40 ms corresponding to the impulse durations of the major acceleration peaks for the two impacts, respectively), rather than the initial or long-term shear moduli, were the most indicative of impact-induced brain strains. These results underscored the need to implement measured brain material properties into an actual head injury model for evaluation. They may also provide guidelines to better characterize brain material properties in future experiments and head injury models. Finally, our finding provided a practical solution to satisfy head injury model validation requirements at both ends of the impact severity spectrum. This would improve the confidence in model simulation performance across a range of time scales relevant to concussion and sub-concussion in the real-world.

Propagation of errors from skull kinematic measurements to finite element tissue responses

Real-time quantification of head impacts using wearable sensors is an appealing approach to assess concussion risk. Traditionally, sensors were evaluated for accurately measuring peak resultant skull accelerations and velocities. With growing interest in utilizing model-estimated tissue responses for injury prediction, it is important to evaluate sensor accuracy in estimating tissue response as well. Here, we quantify how sensor kinematic measurement errors can propagate into tissue response errors. Using previous instrumented mouthguard validation datasets, we found that skull kinematic measurement errors in both magnitude and direction lead to errors in tissue response magnitude and distribution. For molar design instrumented mouthguards susceptible to mandible disturbances, 150-400% error in skull kinematic measurements resulted in 100% error in regional peak tissue response. With an improved incisor design mitigating mandible disturbances, errors in skull kinematics were reduced to <50%, and several tissue response errors were reduced to <10%. Applying 30[Formula: see text] rotations to reference kinematic signals to emulate sensor transformation errors yielded below 10% error in regional peak tissue response; however, up to 20% error was observed in peak tissue response for individual finite elements. These findings demonstrate that kinematic resultant errors result in regional peak tissue response errors, while kinematic directionality errors result in tissue response distribution errors. This highlights the need to account for both kinematic magnitude and direction errors and accurately determine transformations between sensors and the skull.