Head impact accelerations for brain strain-related responses in contact sports: a model-based investigation

On-field Head Acceleration Exposure Measurements Using Instrumented Mouthguards: Multi-stage Screening to Optimize Data Quality

Instrumented mouthguards (iMGs) are widely applied to measure head acceleration event (HAE) exposure in sports. Despite laboratory validation, on-field factors including potential sensor skull-decoupling and spurious recordings limit data accuracy. Video analysis can provide complementary information to verify sensor data but lacks quantitative kinematics reference information and suffers from subjectivity. The purpose of this study was to develop a rigorous multi-stage screening procedure, combining iMG and video as independent measurements, aimed at improving the quality of on-field HAE exposure measurements. We deployed iMGs and gathered video recordings in a complete university men's ice hockey varsity season. We developed a four-stage process that involves independent video and sensor data collection (Stage I), general screening (Stage II), cross verification (Stage III), and coupling verification (Stage IV). Stage I yielded 24,596 iMG acceleration events (AEs) and 17,098 potential video HAEs from all games. Approximately 2.5% of iMG AEs were categorized as cross-verified and coupled iMG HAEs after Stage IV, and less than 1/5 of confirmed or probable video HAEs were cross-verified with iMG data during stage III. From Stage I to IV, we observed lower peak kinematics (median peak linear acceleration from 36.0 to 10.9 g; median peak angular acceleration from 3922 to 942 rad/s2) and reduced high-frequency signals, indicative of potential reduction in kinematic noise. Our study proposes a rigorous process for on-field data screening and provides quantitative evidence of data quality improvements using this process. Ensuring data quality is critical in further investigation of potential brain injury risk using HAE exposure data.

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.

Potential of Soft-Shelled Rugby Headgear to Lower Regional Brain Strain Metrics During Standard Drop Tests

BACKGROUND

The growing concern for player safety in rugby has led to an increased focus on head impacts. Previous laboratory studies have shown that rugby headgear significantly reduces peak linear and rotational accelerations compared to no headgear. However, these metrics may have limited relevance in assessing the effectiveness of headgear in preventing strain-based brain injuries like concussions. This study used an instantaneous deep-learning brain injury model to quantify regional brain strain mitigation of rugby headgear during drop tests. Tests were conducted on flat and angled impact surfaces across different heights, using a Hybrid III headform and neck.

RESULTS

Headgear presence generally reduced the peak rotational velocities, with some headgear outperforming others. However, the effect on peak regional brain strains was less consistent. Of the 5 headgear tested, only the newer models that use open cell foams at densities above 45 kg/m3 consistently reduced the peak strain in the cerebrum, corpus callosum, and brainstem. The 3 conventional headgear that use closed cell foams at or below 45 kg/m3 showed no consistent reduction in the peak strain in the cerebrum, corpus callosum, and brainstem.

CONCLUSIONS

The presence of rugby headgear may be able to reduce the severity of head impact exposure during rugby. However, to understand how these findings relate to brain strain mitigation in the field, further investigation into the relationship between the impact conditions in this study and those encountered during actual gameplay is necessary.

In Silico Investigation of Biomechanical Response of a Human Brain Subjected to Primary Blast

The brain response to the explosion-induced primary blast waves is actively sought. Over the past decade, reasonable progress has been made in the fundamental understanding of blast traumatic brain injury (bTBI) using head surrogates and animal models. Yet, the current understanding of how blast waves interact with human is in nascent stages, primarily due to the lack of data in human. The biomechanical response in human is critically required to faithfully establish the connection to the aforementioned bTBI models. In this work, the biomechanical cascade of the brain under a primary blast has been elucidated using a detailed, full-body human model. The full-body model allowed us to holistically probe short- (<5 ms) and long-term (200 ms) brain responses. The full-body model has been extensively validated against impact loading in the past. We have further validated the head model against blast loading. We have also incorporated the structural anisotropy of the brain white matter. The blast wave transmission, and linear and rotational motion of the head were dominant pathways for the loading of the brain, and these loading paradigms generated distinct biomechanical fields within the brain. Blast transmission and linear motion of the head governed the volumetric response, whereas the rotational motion of the head governed the deviatoric response. Blast induced head rotation alone produced diffuse injury pattern in white matter fiber tracts. The biomechanical response under blast was comparable to the impact event. These insights will augment laboratory and clinical investigations of bTBI and help devise better blast mitigation strategies.

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.

Biomechanical Response of Head Surrogate With and Without the Helmet

Measurements of brain deformations under injurious loading scenarios are actively sought. In this work, we report experimentally measured head kinematics and corresponding dynamic, two-dimensional brain simulant deformations in head surrogates under a blunt impact, with and without a helmet. Head surrogates used in this work consisted of skin, skull, dura, falx, tentorium, and brain stimulants. The head surrogate geometry was based on the global human body models consortium's head model. A base head surrogate consisting of skin-skull-brain was considered. In addition, the response of two other head surrogates, skin-skull-dura-brain, and skin-skull-dura-brain-falx-tentorium, was investigated. Head surrogate response was studied for sagittal and coronal plane rotations for impactor velocities of 1 and 3 m/s. Response of head surrogates was compared against strain measurements in PMHS. The strain pattern in the brain simulant was heterogenous, and peak strains were established within ∼30 ms. The choice of head surrogate affect the spatiotemporal evolution of strain. For no helmet case, peak MPS of ∼50-60% and peak MSS of ∼35-50% were seen in brain simulant corresponding to peak rotational accelerations of ∼5000-7000 rad/s2. Peak head kinematics and peak MPS have been reduced by up to 75% and 45%, respectively, with the conventional helmet and by up to 90% and 85%, respectively, with the helmet with antirotational pads. Overall, these results provide important, new data on brain simulant strains under a variety of loading scenarios-with and without the helmets.

Design and characterization of 3-D printed hydrogel lattices with anisotropic mechanical properties

The goal of this study was to design, fabricate, and characterize hydrogel lattice structures with consistent, controllable, anisotropic mechanical properties. Lattices, based on three unit-cell types (cubic, diamond, and vintile), were printed using stereolithography (SLA) of polyethylene glycol diacrylate (PEGDA). To create structural anisotropy in the lattices, unit cell design files were scaled by a factor of two in one direction in each layer and then printed. The mechanical properties of the scaled lattices were measured in shear and compression and compared to those of the unscaled lattices. Two apparent shear moduli of each lattice were measured by dynamic shear tests in two planes: (1) parallel and (2) perpendicular to the scaling direction, or cell symmetry axis. Three apparent Young's moduli of each lattice were measured by compression in three different directions: (1) the "build" direction or direction of added layers, (2) the scaling direction, and (3) the unscaled direction perpendicular to both scaling and build directions. For shear deformation in unscaled lattices, the apparent shear moduli were similar in the two perpendicular directions. In contrast, scaled lattices exhibit clear differences in apparent shear moduli. In compression of unscaled lattices, apparent Young's moduli were independent of direction in cubic and vintile lattices; in diamond lattices Young's moduli differed in the build direction, but were similar in the other two directions. Scaled lattices in compression exhibited additional differences in apparent Young's moduli in the scaled and unscaled directions. Notably, the effects of scaling on apparent modulus differed between each lattice type (cubic, diamond, or vintile) and deformation mode (shear or compression). Scaling of 3D-printed, hydrogel lattices may be harnessed to create tunable, structures of desired shape, stiffness, and mechanical anisotropy, in both shear and compression.

EFFECTS OF EXERCISE TRAINING COMBINED WITH ACUPUNCTURE AND MOXIBUSTION ON FATIGUE

ABSTRACT Introduction: Muscle fatigue bothers athletes, affecting training level and competitive performance, it also has a great impact on the physical health of athletes, predisposing them to accidents and an early termination of their career. Relieving sports fatigue is the focus of research in the field of sports health nowadays. Objective: Study the effect of acupuncture and moxibustion rehabilitation combined with physical training on sports fatigue. Methods: A controlled experiment was used. The experimental group used acupuncture and moxibustion combined with exercise training, while the control group used acupuncture and moxibustion. The same group of doctors performed the acupuncture and moxibustion treatment according to the actual situation of the patients, and they tested the changes in VAS pain score and PRI pain score. Results: The VAS pain score in the experimental group was 7.88 points before the procedure and 2.96 points after the sixth week of the procedure. The control group score was 7.67 before the start and 5.03 after training. The total PRI pain score in the experimental group was 6.52 points before training and 2.05 points in the sixth week of training. The control group scored 6.66 before the procedure and 3.89 in the sixth week. Conclusion: The combination of training and exercises can achieve a better rehabilitation effect compared to the isolated treatment of acupuncture and moxibustion. Level of evidence II; Therapeutic studies - investigation of treatment outcomes.

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.

Drop Test Kinematics Using Varied Impact Surfaces and Head/Neck Configurations for Rugby Headgear Testing

World Rugby employs a specific drop test method to evaluate headgear performance, but almost all researchers use a different variation of this method. The aim of this study was, therefore, to quantify the differences between variations of the drop testing method using a Hybrid III headform and neck in the following impact setups: (1) headform only, with a flat steel impact surface, approximating the World Rugby method, (2 and 3) headform with and without a neck, respectively, onto a flat MEP pad impact surface, and (4) headform and neck, dropped onto an angled MEP pad impact surface. Each variation was subject to drop heights of 75-600 mm across three orientations (forehead, side, and rear boss). Comparisons were limited to the linear and rotational acceleration and rotational velocity for simplicity. Substantial differences in kinematic profile shape manifested between all drop test variations. Peak accelerations varied highly between variations, but the peak rotational velocities did not. Drop test variation also significantly changed the ratios of the peak kinematics to each other. This information can be compared to kinematic data from field head impacts and could inform more realistic impact testing methods for assessing headgear.

American Football Helmet Effectiveness Against a Strain-Based Concussion Mechanism

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

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

Real-time dynamic simulation remains a significant challenge for spatiotemporal data of high dimension and resolution. In this study, we establish a transformer neural network (TNN) originally developed for natural language processing and a separate convolutional neural network (CNN) to estimate five-dimensional (5D) spatiotemporal brain-skull relative displacement resulting from impact (isotropic spatial resolution of 4 mm with temporal resolution of 1 ms). Sequential training is applied to train (N = 5184 samples) the two neural networks for estimating the complete 5D displacement across a temporal duration of 60 ms. We find that TNN slightly but consistently outperforms CNN in accuracy for both displacement and the resulting voxel-wise four-dimensional (4D) maximum principal strain (e.g., root mean squared error (RMSE) of ~1.0% vs. ~1.6%, with coefficient of determination, R 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.

Head Impact Kinematics and Brain Deformation in Paired Opposing Youth Football Players

Head impact exposure is often quantified using peak resultant kinematics. While kinematics describes the inertial response of the brain to impact, they do not fully capture the dynamic brain response. Strain, a measure of the tissue-level response of the brain, may be a better predictor of injury. In this study, kinematic and strain metrics were compared to contact characteristics in youth football. Players on 2 opposing teams were instrumented with head impact sensors to record impact kinematics. Video was collected to identify contact scenarios involving opposing instrumented players (ie, paired contact scenarios) and code contact characteristics (eg, player role, impact location). A previously validated, high-resolution brain finite element model, the atlas-based brain model, was used to simulate head impacts and calculate strain metrics. Fifty-two paired contact scenarios (n = 105 impacts) were evaluated. Lighter players tended to have greater biomechanical metrics compared to heavier players. Impacts to the top of the helmet were associated with lower strain metrics. Overall, strain was better correlated with rotational kinematics, suggesting these metrics may be better predictors of the tissue-level brain response than linear kinematics. Understanding the effect of contact characteristics on brain strain will inform future efforts to improve sport safety.

Altered Brain Functional Connectivity in the Frontoparietal Network following an Ice Hockey Season

AbstractSustaining sports-related head impacts has been reported to result in neurological changes that potentially lead to later-life neurological disease. Advanced neuroimaging techniques have been used to detect subtle neurological effects resulting from head impacts, even after a single competitive season. The current study used resting-state functional magnetic resonance imaging to assess changes in functional connectivity of the frontoparietal network, a brain network responsible for executive functioning, in collegiate club ice hockey players over one season. Each player was scanned before and after the season and wore accelerometers to measure head impacts at practices and home games throughout the season. We examined pre- to post-season differences in connectivity within the frontoparietal and default mode networks, as well as the relationship between the total number of head impacts sustained and changes in connectivity. We found a significant interaction between network region of interest and time point (p = 0.016), in which connectivity between the left and right posterior parietal cortex seed regions increased over the season (p < 0.01). Number of impacts had a significant effect on frontoparietal network connectivity, such that more impacts were related to greater connectivity differences over the season (p = 0.042). Overall, functional connectivity increased in ice hockey athletes over a season between regions involved in executive functioning, and sensory integration, in particular. Furthermore, those who sustained more impacts had the greatest changes in connectivity. Consistent with prior findings in resting-state sports-related head impact literature, these findings have been suggested to represent brain injury.

A Finite Element Model of Cerebral Vascular Injury for Predicting Microbleeds Location

Finite Element (FE) models of brain mechanics have improved our understanding of the brain response to rapid mechanical loads that produce traumatic brain injuries. However, these models have rarely incorporated vasculature, which limits their ability to predict the response of vessels to head impacts. To address this shortcoming, here we used high-resolution MRI scans to map the venous system anatomy at a submillimetre resolution. We then used this map to develop an FE model of veins and incorporated it in an anatomically detailed FE model of the brain. The model prediction of brain displacement at different locations was compared to controlled experiments on post-mortem human subject heads, yielding over 3,100 displacement curve comparisons, which showed fair to excellent correlation between them. We then used the model to predict the distribution of axial strains and strain rates in the veins of a rugby player who had small blood deposits in his white matter, known as microbleeds, after sustaining a head collision. We hypothesised that the distribution of axial strain and strain rate in veins can predict the pattern of microbleeds. We reconstructed the head collision using video footage and multi-body dynamics modelling and used the predicted head accelerations to load the FE model of vascular injury. The model predicted large axial strains in veins where microbleeds were detected. A region of interest analysis using white matter tracts showed that the tract group with microbleeds had 95th percentile peak axial strain and strain rate of 0.197 and 64.9 s⁻¹ respectively, which were significantly larger than those of the group of tracts without microbleeds (0.163 and 57.0 s⁻¹). This study does not derive a threshold for the onset of microbleeds as it investigated a single case, but it provides evidence for a link between strain and strain rate applied to veins during head impacts and structural damage and allows for future work to generate threshold values. Moreover, our results suggest that the FE model has the potential to be used to predict intracranial vascular injuries after TBI, providing a more objective tool for TBI assessment and improving protection against it.

Impact and workload are dominating on-field data monitoring techniques to track health and well-being of team-sports athletes

Participation in sports has become an essential part of healthy living in today's world. However, injuries can often occur during sports participation. With advancements in sensor technology and data analytics, many sports have turned to technology-aided, data-driven, on-field monitoring techniques to help prevent injuries and plan better player management. This review searched three databases, Web of Science, IEEE, and PubMed, for peer-reviewed articles on on-field data monitoring techniques that are aimed at improving the health and well-being of team-sports athletes. It was found that most on-field data monitoring methods can be categorized as either player workload tracking or physical impact monitoring. Many studies covered during this review attempted to establish correlations between captured physical and physiological data, as well as injury risk. In these studies, workloads are frequently tracked to optimize training and prevent overtraining in addition to overuse injuries, while impacts are most often tracked to detect and investigate traumatic injuries. This review found that current sports monitoring practices often suffer from a lack of standard metrics and definitions. Furthermore, existing data-analysis models are created on data that are limited in both size and diversity. These issues need to be addressed to create ecologically valid approaches in the future.

Head Impact Sensor Triggering Bias Introduced by Linear Acceleration Thresholding

Contact sports players frequently sustain head impacts, most of which are mild impacts exhibiting 10-30 g peak head center-of-gravity (CG) linear acceleration. Wearable head impact sensors are commonly used to measure exposure and typically detect impacts using a linear acceleration threshold. However, linear acceleration across the head can substantially vary during 6-degree-of-freedom motion, leading to triggering biases that depend on sensor location and impact condition. We conducted an analytical investigation with impact characteristics extracted from on-field American football and soccer data. We assumed typical mouthguard sensor locations and evaluated whether simulated multi-directional impacts would trigger recording based on per-axis or resultant acceleration thresholding. Across 1387 impact directions, a 10g peak CG linear acceleration impact would trigger at only 24.7% and 31.8% of directions based on a 10 g per-axis and resultant acceleration threshold, respectively. Anterior impact locations had lower trigger rates and even a 30 g impact would not trigger recording in some directions. Such triggering biases also varied by sensor location and linear-rotational head kinematics coupling. Our results show that linear acceleration-based impact triggering could lead to considerable bias in head impact exposure measurements. We propose a set of recommendations to consider for sensor manufacturers and researchers to mitigate this potential exposure measurement bias.

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.

Head Impact Research Using Inertial Sensors in Sport: A Systematic Review of Methods, Demographics, and Factors Contributing to Exposure

BACKGROUND

The number and magnitude of head impacts have been assessed in-vivo using inertial sensors to characterise the exposure in various sports and to help understand their potential relationship to concussion.

OBJECTIVES

We aimed to provide a comprehensive review of the field of in-vivo sensor acceleration event research in sports via the summary of data collection and processing methods, population demographics and factors contributing to an athlete's exposure to sensor acceleration events.

METHODS

The systematic search resulted in 185 cohort or cross-sectional studies that recorded sensor acceleration events in-vivo during sport participation.

RESULTS

Approximately 5800 participants were studied in 20 sports using 18 devices that included instrumented helmets, headbands, skin patches, mouthguards and earplugs. Female and youth participants were under-represented and ambiguous results were reported for these populations. The number and magnitude of sensor acceleration events were affected by a variety of contributing factors, suggesting sport-specific analyses are needed. For collision sports, being male, being older, and playing in a game (as opposed to a practice), all contributed to being exposed to more sensor acceleration events.

DISCUSSION

Several issues were identified across the various sensor technologies, and efforts should focus on harmonising research methods and improving the accuracy of kinematic measurements and impact classification. While the research is more mature for high-school and collegiate male American football players, it is still in its early stages in many other sports and for female and youth populations. The information reported in the summarised work has improved our understanding of the exposure to sport-related head impacts and has enabled the development of prevention strategies, such as rule changes.

CONCLUSIONS

Head impact research can help improve our understanding of the acute and chronic effects of head impacts on neurological impairments and brain injury. The field is still growing in many sports, but technological improvements and standardisation of processes are needed.

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.

Orientation of neurites influences severity of mechanically induced tau pathology

Chronic traumatic encephalopathy is a neurodegenerative disease associated with repeated traumatic brain injury (TBI). Chronic traumatic encephalopathy is a tauopathy, in which cognitive decline is accompanied by the accumulation of neurofibrillary tangles of the protein tau in patients' brains. We recently found that mechanical force alone can induce tau mislocalization to dendritic spines and loss of synaptic function in in vitro neuronal cultures with random cell organization. However, in the brain, neurons are highly aligned, so here we aimed to determine how neuronal organization influences early-stage tauopathy caused by mechanical injury. Using microfabricated cell culture constructs to control the growth of neurites and an in vitro simulated TBI device to apply controlled mechanical deformation, we found that neuronal orientation with respect to the direction of a uniaxial high-strain-rate stretch injury influences the degree of tau pathology in injured neurons. We found that a mechanical stretch applied parallel to the neurite alignment induces greater mislocalization of tau proteins to dendritic spines than does a stretch with the same strain applied perpendicular to the neurites. Synaptic function, characterized by the amplitude of miniature excitatory postsynaptic currents, was similarly decreased in neurons with neurites aligned parallel to stretch, whereas in neurons aligned perpendicular to stretch, it had little to no functional loss. Experimental injury parameters (strain, strain rate, direction of stretch) were combined with a standard viscoelastic solid model to show that in our in vitro model, neurite work density during stretch correlates with tau mislocalization. These findings suggest that in a TBI, the magnitude of brain deformation is not wholly predictive of neurodegenerative consequences of TBI but that deformation relative to local neuronal architecture and the neurite mechanical energy during injury are better metrics for predicting trauma-induced tauopathy.

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.

The effect of large deformation on Poisson's ratio of brain white matter: An experimental study

A more Accurate description of the mechanical behavior of brain tissue could improve the results of computational models. While most studies have assumed brain tissue as an incompressible material with constant Poisson's ratio of almost 0.5 and constructed their modeling approach according to this assumption, the relationship between this ratio and levels of applied strains has not yet been studied. Since the mechanical response of the tissue is highly sensitive to the value of Poisson's ratio, this study was designed to investigate the characteristics of the Poisson's ratio of brain tissue at different levels of applied strains. Samples were extracted from bovine brain tissue and tested under unconfined compression at strain values of 5%, 10%, and 30%. Using an image processing method, the axial and transverse strains were measured over a 60-s period to calculate the Poisson's ratio for each sample. The results of this study showed that the Poisson's ratio of brain tissue at strain levels of 5% and 10% was close to 0.5, and assuming brain tissue as an incompressible material is a valid assumption at these levels of strain. For samples under 30% compression, this ratio was higher than 0.5, which could suggest that under strains higher than the brain injury threshold (approximately 18%), tissue integrity was impaired. Based on these observations, it could be concluded that for strain levels higher than the injury threshold, brain tissue could not be assumed as an incompressible material, and new material models need to be proposed to predict the material behavior of the tissue. In addition, the results showed that brain tissue under unconfined compression uniformly stretched in the transverse direction, and the bulging in the samples is negligible.

Hyper-viscoelastic mechanical behavior of cranial pia mater in tension

BACKGROUND

Cranial pia mater, the innermost layer of the meninges, protects the central nervous system by tightly wrapping the brain and damping the external impact force to the brain. Accurate experimental data of the mechanical property of the cranial pia mater can enhance the theoretical prediction of traumatic brain injury or the scientific surgery design for brain disease. The aim of this study is to characterize the mechanical behavior of the cranial pia mater.

METHODS

In vitro tensile and stress-relaxation experiments of ovine cranial pia mater specimens were conducted at eight strain rates to characterize the rate-dependent viscoelastic property. The tensile and stress-relaxation experimental data were fitted by an Ogden hyper-viscoelastic model with a strain rate function to describe the mechanical behavior of the cranial pia mater.

FINDINGS

The elastic modulus and the ultimate stress are significantly increased from 5.545 MPa and 0.535 MPa at 0.00167 s^-1 to 18.345 MPa and 2.547 MPa at 0.83 s^-1 (p < .0001), respectively. The initial stress and the long-term stress (300 s) are also increased significantly with the increasing strain rates (p < .0001). A good fit of the experimental data with the Ogden hyper-viscoelastic model incorporated with a strain rate function was achieved (R² > 0.93).

INTERPRETATION

The cranial pia mater exhibits as a rate-dependent hyper-viscoelastic material in the tensile and stress-relaxation experiments. Compared with the brain, the stiffer nature of the cranial pia mater indicates its essential role in brain protection. The rate-dependent constitutive model provides a proper description of the hyper-viscoelastic characteristics of the cranial pia mater in tension and may provide a basic constitutive relationship for numerical simulations of traumatic brain injury.

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

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

Structural Anisotropy vs. Mechanical Anisotropy: The Contribution of Axonal Fibers to the Material Properties of Brain White Matter

Brain's micro-structure plays a critical role in its macro-structure material properties. Since the structural anisotropy in the brain white matter has been introduced due to axonal fibers, considering the direction of axons in the continuum models has been mediated to improve the results of computational simulations. The aim of the current study was to investigate the role of fiber direction in the material properties of brain white matter and compare the mechanical behavior of the anisotropic white matter and the isotropic gray matter. Diffusion tensor imaging (DTI) was employed to detect the direction of axons in white matter samples, and tensile stress-relaxation loads up to 20% strains were applied on bovine gray and white matter samples. In order to calculate the nonlinear and time-dependent properties of white matter and gray matter, a visco-hyperelastic model was used. The results indicated that the mechanical behavior of white matter in two orthogonal directions, parallel and perpendicular to axonal fibers, are significantly different. This difference indicates that brain white matter could be assumed as an anisotropic material and axons have contribution in the mechanical properties. Also, up to 15% strain, white matter samples with axons parallel to the force direction are significantly stiffer than both the gray matter samples and white matter samples with axons perpendicular to the force direction. Moreover, the elastic moduli of white matter samples with axons both parallel and perpendicular to the loading direction and gray matter samples at 15-20% strain are not significantly different. According to these observations, it is suggested that axons have negligible roles in the material properties of white matter when it is loaded in the direction perpendicular to the axon direction. Finally, this observation showed that the anisotropy of brain tissue not only has effects on the elastic behavior, but also has effects on the viscoelastic behavior.

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.

Biomechanical Analysis of Woodpecker Response During Pecking Using a Two-Dimensional Computational Model

Traumatic brain injury (TBI) and chronic traumatic encephalopathy (CTE) due to the impact is a critical health concern. Impact mitigation strategy is a vital design paradigm to reduce the burden of TBI and CTE. In this regard, woodpecker biomimicry continues to attract attention. However, a direct comparison between a woodpecker and human biomechanical responses is lacking. Toward this end, we investigate the biomechanical response of a woodpecker during pecking using a two-dimensional head model. We also analyze the response of concurrent human head model to facilitate direct comparison with woodpecker response. The head models of woodpecker and human were built from medical images, the material properties were adopted from the literature. Both woodpecker and human head models were subjected to head kinematics obtained during pecking and resulting biomechanical response is studied. For the pecking cycle simulated in this work, peak rotational velocity and acceleration were ∼15 rad/s and 7,057 rad/s². These peak values are commensurate with the kinematics threshold values reported in human TBI. Our results show that, for the same input acceleration, the strains and stresses in the woodpecker brain are approximately six times lower than that of the human brain. The stress reduction is mainly attributed to the smaller size of the woodpecker head. The effect of pecking frequency and multiple pecking cycles have also been studied. It is observed that the strains and stresses in the brain are increased by ∼100% as pecking frequency is doubled. During multiple pecking cycle, dwell period of ∼90 ms tend to relax the stresses in the woodpecker brain; however, the amount of relaxation depends on the value of the decay constant. The comparison of biomechanical response against the axonal injury threshold suggests that for peak rotational acceleration of 7,057 rad/s² the maximum principal strain in the brains of woodpecker and human exceed the threshold limit. Acceleration scaling relationship between a woodpecker and equivalent human response is also developed as a function of head size. We obtain a scaling factor, ahaw, of 0.11 for baseline head sizes and a scaling factor of 1.03 as the human head size approaches woodpecker head size.

Tension Strain-Softening and Compression Strain-Stiffening Behavior of Brain White Matter

Brain, the most important component of the central nervous system (CNS), is a soft tissue with a complex structure. Understanding the role of brain tissue microstructure in mechanical properties is essential to have a more profound knowledge of how brain development, disease, and injury occur. While many studies have investigated the mechanical behavior of brain tissue under various loading conditions, there has not been a clear explanation for variation reported for material properties of brain tissue. The current study compares the ex-vivo mechanical properties of brain tissue under two loading modes, namely compression and tension, and aims to explain the differences observed by closely examining the microstructure under loading. We tested bovine brain samples under uniaxial tension and compression loading conditions, and fitted hyperelastic material parameters. At 20% strain, we observed that the shear modulus of brain tissue in compression is about 6 times higher than in tension. In addition, we observed that brain tissue exhibited strain-stiffening in compression and strain-softening in tension. In order to investigate the effect of loading modes on the tissue microstructure, we fixed the samples using a novel method that enabled keeping the samples at the loaded stage during the fixation process. Based on the results of histology, we hypothesize that during compressive loading, the strain-stiffening behavior of the tissue could be attributed to glial cell bodies being pushed against surroundings, contacting each other and resisting compression, while during tension, cell connections are detached and the tissue displays softening behavior.

A knowledge map analysis of brain biomechanics: Current evidence and future directions

Although brain, one of the most complex organs in the mammalian body, has been subjected to many studies from physiological and pathological points of view, there remain significant gaps in the available knowledge regarding its biomechanics. This article reviews the research trends in brain biomechanics with a focus on injury. We used published scientific articles indexed by Web of Science database over the past 40 years and tried to address the gaps that still exist in this field. We analyzed the data using VOSviewer, which is a software tool designed for scientometric studies. The results of this study showed that the response of brain tissue to external forces has been one of the significant research topics among biomechanicians. These studies have addressed the effects of mechanical forces on the brain and mechanisms of traumatic brain injury, as well as characterized changes in tissue behavior under trauma and other neurological diseases to provide new diagnostic and monitoring methods. In this study, some challenges in the field of brain injury biomechanics have been identified and new directions toward understanding the gaps in this field are suggested.

Mitigating the Consequences of Subconcussive Head Injuries

Subconcussive head injury represents a pathophysiology that spans the expertise of both clinical neurology and biomechanical engineering. From both viewpoints, the terms injury and damage, presented without qualifiers, are synonymously taken to mean a tissue alteration that may be recoverable. For clinicians, concussion is evolving from a purely clinical diagnosis to one that requires objective measurement, to be achieved by biomedical engineers. Subconcussive injury is defined as subclinical pathophysiology in which underlying cellular- or tissue-level damage (here, to the brain) is not severe enough to present readily observable symptoms. Our concern is not whether an individual has a (clinically diagnosed) concussion, but rather, how much accumulative damage an individual can tolerate before they will experience long-term deficit(s) in neurological health. This concern leads us to look for the history of damage-inducing events, while evaluating multiple approaches for avoiding injury through reduction or prevention of the associated mechanically induced damage.

A modular impact diverting mechanism for football helmets

To mitigate the injurious effect of the rotational acceleration of the brain, a modular Impact Diverting Mechanism (IDM) has been developed. The IDM can replace stickers (decals) that normally attach to the exterior of a football helmet. The IDM decals reduce friction and catch points between the covered area with the IDM on the outer shell of the helmet and the impacting surface, thereby decreasing rotational acceleration acting on the player's head. A Riddell Speed helmet's exterior was prepared with the IDM and outfitted to a headform equipped with linear accelerometers and gyroscopes. The helmets were tested at an impact velocity of 5.5 m/s at 15°, 30°, and 45° to the vertical: on the front, side, and back of the helmet. Results of 135 impact tests in the lab show that the IDM decal, when compared to helmets without it, reduced the rotational acceleration, rotational velocity, SI, HIC, and RIC ranging from 22% to 77%, 20% to 74%, 13% to 68%, 7% to 68%, 31% to 94%, respectively. Protection against rotational acceleration from oblique impacts is not prioritized in modern football helmets, as evident by current standard helmet testing protocols. This study demonstrates that the inclusion of the IDM decals in football helmets can help reduce the effects of rotational acceleration of the head during oblique impacts.

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.

Brain tissue strains vary with head impact location: A possible explanation for increased concussion risk in struck versus striking football players

BACKGROUND

On-field football helmet impacts over a large range of severities have caused concussions in some players but not in other players. One possible explanation for this variability is the struck player's helmet impact location.

METHODS

We examined the effect of impact location on regional brain tissue strain when input energy was held constant. Laboratory impacts were performed at 12 locations distributed over the helmet and the resulting head kinematics were simulated in two finite element models of the brain: the Simulated Injury Monitor and the Global Human Body Model Consortium brain model.

FINDINGS

Peak kinematics, injury metrics and brain strain varied significantly with impact location. Differences in impact location explained 33 to 37% of the total variance in brain strain for the whole brain and cerebrum, considerably more than the variance explained by impact location for the peak resultant head kinematics (8 to 23%) and slightly more than half of the variance explained by the difference in closing speed (57 to 61%). Both finite element models generated similar strain results, with minor variations for impacts that generated multi-axial rotations, larger variations in brainstem strains for some impact locations and a small bias for the cerebellum.

INTERPRETATION

Based on this experimental and computational simulation study, impact location on the football helmet has a large effect on regional brain tissue strain. We also found that the lowest strains consistently occurred in impacts to the crown and forehead, helmet locations commonly associated with the striking player.

Homogenization of heterogeneous brain tissue under quasi-static loading: a visco-hyperelastic model of a 3D RVE

Researches, in the recent years, reveal the utmost importance of brain tissue assessment regarding its mechanical properties, especially for automatic robotic tools, surgical robots and helmet producing. For this reason, experimental and computational investigation of the brain behavior under different conditions seems crucial. However, experiments do not normally show the distribution of stress and injury in microscopic scale, and due to various factors are costly. Development of micromechanical methods, which could predict the brain behavior more appropriately, could highly be helpful in reducing these costs. This study presents computational analysis of heterogeneous part of the brain tissue under quasi-static loading. Heterogeneity is created by irregular distribution of neurons in a representative volume element (RVE). Considering time-dependent behavior of the tissue, a visco-hyperelastic constitutive model is developed to predict the RVE behavior more realistically. The RVE is studied in different loads and load rates; 1, 2, 3, 10 and 15% strain load are applied at 0.03 and 0.2 s on the RVE as tensile and shear loads. Due to complexity in geometry, self-consistent approximation method is employed to increase the volume fraction of neurons and analyze RVE behavior in various NVFs. The results show increasing the load rate leads to a raise in the maximum stress that indicates the tissue is more vulnerable at higher rates. Moreover, stiffness of the tissue is enhanced in higher NVFs. Additionally, it is found that axons undergo higher stresses; hence, they are more sensitive in accidents which lead to axonal death and would cause TBI and DAI.

Evaluation of Brain Response during Head Impact in Youth Athletes Using an Anatomically Accurate Finite Element Model

During normal participation in football, players are exposed to repetitive subconcussive head impacts, or impacts that do not result in signs and symptoms of concussion. To better understand the effects of repetitive subconcussive impacts, the biomechanics of on-field head impacts and resulting brain deformation need to be well characterized. The current study evaluates local brain response to typical youth football head impacts using the atlas-based brain model (ABM), an anatomically accurate brain finite element (FE) model. Head impact kinematic data were collected from three local youth football teams using the Head Impact Telemetry (HIT) System. The azimuth and elevation angles were used to identify impacts near six locations of interest, and low, moderate, and high acceleration magnitudes (5th, 50th, and 95th percentiles, respectively) were calculated from the grouped impacts for FE simulation. Strain response in the brain was evaluated by examining the range and peak maximum principal strain (MPS) values in each element. A total of 40,538 impacts from 119 individual athletes were analyzed. Impacts to the facemask resulted in 0.18 MPS for the high magnitude impact category. This was 1.5 times greater than the oblique impact location, which resulted in the lowest strain value of 0.12 for high magnitude impacts. Overall, higher strains resulted from a 95th percentile lateral impact (41.0g, 2556 rad/sec²) with two predominant axes of rotation than from a 95th percentile frontal impact (67.6g, 2641 rad/sec²) with a single predominant axis of rotation. These findings highlight the importance of accounting for directional dependence and relative contribution of axes of rotation when evaluating head impact response.

Head Impact Kinematics Estimation With Network of Inertial Measurement Units

Wearable sensors embedded with inertial measurement units have become commonplace for the measurement of head impact biomechanics, but individual systems often suffer from a lack of measurement fidelity. While some researchers have focused on developing highly accurate, single sensor systems, we have taken a parallel approach in investigating optimal estimation techniques with multiple noisy sensors. In this work, we present a sensor network methodology that utilizes multiple skin patch sensors arranged on the head and combines their data to obtain a more accurate estimate than any individual sensor in the network. Our methodology visually localizes subject-specific sensor transformations, and based on rigid body assumptions, applies estimation algorithms to obtain a minimum mean squared error estimate. During mild soccer headers, individual skin patch sensors had over 100% error in peak angular velocity magnitude, angular acceleration magnitude, and linear acceleration magnitude. However, when properly networked using our visual localization and estimation methodology, we obtained kinematic estimates with median errors below 20%. While we demonstrate this methodology with skin patch sensors in mild soccer head impacts, the formulation can be generally applied to any dynamic scenario, such as measurement of cadaver head impact dynamics using arbitrarily placed sensors.

White Matter Anisotropy for Impact Simulation and Response Sampling in Traumatic Brain Injury

Advanced neuroimaging provides new opportunities to enhance head injury models, including the incorporation of white matter (WM) structural anisotropy. Information from high-resolution neuroimaging, however, usually has to be "down-sampled" to match a typically coarse brain mesh. To understand how this mesh-image resolution mismatch affects impact simulation and subsequent response sampling, we compared three competing anisotropy implementations (using either voxels, tractography, or a multiscale submodeling) and two response sampling strategies (element-wise or tractography-based, using brain mesh or neuroimaging for region segmentation, respectively). Using the combination of high resolution options as a baseline, we studied how the choice in each individual category affected the resulting injury metrics. By simulating a recorded loss of consciousness head impact, we found that injury metrics including peak strain and injury susceptibility in the deep WM regions based on fiber strain, but not on maximum principal strain, were sensitive to the anisotropy implementation, response sampling, and region segmentation. Overall, it was recommended to use tractography for anisotropy implementation and response sampling, and to employ neuroimaging for region segmentation, because they led to more accurate injury metrics. Further refining mesh locally via submodeling was unnecessary. Brain strain responses were also parametrically found to be closer to that from minimum fiber reinforcement, consistent with the fact that the majority of WM had a rather high degree of fiber dispersion. Finally, the upgraded Worcester Head Injury Model incorporating WM anisotropy was successfully re-validated against cadaveric impacts and an in vivo head rotation ("good" to "excellent" validation with an average Correlation Analysis score of 0.437 and 0.509, respectively). These investigations may facilitate further continual development of head injury models to better study traumatic brain injury.

Spinal constraint modulates head instantaneous center of rotation and dictates head angular motion

The head is kinematically constrained to the torso through the spine and thus, the spine dictates the amount of output head angular motion expected from an input impact. Here, we investigate the spinal kinematic constraint by analyzing the head instantaneous center of rotation (HICOR) with respect to the torso in head/neck sagittal extension and coronal lateral flexion during mild loads applied to 10 subjects. We found the mean HICOR location was near the C5-C6 intervertebral joint in sagittal extension, and T2-T3 intervertebral joint in coronal lateral flexion. Using the impulse-momentum relationship normalized by subject mass and neck length, we developed a non-dimensional analytical ratio between output angular velocity and input linear impulse as a function of HICOR location. The ratio was 0.65 and 0.50 in sagittal extension and coronal lateral flexion respectively, implying 30% greater angular velocities in sagittal extension given an equivalent impulse. Scaling to subject physiology also predicts larger required impulses given greater subject mass and neck length to achieve equivalent angular velocities, which was observed experimentally. Furthermore, the HICOR has greater motion in sagittal extension than coronal lateral flexion, suggesting the head and spine can be represented with a single inverted pendulum in coronal lateral flexion, but requires a more complex representation in sagittal extension. The upper cervical spine has substantial compliance in sagittal extension, and may be responsible for the complex motion and greater extension angular velocities. In analyzing the HICOR, we can gain intuition regarding the neck's role in dictating head motion during external loading.

Physical Biology of Axonal Damage

Excessive physical impacts to the head have direct implications on the structural integrity at the axonal level. Increasing evidence suggests that tau, an intrinsically disordered protein that stabilizes axonal microtubules, plays a critical role in the physical biology of axonal injury. However, the precise mechanisms of axonal damage remain incompletely understood. Here we propose a biophysical model of the axon to correlate the dynamic behavior of individual tau proteins under external physical forces to the evolution of axonal damage. To propagate damage across the scales, we adopt a consistent three-step strategy: First, we characterize the axonal response to external stretches and stretch rates for varying tau crosslink bond strengths using a discrete axonal damage model. Then, for each combination of stretch rates and bond strengths, we average the axonal force-stretch response of n = 10 discrete simulations, from which we derive and calibrate a homogenized constitutive model. Finally, we embed this homogenized model into a continuum axonal damage model of [1-d]-type in which d is a scalar damage parameter that is driven by the axonal stretch and stretch rate. We demonstrate that axonal damage emerges naturally from the interplay of physical forces and biological crosslinking. Our study reveals an emergent feature of the crosslink dynamics: With increasing loading rate, the axonal failure stretch increases, but axonal damage evolves earlier in time. For a wide range of physical stretch rates, from 0.1 to 10 /s, and biological bond strengths, from 1 to 100 pN, our model predicts a relatively narrow window of critical damage stretch thresholds, from 1.01 to 1.30, which agrees well with experimental observations. Our biophysical damage model can help explain the development and progression of axonal damage across the scales and will provide useful guidelines to identify critical damage level thresholds in response to excessive physical forces.

Estimated Brain Tissue Response Following Impacts Associated With and Without Diagnosed Concussion

Kinematic measurements of head impacts are sensitive to sports concussion, but not highly specific. One potential reason is these measures reflect input conditions only and may have varying degrees of correlation to regional brain tissue deformation. In this study, previously reported head impact data recorded in the field from high school and collegiate football players were analyzed using two finite element head models (FEHM). Forty-five impacts associated with immediately diagnosed concussion were simulated along with 532 control impacts without identified concussion obtained from the same players. For each simulation, intracranial response measures (max principal strain, strain rate, von Mises stress, and pressure) were obtained for the whole brain and within four regions of interest (ROI; cerebrum, cerebellum, brain stem, corpus callosum). All response measures were sensitive to diagnosed concussion; however, large inter-athlete variability was observed and sensitivity strength depended on measure, ROI, and FEHM. Interestingly, peak linear acceleration was more sensitive to diagnosed concussion than all intracranial response measures except pressure. These findings suggest FEHM may provide unique and potentially important information on brain injury mechanisms, but estimations of concussion risk based on individual intracranial response measures evaluated in this study did not improve upon those derived from input kinematics alone.

A brain impact stress analysis using advanced discretization meshless techniques

This work has the objective to compare the mechanical behaviour of a brain impact using an alternative numerical meshless technique. Thus, a discrete geometrical model of a brain was constructed using medical images. This technique allows to achieve a discretization with realistic geometry, allowing to define locally the mechanical properties according to the medical images colour scale. After defining the discrete geometrical model of the brain, the essential and natural boundary conditions were imposed to reproduce a sudden impact force. The analysis was performed using the finite element analysis and the radial point interpolation method, an advanced discretization technique. The results of both techniques are compared. When compared with the finite element analysis, it was verified that meshless methods possess a higher convergence rate and that they are capable of producing smoother variable fields.

Viscoelastic computational modeling of the human head-neck system: Eigenfrequencies and time-dependent analysis

A subject-specific 3-dimensional viscoelastic finite element model of the human head-neck system is presented and investigated based on computed tomography and magnetic resonance biomedical images. Ad hoc imaging processing tools are developed for the reconstruction of the simulation domain geometry and the internal distribution of bone and soft tissues. Material viscoelastic properties are characterized point-wise through an image-based interpolating function used then for assigning the constitutive prescriptions of a heterogenous viscoelastic continuum model. The numerical study is conducted both for modal and time-dependent analyses, compared with similar studies and validated against experimental evidences. Spatiotemporal analyses are performed upon different exponential swept-sine wave-localized stimulations. The modeling approach proposes a generalized, patient-specific investigation of sound wave transmission and attenuation within the human head-neck system comprising skull and brain tissues. Model extensions and applications are finally discussed.

Pilot Findings of Brain Displacements and Deformations during Roller Coaster Rides

With 300,000,000 riders annually, roller coasters are a popular recreational activity. Although the number of roller coaster injuries is relatively low, the precise effect of roller coaster rides on our brains remains unknown. Here we present the quantitative characterization of brain displacements and deformations during roller coaster rides. For two healthy adult male subjects, we recorded head accelerations during three representative rides, and, for comparison, during running and soccer headers. From the recordings, we simulated brain displacements and deformations using rigid body dynamics and finite element analyses. Our findings show that despite having lower linear accelerations than sports head impacts, roller coasters may lead to brain displacements and strains comparable to mild soccer headers. The peak change in angular velocity on the rides was 9.9 rad/sec, which was higher than the 5.6 rad/sec in soccer headers with ball velocities reaching 7 m/sec. Maximum brain surface displacements of 4.0 mm and maximum principal strains of 7.6% were higher than in running and similar to soccer headers, but below the reported average concussion strain. Brain strain rates during roller coaster rides were similar to those in running, and lower than those in soccer headers. Strikingly, on the same ride and at a similar position, the two subjects experienced significantly different head kinematics and brain deformation. These results indicate that head motion and brain deformation during roller coaster rides are highly sensitive to individual subjects. Although our study suggests that roller coaster rides do not present an immediate risk of acute brain injury, their long-term effects require further longitudinal study.