Three-dimensional human head finite-element model validation against two experimental impacts

Quantitative assessment of traumatic brain injury risk in diverse age groups of females: Insights from computational biomechanics

Traumatic Brain Injury (TBI) stands as a multifaceted health concern, exhibiting varying influences across human population. This study delves into the biomechanical complexities of TBI within gender-specific contexts, focusing on females. Our primary objective is to investigate distinctive injury mechanisms and risks associated with females, emphasizing the imperative for tailored investigations within this cohort. By employing Fluid-Structure Interaction (FSI) Analysis, we conducted simulations to quantify biomechanical responses to traumatic forces across diverse age groups of females. The study utilized a scaling technique to create finite element models (FEMs). The young female FEM, based on anthropometric data, showcased a 15 % smaller head geometry compared to the young male FEM. Moreover, while the elderly female FEM closely mirrored the young female FEM in most structural aspects, it showed distinctive features such as brain atrophy and increased cerebrospinal fluid (CSF) layer thickness. Notably, the child female FEM (ages 7-11 years) replicated around 95 % of the young female FEM's geometry. These structural distinctions meticulously captured age-specific variations across our modeled female age groups. It's noteworthy that identical conditions, encompassing impact intensity, loading type, and boundary conditions, were maintained across all FEMs in this biomechanical finite element analysis, ensuring comparative results. The findings unveiled significant variations in frontal and occipital pressures among diverse age groups, highlighting potential age-related discrepancies in TBI susceptibility among females. These variations were primarily linked to differences in anatomical features, including brain volume, CSF thickness, and brain condition, as the same material properties were used in the FEMs. These results were approximately 4.70, 6.33 and 6.43 % in frontal area of brain in diverse age groups of females (young, elderly, and child) respectively compared to young male FEM. Comparing the FEM results between the young female and the elderly female, we observed a decrease in occipital brain pressure at the same point, reducing from 171,993 to 167,793 Pa, marking an approximate 2.5 % decrease. While typically the elderly exhibit greater brain vulnerability compared to the young, our findings showcase a reduction in brain pressure. Notably, upon assessing the relative movement between the brain and the skull at the point located in occipital area, we observed greater relative movement in the elderly (1.8 mm) compared to the young female (1.04 mm). Therefore, brain atrophy increases the range of motion of the brain within the cranial space. The study underscores the critical necessity for nuanced TBI risk assessment tailored to age and gender, emphasizing the importance of age-specific protective strategies in managing TBIs across diverse demographics. Future research employing individual modeling techniques and exploring a wider age spectrum holds promise in refining our understanding of TBI mechanisms and adopting targeted approaches to mitigate TBI in diverse groups.

Revealing the role of material properties in impact-related injuries: Investigating the influence of brain and skull density variations on head injury severity

Traumatic brain injuries (TBI) resulting from head impacts are a major public health concern, which prompted our research to investigate the complex relationship between the material properties of brain tissue and the severity of TBI. The goal of this research is to investigate how variations in brain and skull density influence the vulnerability of brain tissue to traumatic injury, thereby enhancing our understanding of injury mechanism. To achieve this goal, we employed a well-validated finite element head model (FEHM). The current investigation was divided into two phases: in the first one, three distinct brain viscoelastic materials that had been utilized in prior studies were analyzed. The review of the properties of these three materials has been meticulous, encompassing both the spectrum of mechanical properties and the behaviors that are relevant to the way in which brain tissue reacts to traumatic loading conditions. In the second phase, the material properties of both the brain and skull tissue, alongside the impact conditions, were held constant. After this step, the focus was directed towards the variation of density in the brain and skull, which was consistent with the results obtained from previous experimental investigations, in order to determine the precise impact of these variations in density. This approach allowed a more profound comprehension of the impacts that density had on the simulation results. In the first phase, Material No. 2 exhibited the highest maximum first principal strain value in the frontal region (ε max = 15.41 % ), indicating lower stiffness to instantaneous deformation. This characteristic suggests that Material No. 2 may deform more extensively upon impact, potentially increasing the risk of injury due to its viscoelastic behavior. In contrast, Material No. 1, with a lower maximum first principal strain in the frontal region (ε max = 7.87 % ), displayed greater stiffness to instantaneous deformation, potentially reducing the risk of brain injury upon head impact. The second phase provided quantitative findings revealing a proportional relationship between brain tissue density and the pressures experienced by the brain. A 2 % increase in brain tissue density corresponded to approximately a 1 % increase in pressure on the brain tissue. Similarly, changes in skull density exhibited a similar quantitative relationship, with a 6 % increase in skull density leading to a 2.5 % increase in brain pressure. This preliminary approximate ratio of 2 to 1 between brain and skull density variations provides an initial quantitative framework for assessing the impact of density changes on brain vulnerability. These findings have several implications for the development of protective measures and injury prevention strategies, particularly in contexts where head trauma is a major issue.

An atlas of the heterogeneous viscoelastic brain with local power-law attenuation synthesised using Prony-series

This review addresses the acute need to acknowledge the mechanical heterogeneity of brain matter and to accurately calibrate its local viscoelastic material properties accordingly. Specifically, it is important to compile the existing and disparate literature on attenuation power-laws and dispersion to make progress in wave physics of brain matter, a field of research that has the potential to explain the mechanisms at play in diffuse axonal injury and mild traumatic brain injury in general. Currently, viscous effects in the brain are modelled using Prony-series, i.e., a sum of decaying exponentials at different relaxation times. Here we collect and synthesise the Prony-series coefficients appearing in the literature for twelve regions: brainstem, basal ganglia, cerebellum, corona radiata, corpus callosum, cortex, dentate gyrus, hippocampus, thalamus, grey matter, white matter, homogeneous brain, and for eight different mammals: pig, rat, human, mouse, cow, sheep, monkey and dog. Using this data, we compute the fractional-exponent attenuation power-laws for different tissues of the brain, the corresponding dispersion laws resulting from causality, and the averaged Prony-series coefficients. STATEMENT OF SIGNIFICANCE: Traumatic brain injuries are considered a silent epidemic and finite element methods (FEMs) are used in modelling brain deformation, requiring access to viscoelastic properties of brain. To the best of our knowledge, this work presents 1) the first multi-frequency viscoelastic atlas of the heterogeneous brain, 2) the first review focusing on viscoelastic modelling in both FEMs and experimental works, 3) the first attempt to conglomerate the disparate existing literature on the viscoelastic modelling of the brain and 4) the largest collection of viscoelastic parameters for the brain (212 different Prony-series spanning 12 different tissues and 8 different animal surrogates). Furthermore, this work presents the first brain atlas of attenuation power-laws essential for modelling shear waves in brain.

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

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

Development of detailed finite element models for in silico analyses of brain impact dynamics

BACKGROUND AND OBJECTIVE

In the last few decades, several studies have been performed to investigate traumatic brain injuries (TBIs) and to understand the biomechanical response of brain tissues, by using experimental and computational approaches. As part of computational approaches, human head finite element (FE) models show to be important tools in the analysis of TBIs, making it possible to estimate local mechanical effects on brain tissue for different accident scenarios. The present study aims to contribute to the computational approach by means of the development of three advanced FE head models for accurately describing the head tissue dynamics, the first step to predict TBIs.

METHODS

We have developed three detailed FE models of human heads from magnetic resonance images of three volunteers: an adult female (32 yrs), an adult male (35 yrs), and a young male (16 yrs). These models have been validated against experimental data of post mortem human subjects (PMHS) tests available in the literature. Brain tissue displacements relative to the skull, hydrostatic intracranial pressure, and head acceleration have been used as the parameters to compare the model response with the experimental response for validation. The software CORAplus (CORrelation and Analysis) has been adopted to evaluate the bio-fidelity level of FE models.

RESULTS

Numerical results from the three models agree with experimental data. FE models presented in this study show a good bio-fidelity for hydrostatic pressure (CORA score of 0.776) and a fair bio-fidelity brain tissue displacements relative to the skull (CORA score of 0.443 and 0.535). The comparison among numerical simulations carried out with the three models shows negligible differences in the mechanical state of brain tissue due to the different morphometry of the heads, when the same acceleration history is considered.

CONCLUSIONS

The three FE models, thanks to their accurate description of anatomical morphology and to their bio-fidelity, can be useful tools to investigate brain mechanics due to different impact scenarios. Therefore, they can be used for different purposes, such as the investigation of the correlation between head acceleration and tissue damage, or the effectiveness of helmet designs. This work does not address the issue to define injury thresholds for the proposed models.

The Effect of 3D Whole, Major, and Small Vasculature On Mouse Brain Strain Under Both Diffuse and Focal Brain Injury Loading

Blood vessels are much stiffer than brain parenchyma and their effects in finite element (FE) brain models need to be investigated. Despite the publication of some comprehensive three-dimensional (3D) brain vasculature models, no mechanical model exists for the mouse brain vasculature. Moreover, how the vasculature affects the mechanical behavior of brain tissue remains controversial. Therefore, we developed FE mouse brain models with detailed 3D vasculature to investigate the effect of the vasculature on brain strains under both diffuse (closed-head impact) and focal injury (controlled cortical impact (CCI)) loading, two commonly laboratory models of traumatic brain injury. The effect of the vasculature was examined by comparing maximum principal strain in mouse brain FE models with and without the vasculature. On average, modeling comprehensive vasculature under diffuse injury loading reduced average brain strain predictions by 32% with non-linear elastic properties. Nearly three-fourths of the 32% strain reduction was attributable to the effects of the major branches of the vasculature. Meanwhile, during focal open-skull CCI injury loading, the contribution of the vasculature was limited, producing a less than 5% reduction in all cases. Overall, the vasculature, especially the major branches, increased the load-bearing capacity of the brain FE model and thus reduced brain strain predictions.

Sudden death after facial impacts: Is the brainstem involved?

Three deaths following facial impacts in the presence of witnesses and resulting in brain lesions that were visualized only on pathological examination were studied at the forensic medicine institute of Marseille. Craniofacial impacts, even of low intensity, received during brawls may be associated with brain lesions ranging from a simple knock-out to fatal injuries. In criminal cases that are brought to court, even by autopsy it is still difficult to establish a direct link between the violence of the impact and the injuries that resulted in death. During a facial impact, the head undergoes a movement of violent forced hyperextension. Death may thus be secondary to the transmission of forces to the brain, either by a mechanism involving nerve conduction that may be termed a reflex mechanism (for example by vagal hyperstimulation) or by injury to the central nervous system (axonal damage). In such situations, autopsy does not make it possible to determine the cause of death, but only to suspect it in a context of voluntary violence in the presence of witnesses, with or without violent injury observed on external examination or on superficial incisions to determine the extent of bruises or hematoma. Systemic and comprehensive investigation involving pathology and toxicology is essential in any medicolegal case for positive interpretation and discrimination of other causes of death.

On the Sensitivity Analysis of Porous Finite Element Models for Cerebral Perfusion Estimation

Computational physiological models are promising tools to enhance the design of clinical trials and to assist in decision making. Organ-scale haemodynamic models are gaining popularity to evaluate perfusion in a virtual environment both in healthy and diseased patients. Recently, the principles of verification, validation, and uncertainty quantification of such physiological models have been laid down to ensure safe applications of engineering software in the medical device industry. The present study sets out to establish guidelines for the usage of a three-dimensional steady state porous cerebral perfusion model of the human brain following principles detailed in the verification and validation (V&V 40) standard of the American Society of Mechanical Engineers. The model relies on the finite element method and has been developed specifically to estimate how brain perfusion is altered in ischaemic stroke patients before, during, and after treatments. Simulations are compared with exact analytical solutions and a thorough sensitivity analysis is presented covering every numerical and physiological model parameter. The results suggest that such porous models can approximate blood pressure and perfusion distributions reliably even on a coarse grid with first order elements. On the other hand, higher order elements are essential to mitigate errors in volumetric blood flow rate estimation through cortical surface regions. Matching the volumetric flow rate corresponding to major cerebral arteries is identified as a validation milestone. It is found that inlet velocity boundary conditions are hard to obtain and that constant pressure inlet boundary conditions are feasible alternatives. A one-dimensional model is presented which can serve as a computationally inexpensive replacement of the three-dimensional brain model to ease parameter optimisation, sensitivity analyses and uncertainty quantification. The findings of the present study can be generalised to organ-scale porous perfusion models. The results increase the applicability of computational tools regarding treatment development for stroke and other cerebrovascular conditions.

Development and Application of Digital Human Models in the Field of Vehicle Collisions: A Review

In the human-vehicle-road system of collisions, the human is the most important factor, and digital human models (DHMs) are developed with the aim of preventing or at least reducing human injury. Because most of the relevant literature is focused mainly on collisions in traffic accidents (TAs), only some of the literature reviewed in this paper involves research results on other aspects of collisions. In this review, based on the background of DHMs and the application of DHMs regarding human injury biomechanics in collisions field, research results regarding the development of DHMs are described, the methods for verifying such models are introduced, and the application of the research results is discussed based on the aspect of human injury biomechanics. From the research literature, the development and validation of DHMs and their application in human injury biomechanics are summarized, and future research trends are proposed and discussed.

The importance of modeling the human cerebral vasculature in blunt trauma

BACKGROUND

Multiple studies describing human head finite element (FE) models have established the importance of including the major cerebral vasculature to improve the accuracy of the model predictions. However, a more detailed network of cerebral vasculature, including the major veins and arteries as well as their branch vessels, can further enhance the model-predicted biomechanical responses and help identify correlates to observed blunt-induced brain injury.

METHODS

We used an anatomically accurate three-dimensional geometry of a 50th percentile U.S. male head that included the skin, eyes, sinuses, spine, skull, brain, meninges, and a detailed network of cerebral vasculature to develop a high-fidelity model. We performed blunt trauma simulations and determined the intracranial pressure (ICP), the relative displacement (RD), the von Mises stress, and the maximum principal strain. We validated our detailed-vasculature model by comparing the model-predicted ICP and RD values with experimental measurements. To quantify the influence of including a more comprehensive network of brain vessels, we compared the biomechanical responses of our detailed-vasculature model with those of a reduced-vasculature model and a no-vasculature model.

RESULTS

For an inclined frontal impact, the predicted ICP matched well with the experimental results in the fossa, frontal, parietal, and occipital lobes, with peak-pressure differences ranging from 2.4% to 9.4%. For a normal frontal impact, the predicted ICP matched the experimental results in the frontal lobe and lateral ventricle, with peak-pressure discrepancies equivalent to 1.9% and 22.3%, respectively. For an offset parietal impact, the model-predicted RD matched well with the experimental measurements, with peak RD differences of 27% and 24% in the right and left cerebral hemispheres, respectively. Incorporating the detailed cerebral vasculature did not influence the ICP but redistributed the brain-tissue stresses and strains by as much as 30%. In addition, our detailed-vasculature model predicted strain reductions by as much as 28% when compared to current reduced-vasculature FE models that only include the major cerebral vessels.

CONCLUSIONS

Our study highlights the importance of including a detailed representation of the cerebral vasculature in FE models to more accurately estimate the biomechanical responses of the human brain to blunt impact.

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.

Finite element analysis to determine the cause of ring fractures in a motorcyclist's head

The finite element (FE) method can potentially help in reconstructing skull fracture biomechanisms, enabling differentiation of the injury patterns caused by traffic accidents. This study aims to (1) reconstruct a motorcycle driver-car accident case using the total human model for safety and FE simulations; and (2) analyze the biomechanisms of fatal ring fractures in the motorcyclist's skull base to determine if the fatal craniocerebral injuries were caused by a fall onto the highway after hitting a pedestrian or by the subsequent impact of a car. We simulated a series of loading scenarios of falls onto the road and impacts by a car, with and without a helmet being used. We reconstructed the injury processes and compared the biomechanics results to the skull tolerance limit. For the scenario of falling with a helmet, the Von-Mises stress around the foramen magnum indicated ring fractures with a slight fracture at the impact site, consistent with that detected in a traditional forensic pathology autopsy. Moreover, we found that a helmet can significantly protect the skull by controlling the increase in stress around the impact site. However, it has very little effect on the skull base, neck, or cervical spine. We determined that the characteristic ring fracture was most probably caused by the fall onto the highway. Thus, the subsequent car accident did not contribute to the motorcyclist's death. Our study demonstrates that the FE model and method can explore injury biomechanisms, assisting in the identification of injury patterns in forensic practices.

Study of cerebrospinal injuries by force transmission secondary to mandibular impacts using a finite element model

Brain and cervical injuries are often described after major facial impacts but rarely after low-intensity mandibular impacts. Force transmission to the brain and spinal cord from a mandibular impact such as a punch was evaluated by the creation and validation of a complete finite element model of the head and neck. Anteroposterior uppercut impacts on the jaw were associated with considerable extension and strong stresses at the junction of the brainstem and spinal cord. Hook punch impacts transmitted forces directly to the brainstem and the spinal cord without extension of the spinal cord. Deaths after this type of blow with no observed histological lesions may be related to excessive stressing of the brainstem, through which pass the sensory-motor pathways and the vagus nerve and which is the regulatory center of the major vegetative functions. Biological parameters are different in each individual, and by using digital modeling they can be modulated at will (jaw shape, dentition…) for a realistic approach to forensic applications.

Intracranial pressure-based validation and analysis of traumatic brain injury using a new three-dimensional finite element human head model

Traumatic brain injuries are life-threatening injuries that can lead to long-term incapacitation and death. Over the years, numerous finite element human head models have been developed to understand the injury mechanisms of traumatic brain injuries. Many of these models are erroneous and used ellipsoidal or spherical geometries to represent brain. This work is focused on the development of high-quality, comprehensive three-dimensional finite element human head model with accurate representation of cerebral sulci and gyri structures in order to study traumatic brain injury mechanisms. Present geometry, predicated on magnetic resonance imaging data consist of three rudimentary components, that is, skull, cerebrospinal fluid with the ventricular system, and the soft tissues comprising the cerebrum, cerebellum, and brain stem. The brain is modeled as a hyperviscoelastic material. Meshed model with 10 nodes modified tetrahedral type element (C3D10M) is validated against two cadaver-based impact experiments by comparing the intracranial pressures at different locations of the head. Our results indicate a better agreement with cadaver results, specifically for the case of frontal and parietal intracranial pressure values. Existing literature focuses mostly on intracranial pressure validation, while the effects of von Mises stress on brain injury are not analyzed in detail. In this work, a detailed interpretation of neurological damage resulting from impact injury is performed by analyzing von Mises stress and intracranial pressure distribution across numerous segments of the brain. A reasonably good correlation with experimental data signifies the robustness of the model for predicting injury mechanisms based on clinical predictions of injury tolerance criteria.

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

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

On the Development of Interspecies Traumatic Brain Injury Correspondence Rules

Traumatic brain injury analysis in human is exceedingly difficult due to the methods in which data can be collected, thus many researchers commonly implement animal surrogates. However, use of these surrogates is costly and restricted by ethical concerns and test logistics. Computational models and simulations do not have these constraints and can produce significant amounts of data in relatively short periods. This paper shows the development of a human head and neck model and a full body porcine model. Both models are developed from high-resolution CT and MRI scans and the latest low-to-high strain rate mechanical data available in the literature to represent tissue component material behaviors. Both models are validated against experiments from the literature and used to complete an initial interspecies correspondence rule development study for blast overpressure effects. The results indicate the similarities in the way injury develops in the pig brain and human brain but these similarities occur at very different insult levels. These results are extended by a study, which shows that blast peak pressure is the driving factor in injury prediction and, depending on the injury metric used, significantly different injuries could be predicted.

Development and validation of an optimized finite element model of the human orbit

INTRODUCTION

The authors' main purpose was to develop a detailed finite element model (FEM) of the human orbit and to validate it by analyzing its behavior under the stress of blunt traumas.

MATERIALS AND METHODS

A pre-existing 3D FEM of a human head was modified and used in this study. Modifications took into account preliminary research carried out on PubMed database. Data from a CT scan of the head were computed with Mimics^® software to re-create the skull geometry. The mesh production, the model's properties and the simulations of blunt orbital traumas were conducted on Hyperworks^® software.

RESULTS

The resulting 3D FEM was composed of 640 000 elements and was used to perform blunt trauma simulations on an intact orbit. A total of 27 tests were simulated. Fifteen tests were realized with a metallic cylinder impactor; 12 tests simulated a hit by a closed fist. In all the tests conducted (27/27), the orbital floor was fractured. Fracture patterns were similar to those found in real clinical situations according to the buckling and hydraulic theories of orbital floor fractures.

DISCUSSION

The similitude between the fracture patterns produced on the model and those observed in vivo allows for a validation of the model. This model constitutes, at the authors knowledge, the most sophisticated one ever developed.

Assessment of mechanical properties of human head tissues for trauma modelling

Many discrepancies are found in the literature regarding the damage and constitutive models for head tissues as well as the values of the constants involved in the constitutive equations. Their proper definition is required for consistent numerical model performance when predicting human head behaviour, and hence skull fracture and brain damage. The objective of this research is to perform a critical review of constitutive models and damage indicators describing human head tissue response under impact loading. A 3D finite element human head model has been generated by using computed tomography images, which has been validated through the comparison to experimental data in the literature. The threshold values of the skull and the scalp that lead to fracture have been analysed. We conclude that (1) compact bone properties are critical in skull fracture, (2) the elastic constants of the cerebrospinal fluid affect the intracranial pressure distribution, and (3) the consideration of brain tissue as a nearly incompressible solid with a high (but not complete) water content offers pressure responses consistent with the experimental data.

Do blast induced skull flexures result in axonal deformation?

Subject-specific computer models (male and female) of the human head were used to investigate the possible axonal deformation resulting from the primary phase blast-induced skull flexures. The corresponding axonal tractography was explicitly incorporated into these finite element models using a recently developed technique based on the embedded finite element method. These models were subjected to extensive verification against experimental studies which examined their pressure and displacement response under a wide range of loading conditions. Once verified, a parametric study was developed to investigate the axonal deformation for a wide range of loading overpressures and directions as well as varying cerebrospinal fluid (CSF) material models. This study focuses on early times during a blast event, just as the shock transverses the skull (< 5 milliseconds). Corresponding boundary conditions were applied to eliminate the rotation effects and the resulting axonal deformation. A total of 138 simulations were developed– 128 simulations for studying the different loading scenarios and 10 simulations for studying the effects of CSF material model variance–leading to a total of 10,702 simulation core hours. Extreme strains and strain rates along each of the fiber tracts in each of these scenarios were documented and presented here. The results suggest that the blast-induced skull flexures result in strain rates as high as 150–378 s^-1. These high-strain rates of the axonal fiber tracts, caused by flexural displacement of the skull, could lead to a rate dependent micro-structural axonal damage, as pointed by other researchers.

BIOMECHANICAL ANALYSIS OF OCCUPANT’S BRAIN RESPONSE AND INJURY IN VEHICLE INTERIOR SECOND IMPACT UTILIZING A REFINED HEAD FINITE ELEMENT MODEL

In vehicle side collisions, traumatic brain injury caused by the impact between occupant’s head and the interior parts of A or B pillar is a major reason of death and disability. In order to analyze the biomechanical response and injury mechanism of occupant’s brain in side collisions, a refined finite element head model representing the 50th percentile Chinese male was developed. Its improvements of biofidelity comparing to the original head model were illustrated through model simulation against the same post mortem human subjects test. Based on the refined head model, the brain biomechanical responses and injuries in the side impact with interior parts of A pillar and B pillar were analyzed according to FMVSS 201U, and the influences of different impact locations and directions were investigated. The results showed that the brain tissues on impact side sustained positive pressure and those on the opposite side experienced negative pressure. The transmission of pressure wave was easy to cause brain concussion and other diffuse brain injuries. The intracranial pressure distribution exhibited a typical pattern of contrecoup injury. The extreme stress concentration in the junction area of the cerebrum, cerebellum and brain stem could lead to focal injury such as brain contusion and laceration. Moreover, the impact injury of A pillar was more serious than that of B pillar, which was consistent with the traffic injury statistics that the head injury in oblique side collisions was more serious than that of vertical side collisions. Therefore, the interior parts of A pillar should be designed to absorb more energy than those of B pillar under the same conditions. In addition, the severity of brain injury is more sensitive to the variation of the horizontal angle than that of the vertical angle. Both the peak values of the occipital fossa pressure in effect simulations of the horizontal and vertical angles were three to four times of the peak values of the forehead pressure. When the impact horizontal angle was up to 255[Formula: see text], or the vertical angle was up to 45[Formula: see text], the head HIC(d) values would be up to 1320.45 and 1101.06, respectively, which indicated a AIS 3[Formula: see text] injury risk of the head.

Understanding how a sport-helmet protects the head from closed injury by virtual impact tests

Understanding how a helmet protects the head, especially the soft brain tissues, is the prerequisite for improving helmet design. Intracranial pressure and stresses/strains in the brain tissues are the direct indicators of traumatic brain injury and they can be used to measure helmet performance. In this study, the effects of helmet design parameters such as the helmet shell stiffness, liner compliance and thickness on the brain injury indicators were investigated by virtual impact tests. A finite element head model (FEHM) was first constructed from medical images; a personally-fitted helmet made of composite material and foam was virtually prototyped using geometric information extracted from the FEHM; a helmet-head finite element model was then assembled. Virtual impact tests were conducted using the resulting helmet-head model. The obtained results suggested that, if the helmet shell already has adequate strength to resist excessive deformation and fracture, further increasing shell stiffness and strength would not considerably reduce intracranial pressure and brain strains; to reach the maximum protection with the available materials, the key is to effectively use the second stage in the stress-strain history of the liner foam material.

An improved finite element modeling of the cerebrospinal fluid layer in the head impact analysis

The finite element (FE) method has been widely used to investigate the mechanism of traumatic brain injuries (TBIs), because it is technically difficult to quantify the responses of the brain tissues to the impact in experiments. One of technical challenges to build a FE model of a human head is the modeling of the cerebrospinal fluid (CSF) of the brain. In the current study, we propose to use membrane elements to construct the CSF layer. Using the proposed approach, we demonstrate that a head model can be built by using existing meshes available in commercial databases, without using any advanced meshing software tool, and with the sole use of native functions of the FE package Abaqus. The calculated time histories of the intracranial pressures at frontal, posterior fossa, parietal, and occipital positions agree well with the experimental data and the simulations in the literature, indicating that the physical effects of the CSF layer have been accounted for in the proposed modeling approach. The proposed modeling approach would be useful for bioengineers to solve practical problems.

BIOMECHANICS OF PORCINE BRAIN TISSUE UNDER FINITE COMPRESSION

While the developments in finite element method have made it possible to simulate the complex problems in biomechanics, building an accurate constitutive model of living tissues is a major factor in getting reliable finite element analysis (FEA) results. In this study, a set of experiments were performed to test the properties of porcine brain tissue under unconfined uniaxial compression with 20[Formula: see text]s hold time at varied strain levels (10–50%) and strain rates (0.1–1[Formula: see text]s[Formula: see text]). A novel method was developed to build a quasi-linear viscoelasticity (QLV) model. The elastic function and the relaxation function were calculated separately. An Odgen model was adopted to characterize the elastic behavior, while the relaxation response was modeled at six decay rates. A standard to choose the visco parameters was discussed and carried out. The disparity of parameter values in previous models was discussed and explained. It is suggested that this model should be used in the tested loading conditions (relaxation time, strain rate, strain levels).

Finite element simulations of the head-brain responses to the top impacts of a construction helmet: Effects of the neck and body mass

Traumatic brain injuries are among the most common severely disabling injuries in the United States. Construction helmets are considered essential personal protective equipment for reducing traumatic brain injury risks at work sites. In this study, we proposed a practical finite element modeling approach that would be suitable for engineers to optimize construction helmet design. The finite element model includes all essential anatomical structures of a human head (i.e. skin, scalp, skull, cerebrospinal fluid, brain, medulla, spinal cord, cervical vertebrae, and discs) and all major engineering components of a construction helmet (i.e. shell and suspension system). The head finite element model has been calibrated using the experimental data in the literature. It is technically difficult to precisely account for the effects of the neck and body mass on the dynamic responses, because the finite element model does not include the entire human body. An approximation approach has been developed to account for the effects of the neck and body mass on the dynamic responses of the head-brain. Using the proposed model, we have calculated the responses of the head-brain during a top impact when wearing a construction helmet. The proposed modeling approach would provide a tool to improve the helmet design on a biomechanical basis.

A multi-body dynamics study on a weight-drop test of rat brain injury

Traumatic brain injury (TBI), induced by impact of an object with the head, is a major health problem worldwide. Rats are a well-established animal analogue for study of TBI and the weight-drop impact-acceleration (WDIA) method is a well-established model in rats for creating diffuse TBI, the most common form of TBI seen in humans. However, little is known of the biomechanics of the WDIA method and, to address this, we have developed a four-degrees-of-freedom multi-body mass-spring-damper model for the WDIA test in rats. An analytical expression of the maximum skull acceleration, one of the important head injury predictor, was derived and it shows that the maximum skull acceleration is proportional to the impact velocity but independent of the impactor mass. Furthermore, a dimensional analysis disclosed that the maximum force on the brain and maximum relative displacement between brain and skull are also linearly proportional to impact velocity. Additionally, the effects of the impactor mass were examined through a parametric study from the developed multi-body dynamics model. It was found that increasing impactor mass increased these two brain injury predictors.

Effect of helmet liner systems and impact directions on severity of head injuries sustained in ballistic impacts: a finite element (FE) study

The current study aims to investigate the effectiveness of two different designs of helmet interior cushion, (Helmet 1: strap-netting; Helmet 2: Oregon Aero foam-padding), and the effect of the impact directions on the helmeted head during ballistic impact. Series of ballistic impact simulations (frontal, lateral, rear, and top) of a full-metal-jacketed bullet were performed on a validated finite element head model equipped with the two helmets, to assess the severity of head injuries sustained in ballistic impacts using both head kinematics and biomechanical metrics. Benchmarking with experimental ventricular and intracranial pressures showed that there is good agreement between the simulations and experiments. In terms of extracranial injuries, top impact had the highest skull stress, still without fracturing the skull. In regard to intracranial injuries, both the lateral and rear impacts generally gave the highest principal strains as well as highest shear strains, which exceed the injury thresholds. Off-cushion impacts were found to be at higher risk of intracranial injuries. The study also showed that the Oregon Aero foam pads helped to reduce impact forces. It also suggested that more padding inserts of smaller size may offer better protection. This provides some insights on future's helmet design against ballistic threats.

Biofidelic white matter heterogeneity decreases computational model predictions of white matter strains during rapid head rotations

The finite element (FE) brain model is used increasingly as a design tool for developing technology to mitigate traumatic brain injury. We developed an ultra high-definition FE brain model (>4 million elements) from CT and MRI scans of a 2-month-old pre-adolescent piglet brain, and simulated rapid head rotations. Strain distributions in the thalamus, coronal radiata, corpus callosum, cerebral cortex gray matter, brainstem and cerebellum were evaluated to determine the influence of employing homogeneous brain moduli, or distinct experimentally derived gray and white matter property representations, where some white matter regions are stiffer and others less stiff than gray matter. We find that constitutive heterogeneity significantly lowers white matter deformations in all regions compared with homogeneous properties, and should be incorporated in FE model injury prediction.

On the feasibility of life-saving locomotive bumpers

Motivated by the thousands of pedestrians killed each year in train impacts, this paper investigates the life-saving capability of four high-level locomotive bumper concepts. The head motions produced by the four concepts are modeled as one or two square acceleration pulses and are analyzed using the Head Injury Criterion (HIC). Surprisingly, the analyses show that all four concepts can achieve HIC values of less than 200 for an impact with a locomotive traveling at 100 km/h. Two of the concepts eject the pedestrian trackside with at a velocity of roughly 40 km/h and the risk of ground-impact injury is discussed in the context of related automobile accident data. The computed bumper lengths are a fraction of the overall length of a locomotive and are thus feasible for practical implementation. One concept involves an oblique impact and the potential for rotational head injury is analyzed. This basic feasibility research motivates future investigations into the detailed design of bumper shapes, multi-body pedestrian simulations, and finite-element injury models.

Head injury predictors in sports trauma--a state-of-the-art review

Head injuries occur in a great variety of sports. Many of these have been associated with neurological injuries, affecting the central nervous system. Some examples are motorsports, cycling, skiing, horse riding, mountaineering and most contact sports such as football, ice and field hockey, soccer, lacrosse, etc. The outcome of head impacts in these sports can be very severe. The worst-case scenarios of permanent disability or even death are possibilities. Over recent decades, many In recent decades, a great number of head injury criteria and respective thresholds have been proposed. However, the available information is much dispersed and a consensus has still not been achieved regarding the best injury criteria or even their thresholds. This review paper gives a thorough overview of the work carried out by the scientific community in the field of impact biomechanics about head injuries sustained during sports activity. The main goal is to review the head injury criteria, as well as their thresholds. Several are reviewed, from the predictors based on kinematics to the ones based on human tissue thresholds. In this work, we start to briefly introduce the head injuries and their mechanisms commonly seen as a result of head trauma in sports. Then, we present and summarize the head injury criteria and their respective thresholds.

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

Traumatic brain injury is a leading cause of disability and injury-related death. To enhance our ability to prevent such injuries, brain response can be studied using validated finite element (FE) models. In the current study, a high-resolution, anatomically accurate FE model was developed from the International Consortium for Brain Mapping brain atlas. Due to wide variation in published brain material parameters, optimal brain properties were identified using a technique called Latin hypercube sampling, which optimized material properties against three experimental cadaver tests to achieve ideal biomechanics. Additionally, falx pretension and thickness were varied in a lateral impact variation. The atlas-based brain model (ABM) was subjected to the boundary conditions from three high-rate experimental cadaver tests with different material parameter combinations. Local displacements, determined experimentally using neutral density targets, were compared to displacements predicted by the ABM at the same locations. Error between the observed and predicted displacements was quantified using CORrelation and Analysis (CORA), an objective signal rating method that evaluates the correlation of two curves. An average CORA score was computed for each variation and maximized to identify the optimal combination of parameters. The strongest relationships between CORA and material parameters were observed for the shear parameters. Using properties obtained through the described multiobjective optimization, the ABM was validated in three impact configurations and shows good agreement with experimental data. The final model developed in this study consists of optimized brain material properties and was validated in three cadaver impacts against local brain displacement data.

Variation of bone layer thicknesses and trabecular volume fraction in the adult male human calvarium

The human calvarium is a sandwich structure with two dense layers of cortical bone separated by porous cancellous bone. The variation of the three dimensional geometry, including the layer thicknesses and the volume fraction of the cancellous layer across the population, is unavailable in the current literature. This information is of particular importance to mathematical models of the human head used to simulate mechanical response. Although the target geometry for these models is the median geometry of the population, the best attempt so far has been the scaling of a unique geometry based on a few median anthropometric measurements of the head. However, this method does not represent the median geometry. This paper reports the average three dimensional geometry of the calvarium from X-ray computed tomography (CT) imaging and layer thickness and trabecular volume fraction from micro CT (μCT) imaging of ten adult male post-mortem human surrogates (PMHS). Skull bone samples have been obtained and μCT imaging was done at a resolution of 30 μm. Monte Carlo simulation was done to estimate the variance in these measurements due to the uncertainty in image segmentation. The layer thickness data has been averaged over areas of 5mm(2). The outer cortical layer was found to be significantly (p < 0.01; Student's t test) thicker than the inner layer (median of thickness ratio 1.68). Although there was significant location to location difference in all the layer thicknesses and volume fraction measurements, there was no trend. Average distribution and the variance of these metrics on the calvarium have been shown. The findings have been reported as colormaps on a 2D projection of the cranial vault.

EFFECT OF CEREBROSPINAL FLUID MODELED WITH DIFFERENT MATERIAL PROPERTIES ON A HUMAN FINITE ELEMENT HEAD MODEL

The aim of this study was to enhance head-injury prediction, this paper investigated the behavior of cerebrospinal fluid (CSF) in finite element (FE) modeling. Nine different material properties selected according to material definitions and property values were used to represent CSF in FE head models. To evaluate the influence of CSF material parameters on brain mechanical responses, the models were validated against available cadaver experiment data. Results showed that coup pressure increased whereas contrecoup pressure decreased when the head sustained purely translational impact with increased bulk modulus when CSF was modeled as fluid. However, with increased bulk modulus, coup pressure, contrecoup pressure and relative skull-brain motions decreased under rotational impact. When CSF was assumed to be an elastic material, coup pressure increased whereas contrecoup pressure decreased with increased elastic modulus when the head was subjected to purely translational impact. However, the variation trend was not obvious during head rotation. Results also indicated that when subjected to brain strain and von Mises stress, the model was prone to underestimate brain injury when CSF was modeled as an elastic material, especially during purely translational impact to the head. The model with CSF as fluid reduced the strain rate of brain, which was more likely to be realistic than the model with CSF as a viscoelastic material. These findings suggested that significantly higher values of the bulk modulus of CSF modeled as fluid were needed to predict intracranial dynamic responses and brain injury during head impact.

Investigation of the relationship between facial injuries and traumatic brain injuries using a realistic subject-specific finite element head model

In spite of anatomic proximity of the facial skeleton and cranium, there is lack of information in the literature regarding the relationship between facial and brain injuries. This study aims to correlate brain injuries with facial injuries using finite element method (FEM). Nine common impact scenarios of facial injuries are simulated with their individual stress wave propagation paths in the facial skeleton and the intracranial brain. Fractures of cranio-facial bones and intracranial injuries are evaluated based on the tolerance limits of the biomechanical parameters. General trend of maximum intracranial biomechanical parameters found in nasal bone and zygomaticomaxillary impacts indicates that severity of brain injury is highly associated with the proximity of location of impact to the brain. It is hypothesized that the midface is capable of absorbing considerable energy and protecting the brain from impact. The nasal cartilages dissipate the impact energy in the form of large scale deformation and fracture, with the vomer-ethmoid diverging stress to the "crumpling zone" of air-filled sphenoid and ethmoidal sinuses; in its most natural manner, the face protects the brain. This numerical study hopes to provide surgeons some insight in what possible brain injuries to be expected in various scenarios of facial trauma and to help in better diagnosis of unsuspected brain injury, thereby resulting in decreasing the morbidity and mortality associated with facial trauma.

A QCT-Based Nonsegmentation Finite Element Head Model for Studying Traumatic Brain Injury

In the existing finite element head models (FEHMs) that are constructed from medical images, head tissues are usually segmented into a number of components according to the interior anatomical structure of the head. Each component is represented by a homogenous material model. There are a number of disadvantages in the segmentation-based finite element head models. Therefore, we developed a nonsegmentation finite element head model with pointwise-heterogeneous material properties and corroborated it by available experiment data. From the obtained results, it was found that although intracranial pressures predicted by the existing (piecewise-homogeneous) and the proposed (pointwise-heterogeneous) FEHM are very similar to each other, strain/stress levels in the head tissues are very different. The maximum peak strains/stresses predicted by the proposed FEHM are much higher than those by the existing FEHM, indicating that piecewise-homogeneous FEHM may have underestimated the stress/strain level induced by impact and thus may be inaccurate in predicting traumatic brain injuries.

Development of a finite element head model for the study of impact head injury

This study is aimed at developing a high quality, validated finite element (FE) human head model for traumatic brain injuries (TBI) prediction and prevention during vehicle collisions. The geometry of the FE model was based on computed tomography (CT) and magnetic resonance imaging (MRI) scans of a volunteer close to the anthropometry of a 50th percentile male. The material and structural properties were selected based on a synthesis of current knowledge of the constitutive models for each tissue. The cerebrospinal fluid (CSF) was simulated explicitly as a hydrostatic fluid by using a surface-based fluid modeling method. The model was validated in the loading condition observed in frontal impact vehicle collision. These validations include the intracranial pressure (ICP), brain motion, impact force and intracranial acceleration response, maximum von Mises stress in the brain, and maximum principal stress in the skull. Overall results obtained in the validation indicated improved biofidelity relative to previous FE models, and the change in the maximum von Mises in the brain is mainly caused by the improvement of the CSF simulation. The model may be used for improving the current injury criteria of the brain and anthropometric test devices.

A computational study on brain tissue under blast: primary and tertiary blast injuries

In this paper, a biomechanical study of a human head model exposed to blast shock waves followed by a blunt impact with the surface of the enclosing walls of a confined space is carried out. Under blast, the head may experience primary blast injury (PBI) due to exposure to the shockwaves and tertiary blast injury (TeBI) due to a possible blunt impact. We examine the brain response data in a deformable finite element head model in terms of the inflicted stress/pressure, velocity, and acceleration on the brain for several blast scenarios with different intensities. The data will be compared for open space and confined spaces. Following the initial impact of the shock front in the confined space, one can see the fluctuations in biomechanical data due to wave reflections. Although the severity of the PBI and TeBI is dependent on the situation, for the cases studied here, PBI is considerably more pronounced than TeBI in confined spaces.

A modified controlled cortical impact technique to model mild traumatic brain injury mechanics in mice

For the past 25 years, controlled cortical impact (CCI) has been a useful tool in traumatic brain injury (TBI) research, creating injury patterns that includes primary contusion, neuronal loss, and traumatic axonal damage. However, when CCI was first developed, very little was known on the underlying biomechanics of mild TBI. This paper uses information generated from recent computational models of mild TBI in humans to alter CCI and better reflect the biomechanical conditions of mild TBI. Using a finite element model of CCI in the mouse, we adjusted three primary features of CCI: the speed of the impact to achieve strain rates within the range associated with mild TBI, the shape, and material of the impounder to minimize strain concentrations in the brain, and the impact depth to control the peak deformation that occurred in the cortex and hippocampus. For these modified cortical impact conditions, we observed peak strains and strain rates throughout the brain were significantly reduced and consistent with estimated strains and strain rates observed in human mild TBI. We saw breakdown of the blood–brain barrier but no primary hemorrhage. Moreover, neuronal degeneration, axonal injury, and both astrocytic and microglia reactivity were observed up to 8 days after injury. Significant deficits in rotarod performance appeared early after injury, but we observed no impairment in spatial object recognition or contextual fear conditioning response 5 and 8 days after injury, respectively. Together, these data show that simulating the biomechanical conditions of mild TBI with a modified cortical impact technique produces regions of cellular reactivity and neuronal loss that coincide with only a transient behavioral impairment.

Development and validation of two subject-specific finite element models of human head against three cadaveric experiments

Head injury, being one of the main causes of death or permanent disability, continues to remain a major health problem with significant socioeconomic costs. Numerical simulations using the FEM offer a cost-effective method and alternative to experimental methods in the biomechanical studies of head injury. The present study aimed to develop two realistic subject-specific FEMs of the human head with detailed anatomical features from medical images (Model 1: without soft tissue and Model 2: with soft tissue and differentiation of white and gray matters) and to validate them against the intracranial pressure (ICP) and relative intracranial motion data of the three cadaver experimental tests. In general, both the simulated results were in reasonably good agreement with the experimental measured ICP and relative displacements, despite slight discrepancy in a few neutral density targets markers. Sensitivity analysis showed some variations in the brain's relative motion to the material properties or marker's location. The addition of soft tissue in Model 2 helped to damp out the oscillations of the model response. It was also found that, despite the fundamental anatomical differences between the two models, there existed little evident differences in the predicted ICP and relative displacements of the two models. This indicated that the advancements on the details of the extracranial features would not improve the model's predicting capabilities of brain injury.