Development of a Single-Degree-of-Freedom Mechanical Model for Predicting Strain-Based Brain Injury Responses

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 study on the application of diffuse axonal multi-axis general evaluation for brain injury assessment in small overlap barrier crash test

PURPOSE

Head injury criterion (HIC) companied by a rotation-based metric was widely believed to be helpful for head injury prediction in road traffic accidents. Recently, the Euro-New Car Assessment Program utilized a newly developed metric called diffuse axonal multi-axis general evaluation (DAMAGE) to explain test device for human occupant restraint (THOR) head injury, which demonstrated excellent ability in capturing concussions and diffuse axonal injuries. However, there is still a lack of comprehensive understanding regarding the effectiveness of using DAMAGE for Hybrid Ⅲ 50th percentile male dummy (H50th) head injury assessment. The objective of this study is to determine whether the DAMAGE could capture the risk of H50th brain injury during small overlap barrier tests.

METHODS

To achieve this objective, a total of 24 vehicle crash loading curves were collected as input data for the multi-body simulation. Two commercially available mathematical dynamic models, namely H50th and THOR, were utilized to investigate the differences in head injury response. Subsequently, a decision method known as simple additive weighting was employed to establish a comprehensive brain injury metric by incorporating the weighted HIC and either DAMAGE or brain injury criterion. Furthermore, 35 sets of vehicle crash test data were used to analyze these brain injury metrics.

RESULTS

The rotational displacement of the THOR head is significantly greater than that of the H50th head. The maximum linear and rotational head accelerations experienced by H50th and THOR models were (544.6 ± 341.7) m/s2, (2468.2 ± 1309.4) rad/s2 and (715.2 ± 332.8) m/s2, (3778.7 ± 1660.6) rad/s2, respectively. Under the same loading condition during small overlap barrier (SOB) tests, THOR exhibits a higher risk of head injury compared to the H50th model. It was observed that the overall head injury response during the small overlap left test condition is greater than that during the small overlap right test. Additionally, an equation was formulated to establish the necessary relationship between the DAMAGE values of THOR and H50th.

CONCLUSION

If H50th rather than THOR is employed as an evaluation tool in SOB crash tests, newly designed vehicles are more likely to achieve superior performance scores. According to the current injury curve for DAMAGE and brain injury criterion, it is highly recommended that HIC along with DAMAGE was prioritized for brain injury assessment in SOB tests.

Characteristics of head frequency response in blunt impacts: a biomechanical modeling study

Existing evaluation criteria for head impact injuries are typically based on time-domain features, and less attention has been paid to head frequency responses for head impact injury assessment. The purpose of the current study is, therefore, to understand the characteristics of human body head frequency response in blunt impacts via finite element (FE) modeling and the wavelet packet analysis method. FE simulation results show that head frequency response in blunt impacts could be affected by the impact boundary condition. The head energy peak and its frequency increase with the increase in impact; a stiffer impact block is associated with a higher head energy peak, and a bigger impact block could result in a high proportion of the energy peak. Regression analysis indicates that only the head energy peak has a high correlation with exiting head injury criteria, which implies that the amplitude-frequency aggregation characteristic but not the frequency itself of the head acceleration response has predictability for head impact injury in blunt impacts. The findings of the current study may provide additional criteria for head impact injury evaluation and new ideas for head impact injury protection.

Strain Response of an Anatomically Accurate Nonhuman Primate Finite Element Brain Model Under Sagittal Loading

INTRODUCTION

Prevention and treatment of traumatic brain injuries is critical to preserving soldier brain health. Laboratory studies are commonly used to reproduce injuries, understand injury mechanisms, and develop tolerance limits; however, this approach has limitations for studying brain injury, which requires a physiological response. The nonhuman primate (NHP) has been used as an effective model for investigating brain injury for many years. Prior research using the NHP provides a valuable resource to leverage using modern analysis and modeling techniques to improve our understanding of brain injury. The objectives of the present study are to develop an anatomically accurate finite element model of the NHP and determine regional brain responses using previously collected NHP data.

MATERIALS AND METHODS

The finite element model was developed using a neuroimaging-based anatomical atlas of the rhesus macaque that includes both cortical and subcortical structures. Head kinematic data from 10 sagittal NHP experiments, four +Gx (rearward) and six -Gx (frontal), were used to test model stability and obtain brain strain responses from multiple severities and vectors.

RESULTS

For +Gx tests, the whole-brain cumulative strain damage measure exceeding a strain threshold of 0.15 (CSDM15) ranged from 0.28 to 0.89, and 95th percentile of the whole-brain maximum principal strain (MPS95) ranged from 0.21 to 0.59. For -Gx tests, whole-brain CSDM15 ranged from 0.02 to 0.66, and whole-brain MPS95 ranged from 0.08 to 0.39.

CONCLUSIONS

Recognizing that NHPs are the closest surrogate to humans combined with the limitations of conducting brain injury research in the laboratory, a detailed anatomically accurate finite element model of an NHP was developed and exercised using previously collected data from the Naval Biodynamics Laboratory. The presently developed model can be used to conduct additional analyses to act as pilot data for the design of newer experiments with statistical power because of the sensitivity and resources needed to conduct experiments with NHPs.

Machine-learning-based head impact subtyping based on the spectral densities of the measurable head kinematics

BACKGROUND

Traumatic brain injury can be caused by head impacts, but many brain injury risk estimation models are not equally accurate across the variety of impacts that patients may undergo, and the characteristics of different types of impacts are not well studied. We investigated the spectral characteristics of different head impact types with kinematics classification.

METHODS

Data were analyzed from 3262 head impacts from lab reconstruction, American football, mixed martial arts, and publicly available car crash data. A random forest classifier with spectral densities of linear acceleration and angular velocity was built to classify head impact types (e.g., football, car crash, mixed martial arts). To test the classifier robustness, another 271 lab-reconstructed impacts were obtained from 5 other instrumented mouthguards. Finally, with the classifier, type-specific, nearest-neighbor regression models were built for brain strain.

RESULTS

The classifier reached a median accuracy of 96% over 1000 random partitions of training and test sets. The most important features in the classification included both low- and high-frequency features, both linear acceleration features and angular velocity features. Different head impact types had different distributions of spectral densities in low- and high-frequency ranges (e.g., the spectral densities of mixed martial arts impacts were higher in the high-frequency range than in the low-frequency range). The type-specific regression showed a generally higher R2 value than baseline models without classification.

CONCLUSION

The machine-learning-based classifier enables a better understanding of the impact kinematics spectral density in different sports, and it can be applied to evaluate the quality of impact-simulation systems and on-field data augmentation.

Consensus Head Acceleration Measurement Practices (CHAMP): Laboratory Validation of Wearable Head Kinematic Devices

Wearable devices are increasingly used to measure real-world head impacts and study brain injury mechanisms. These devices must undergo validation testing to ensure they provide reliable and accurate information for head impact sensing, and controlled laboratory testing should be the first step of validation. Past validation studies have applied varying methodologies, and some devices have been deployed for on-field use without validation. This paper presents best practices recommendations for validating wearable head kinematic devices in the laboratory, with the goal of standardizing validation test methods and data reporting. Key considerations, recommended approaches, and specific considerations were developed for four main aspects of laboratory validation, including surrogate selection, test conditions, data collection, and data analysis. Recommendations were generated by a group with expertise in head kinematic sensing and laboratory validation methods and reviewed by a larger group to achieve consensus on best practices. We recommend that these best practices are followed by manufacturers, users, and reviewers to conduct and/or review laboratory validation of wearable devices, which is a minimum initial step prior to on-field validation and deployment. We anticipate that the best practices recommendations will lead to more rigorous validation of wearable head kinematic devices and higher accuracy in head impact data, which can subsequently advance brain injury research and management.

Piecewise Multivariate Linearity Between Kinematic Features and Cumulative Strain Damage Measure (CSDM) Across Different Types of Head Impacts

In a previous study, we found that the relationship between brain strain and kinematic features cannot be described by a generalized linear model across different types of head impacts. In this study, we investigate if such a linear relationship exists when partitioning head impacts using a data-driven approach. We applied the K-means clustering method to partition 3161 impacts from various sources including simulation, college football, mixed martial arts, and car crashes. We found piecewise multivariate linearity between the cumulative strain damage (CSDM; assessed at the threshold of 0.15) and head kinematic features. Compared with the linear regression models without partition and the partition according to the types of head impacts, K-means-based data-driven partition showed significantly higher CSDM regression accuracy, which suggested the presence of piecewise multivariate linearity across types of head impacts. Additionally, we compared the piecewise linearity with the partitions based on individual features used in clustering. We found that the partition with maximum angular acceleration magnitude at 4706 rad/s² led to the highest piecewise linearity. This study may contribute to an improved method for the rapid prediction of CSDM in the future.

Concussion biomechanics, head acceleration exposure and brain injury criteria in sport: a review

There are mounting concerns surrounding the risk of neurodegenerative diseases and complications associated with concussion incidence and repetitive head acceleration events (HAE) in sport. The aim of this review is to provide an overview of concussion biomechanics, head acceleration exposure and brain injury criteria in sport. Rotational head motion appears to be the primary contributor to brain injury risk due to the unique mechanical properties of the brain and its location within the body. There is a growing evidence base of different biomechanical brain injury mechanisms, including those involving repetitive HAE. Historically, many studies on concussion biomechanics, head acceleration exposure and brain injury criteria in sport have been limited by validity of the biomechanical approaches undertaken. Biomechanical approaches such as instrumented mouthguards and subject-specific finite element (FE) brain models provide a unique opportunity to develop greater brain injury criteria and aid in on-field athlete removal. Implementing these approaches on a large-scale can gain insight into potential risk factors within sports and certain athletes/cohorts who sustain a greater number and/or severity of HAE throughout their playing career. These findings could play a key role in the development of concussion prevention strategies and techniques that mitigate the severity of HAE in sport.

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.

Predictive Factors of Kinematics in Traumatic Brain Injury from Head Impacts Based on Statistical Interpretation

Brain tissue deformation resulting from head impacts is primarily caused by rotation and can lead to traumatic brain injury. To quantify brain injury risk based on measurements of kinematics on the head, finite element (FE) models and various brain injury criteria based on different factors of these kinematics have been developed, but the contribution of different kinematic factors has not been comprehensively analyzed across different types of head impacts in a data-driven manner. To better design brain injury criteria, the predictive power of rotational kinematics factors, which are different in (1) the derivative order (angular velocity, angular acceleration, angular jerk), (2) the direction and (3) the power (e.g., square-rooted, squared, cubic) of the angular velocity, were analyzed based on different datasets including laboratory impacts, American football, mixed martial arts (MMA), NHTSA automobile crashworthiness tests and NASCAR crash events. Ordinary least squares regressions were built from kinematics factors to the 95% maximum principal strain (MPS95), and we compared zero-order correlation coefficients, structure coefficients, commonality analysis, and dominance analysis. The angular acceleration, the magnitude and the first power factors showed the highest predictive power for the majority of impacts including laboratory impacts, American football impacts, with few exceptions (angular velocity for MMA and NASCAR impacts). The predictive power of rotational kinematics about three directions (x: posterior-to-anterior, y: left-to-right, z: superior-to-inferior) of kinematics varied with different sports and types of head impacts.

Time Window of Head Impact Kinematics Measurement for Calculation of Brain Strain and Strain Rate in American Football

Wearable devices have been shown to effectively measure the head's movement during impacts in sports like American football. When a head impact occurs, the device is triggered to collect and save the kinematic measurements during a predefined time window. Then, based on the collected kinematics, finite element (FE) head models can calculate brain strain and strain rate, which are used to evaluate the risk of mild traumatic brain injury. To find a time window that can provide a sufficient duration of kinematics for FE analysis, we investigated 118 on-field video-confirmed football head impacts collected by the Stanford Instrumented Mouthguard. The simulation results based on the kinematics truncated to a shorter time window were compared with the original to determine the minimum time window needed for football. Because the individual differences in brain geometry influence these calculations, we included six representative brain geometries and found that larger brains need a longer time window of kinematics for accurate calculation. Among the different sizes of brains, a pre-trigger time of 40 ms and a post-trigger time of 70 ms were found to yield calculations of brain strain and strain rate that were not significantly different from calculations using the original 200 ms time window recorded by the mouthguard. Therefore, approximately 110 ms is recommended for complete modeling of impacts for football.

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.

The relationship between brain injury criteria and brain strain across different types of head impacts can be different

Multiple brain injury criteria (BIC) are developed to quickly quantify brain injury risks after head impacts. These BIC originated from different head impact types (e.g. sports and car crashes) are widely used in risk evaluation. However, the accuracy of using the BIC on brain injury risk estimation across head impact types has not been evaluated. Physiologically, brain strain is often considered the key parameter of brain injury. To evaluate the BIC's risk estimation accuracy across five datasets comprising different head impact types, linear regression was used to model 95% maximum principal strain, 95% maximum principal strain at the corpus callosum and cumulative strain damage (15%) on 18 BIC. The results show significantly different relationships between BIC and brain strain across datasets, indicating the same BIC value may suggest different brain strain across head impact types. The accuracy of brain strain regression is generally decreasing if the BIC regression models are fitted on a dataset with a different type of head impact rather than on the dataset with the same type. Given this finding, this study raises concerns for applying BIC to estimate the brain injury risks for head impacts different from the head impacts on which the BIC was developed.

Rapid Estimation of Entire Brain Strain Using Deep Learning Models

OBJECTIVE

Many recent studies have suggested that brain deformation resulting from a head impact is linked to the corresponding clinical outcome, such as mild traumatic brain injury (mTBI). Even though several finite element (FE) head models have been developed and validated to calculate brain deformation based on impact kinematics, the clinical application of these FE head models is limited due to the time-consuming nature of FE simulations. This work aims to accelerate the process of brain deformation calculation and thus improve the potential for clinical applications.

METHODS

We propose a deep learning head model with a five-layer deep neural network and feature engineering, and trained and tested the model on 2511 total head impacts from a combination of head model simulations and on-field college football and mixed martial arts impacts.

RESULTS

The proposed deep learning head model can calculate the maximum principal strain (Green Lagrange) for every element in the entire brain in less than 0.001s with an average root mean squared error of 0.022, and with a standard deviation of 0.001 over twenty repeats with random data partition and model initialization.

CONCLUSION

Trained and tested using the dataset of 2511 head impacts, this model can be applied to various sports in the calculation of brain strain with accuracy, and its applicability can even further be extended by incorporating data from other types of head impacts.

SIGNIFICANCE

In addition to the potential clinical application in real-time brain deformation monitoring, this model will help researchers estimate the brain strain from a large number of head impacts more efficiently than using FE models.

A Computational Study of Liquid Shock Absorption for Prevention of Traumatic Brain Injury

Mild traumatic brain injury (mTBI), more colloquially known as concussion, is common in contact sports such as American football, leading to increased scrutiny of head protective gear. Standardized laboratory impact testing, such as the yearly National Football League (NFL) helmet test, is used to rank the protective performance of football helmets, motivating new technologies to improve the safety of helmets relative to existing equipment. In this work, we hypothesized that a helmet which transmits a nearly constant minimum force will result in a reduced risk of mTBI. To evaluate the plausibility of this hypothesis, we first show that the optimal force transmitted to the head, in a reduced order model of the brain, is in fact a constant force profile. To simulate the effects of a constant force within a helmet, we conceptualize a fluid-based shock absorber system for use within a football helmet. We integrate this system within a computational helmet model and simulate its performance on the standard NFL helmet test impact conditions. The simulated helmet is compared with other helmet designs with different technologies. Computer simulations of head impacts with liquid shock absorption predict that, at the highest impact speed (9.3 m/s), the average brain tissue strain is reduced by 27.6% ± 9.3 compared to existing helmet padding when tested on the NFL helmet protocol. This simulation-based study puts forth a target benchmark for the future design of physical manifestations of this technology.

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

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

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

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.

Effect of foam densification and impact velocity on the performance of a football helmet using computational modeling

The NFL recently released validated helmet-impact models to study the performance of currently used helmets. This study used the model of a Riddell Speed Classic helmet to determine the influence of the properties of protective foam padding on acceleration and deformation at two common impact locations to cause concussions. The performance of the helmet was measured before and after manipulating the material properties of the protective foam liner material using FEA software. The densification strain was adjusted by using the scale factor tool in LS-DYNA to create four material categories - soft, standard, stiff, and rigid. The helmet was tested under side and rear impacts using the four material properties at 2.0, 5.5, 7.4, 9.3 and 12.3 m/s impact speeds using the NOCSAE linear impactor model. This study suggests that the standard foam material compresses to a range that could be considered to have "bottomed out" at impact speeds at 5.5 m/s for side impacts. Despite testing a wide range of material properties, the measured accelerations did not vary dramatically across material properties. Rather, impact speed played the dominant role on measured acceleration. This is the first study to demonstrate how open-source impact models can be used to run a design of experiments and investigate the role between different materials used inside a helmet and football helmet performance.

Mechanisms and variances of rotation-induced brain injury: a parametric investigation between head kinematics and brain strain

There lacks a comprehensive understanding of the correlation between head kinematics and brain strain especially deep-brain strain, partially resulting the deficiency of understanding brain injury mechanisms and the difficulty of choosing appropriate brain injury metrics. Hence, we simulated 76 impacts that were focused on concussion-relevant rotational kinematics and evaluated cumulative strain damage measure (CSDM) and average strain that could represent brain strain distribution. For the whole brain, axial rotation induced the highest CSDM, while lateral bending produced the lowest CSDM. However, for the deep-brain components, lateral bending produced the highest CSDM to the corpus callosum and thalamus. We further confirmed that brain strain was mainly produced by rotational kinematics, for which the effect of rotational deceleration could not be ignored with the deceleration influencing CSDM20 up to 27%. Our data supported that peak rotational velocity correlated to brain strain with an average R² of 0.77 across various impact directions and different shapes of loading curves. The correlation between peak rotational velocity and brain strain reached to an average R² of 0.99 for each specific impact direction. Our results supported using direction-specific peak rotation velocity for predicting strain-related brain injury. Additionally, we highlighted the importance of investigating whole-brain and deep-brain strain, as well as considering rotational deceleration.

Multi-Directional Dynamic Model for Traumatic Brain Injury Detection

Given the worldwide adverse impact of traumatic brain injury (TBI) on the human population, its diagnosis and prediction are of utmost importance. Historically, many studies have focused on associating head kinematics to brain injury risk. Recently, there has been a push toward using computationally expensive finite element (FE) models of the brain to create tissue deformation metrics of brain injury. Here, we develop a new brain injury metric, the brain angle metric (BAM), based on the dynamics of a 3 degree-of-freedom lumped parameter brain model. The brain model is built based on the measured natural frequencies of an FE brain model simulated with live human impact data. We show that it can be used to rapidly estimate peak brain strains experienced during head rotational accelerations that cause mild TBI. In our data set, the simplified model correlates with peak principal FE strain (R2 = 0.82). Further, coronal and axial brain model displacement correlated with fiber-oriented peak strain in the corpus callosum (R2 = 0.77). Our proposed injury metric BAM uses the maximum angle predicted by our brain model and is compared against a number of existing rotational and translational kinematic injury metrics on a data set of head kinematics from 27 clinically diagnosed injuries and 887 non-injuries. We found that BAM performed comparably to peak angular acceleration, translational acceleration, and angular velocity in classifying injury and non-injury events. Metrics that separated time traces into their directional components had improved model deviance compare with those that combined components into a single time trace magnitude. Our brain model can be used in future work to rapidly approximate the peak strain resulting from mild to moderate head impacts and to quickly assess brain injury risk.

Head Rotational Kinematics, Tissue Deformations, and Their Relationships to the Acute Traumatic Axonal Injury

Head rotational kinematics and tissue deformation metrics obtained from finite element models (FEM) have the potential to be used as traumatic axonal injury (TAI) assessment criteria and headgear evaluation standards. These metrics have been used to predict the likelihood of TAI occurrence; however, their ability in the assessment of the extent of TAI has not been explored. In this study, a pig model of TAI was used to examine a wide range of head loading conditions in two directions. The extent of TAI was quantified through histopathology and correlated to the FEM-derived tissue deformations and the head rotational kinematics. Peak angular acceleration and maximum strain rate of axonal fiber and brain tissue showed relatively good correlation to the volume of axonal injury, with similar correlation trends for both directions separately or combined. These rotational kinematics and tissue deformations can estimate the extent of acute TAI. The relationships between the head kinematics and the tissue strain, strain rate, and strain times strain rate were determined over the experimental range examined herein, and beyond that through parametric simulations. These relationships demonstrate that peak angular velocity and acceleration affect the underlying tissue deformations and the knowledge of both help to predict TAI risk. These relationships were combined with the injury thresholds, extracted from the TAI risk curves, and the kinematic-based risk curves representing overall axonal and brain tissue strain and strain rate were determined for predicting TAI. After scaling to humans, these curves can be used for real-time TAI assessment.

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.

An Injury Risk Function for the Leg, Foot, and Ankle Exposed to Axial Impact Loading Using Force and Impulse

Most injury risk functions (IRFs) for dynamic axial loading of the leg have been targeted toward automotive applications such as predicting injury caused by intrusion into the occupant compartment from frontal collisions. Recent focus on leg injuries in the military has led to questions about the applicability of these IRFs shorter duration, higher amplitude loading associated with underbody blast (UBB). To investigate these questions, data were collected from seven separate test series that subjected post-mortem human legs to axial impact. A force and impulse-based Weibull survival model was developed from these studies to estimate fracture risk. Specimen age was included as a covariate to reduce variance and improve survival model fit. The injury criterion estimated 50% risk of injury for a leg exposed to 13 N s of impulse at peak force and 8.07 kN of force for force durations less than and greater than half the natural period of the leg, respectively. A supplemental statistical analysis estimated that the proposed IRF improves injury prediction accuracy by more than 9% compared to the predictions from automobile-based risk functions developed for automotive intrusion. The proposed leg IRF not only improves injury prediction for higher rate conditions but also provides a single injury prediction tool for an expanded range of load durations ranging from 5 to 90 ms, which spans both automotive and military loading environments.

Development of a Second-Order System for Rapid Estimation of Maximum Brain Strain

Diffuse brain injuries are assessed with deformation-based criteria that utilize metrics based on rotational head kinematics to estimate brain injury severity. Although numerous metrics have been proposed, many are based on empirically-derived models that use peak kinematics, which often limit their applicability to a narrow range of head impact conditions. However, over a broad range of impact conditions, brain deformation response to rotational head motion behaves similarly to a second-order mechanical system, which utilizes the full kinematic time history of a head impact. This study describes a new brain injury metric called Diffuse Axonal Multi-Axis General Evaluation (DAMAGE). DAMAGE is based on the equations of motion of a three-degree-of-freedom, coupled 2nd-order system, and predicts maximum brain strain using the directionally dependent angular acceleration time-histories from a head impact. Parameters for the effective mass, stiffness, and damping were determined using simplified rotational pulses which were applied multiaxially to a 50th percentile adult human male finite element model. DAMAGE was then validated with a separate database of 1747 head impacts including helmet, crash, and sled tests and human volunteer responses. Relative to existing rotational brain injury metrics that were evaluated in this study, DAMAGE was found to be the best predictor of maximum brain strain.

Development of a Metric for Predicting Brain Strain Responses Using Head Kinematics

Diffuse brain injuries are caused by excessive brain deformation generated primarily by rapid rotational head motion. Metrics that describe the severity of brain injury based on head motion often do not represent the governing physics of brain deformation, rendering them ineffective over a broad range of head impact conditions. This study develops a brain injury metric based on the response of a second-order mechanical system, and relates rotational head kinematics to strain-based brain injury metrics: maximum principal strain (MPS) and cumulative strain damage measure (CSDM). This new metric, universal brain injury criterion (UBrIC), is applicable over a broad range of kinematics encountered in automotive crash and sports. Efficacy of UBrIC was demonstrated by comparing it to MPS and CSDM predicted in 1600 head impacts using two different finite element (FE) brain models. Relative to existing metrics, UBrIC had the highest correlation with the FE models, and performed better in most impact conditions. While UBrIC provides a reliable measurement for brain injury assessment in a broad range of head impact conditions, and can inform helmet and countermeasure design, an injury risk function was not incorporated into its current formulation until validated strain-based risk functions can be developed and verified against human injury data.

Sport-Related Concussion Alters Indices of Dynamic Cerebral Autoregulation

Sport-related concussion is known to affect a variety of brain functions. However, the impact of this brain injury on cerebral autoregulation (CA) is poorly understood. Thus, the goal of the current study was to determine the acute and cumulative effects of sport-related concussion on indices of dynamic CA. Toward this end, 179 elite, junior-level (age 19.6 ± 1.5 years) contact sport (ice hockey, American football) athletes were recruited for preseason testing, 42 with zero prior concussions and 31 with three or more previous concussions. Eighteen athletes sustained a concussion during that competitive season and completed follow-up testing at 72 h, 2 weeks, and 1 month post injury. Beat-by-beat arterial blood pressure (BP) and middle cerebral artery blood velocity (MCAv) were recorded using finger photoplethysmography and transcranial Doppler ultrasound, respectively. Five minutes of repetitive squat-stand maneuvers induced BP oscillations at 0.05 and 0.10 Hz (20- and 10-s cycles, respectively). The BP-MCAv relationship was quantified using transfer function analysis to estimate Coherence (correlation), Gain (amplitude ratio), and Phase (timing offset). At a group level, repeated-measures ANOVA indicated that 0.10 Hz Phase was significantly reduced following an acute concussion, compared to preseason, by 23% (-0.136 ± 0.033 rads) at 72 h and by 18% (-0.105 ± 0.029 rads) at 2 weeks post injury, indicating impaired autoregulatory functioning; recovery to preseason values occurred by 1 month. Athletes were cleared to return to competition after a median of 14 days (range 7-35), implying that physiologic dysfunction persisted beyond clinical recovery in many cases. When comparing dynamic pressure buffering between athletes with zero prior concussions and those with three or more, no differences were observed. Sustaining an acute sport-related concussion induces transient impairments in the capabilities of the cerebrovascular pressure-buffering system that may persist beyond 2 weeks and may be due to a period of autonomic dysregulation. Athletes with a history of three or more concussions did not exhibit impairments relative to those with zero prior concussions, suggesting recovery of function over time. Findings from this study support the potential need to consider physiological recovery in deciding when patients should return to play following a concussion.

A Novel Method for Quantifying Human In Situ Whole Brain Deformation under Rotational Loading Using Sonomicrometry

Traumatic brain injuries (TBI) are one of the least understood injuries to the body. Finite element (FE) models of the brain have been crucial for understanding concussion and for developing injury mitigation systems; however, the experimental brain deformation data currently used to validate these models are limited. The objective of this study was to develop a methodology for the investigation of in situ three-dimensional brain deformation during pure rotational loading of the head, using sonomicrometry. Sonomicrometry uses ultrasonic pulses to measure the dynamic distances between piezoelectric crystals implanted in any sound-transmitting media. A human cadaveric head-neck specimen was acquired 14 h postmortem and was instrumented with an array of 32 small sonomicrometry crystals embedded in the head: 24 crystals were implanted in the brain, and 8 were fixed to the inner skull. A dynamic rotation was then applied to the head using a closed-loop controlled test device. Four pulses with different severity levels were applied around three orthogonal anatomical axes of rotation. A repeated test of the highest severity rotation was conducted in each axis to assess repeatability. All tests were completed within 56 h postmortem. Overall, the combined experimental and sonomicrometry methods were demonstrated to reliably and repeatedly capture three-dimensional dynamic deformation of an intact human brain. These methods provide a framework for using sonomicrometry to acquire multidimensional experimental data required for FE model development and validation, and will lend insight into the deformations sustained by the brain during impact.